Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://qi.phys.msu.ru/papers/2013-njp-15-073016.pdf Äàòà èçìåíåíèÿ: Wed Oct 2 18:28:57 2013 Äàòà èíäåêñèðîâàíèÿ: Thu Feb 27 20:33:50 2014 Êîäèðîâêà:

Home

Search

Collections

Journals

About

Contact us

My IOPscience

Schmidt-like coherent mode decomposition and spatial intensity correlations of thermal light

This content has been downloaded from IOPscience. Please scroll down to see the full text. 2013 New J. Phys. 15 073016 (http://iopscience.iop.org/1367-2630/15/7/073016) View the table of contents for this issue, or go to the journal homepage for more

Download details: IP Address: 93.180.54.115 This content was downloaded on 02/10/2013 at 15:28

Please note that terms and conditions apply.


Schmidt-like coherent mode decomposition and spatial intensity correlations of thermal light
I B Bobrov, S S Straupe1 , E V Kovlakov and S P Kulik Faculty of Physics, M V Lomonosov Moscow State University, Moscow, Russia E-mail: straups@yandex.ru
New Journal of Physics 15 (2013) 073016 (10pp)

Received 16 January 2013 Published 5 July 2013 Online at http://www.njp.org/
doi:10.1088/1367-2630/15/7/073016

Abstract. We experimentally study the properties of coherent mode decomposition for the intensity correlation function of quasi-thermal light. We use the technique of spatial mode selection developed for studying the transverse entanglement of photon pairs, and show that it can be extended to characterize classical spatial correlations. Our results demonstrate the existence of a unique, for a given thermal source, basis of coherent modes, correlated in a way much resembling the Schmidt modes of spatially entangled photons. Contents

1. Introduction 2. Coherent mode decomposition for quasi-thermal light 3. Experimental realization 4. Discussion and conclusion Acknowledgments References
1. Introduction

1 2 5 8 9 10

Entanglement in spatial degrees of freedom of quantum light, such as pairs of photons generated in spontaneous parametric down-conversion (SPDC), is currently an object of active research.
1

Author to whom any correspondence should be addressed.

Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any fur ther distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.

New Journal of Physics 15 (2013) 073016 1367-2630/13/073016+10$33.00

© IOP Publishing Ltd and Deutsche Physikalische Gesellschaft


2 These studies are motivated by applications in quantum information science, where the infinitedimensional Hilbert space of the spatial states of photons offers attractive capabilities for high-dimensional quantum state engineering. The effective dimensionality of accessible state space is closely related to the degree of spatial entanglement for bipartite states, making its experimental quantification an important task. One of the most successful approaches makes use of a Schmidt decomposition. As stated by Law and Eberly [1] in their seminal work, a unique basis of coherent modes {|u k , |vk } exists, such that a spatial state of biphoton pairs | 12 = dx1 dx2 (x1 , x2 )a (x1 )a (x2 ) |vac can be decomposed as | = k k |u k |vk , with these so-called Schmidt modes being eigenvectors of reduced single-particle density matrices and k the corresponding eigenvalues. A remarkable feature of this decomposition is its single-sum form, implying perfect one-to-one correlations between Schmidt modes. These correlations have been studied in several recent works [2­4], and the number of significant terms in the decomposition is routinely used as an entanglement quantifier [1, 5­12]. A natural question is can a similar approach be used to study spatial correlations of other origin, for example correlations in spatially multi-mode classical light? The form of Schmidt decomposition, and the way in which explicit expressions for Schmidt modes are derived in the case of SPDC biphotons [1, 13], resembles the well-known notion of coherent mode decomposition for spatial correlation functions in classical statistical optics [14­16]. Here we show that this similarity does not reside only in the form of mathematical equations, but in underlying physics as well. Namely, for a given source of light with quasi-thermal statistics there exists a unique physically distinguished set of coherent modes, correlated in much the same way as Schmidt modes of SPDC. Similar experimental settings would provide similar correlations, with only a quantitative difference in visibility. The situation here resembles that with ghost imaging, which was for some time believed to be possible only due to entanglement [17] and thus only with quantum light. It is, however, well understood now that spatial correlations required to obtain ghost images should not necessarily be of quantum origin, and a purely classical thermal light may be used as well [18, 19]. In this work we try to push this analogy further and show that a framework of Schmidt decomposition may also be generalized to classical states of light. This paper is organized as follows. In section 2 we review some properties of coherent mode decomposition for a Gaussian Schell-model of quasi-thermal light and derive the expressions for measured correlation functions, experimental realization is described in section 3, and we end with a discussion of the results and the conclusion in section 4.
2. Coherent mode decomposition for quasi-thermal light

The coherence properties of classical light are described by the coherence function G (1) (r1 , r2 , ) = E (r1 , ) E (r2 , ) , where we neglect the polarization and use scalar complex amplitudes E (r, ) to describe the electromagnetic field. In the following, we will consider the monochromatic case and omit the frequency dependence for simplicity. One of the results of classical coherence theory states that it can be decomposed in series of coherent mode functions as G (1) (r1 , r2 ) =
n n n (r1 )n (r2 ),

(1)

New Journal of Physics 15 (2013) 073016 (http://www.njp.org/)


3 where n (r) are eigenfunctions of an integral operator with a kernel G (1) (r1 , r2 ) and n are corresponding eigenvalues [16]. Higher-order correlation functions, such as the intensity correlation function G (2) (r1 , r2 ) = I (r1 ) I (r2 ) , may be accessed via Hanbury Brown­Twiss (HBT) interferometry, which is exactly the experimental setting commonly used for studying the spatial entanglement of photon pairs. In the biphoton case, the fourth-order correlation function G (2) (r1 , r2 ) is measured by counting coincidences between photocounts of two single-photon detectors positioned at points r1 and r2 , respectively: Rc a (r1 )a (r2 )a (r1 )a (r2 ) = | (r1 , r2 )|2 with a (r) being the photon annihilation operator and (r1 , r2 ) the biphoton amplitude. A direct analogue of coherent mode decomposition (1) for the biphoton amplitude is exactly the Schmidt decomposition. To see whether this analogy goes beyond simply similar mathematical expressions, let us consider the behavior of a classical light field in the HBT scheme with spatial mode filters in the arms of the interferometer. Such filters are usually composed of phase holograms followed by singlemode fibers (SMFs) and bucket detectors as shown in figure 1. The main feature of Schmidt decomposition--its single-sum form--manifests itself in the absence of coincidences between photocounts in such a setup when appropriate mode filters are used to select orthogonal Schmidt modes in different arms of the setup [2]. We will show below that similar effects may be observed with classical quasi-thermal light. For quasi-thermal light, the intensity correlation function G (2) (r1 , r2 ) can be expressed in terms of the second-order correlation function G (1) (r1 , r2 ) using the Siegert relation: G (2) (r1 , r2 ) = I (r1 ) I (r2 ) + |G (1) (r1 , r2 )|2 . (2) If the spatial mode filters in the two arms of the HBT scheme are described by impulse response functions h 1 (r1 , r1 ) and h 2 (r2 , r2 ), where coordinate r corresponds to the source plane and coordinate r to the detection plane, the transformed intensity correlation function takes the form [18, 20] G (2) (r1 , r2 ) = I (r1 ) I (r2 )
2

+

dr1

dr2 h (r1 , r1 )h 2 (r2 , r2 ) G (r1 , r2 ) ,
1 2

(1)

(3)

where I (r1,2 ) = dr1,2 h 1,2 (r1,2 , r1,2 ) E (r1,2 ) is the intensity distribution in the detection plane. One can always choose the mode filters to satisfy the relation: dr1 h
(m ) 1

(r1 , r1 )n (r1 ) e

-

(r1 )2 2 wf (r2 )2 2 wf

n m , n k ,

(4)

dr1 h (k ) (r2 , r2 )n (r2 ) e 2

-

where wf is the waist of the fundamental Gaussian mode of the fiber. This choice of propagators corresponds to `projection' on modes m (r1 ) and k (r2 ); we will call the transformation corresponding to such propagators an eigenmode filter. Substituting (1) and (4) into (3) and taking into account that bucket detectors placed after the SMF integrate the signal over the whole detection plane, we obtain the following expression for coincidence to single counts ratio: ( Rcm ,k ) (2 g(m),k ) = 1 + m k , (5) (m ) (k ) Rs Rs
New Journal of Physics 15 (2013) 073016 (http://www.njp.org/)


4 where
(2 g(m),k ) =

dr1 dr2 G (2) (r1 , r2 ) dr1 I (r1 ) dr2 I (r2 )

(6)

is the normalized second-order correlation function. It is clear from (5) that spatial correlations between appropriately chosen coherent modes of thermal light show the same correlation features as Schmidt modes of entangled photons--namely, they are only pairwise correlated. The difference turns out to be rather quantitative than qualitative, similar to the situation (2 with thermal and quantum ghost imaging--the achievable visibility defined as V = (g(m),m ) - (2 (2 (2 g(m),n ) )/(g(m),m ) + g(m),n) ), m = n cannot exceed 1/3 for thermal light, while it can reach values arbitrarily close to unity for biphotons. A simple expression (5) is only valid for eigenmode filters; any other mode filters ( h 1m2,k ) (r1,2 , r1,2 ) = f m ,k (r1,2 ) f m ,k (r1,2 ) will give rise to a non-unity value of the normalized , (2 intensity correlation function for m = k . Indeed, the expression for g(m),k ) in the case of arbitrary filters has the following form:
(2 g(m),k )

= 1+

i, j

i j (i(m ) ) i(k ) (
i

i |i(m ) |2

j



(m ) (k ) j ) j (k ) 2 j | j |

,

(7)

where i(m ) = dr f m (r)i (r) are coefficients in expansion of mode functions f m (r) in the basis of coherent eigenmodes. This expression reduces to (5) only if i(m ) = m i , which is the case of eigenmode filters. This fact is a unique and, as we will show below, experimentally verifiable property of coherent eigenmodes for thermal light. Let us now derive the explicit expressions for coherent modes constituting decomposition (1) for a simple model of partially spatially coherent thermal light. Following [14] we describe light from the quasi-thermal source with a Gaussian Schell-model for the second-order coherence function: G (1) (x1 , x2 ) = I (x1 ) I (x2 )µ(x1 - x2 ),
2 µ(x1 - x2 ) = exp[-(x1 - x2 )2 /2µ ]

(8) (9)

where spectral intensity distribution and degree of spatial coherence are taken in the form
2 I (x1,2 ) = A exp[-x1,2 /2 I2 ],

with x1 and x2 being the transverse coordinates, A a positive normalization constant, I beam waist and µ its coherence radius. For the sake of simplicity, we will first consider a one-dimensional source; generalization to an experimentally relevant two-dimensional case is straightforward. In the case of partially spatial coherent light described by a Gaussian Schell-model, one can find the eigenfunctions in a closed analytical form [14, 16]: 2c n (x ) = n = A where a= 1 , 4 I2 b= 1 , 2 4µ c= a 2 + 2a b (12)
1/4



1 2n n !

Hn (x 2c) exp(-c x 2 ),
n

(10) (11)

a+b+c

1/2

b a+b+c

,

New Journal of Physics 15 (2013) 073016 (http://www.njp.org/)


5 and Hn (x ) are Hermite polynomials. Eigenvalues in this coherent mode decomposition are determined by the ratio of coherence radius to beam waist = µ : I
n 1 n (13) = . 0 2 /2 + 1 + [( /2)2 + 1]1/2 In this simple but experimentally relevant model, the eigenmodes turn out to be familiar Hermite­Gaussian (HG) modes, which were also found to be Schmidt modes of SPDC biphotons under a double-Gaussian approximation for the biphoton amplitude [2, 3, 13]. We note that similar decomposition may be carried out in polar coordinates giving rise to coherent modes with orbital angular momentum (OAM) [21­23]. The spectrum of OAM eigenvalues in this decomposition was measured experimentally in [24] with the interferometric technique.

3. Experimental realization

The experimental setup is shown in figure 1. We used a well-known technique to create a pseudothermal light source with controlled spatial coherence [25]. In this scheme, a beam of a He­Ne laser with wavelength = 632.8 nm is scattered on a slowly rotating grounded-glass disc. It is important to make the characteristic time of variation of the scattering pattern larger than the detection time det 2 ns, defined by the coincidence circuit time window. The disc was located at the focus of an effective lens consisting of two microscope objectives O1 and O2, forming a collimated beam. The spatial coherence properties of the beam are controlled by iris apertures S and S2 , determining the coherence radius µ and the beam waist I , respectively. The detection part of the setup is an HBT scheme, consisting of a 50/50 non-polarizing beam splitter; spatial mode filters--a reflective liquid-crystal-on-silicon (LCoS) spatial light modulator (SLM) in the transmitted channel and a glass step mask in the reflected one, followed by SMF1,2 placed in the focal planes of 8â microscope objectives C1,2 ; electrical signals from two singlephoton detectors D1,2 (Perkin­Elmer) were fed to a coincidence circuit. The model of SLM used was limited to 0.8 phase shifts only, forcing us to use double reflection as shown in figure 1. Polarizers P2,3 were used to specify the correct polarization required for operation of the polarization-sensitive SLM matrix and to remove the undiffracted components, while polarizer P1 was used to control the overall intensity. The propagator for each arm of the HBT scheme described above is (r )2 rr k h (r, r ) exp - 2 â dr exp -i f wf (r - rf )2 â exp[i arg( Hm ( 2c x ) Hn ( 2c y ))] â exp - , (14) wf2 where k = 2 / is the wave number, f the focal length of the microscope objectives C1,2 , wf the waist of the fundamental gaussian mode of the fiber and rf the fiber tip displacement in the focal plane. This propagator defines the detection modes, which for our purposes should be matched with the eigenmodes given by (10). The mode matching is accomplished by choosing the focal length f = the eigenmodes:
2 k wf 4c

; in this case the detection mode waist coincides with the waists of 2 . 2

k 2 wf2 1 =c= 2 4f 4

2 I

1+

New Journal of Physics 15 (2013) 073016 (http://www.njp.org/)


6

Figure 1. Experimental setup (see text for a detailed description).

The mode filters were tested with a coherent single-mode laser beam. The sizes and positions of the phase masks relative to the beam were determined experimentally following the operational criterion of minimizing the counting rate for higher-order modes in the central position of the SMF2 fiber. k 2 w2 Satisfying the condition 4 f 2f = c is essential for obtaining the desired matching of the detected modes and the eigenmodes of the coherence function. In the actual experiment the detection mode waist was fixed, while the mode matching was done by changing the parameters of the input beam. The widths of the apertures S and S2 determine the coherence radius µ and the beam waist I , respectivetly, and thus define the eigenmode width c through (12). One can (2 see from (7) that the value of the normalized intensity coherence function g(m),0) (rf = 0) = 1 for m = 0 only if the mode-matching condition is satisfied; otherwise it is always larger than 1. This fact was used as an operational criterion of good mode matching. The corresponding beam waist and coherence radius were measured as I = (2.3 ± 0.1) mm and µ = (0.57 ± 0.02) mm, corresponding to the value of = 0.24 ± 0.02. With the mode-matching conditions satisfied, the spatial mode filters in both arms of the setup act as projectors on the eigenmodes of coherent mode decomposition (1). In this case, the single count rate of a detector in one of the arms is proportional to the partial intensity of a constituent mode selected by the filter: I
(m )



dr

drh

(m ) 1,2

2

(r, r , ) E (r) m .

(15)

It provides a way of measuring the eigenvalues in (1) experimentally. The distribution of eigenvalues, normalized to 0 , is shown in figure 2. We find excellent agreement with theoretical predictions obtained from (13) using the independently measured value = 0.24 ± 0.02. The spatial shape of modes (10) selected by the corresponding mask in the transmitted channel is revealed in the dependence of intensity correlation function on fiber displacement in the reflected channel. Indeed, if a phase mask corresponding to the mth eigenmode is installed
New Journal of Physics 15 (2013) 073016 (http://www.njp.org/)


7

Figure 2. Experimentally measured eigenvalues of coherent mode decom-

position. (a) Two-dimensional distribution of normalized eigenvalues m n /0 (indexes m and n correspond to HGm 0 and HG0n modes, respectively). (b) Onedimensional section (n = 0) of the distribution. Blue bars are experimental data and black circles are values calculated from (13) using the independently measured value = 0.24 ± 0.02. Error bars correspond to error in the measured value of and error bars for experimental points are smaller than point markers. in the transmitted channel with no mask in the reflected channel, the general expression (7) for a normalized intensity correlation function reduces to the following form:
(2 g(m),0) (rf ) = 1 +

m |m (rf )|2 , k |k (rf )|2
k

(16)

where rf is the position of the SMF in the reflected channel and k (rf ) = drk (r)0 (r - rf ). Experimental data for HGm 0 modes are shown in figure 3(b). It is clear that the selected (2 modes are indeed orthogonal in the sense that g(m),0) (0) = 1, m = 0, consistently with (5). Dependence on the fiber displacement is also in qualitative agreement with (16). However, for (2 large displacements we observe anomalously high values of g(m),0) (rf ) for large m . This is most probably an artifact of our method for mode selection based on using phase-only step masks, inevitably causing higher-order mode contributions and compromising mode purity for large displacements. Figure 3(a) shows the results of numerical simulations for the imperfect masks, showing reasonable agreement with experimental data. For zero displacement and low-order modes, we find the detected modes to be close to Hermite­Gaussian eigenmodes, and thus they are pairwise correlated, as illustrated by figure 4. (2 We expect g(m),n ) (rf = 0) to behave in accordance with expression (5), i.e. to be essentially unity for modes of different order and larger than unity only for modes of the same order. Experimental results of figure 4 confirm this statement, although some amount of anomalous correlation is observed for modes of high order. These correlations are more pronounced in figure 4(b), which may be partially explained by the much poorer quality of the glass phase masks used in the reflected channel. Nevertheless, one can clearly observe the good agreement with the expected correlation properties of quasi-Schmidt coherent modes. The main source of discrepancies between our experimental results and theoretical expectations for ideal projective measurements in a Schmidt-like basis is the use of non-ideal mode filters. Technical limitations of the SLM used do not allow us to use it for amplitude
New Journal of Physics 15 (2013) 073016 (http://www.njp.org/)


8

Figure 3. Spatial structure of eigenmodes: dependence of normalized g(

function on fiber displacement in the reflected channel. (a) Numerical simulation and (b) experimental results. Unity values for zero displacement correspond to the absence of correlations for modes of different order, predicted by theory.

(2) m ,0)

Figure 4. Mode cross-correlations: (a) g(

- 1 values measured with HGm 0 modes selected in the transmitted arm and the HG00 mode in the reflected one; (2 (b) g(m),1) - 1 values, corresponding to HGm 0 modes filtered in the transmitted arm, and the HG01 mode selected in the reflected arm.

(2) m ,0)

modulation, which is necessary for perfect transformations of Hermite­Gaussian modes. Nevertheless, our results are very similar to what is expected for exact Hermite­Gaussian filters.
4. Discussion and conclusion

In conclusion, we have investigated spatial intensity correlations in partially spatially coherent thermal light using the technique of projective measurements originally developed to study the spatial entanglement of photon pairs. In particular, our detection scheme is similar to those used in experiments on orbital angular momentum entanglement [27] and essentially the same as was used by authors to study approximate Schmidt decomposition for SPDC biphotons [2]. We have found much similarity in the properties of coherent mode decomposition for the intensity correlation function of classical thermal light and Schmidt decomposition for biphoton
New Journal of Physics 15 (2013) 073016 (http://www.njp.org/)


9 amplitude. Although the expressions for eigenmodes and eigenvalues are different, they are still governed by a single experimental parameter = µ / I --the ratio of coherence radius to the width of intensity distribution. A similar quantity for entangled photons is known as the Fedorov ratio R = xc / xs --the ratio of widths of coincidence (conditional) and single count (single particle) distributions in either the coordinate or momentum space of photons [28]. This quantity is shown to be equal to Schmidt number K , which is the effective number of nonzero eigenvalues i in the Schmidt decomposition, defined as K = i i-2 . Thus for a pure two-photon state it can serve as a measure of spatial entanglement. It is clear that this quantity has a purely classical analogue--optical etendue [6]. Results of this work show that the analogy between spatial entanglement and classical spatial correlations goes further. In purely operational terms the qualitative similarity is almost complete, with the difference being only quantitative--excess of g (2) over unity is limited to 1, while it can be unlimited in the case of entangled photons. One can consider the described experiments as a special case of thermal light ghost interference/imaging with mode-sensitive detection. In this context, our work adds new arguments to the discussion of quantum and classical features of spatial correlations underlying this method. Coherent mode decomposition for the intensity correlation function inherits most of the properties of its `quantum' counterpart. Most importantly, it defines a naturally preferred set of modes that are both -correlated and have maximal partial intensities. Recently, an essentially similar concept of `optical eigenmodes' was introduced by De Luca et al [29] in the context of compressive imaging. These authors used coherent mode decomposition to reconstruct the optical field passing through a complex amplitude and/or phase mask, and argued that coherent eigenmodes are the best choice of modes to fulfill this task. Our work demonstrates that in the case of spatially incoherent light with strong intensity correlations (such as quasi-thermal light), eigenmodes also exhibit stronger correlations than any other coherent spatial mode. This means that they are to be used for a `ghost' generalization of eigenmode imaging. Our experiment is a proof-of-principle one, and treats only the simplest case of a freely propagating Gaussian Schell-model beam; that is why eigenmodes are Hermite­Gaussian in our case. For a more complex field distribution, the eigenmodes may be different but the general framework stays the same, making our experiments an important step toward compressive ghost imaging with thermal light. From a more fundamental point of view, the spatial mode-filtering technique used in this type of experiment may serve as an additional tool to study the differences between classical and non-classical spatial correlations. The differences here may go further than obvious reduction of visibility. The statistics of mode-filtered light is determined by both the initial photon statistics and the spatial structure of detection modes relative to the (quasi)-Schmidt eigenmodes of the input field. The presented technique may contribute to a deeper understanding of what is classical and what is really quantum in spatial correlations of multi-mode light.
Acknowledgments

We are grateful to A A Kalinkin for considerable help at an early stage of the experiment. This work was funded by the European Union Seventh Framework Programme under grant agreement no. 308803 (project BRISQ2), the European Sixth Framework Programme under grant agreement no. 175 (project NANOQUINT), RFBR grant 12-02-31041 and the Federal Program of the Russian Ministry of Education and Science (grants 8393 and 11.519.11.4009). SSS and IBB are grateful to the `Dynasty' foundation for financial support.
New Journal of Physics 15 (2013) 073016 (http://www.njp.org/)


10
References
[1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] Law C K and Eberly J H 2004 Phys. Rev. Lett. 92 127903 Straupe S et al 2011 Phys. Rev. A 83 060302 Miatto F M, Brougham T and Yao A M 2012 Eur. Phys. J. D 66 183 Miatto F M et al 2012 Eur. Phys. J. D 66 263 Torres J P, Alexandrescu A and Torner L 2003 Phys. Rev. A 68 050301 van Exter M P et al 2006 Phys. Rev. A 74 012309 Fedorov M V et al 2007 Phys. Rev. Lett. 99 063901 Peeters W H, Verstegen E J K and van Exter M P 2007 Phys. Rev. A 76 042302 Pors J B et al 2008 Phys. Rev. Lett. 101 120502 Miatto F M, Yao A M and Barnett S M 2011 Phys. Rev. A 83 033816 Miatto F M et al 2012 Eur. Phys. J. D 66 178 ¨ Salakhutdinov V D, Eliel E R and Loffler W 2012 Phys. Rev. Lett. 108 173604 Fedorov M V, Mikhailova Yu M and Volkov P A 2009 J. Phys. B: At. Mol. Opt. Phys. 42 175503 Starikov A and Wolf E 1982 J. Opt. Soc. Am. 72 923 Agarwal G S and Wolf E J. Mod. Opt. 40 1489 Mandel L and Wolf E 1995 Optical Coherence and Quantum Optics (Cambridge: Cambridge University Press) Pittman T, Shih Y, Strekalov D and Sergienko A 1995 Phys. Rev. A 52 3429 Gatti A, Brambilla E, Bache M and Lugiato L A 2004 Phys. Rev. Lett. 93 093602 Erkmen I and Shapiro H 2008 Unified theory of ghost imaging with Gaussian-state light Phys. Rev. A 77 043809 Bache M et al 2006 Phys. Rev. A 73 053802 Simon R and Mukunda N 1993 J. Opt. Soc. Am. A 10 95 Serna J and Movilla J M 2001 Opt. Lett. 26 405 Liu Y-D et al 2008 Opt. Commun. 281 1968 Di Lorenzo, Pires H, Woudenberg J and van Exter M P 2010 Opt. Lett. 35 889 Martienssen W and Spiller K 1964 Am. J. Phys. 32 919 Klyshko D N 1988 Two-photon light: influence of filtration and a new possible EPR experiment Phys. Lett. A 128 133­7 Molina-Terriza G, Torres J P and Torner L 2007 Nature Phys. 3 305 Fedorov M V et al 2005 Phys. Rev. A 72 032110 De Luca A C, Kosmeier S, Dholakia K and Mazilu M 2011 Phys. Rev. A 84 021803

New Journal of Physics 15 (2013) 073016 (http://www.njp.org/)