Quantum correlations across two octaves from combined up and down conversion
Abstract
We propose and analyse a cascaded optical parametric system which involves three interacting modes across two octaves of frequency difference. Our system, combining degenerate optical parametric oscillation (OPO) with second harmonic generation (SHG), promises to be a useful source of squeezed and entangled light at three differing frequencies. We show how changes in damping rates and the ratio of the two concurrent nonlinearities affect the quantum correlations in the output fields. We analyse the threshold behaviour, showing how the normal OPO threshold is changed by the addition of the SHG interactions. We also find that the inclusion of the OPO interaction removes the self-pulsing behaviour found in normal SHG. Finally, we show how the Einstein-Podolsky-Rosen correlations can be controlled by the injection of a coherent seed field at the lower frequency.
I Introduction
The theory of the interaction of light fields at one frequency with nonlinear materials to produce fields at different frequencies goes back at least to Armstrong et al. and their seminal work which included downconversion and second and third harmonic generation Armstrong . Since the publication of that work, the optical parametric oscillator (OPO) in both its degenerate and nondegenerate forms Giordmaine ; Nassau ; Yariv has become a standard workhorse for quantum optics and quantum information, especially with respect to the Einstein Podolsky Rosen paradox RMPMargaret . The related process of intracavity second harmonic generation (SHG) has also long been known to produce quantum states of the optical field SHGPereira .
In the degenerate OPO, any entanglement will necessarily be across one octave, with the same being true of SHG sumdiff ; PingKoy . In this work we combine these two processes in either a cascaded or concurrent manner, to produce entangled beams and states exhibiting EPR steering across two octaves of frequency difference. Such a difference in frequencies has previously been predicted for a system which cascades two SHG processes to produce entangled outputs at three different frequencies, with both bipartite 4HG and tripartite correlations 4HGtri . The three level system we analyse here differs essentially only in the choice of cavity field which is externally pumped. In these previous two octave systems, this was the field at the lowest frequency. In this work it is the field at the intermediate frequency which is pumped. Just as with the normal OPO and SHG processes, this small change leads to markedly different behaviours.
The system we analyse has the potential to provide enhanced flexibility for quantum interfaces between light and atomic ensembles, quantum state engineering, multiplexing in quantum communications multiplex , the entanglement of atomic ensembles, and quantum teleportation Hammerer . The availability of entanglement and EPR-steering over such a large frequency range will bring further flexibility to the linking of quantum processes at different wavelengths, for example the telecommunications frequencies and atomic systems used in quantum information processing, particularly with regard to quantum memory Julsgaard .
In this article we first provide the Hamiltonian, then develop the equations of motion in the positive-P representation P+ . These equations are then solved numerically to find the time evolution of the intracavity fields. We check the full quantum numerical results against those found analytically for the classical steady states, finding that these agree in most parameter regimes. One regime where they do not agree is that in which the classical solutions exhibit self pulsing behaviour. In other regimes we use the steady state solutions for a linearised fluctuation analysis. This allows us to find the oscillation threshold, which is changed from that in the standard OPO. Using the standard input-output relations mjc , we are able to calculate the expressions for squeezing and both bipartite and tripartite EPR steering and inseparability in the output modes. In cases where the output expressions are rather simple, we give these analytically. In other cases the results are produced graphically. We look at the effects of changing the ratio of the two nonlinearities and the cavity damping rates. Finally we examine the effects of an injected signal at the lowest frequency. The range of interesting quantum states found suggests that this system shows promise for emerging quantum technological applications.
Ii Hamiltonian and equations of motion
The system we investigate here uses two nonlinear interactions within the same pumped optical cavity which is resonant for all three frequencies of interest. These could be either two crystals or one customised dielectric Zhu which converts the input field via both up and down conversion. The three interacting electromagnetic fields are the central externally pumped field at frequency , and two others at and . The field at interacts via a nonlinearity represented by to produce a downconverted field at , where . It also interacts via the nonlinearity represented by to produce an upconverted field at . This field is therefore the fourth harmonic of , with the interacting fields spanning two octaves of frequency difference.
The low frequency field at , is represented by the bosonic operator . The second harmonic, at , which will be externally pumped, is represented by , and the fourth harmonic, at , is represented by . The unitary interaction Hamiltonian in a rotating frame is then written as
(1) |
Since we are analysing the intracavity configuration, we also have the pumping Hamiltonian,
(2) |
where represents an external pumping field which is usually taken as coherent, although this is not necessary Liz . The damping of the cavity into a zero temperature Markovian reservoir is described by the Lindblad superoperator
(3) |
where is the system density matrix and is the cavity loss rate at . In this work we will treat all three optical fields as being at resonance with the optical cavity. While including detuning is possible, this makes analytical results very difficult to obtain, so we will stick to the simplest case here. In general, any detuning acts to degrade the correlations used to measure squeezing and entanglement in a system Granja .
In order to analyse this system, we will use the well known and exact quantum phase space method, the positive-P representation P+ , which allows us to readily calculate any time-normally-ordered operator moments. Following the usual procedures DFW , we derive equations of motion in the positive-P representation P+ ,
(4) |
It should be noted that these have the same form in either Itô or Stratonovich calculus SMCrispin . In the above, the complex variable pairs correspond to the operator pairs in the sense that stochastic averages of products converge to normally-ordered operator expectation values, e.g. . The are Gaussian noise terms with the properties and . Although there can be divergence problems with the positive-P representation, it is known to be accurate where it converges, which is the case with all results presented here.
Iii Steady-state and threshold properties
In order to obtain analytical steady-state results for the intracavity intensities and amplitudes, we solve the semi-classical equivalents of Eq. 4, simply obtained by removing the noise terms. The results thus obtained can be checked against stochastic integration of the full equations. This procedure also allows us to calculate the threshold pumping value at which the downconversion process begins to produce non-zero amplitudes in the low frequency mode. This threshold behaviour is well known from the theory of the optical parametric oscillator (OPO) DMW ; Arabe . A stability analysis of the system allows the threshold pumping amplitude to be calculated as
(5) |
We immediately see that this is higher than the threshold for isolated downconversion, where the threshold is . The increased pump power is required because the upconversion process to produce the mode at also depletes the pump in our system.
The steady state amplitudes for the three modes can be found in the two different cases:
(i) below threshold ,
(6) |
where
(7) |
and
(ii) above threshold ,
(8) |
As with the standard OPO, the system exhibits similar behaviour to a second-order phase transition at . When the pumping is above threshold, the below-threshold solution for the fundamental frequency field becomes unstable and the system moves onto a new stable branch witht two solutions of the the fundamental field having equal amplitude and opposite phase. The steady amplitudes of the central frequency and higher frequency modes have opposite phases whether the system is running below or above threshold. What is noticeable is that the steady state solutions above threshold for and have no dependence on the pump power. Once the cavity is being pumped above the oscillation threshold, these two fields do not change with changes in the pumping.
The time development of the intensities above threshold is shown in Fig. 1 in the fully quantum picture with the positive-P equations integrated over stochastic trajectories. With we see that the analytical steady-state values, plotted as dashed lines, are in good agreement with the quantum solutions.
It is also well known that in normal second harmonic generation (SHG) there is a pumping threshold above which the output intensities exhibit a periodic pulsing behaviour pulse ; Bache . In the present case the classical behaviour of the system is similar and a hard mode transition can be found above which self-pulsing occurs. However, this does not survive the full quantum treatment, with the oscillations disappearing completely. A less pronounced damping of self-pulsing oscillations has recently been found in a full quantum treatment of other cascaded systems 3HG ; 4HG and shows the dangers of relying on classical analyses of quantum optical systems. The canonical method to calculate self pulsing in SHG is to numerically integrate the classical equations with a small complex seed in one or both the modes. Without this seed, the self-pulsing is not found, although it appears with integration of the positive-P equations without needing any seed at all. For our system, small complex seeds in the initial condition of the classical simulations gives self-pulsing, as shown in Fig. 2. On the other hand, the quantum solution diverges from this at short times, to enter a steady state with a much lower average value. The reason for this is that the classical solutions stay on the unstable branch of the solutions for , remaining at zero. The classical solution is unphysical. In the quantum case, spontaneous downconversion early in the evolution leads to stimulated downconversion and the steady state remains on the stable branch. A small injected signal in the classical integration will also push the solutions onto the stable branch, and in this case self-pulsing is found neither classically nor quantum mechanically.
Iv Ornstein-Uhlenbeck analysis and fluctuation spectra
When nonlinear optical media are held inside a pumped optical cavity, the accessible observables are usually the output spectral correlations, which are accessible using homodyne measurement techniques mjc . These are readily calculated in the steady state by treating the system as an Ornstein-Uhlenbeck process SMCrispin . In order to do this, we begin by expanding the positive-P variables into their steady-state expectation values plus delta-correlated Gaussian fluctuation terms, e.g.
(9) |
Given that we can calculate the , we may now write the equations of motion for the fluctuation terms. The resulting equations are written for the vector of fluctuation terms as
(10) |
where is the drift matrix containing the steady-state solution, is found from the factorisation of the drift matrix of the original Fokker-Planck equation, , with the steady-state values substituted in, and is a vector of Wiener increments. As long as the matrix has no eigenvalues with negative real parts, this method may be used to calculate the intracavity spectra via
(11) |
from which the output spectra are calculated using the standard input-output relations mjc .
In this case, is found as
(12) |
and is a matrix with on the diagonal. In the above, the should be read as their steady-state mean values, so that , for example. These are now complex numbers that are the averages of the positive-P stochastic variables. Because we have parametrised our system using , the frequency is in units of . is now in terms of quadratic products of the fluctuation operators such as and .
Since quadrature properties are what is measured by homodyne detection, we define the amplitude and phase quadrature operators as
(13) |
We note here that other definitions are sometimes used in the literature and that this changes the numerical value of the Heisenberg uncertainty principle. Our choice gives and means that squeezing in a particular quadrature exists whenever its variance is found to be less than 1.
To express the fluctuation expressions in terms of the canonical quadratures, we calculate
(14) |
where is the block diagonal matrix constructed from
(15) |
gives us the products from which we construct the output variances and covariances for modes and as,
(16) |
in which the variances and covariances are defined as and
V Steady state bipartite correlations
The squeezing in the amplitude and phase quadrature for the three different modes can be calculated analytically following from Eq. 6 and Eq. 14. Since the fundamental mode has a mean amplitude of zero below threshold, this simplifies the drift matrix and we can derive the below threshold output squeezing spectra as
(17) |
where
(18) |
and The spectral variances of the squeezed quadratures are shown in Fig. 3, for . We note here that all spectra shown are symmetric about zero frequency. What we notice is that the quadratures which exhibit squeezing are those we expect from parametric downconversion, with being squeezed, and from second harmonic generation, with both and being squeezed.
Above threshold the analytical expressions for the output squeezing are quite lengthy, mainly due to that fact that the low frequency mode now has a non-zero solution. We will not give these here, but will illustrate the results in Fig. 4, for . We see that the same quadratures are squeezed as below threshold, but that the degree of squeezing has been reduced.
The next question we raise is whether any of the possible bipartitions will exhibit the Einstein-Podolsky-Rosen paradox EPR , now commonly known as EPR steering Erwin ; Jonesteer . In the continuous variable case, this is usually measured using the Reid inequalities for the inferred variances EPRMDR ; ZYOu . This is written for the output spectral variances as
(19) |
where
(20) |
In the language of EPR-steering, shows that mode can be steered by measurements of mode . In some cases asymmetric steering is possible, where while . The question as to whether this was possible was first raised by Wiseman et al. Wiseman , and answered in the affirmative for Gaussian measurements by Olsen and Bradley SFG , Midgley et al. sapatona , and Händchen et al. Handchen . It has since been shown that asymmetric steering is generally possible Bowles , without any restriction on measurements. Because EPR steerable states are a strict subset of the entangled states, both symmetric and asymmetric steering demonstrate that the two modes concerned are fully bipartite entangled. We will therefore use the Reid inequalities to demonstrate both EPR steering and bipartite entanglement.
We obtain the below threshold covariances between each pair of modes as
(21) |
Since the covariances between modes 1 and 2 and modes 1 and 3 are zero, we can easily find four of the possible EPR correlations as
(22) |
It is obvious that none of these bipartitions can exhibit EPR steering below threshold, due to to Heisenberg Uncertainty Principal. An interesting result is that, although and are products of variances for different modes, they have equal values, with neither falling below one. This is not the case above threshold, where these two are no longer equal.
The case for modes 2 and 3, however, is different. A complicated analytical expression tells us that , so that any EPR steering here is completely symmetric. The result for the same parameters as in Fig. 3 is shown in Fig. 5. We see that the Reid inequalities are violated over a range near zero frequency, meaning that modes 2 and 3 are genuinely bipartite entangled.
Above threshold, the analytical expressions for all bipartitions become extremely complicated, and are best represented graphically. We will begin with and all cavity loss rates being equal, showing the effects of varying these later in the article. We find that modes 1 and 2 exhibit symmetric EPR steering over a broad range, while 1 and 3 exhibit completely asymmetric EPR steering over a narrower range of frequencies. The two higher frequency modes, which exhibit EPR steering below threshold, lose this property completely as the solution for moves onto the stable branch where it has non-zero amplitude. In terms of entanglement and EPR steering properties, the system changes completely at threshold.
We find that the symmetry or asymmetry of the EPR steering between the output modes above threshold can be simply controlled by the ratio of loss rates and the ratio of nonlinearities. Firstly, in Fig. 7, we show the results of a loss rate for the middle frequency which is one tenth of that for the other two, i.e. . Whereas modes 1 and 2 exhibited symmetric steering for equal loss rates, their steering is now asymmetric. The opposite has happened with modes 1 and 3, with their steering now being symmetric. The symmetry properties of the EPR steering can be controlled by adjusting the cavity loss rates, as was also found with intracavity second harmonic generation SHGEPR .
Changing the ratio also has an effect on the EPR steering properties above threshold. We can see in Fig. 8 that this can result in asymmetric steering in the bipartition of modes 1 and 2, with this swapping over at a certain frequency. Below , mode 1 can steer mode 2, while above this frequency there is a small violation of the inequality by . The pairing of 1 and 3 exhibits both symmetric and asymmetric EPR steering as the measurement frequency changes. We did not find any any steering involving the pair of fields at and , for the whole parameter range investigated with this ratio of the nonlinearities.
Vi Tripartite correlations
There are several methods of detecting tripartite inseparability and entanglement, with one common technique being based on inequalities developed by van Loock and Furusawa (vLF) vLF . These have proven useful for other cascaded systems AxMuzzJPB ; AxMuzz . The spectral inequalities we will use here are the set
(23) |
the violation of any one of which is sufficient to prove bipartite inseparability. Following the work of Teh and Reid Teh&Reid , any one of these less than demonstrates genuine tripartite entanglement, while one of them less than demonstrates genuine tripartite EPR steering. We did not find a violation of these inequalities below threshold. Above threshold we found that some, but not all, of the set of inequalities are violated for particular parameter regimes, as shown in Fig. 9, where we have divided the values of by four so as to be directly comparable with the tripartite EPR steering inequality to be described below. This value of demonstrates tripartite inseparability for the system.
With our three mode system, investigating tripartite EPR-steering is also of interest. It has been shown by Wang et al. Wang that, in a multipartite system, the steering of a given quantum mode is allowed when not less than half of the total number of modes take part in the steering group. In a tripartite system, this means that measurements on two of the modes are needed to steer the third. In order to quantify this, we will use the correlation functions developed by Olsen, Bradley, and Reid OBR . With spectral tripartite inferred variances defined as
(24) |
we define
(25) |
so that a value of less than one means that there is an inferred violation of the Heisenberg uncertainty principal and mode can be steered by the combined forces of modes and . According to the work of He and Reid HeReid , genuine tripartite steering is demonstrated whenever
(26) |
We did not find genuine tripartite steering for this system. As shown in Fig. 9, we found that modes 1 and 2 could combine for some parameters to steer mode 3. We investigated a wide parameter regime numerically, but did not find any for which more than one of the modes could be steered by the remaining pair simultaneously.
Vii An injected signal at the lower frequency
It is also possible to pump one of the cavity modes other than that at . The process of optical parametric downconversion with an injected signal has been experimentally and theoretically studied in some depth Bjorkholm ; Haub ; Hovde ; Plusquellic , with the injected signal often used for frequency stabilisation. An injected signal has also been shown to have a strong effect on any quantum correlations kaled , both changing the quadratures where squeezing is found and allowing for control of the asymmetry of EPR steering signal . For these reasons, we will examine here the effects of injecting a coherent signal at . Theoretically, this involves another term in the pumping Hamiltonian, so that
(27) |
where is the pumping Hamiltonian with injected signal. This change means that the equations of motion for and will have and added to them.
The immediate effect of an injected signal is to change the threshold properties of the system, with the low frequency mode developing a steady-state non-zero amplitude for all finite values of . There is now no critical pump value for , with the solutions remaining on the stable branch for all pumping values. The injected signal has an even more dramatic effect on the EPR steering properties of the system. As seen above in Fig. 5 the only two modes exhibiting EPR steering below threshold without injected signal were modes 2 and 3. With injected signal, the EPR steering of this bipartition soon vanishes as the signal is increased, which can be seen on the left hand side of Fig. 10, which shows the results for steerable bipartitions as the amplitude of the injected signal is increased. The quantities plotted are the minimum values of the Reid EPR correlations across all frequencies, ( numerically), so that a value of one means that the values near the carrier frequency can actually be larger than one. The addition of even a small injected signal (by comparison with ) has a dramatic effect on the and bipartitions, These become highly steerable for small injection and then less so as is increased. While exhibits symmetric steering, is totally asymmetric for these parameters, with across the whole range shown. The steerability of disappears on the same sort of scale of injection with which the others increase.
The spectral values of the Reid EPR correlations for the bipartitions which exhibit steering for similar parameters as in Fig. 10, but at a fixed , are shown in Fig. 11. The asymmetry of the EPR steering demonstrated by and is clearly shown. Nevertheless, this result shows that modes 1 and 3 are entangled across two octaves of frequency difference, and that this system is therefore a potentially important resource for any quantum processes linking resources over a large bandwidth. The injection of the coherent signal allows for a simple means of control over the entanglement properties of the system.
Viii Conclusion
In conclusion, the proposed system is a good candidate for novel quantum technologies which need squeezed and entangled optical states spanning a wide range of frequencies. With a single cavity input field it produces three output fields which are quadrature squeezed and different pairs of modes which are EPR steerable, with selection of the desired pairs being possible either by increasing the pump power or by injected signal. The quantum correlations of interest change depending on whether the system is being operated above or below the oscillation threshold, with good EPR steering being available in both regimes. The tripartite entanglement inequalities are only violated above threshold, where the lowest frequency mode develops a non-zero mean amplitude. An injected signal at the lowest frequency removes the threshold altogether and can provide either symmetric EPR steering across one octave or asymmetric EPR steering across two octaves. The flexibility and easy controllability of this system make it an attractive candidate for experimental investigation and future technological use.
Acknowledgements.
J.Y. Li was supported by the National Natural Science Foundation of China (Grant No. 11704287).References
- (1) J.A. Armstrong, N. Bloembergen, J. Ducuing, and P.S. Pershan, Phys. Rev. 127, 1918 (1962).
- (2) J.A. Giordmaine and R.C. Miller, Phys. Rev. Lett. 14, 973 (1965).
- (3) K. Nassau and H.J. Levinstein, Appl. Phys. Lett. 7, 69 (1965).
- (4) A. Yariv and W.H. Louisell, IEEE J. Quantum Elect. QE2, 418 (1966).
- (5) M.D. Reid, P.D. Drummond, W.P. Bowen, E.G. Cavalcanti, P.K. Lam, H.A. Bachor, U.L. Andersen, and G. Leuchs, Rev. Mod. Phys. 81, 1727 (2009).
- (6) S.F. Pereira, Min Xiao, H.J. Kimble, and J.L. Hall, Phys. Rev. A38, 4931 (1988).
- (7) M.K. Olsen and R.J. Horowicz, Opt. Commun. 168, 135 (1999).
- (8) N.B. Grosse, S. Assad, M. Mehmet, R. Schnabel, T. Symul, and P.K. Lam, Phys. Rev. Lett. 100, 243601 (2008).
- (9) M.K. Olsen, arXiv:1707.02537.
- (10) M.K. Olsen, Opt. Commun. in press, https://doi.org/10.1016/j.optcom.2017.09.090
- (11) C. Baune, J. Gniesmer, S. Kocsis, C.E. Vollmer, P. Zell, J. Fiurasek, and R. Schnabel, Phys. Rev. A93, 010302 (2016).
- (12) K. Hammerer, A.S. Sørensen, and E.S. Polzik, Rev. Mod. Phys. 82, 1041 (2010).
- (13) B. Julsgaard, J. Sherson, I. Cirac, L. Fiurás̆ek, and E.S. Polzik, Nature 432, 482 (2004).
- (14) P.D. Drummond and C.W. Gardiner, J. Phys. A 13, 2353 (1980).
- (15) C.W. Gardiner and M.J. Collett, Phys. Rev. A31, 3761 (1985).
- (16) S. Zhu, Y. Zhu, and N. Ming, Science 278, 843 (1997).
- (17) E. Marcellina, J .F.Corney, and M. K. Olsen Opt. Commun. 309, 9 (2013).
- (18) M.K. Olsen, S.C.G. Granja, and R.J. Horowicz, Opt. Commun. 165, 293 (1999).
- (19) D.F. Walls and G.J. Milburn, Quantum Optics (Springer, Berlin, 1995).
- (20) C.W. Gardiner, Handbook of Stochastic Methods, (Springer, Berlin, 1985).
- (21) P.D. Drummond, K.J. McNeil, and D.F. Walls, Optica Acta 28, 211 (1981).
- (22) S. Chaturvedi, K. Dechoum, and P.D. Drummond, Phys. Rev. A65, 033805 (2002).
- (23) K.J. McNeil, P.D. Drummond, and D.F. Walls, Opt. Commun. 27, 292 (1978).
- (24) M. Bache, P. Lodahl, A.V. Mamaev, M. Marcus, and M. Saffman, Phys. Rev. A65, 033811 (2002).
- (25) M.K. Olsen, arXiv:1706.05174.
- (26) A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. 47, 777 (1937).
- (27) E. Schrödinger, Proc. Cam. Philos. Soc. 31, 555 (1935).
- (28) S.J. Jones, H.M. Wiseman, and A.C. Doherty, Phys. Rev. A76, 052116 (2007).
- (29) M.D. Reid, Phys. Rev. A40, 913 (1989).
- (30) Z.Y. Ou, S.F. Pereira, H.J. Kimble, and K.C. Peng, Phys. Rev. Lett. 68, 3663 (1992).
- (31) H.M. Wiseman, S.J. Jones, and A.C. Doherty, Phys. Rev. Lett. 98, 140402 (2007).
- (32) M.K. Olsen and A.S. Bradley, Phys. Rev. A77, 023813 (2008).
- (33) S.L.W Midgley, A.J. Ferris, and M.K. Olsen, Phys. Rev. A81, 022101 (2010).
- (34) V. Händchen, T. Eberle, S. Steinlechner, A. Samblowski, T. Franz, R.F. Werner and R. Schnabel, Nat. Photon. 6, 596 (2012).
- (35) J. Bowles, T. Vértesi, M.T. Quintino, and N. Brunner, Phys. Rev. Lett. 112, 200402 (2014).
- (36) M.K. Olsen, Phys. Rev. A88, 051802 (2013).
- (37) P. van Loock and A. Furusawa, Phys. Rev. A67, (2003) 052315.
- (38) M.K. Olsen and A.S. Bradley, J. Phys. B 39, 127 (2005).
- (39) M.K. Olsen and A.S. Bradley, Phys. Rev. A74 063809 (2006).
- (40) R.Y. Teh and M.D. Reid, Phys. Rev. A90, 062337 (2014).
- (41) M. Wang, Y. Xiang, Q. He, Q. Gong, Phys. Rev. A91, 012112 (2015).
- (42) M.K. Olsen, A.S. Bradley, M.D. Reid, J. Phys. B 39, 2515 (2006).
- (43) Q.Y. He and M.D. Reid, Phys. Rev. Lett. 111, 250403 (2013).
- (44) J.E. Bjorkholm and H.G. Danielmeyer, Appl. Phys. Lett. 15, 171 (1969).
- (45) J.G. Haub, M.J. Johnson, B.J. Orr, and R. Wallenstein, Appl. Phys. Lett. 58, 1718 (1991).
- (46) D.C. Hovde, J.H. Timmermans, G. Scoles, and K.K. Lehmann, Opt. Commun. 86, 294 (1991).
- (47) D.F. Plusquellic, O. Votava, and D.J. Nesbitt, Appl. Optics 35, 1464 (1996).
- (48) M.K. Olsen, K. Dechoum, and L.I. Plimak, Opt. Commun. 223, 123 (2003).
- (49) M.K. Olsen, Phys. Rev. Lett. 119, 160501 (2017).