Ghost Imaging with Twisted Gaussian Schell ... - Semantic Scholar

University of Miami

Scholarly Repository Physics Articles and Papers

Physics

2-1-2009

Ghost Imaging with Twisted Gaussian SchellModel Beam Yangjian Cai Qiang Lin Olga Korotkova University of Miami, [email protected]

Recommended Citation Cai, Yangjian; Lin, Qiang; and Korotkova, Olga, "Ghost Imaging with Twisted Gaussian Schell-Model Beam" (2009). Physics Articles and Papers. Paper 10. http://scholarlyrepository.miami.edu/physics_articles/10

This Article is brought to you for free and open access by the Physics at Scholarly Repository. It has been accepted for inclusion in Physics Articles and Papers by an authorized administrator of Scholarly Repository. For more information, please contact [email protected].

Ghost imaging with twisted Gaussian Schellmodel beam Yangjian Cai1 , Qiang Lin2 and Olga Korotkova3 1. School of Physical Science and Technology, Suzhou University, Suzhou 215006, China 2.Institute of Optics, Department of Physics, Zhejiang University, Hangzhou 310027, China 3. Department of Physics, University of Miami, Coral Gables, Florida 33146, USA *Corresponding author: [email protected]

Abstract: Based on the classical optical coherence theory, ghost imaging with twisted Gaussian Schell-model (GSM) beams is analyzed. It is found that the twist phase of the GSM beam has strong influence on ghost imaging. As the absolute value of the twist factor increases, the ghost image disappears gradually, but its visibility increases. This phenomenon is caused by the fact that the twist phase enhances the transverse spatial coherence of the twisted GSM beam on propagation. ”2009 Optical Society of America OCIS codes: (030.1640) Coherence; (030.1670) Coherent optical effects; (110.2990) Image formation theory

References and links 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14. 15. 16. 17. 18. 19.

L. Mandel and E. Wolf, Optical Coherence and Quantum Optics (Cambridge U. Press,1995). J. T. Foley and M. S. Zubairy, “The directionality of Gaussian Schell-model beams,” Opt. Commun. 26, 297-300 (1976). E. Collett and E. Wolf, "Is complete spatial coherence necessary for the generation of highly directional light beams?," Opt. Lett. 2, 27-29 (1978) E. Wolf and E. Collett, “Partially coherent sources which produce same far-field intensity distribution as a laser,” Opt. Commun. 25, 293-296 (1978). F. Gori, “Collet-Wolf sources and multimode lasers,” Opt. Commun. 34, 301-305 (1978). P. de Santis, F. Gori, G. Guattari, and C. Palma. “An example of a Collett-Wolf source,” Opt. Commun. 29, 256-260 (1979). F. Wang and Y. Cai, “Experimental observation of fractional Fourier transform for a partially coherent optical beam with Gaussian statistics,” J. Opt. Soc. Am. A 24, 1937-1944 (2007). A. T. Friberg, and R. J. Sudol. Propagation parameters of Gaussian Schell-model beams. Opt. Commun 41, 383-387 (1982). E. Tervonen, A. T. Friberg, and J. Turunen. Gaussian Schell-model beams generated with synthetic acousto-optic holograms. J. Opt. Soc. Am. A. 9, 796-803 (1992). R. Simon and N. Mukunda, “Twisted Gaussian Schell-model beams,” J. Opt. Soc. Am. A 10, 95-109 (1993). A. T. Friberg, E. Tervonen, and J. Turunen, “Interpretation and experimental demonstration of twisted Gaussian Schell-model beams,” J. Opt. Soc. Am. A 11, 1818-1826 (1994). R. Simon and N. Mukunda, “Twist phase in Gaussian-beam optics,” J. Opt. Soc. Am. A 15, 2373-2382 (1998). D. Ambrosini, V. Bagini, F. Gori, and M. Santarsiero, “Twisted Gaussian Schell-model beams: a superposition model,” J. Mod. Opt. 41, 1391-1399 (1994). R. Simon, K. Sundar, and N. Mukunda, “Twisted Gaussian Schell-model beams. I. Symmetry structure and normal-mode spectrum,” J. Opt. Soc. Am. A 10, 2008-2016 (1993). K. Sundar, R. Simon, and N. Mukunda, “Twisted Gaussian Schell-model beams. II. Spectrum analysis and propagation characteristics,” J. Opt. Soc. Am. A 10, 2017-2023 (1993). R. Simon, A. T. Friberg, and E. Wolf, “Transfer of radiance by twisted Gaussian Schell-model beams in paraxial systems,” Pure Appl. Opt. 5, 331-343 (1996). M. J. Bastiaans. Wigner-distribution function applied to twisted Gaussian light propagating in first-order optical systems. J. Opt. Soc. Am. A 17, 2475-2480 (2000). Q. Lin and Y. Cai, “Tensor ABCD law for partially coherent twisted anisotropic Gaussian-Schell model beams,” Opt. Lett. 27, 216-218 (2002). R. Simon and N. Mukunda, “Gaussian Schell-model beams and general shape invariance” J. Opt. Soc. Am. A 16, 2465-2475 (1999).

#102691 - $15.00 USD

(C) 2009 OSA

Received 13 Oct 2008; revised 29 Nov 2008; accepted 2 Jan 2009; published 5 Feb 2009

16 February 2009 / Vol. 17, No. 4 / OPTICS EXPRESS 2453

20. Q. Lin and Y. Cai, “Fractional Fourier transform for partially coherent Gaussian-Schell model beams,” Opt. Lett. 27, 1672-1674 (2002). 21. S. A. Ponomarenko, “Twisted Gaussian Schell-mode solitons,” Phys. Rev. E 64, 036618 (2001). 22. Y. Cai and L. Hu, “Propagation of partially coherent twisted anisotropic Gaussian Schell-model beams through an apertured astigmatic optical system,” Opt. Lett. 31, 685-687 (2006). 23. Y. Cai and S. He, “Propagation of a partially coherent twisted anisotropic Gaussian Schell-model beam in a turbulent atmosphere,” Appl. Phys. Lett. 89, 041117 (2006). 24. J. Serna and J. M. Movilla, “Orbital angular momentum of partially coherent beams,” Opt. Lett. 26, 405406 (2001). 25. Y. Cai and U. Peschel, “Second-harmonic generation by an astigmatic partially coherent beam,” Opt. Express 15, 15480-15492 (2007). 26. T. B. Pittman, Y. H. Shih, D. V. Strekalov, and A. V. Sergienko, “Optical imaging by means of two-photon quantum entanglement,” Phys. Rev. A 52, R3429 (1995). 27. G. A. Barbosa, “Quantum images in double-slit experiments with spontaneous down-conversion light,” Phys. Rev. A 54, 4473-4478 (1996). 28. A. Gatti, E. Brambilla, L. A. Lugiato and M. I. Kolobov, “Quantum entangled images,” Phys. Rev. Lett. 83, 1763-1766 (1999). 29. B. E. A. Saleh, A. F. Abouraddy, A. V. Sergienko, and M. C. Teich, “Duality between partial coherence and partial entanglement,” Phys. Rev. A 62, 043816 (2000). 30. A. F. Abouraddy, B. E. A. Saleh, A.V. Sergienko, M. C. Teich, “Quantum holography,” Opt Express 9, 498-505 (2001). 31. A. F. Abouraddy, B. E. A. Saleh, A. V. Sergienko, and M. C. Teich, “Entangled-photon Fourier optics,” J. Opt. Soc. Am. B 19, 1174-1184 (2002). 32. A. Gatti, E. Brambilla, and L. A. Lugiato, “Entangled imaging and wave-particle duality: from the microscopic to the macroscopic realm,” Phys. Rev. Lett. 90, 133603 (2003). 33. A. F. Abouraddy, B. E. A. Saleh, A. V. Sergienko, and M. C. Teich, “Role of entanglement in two-photon imaging,” Phys. Rev. Lett. 87, 123602 (2001). 34. R. S. Bennink, S. J. Bentley, and R. W. Boyd, “Two-photon coincidence imaging with a classical source,” Phys. Rev. Lett. 89, 113601 (2002). 35. R. S. Bennink, S. J. Bentley, and R. W. Boyd, “Quantum and classical coincidence imaging,” Phys. Rev. Lett. 92, 033601 (2004). 36. A. Gatti, E. Brambilla, M. Bache, and L. A. Lugiato, “Ghost imaging with thermal light: comparing entanglement and classical correlation,” Phys. Rev. Lett. 93, 093602 (2004). 37. A. Gatti, E. Brambilla, M. Bache, and L. A. Lugiato, “Correlated imaging, quantum and classical,” Phys. Rev. A 70, 013802 (2004). 38. J. Cheng and S. Han, "Incoherent coincidence imaging and its applicability in x-ray diffraction," Phys. Rev. Lett. 92, 093903 (2004). 39. D. Cao, J. Xiong and K. Wang, “Geometrical optics in correlated imaging systems,” Phys. Rev. A 71, 013801 (2005). 40. Y. Cai and S. Zhu, “Ghost imaging with incoherent and partially coherent light radiation,” Phys. Rev. E 71, 056607 (2005). 41. F. Ferri, D. Magatti, A. Gatti, M. Bache, E. Brambilla, and L. A. Lugiato, “High-resolution ghost image and ghost diffraction experiments with thermal light,” Phys. Rev. Lett. 94, 183602 (2005). 42. A. Valencia, G. Scarcelli, M. D'Angelo, and Y. Shih, “Two-photon imaging with thermal light,” Phys. Rev. Lett. 94, 063601 (2005). 43. A. Gatti, M. Bache, D. Magatti, E. Brambilla, F. Ferri, and L. A. Lugiato, “Coherent imaging with pseudothermal incoherent light,” J. Mod. Opt. 53, 739-760 (2006). 44. D. Zhang, Y. Zhai, L. Wu, and X. Chen, “Correlated two-photon imaging with true thermal light,” Opt. Lett. 30, 2354-2356 (2005). 45. G. Scarcelli, V. Berardi, and Y. Shih, “Phase-conjugate mirror via two-photon thermal light imaging,” Appl. Phys. Lett. 88, 061106 (2006). 46. L. Basano and P. Ottonello, “Experiment in lensless ghost imaging with thermal light,” Appl. Phys. Lett. 89, 091109 (2006). 47. M. Bache, D. Magatti, F. Ferri, A. Gatti, E. Brambilla, and L. A. Lugiato, “Coherent imaging of a pure phase object with classical incoherent light,” Phys. Rev. A 73, 053802 (2006) 48. R. Borghi, F. Gori, and M. Santarsiero, ”Phase and amplitude retrieval in ghost diffraction from fieldcorrelation measurements,” Phys. Rev. Lett. 96, 183901 (2006). 49. M. Zhang, Q. Wei, X. Shen, Y. Liu, H. Liu, J. Cheng, and S. Han, “Lensless Fourier-transform ghost imaging with classical incoherent light,” Phys. Rev. A 75, 021803 (2007). 50. Y. Bai, H. Liu, and S. Han, “Transmission area and correlated imaging,” Opt. Express 15, 6062-6068 (2007). 51. J. Chen and S. Han, “Classical correlated imaging from the perspective of coherent-mode representation,” Phys. Rev. A 76, 023824 (2007).

#102691 - $15.00 USD

(C) 2009 OSA

Received 13 Oct 2008; revised 29 Nov 2008; accepted 2 Jan 2009; published 5 Feb 2009

16 February 2009 / Vol. 17, No. 4 / OPTICS EXPRESS 2454

52. S. Crosby, S. Castelletto, C. Aruldoss, R. E. Scholten, and A. Roberts, “Modelling of classical ghost images obtained using scattering light,” New J. Phys. 9, 285 (2007). 53. V. Torres-Company, H. Lajunen, J. Lancis, and A. T. Friberg, “Ghost interference with classical partially coherent light pulses,” Phys. Rev. A 77, 043811(2008). 54. D. Cao, J. Xiong, S. Zhang, L. Lin, L. Gao, and K. Wang, “Enhancing visibility and resolution in Nth-order intensity correlation of thermal light,” Appl. Phys. Lett. 92, 201102 (2008). 55. H. Liu and S. Han, “Spatial longitudinal coherence length of a thermal source and its influence on lensless ghost imaging,” Opt. Lett. 33, 824-826 (2008). 56. F. Ferri, D. Magatti, V. G. Sala, and A. Gatti, “Longitudinal coherence in thermal ghost imaging,” Appl. Phys. Lett. 92, 261109 (2008). 57. A. Dogariu and E. Wolf, “Coherence theory of pairs of correlated wave fields,” J. Mod. Opt. 50, 1791-1796 (2003). 58. M. Santarsiero and R. Borghi, “Correspondence between super-Gaussian and flattened Gaussian beams,” J. Opt. Soc. Am. A 16, 188-190 (1999). 59. F. Gori, “Flattened Gaussian beams,” Opt. Commun. 107, 335-341 (1994). 60. Y. Li, “Light beam with flat-topped profiles,” Opt. Lett. 27, 1007-1009 (2002). 61. A. A. Tovar, “Propagation of flat-topped multi-Gaussian laser beams,” J. Opt. Soc. Am. A 18, 1897-1904 (2001). 62. F. Wang, Y. Cai, and S. He, “Experimental observation of coincidence fractional Fourier transform with a partially coherent beam,” Opt. Express 14, 6999-7004 (2006).

1. Introduction Conventionally partially coherent beams are characterized by a second-order correlation function of the electric field which depends on two spatial arguments and the oscillation frequency and being called the cross-spectral density [1]. Gaussian Schell-model (GSM) beams are partially coherent beams whose spectral density, i.e. the cross-spectral density evaluated at a single argument and spectral degree of coherence, i.e. a suitably normalized version of the cross-spectral density have Gaussian shapes [2-9]. By scattering a coherent laser beam from a rotating grounded glass, then transforming the spectral density distribution of the scattered light into Gaussian profile with a Gaussian amplitude filter, a GSM beam can be generated [6, 7]. GSM beams can also be generated with specially synthesized rough surfaces, spatial light modulators and synthetic acousto-optic holograms (c.f. [9]). A more general partially coherent beam can possess a twist phase, which differs in many respects from the customary quadratic phase factor. Simon and Mukunda first introduced the partially coherent twisted Gaussian Schell-model (GSM) beam [10] opening up “a new dimension” in the area of partially coherent fields [10, 11]. Unlike the usual phase curvature, the twist phase is bounded in strength due to the fact that the cross-spectral density function must be non-negative definite and, moreover, it is absent in a coherent Gaussian beam [10]. The twist phase has an intrinsic chiral or handedness property and is responsible for the rotation of the beam spot on propagation [10-12]. An essential aspect of the twist phase is that it is intrinsically two-dimensional, and it cannot be separated into a sum of one-dimensional contributions [10-12]. Experimental observation of the twisted GSM beam has been reported in Ref. [11]. Later the studies relating to superposition [13], coherent-mode decomposition [14, 15] and the analysis of the transfer of radiance [16] of the twisted GSM beams have been carried out. The conventional method for treating the propagation of twisted GSM beams is Wignerdistribution function [10, 17]. Lin and Cai have introduced a convenient alternative tensor method for treating propagation of GSM and twisted GSM beams [18]. Paraxial propagation of the twisted GSM beams through free space, paraxial optical systems and some random linear media, have been then studied in details [16-23]. Dependence of the orbital angular momentum of a partially coherent beam on its twist phase was revealed in Ref. [24]. More recently, the influence of the twist phase on the second-harmonic generation by a partially coherent beam has been investigated [25]. Ghost imaging, also known as quantum or correlated, or coincidence or two-photon imaging, has been studied extensively in recent years [26-56]. It provides a technique for nonlocal image production of an object through the measurement of the fourth-order

#102691 - $15.00 USD

(C) 2009 OSA

Received 13 Oct 2008; revised 29 Nov 2008; accepted 2 Jan 2009; published 5 Feb 2009

16 February 2009 / Vol. 17, No. 4 / OPTICS EXPRESS 2455

correlation function, i.e., intensity correlation function, calculated for two optical paths with the object placed in one of them (see Fig.1). Ghost imaging was first demonstrated with entangled photon pairs generated by spontaneous parametric down conversion in 1995 [26]. Since then many relating experimental results were obtained and, soon, ghost imaging has found wide applications in quantum metrology, lithography and holography [17-33]. Recently, both theoretical and experimental results have shown that ghost imaging can be realized with classical momentum-correlated laser pulses [34, 35] and partially coherent light (i.e., pseudo-thermal light) [36-56]. In particular, it has been found that ghost imaging with partially coherent light is highly dependent on coherence properties of the source [36-56]. The visibility (typical quality criterion) of the image is generally higher in the quantum case than in the classical, however, due to the fact that ghost imaging with classical partially coherent light is much simpler realizable, it has its own advantage for practical applications [6, 7]. In any of the previous theoretical investigations on ghost imaging with partially coherent light, the twist phase of the cross-spectral density of the light source has not been taken into consideration. Therefore, since in the absence of the twist phase the two-dimensional ghost imaging problem can be regarded as a product of two one-dimensional ghost imaging problems, only one-dimensional ghost imaging has been previously studied. If the partially coherent light source possesses a twist phase, two-dimensional ghost imaging problem cannot longer be regarded as a product of two one-dimensional ghost imaging problems, and one, therefore, has to study ghost imaging in the two-dimensional setting. In this paper, we study the ghost imaging with partially coherent twisted GSM beam in two-dimensions with the help of the tensor method. Thus for the first time, the twist phase of the cross-spectral density of the light source is taken into account for ghost imaging. The influence of the twist phase on ghost imaging is studied in detail. It is found that the ghost image disappears gradually but its visibility increases as the absolute value of the twist factor increases. This phenomenon is caused by the fact that the twist phase enhances the transverse spatial coherence of the twisted GSM beam on propagation. 2. Second-order correlation function of partially coherent twisted GSM beam We begin by a brief review of the statistical description of the twisted GSM beams using the second-order coherence theory in the space-frequency domain [1] and the tensor method for treatment of interaction of the beam with paraxial optical systems. According to Refs. [13-15], the cross-spectral density of a twisted GSM beam in source plane z = 0 is expressed as follows ª r 2  r 2  r  r 2 ik P º T 0 exp «  1 2 2  1 22   r1  r2 J  r1  r2 » , (1)  2V g 0 2 «¬ 4V I 0 »¼ where frequency dependence was omitted for conciseness. In Eq. (1) U is a scalar field, k 2S / O is the wave number with O being the wavelength of light field, r1 and r2 are two W  r1 , r2

U (r1 )U * (r2 )

arbitrary position vectors in the source plane, V I 0 and V g 0 denote the transverse beam width and spectral coherence width, respectively, P 0 is a scalar real-valued twist factor with the 1

dimension of an inverse distance, limited by the double inequality 0 d P02 d ª k 2V g4 0 º due to ¬ ¼ the non-negative definiteness requirement of Eq. (1). We note that in the coherent limit, V g 0 o f , the twist factor P0 disappears. J is an anti-symmetric matrix, viz.,

§ 0 1· (2) ¨ 1 0 ¸ . © ¹ Under condition P0 0 the cross-spectral density in Eq. (1) reduces to the one of the conventional GSM beam [1-9]. Due to the existence of the term J

#102691 - $15.00 USD

(C) 2009 OSA

Received 13 Oct 2008; revised 29 Nov 2008; accepted 2 Jan 2009; published 5 Feb 2009

16 February 2009 / Vol. 17, No. 4 / OPTICS EXPRESS 2456

 r1  r2 J  r1  r2 T

x1 y2  x2 y1 in the right side of Eq. (1), the two-dimensional cross-spectral density cannot be split in a product of two one-dimensional cross-spectral densities. Thus, we must treat the propagation of the twisted GSM beams using full two-dimensional versions of their cross-spectral densities. We note here that one should not confuse the twist factor of the GSM beam with its wave-front radius of curvature (i.e., the usual quadratic phase factor): while the former is present in Eq. (1) the later is absent. We also should point out that the modulus of the spectral degree of coherence ( W  r1 , r2 / W  r1 , r1 W  r2 , r2 ) of a twisted GSM

beam in the source plane is determined solely by parameter V g 0 , and is independent of parameters V I 0 and P0 . Equation (1) can alternatively be expressed in the following tensor form [18]

§ ik · (3) W  r exp ¨  r T M 0 1r ¸ , 2 © ¹ y2 , and M 0 1 is a 4 u 4 matrix, called partially coherent

where r T  r1T r2T  x1 y1 x2 complex curvature tensor, takes the form

M 01

§§ i i ·  2 ¸I ¨¨  2 ¨ ¸ 2 V V k k ¨© I0 g0 ¹ ¨ i ¨ T ¨ kV 2 I  P0 J g0 ©

i kV

2 g0

§ i ¨¨  2 V k 2 I0 ©

· ¸ ¸ ¸, i · ¸  2 ¸I¸ kV g 0 ¹¸ ¹

I  P0 J

(4)

I being a 2 u 2 unit matrix. Following Ref. [18], after propagation through a general astigmatic ABCD optical system, the cross-spectral density of the twisted GSM beam in the output plane is expressed as follows W  u

1 § ik · exp ¨  u T M 1u ¸ , 1/ 2 1 2   © ¹ ª¬ Det( A  BM 0 ) º¼

where Det stands for the determinant of a matrix, u T

u

T 1

(5)

uT2 with u1 and u 2 being two

arbitrary position vectors in the output plane, M 1 is the partially coherent complex curvature tensor in the output plane, and is related with M 0 1 through the following tensor ABCD law M 1

  C  DM  A  BM 1

0

1

0

1

,

(6)

 , B,C   and D  have the forms Here matrices A  A

§ A 0I ·  § B 0I ·  § C 0I ·  § D 0I · ,B ¨ ,C ¨ ,D ¨ , ¨ *¸ *¸ *¸ *¸ © 0I A ¹ © 0I  B ¹ © 0I C ¹ © 0I D ¹

(7)

with A, B,C and D being the 2 u 2 sub-matrices of the astigmatic optical system. Symbol “*” in Eq. (7) which denotes the complex conjugate is required for a general optical system in the presence of loss or gain, although it does not appear in Eq. (13) of Ref. [18]. 3. Ghost imaging with partially coherent twisted GSM beam

We will now use the twisted GSM beam as the illumination field in the typical lensless ghost imaging setup shown in Fig. 1 [45, 46, 55, 56]. A beam generated by a twisted GSM source is split into two beams which are sent to paths 1 and 2. Along path 1 the beam propagates first to an object with transmision function H  v1 and then to single-photon detector D1 (located

#102691 - $15.00 USD

(C) 2009 OSA

Received 13 Oct 2008; revised 29 Nov 2008; accepted 2 Jan 2009; published 5 Feb 2009

16 February 2009 / Vol. 17, No. 4 / OPTICS EXPRESS 2457

at u1 0 ). Along path 2 the beam propagates to bucket detector D 2 . We assume that the distance from the light source to the object is the same as the distance from the light source to deterctor D2 [45, 46, 55, 56]. Finally, the output signals from detectors D1 and D2 are sent to an electronic coincidence circuit in order to measure the coincidence counting rate [i.e., the fourth-order correlation G  2  u1 , u 2 ] of two detected signals.

Twisted GSM source

z1

z

D1

Beam splitter

z

v1

Path two

r

u1

Object

Path one

Coincidence circuit

D2

u2

Fig. 1. Scheme for lensless ghost imaging with a twisted GSM beam

With the help of the moment theorem for the complex Gaussian random fields in classical optical coherence theory, the fourth-order correlation (i.e., intensity correlation) function of the fields incident onto detectors D1 and D 2 can be written as [1, 29, 35-53] G (2)  u1

0, u 2  I (u1

0) I (u 2 ) !  I (u1

0) ! I (u 2 ) !  W  u1

0, u 2 , 2

(8)

where  I (u1 0) ! and  I (u 2 ) ! denote the spectral densities of the beams on the surfaces of detectors D1 and D 2 , respectively and, as before, W  u1 0, u 2 is the mutual cross-spectral density of fields at two detectors [57]. The term W  u1 0, u 2 2 denotes, in fact, the intensity fluctuation correlation (  I (u1 0) I (u2 ) !   I (u1 0) ! I (u2 ) ! ), which contains the image of the object. The term  I (u1 0) ! I (u 2 ) ! represents the background component and should be, therefore, subtracted from the intensity correlation term. Within the validity of the paraxial approximation  I (u1 0) ! ,  I (u 2 ) ! and

W  u1

0, u 2 at the two detectors are given by the expressions W  u1 , u 2  I (ui ) !

f

f

f

f

³ ³ ³ ³ U (r )U (r ) h  u , r h  u , r dr dr , ³ ³ ³ ³ h  u , r h  u , r U  r U  r dr dr , i f

*

1

f f f f f f f i

f f f f

i

2

i

1

1

1

* 2

1

2

2

1

i

2

1

2

1

2

(9)

2

1, 2

(10)

where U  r1 U  r2 is the second-order correlation function of the twisted GSM field in the source plane given by Eq. (1) or Eq. (3); h1  u1 , r1 , h2  u 2 , r2 are the response functions of the two optical paths having the forms ª ik ik 2 2º exp «  r1  v1   v1  u1 »H ( v1 )dv1 , z z 2 2 ¬ ¼ 1 i 2º ª ik h2 (u 2 , r2 )  exp «   r2  u 2 » , Oz ¬ 2z ¼

h1 (u1 , r1 )

1

f

f

O 2 zz1 ³f ³f

(11) (12)

Substituting from Eqs. (3), (11) and (12) into Eq. (9), we obtain (after some operation) the formula #102691 - $15.00 USD

(C) 2009 OSA

Received 13 Oct 2008; revised 29 Nov 2008; accepted 2 Jan 2009; published 5 Feb 2009

16 February 2009 / Vol. 17, No. 4 / OPTICS EXPRESS 2458

W  u1

0, u 2 =

where v Tu

f f f f f f i ik ª ik º exp «  r T M 01r   r T B 1r T  2r T B 1 v u  v Tu B 1 v u » 1/ 2 ³f ³f ³f ³f ³f ³f 2 2  ¬ ¼ ª º O z1 ¬ Det  B ¼ ª ik 2 º v1 » H ( v1 ) drdv1 , (13) u exp «  ¬ 2 z1 ¼ 3

 v1

u 2 , dr

dx1dy1dx2 dy2 , and

dr1dr2

 § zI 0 I · . B ¨ ¸ © 0I  zI ¹

(14)

After vector integration over r , Eq. (13) reduces to f f ª ik i ik 2 º W  u1 0, u 2 = exp «  v Tu M 1 v u  v1 »H ( v1 )dv1 , 1/ 2 ³f ³f 1  z1 ¼ 2 2  º ¬ O z1 ª¬ Det  I  BM 0 ¼ where

(15)

I denotes the 4 u 4 unit matrix, and

M 1

M

0

1 § M11 ¨ 1 T ¨  M12 ©

1  B

1 · M12 ¸. 1 ¸ M 22 ¹

(16)

and the first part of the exponential term under the integral sign, on the right side of Eq. (15), can be expanded as ik § ik · § ik · 1 1 exp ¨  v Tu M 1 v u ¸ exp ¨  v1T M11 v1  ikv1T M12 u 2  uT2 M 221u 2 ¸ . 2 2 2 © ¹ © ¹

(17)

After performing several algebraic operations, we find that Eq. (15) becomes W  u1

with B1

0, u 2

i

§ ik · exp ¨  uT2 M 221u 2 ¸ © 2 ¹

 ª Det  I  BM º¼ ¬ f f ik ª ik º 1 1 v1  v1T B11 v1  ikv1T M12 u 2 »H ( v1 )dv1 , u ³ ³ exp «  v1T M11 f f 2 ¬ 2 ¼

O ª¬ Det  B1 º¼

1/ 2

1 0

1/ 2

(18)

z1I .

Fig. 2. 3D-distribution and corresponding contour graph of the square of the transmission function ( H ( v ) 2 ) of a soft circular aperture versus v and v with a1 1 1x 1y

0.1mm and N=30

For the convenience of calculations and analysis, we assume that the object in path one is a soft circular aperture, whose transmission function can be expressed as a finite sum of Gaussian modes H ( v1 )

where C

0

N

¦ n 1

 1

n 1

N

1 C0

N

¦ n 1

 1 N

n 1

§ nv12 · §N· ¨ ¸ exp ¨  2 ¸. ©n¹ © a1 ¹

(19)

§ N · is a normalization factor. Note that there are many different models ¨ ¸ ©n¹

for describing the transmission function of a soft circular aperture, such as super-Gaussian

#102691 - $15.00 USD

(C) 2009 OSA

Received 13 Oct 2008; revised 29 Nov 2008; accepted 2 Jan 2009; published 5 Feb 2009

16 February 2009 / Vol. 17, No. 4 / OPTICS EXPRESS 2459

model [58], flattened-Gaussian model [59], flat-topped model [60] and flat-topped multiGaussian model [61]. We have chosen the flat-topped model in this paper just for mathematical convenience. Figure 1 shows the 3D-distribution and the corresponding contour graph of the square of the transmission function ( H ( v1 ) 2 ) of a soft circular aperture versus v1x and v1y with a1

0.1mm and N=30. It may be readily demonstrated, at least

numerically, that the width of the aperture increases with increase in a1 . We note that in the ghost imaging experiments, there is no special restriction on the transmission function of the object [26-56]. In other words, the ghost imaging method works similarly for all objects with soft boundaries, for sufficiently large number N. Equation (19) can alternatively be expressed in the following tensor form as 1 C0

H ( v1 )

N

¦

 1

n 1

N

n 1

§N· § ik T · ¨ ¸ exp ¨  v1 H n v1 ¸, n 2 © ¹ © ¹

(20)

where Hn

2n I. ika12

(21)

Substituting from Eq. (20) into Eq. (18), we obtain, after integration over v1 , the expression W  u1

0, u 2

 1

n 1

§N· 1 ¨ ¸ 1/ 2 N © n ¹ ª Det  B M 1  B H  I Det  I  BM  1 º n 1 1 11 1 n 1 ¬ ¼ 1 ik § ik · 1 T 1 1 M11 u2 ¸ . u exp ¨  uT2 M 221u 2  uT2  M12  B11  H n M12  2 2 © ¹ 1 C0

N

¦

(22)

In a similar way, we obtain the following expressions for the intensities of beams at two detectors 1 N N  1 ¦¦ C02 n 1 m 1 N 2

n m

 I (u1

0) ! =

§ N ·§ N · 1 , ¨ ¸¨ ¸ 1/ 2  1  I º © n ¹ © m ¹ ª Det  B 1M 1  B 1H nm  I Det  BM 0 ¬ ¼

 I (u 2 ) ! =

where u T2

u

T 2

1 ª ik exp «  u T2 M 1u 2 1/ 2 1   ¬ 2 ª Det  I  BM 0 º ¬ ¼

º »¼ ,

(23)

(24)

 is given by Eq. (14), M 1 is given by Eq. (16), B  and H are uT2 , B 1 nm

expressed as follows B 1

§ z1I 0I · ¨ ¸, © 0I  z1I ¹

H nm

2 § nI 0I · ¨ ¸. ika12 © 0I mI ¹

(25)

4. Numerical example In this section, we calculate numerically the ghost image formed with a typical twisted GSM beam using the formulae derived in section 3.

#102691 - $15.00 USD

(C) 2009 OSA

Received 13 Oct 2008; revised 29 Nov 2008; accepted 2 Jan 2009; published 5 Feb 2009

16 February 2009 / Vol. 17, No. 4 / OPTICS EXPRESS 2460

Fig. 3. Normalized ghost image W  u 1

0, u 2 / W  u1 2

0, u 2

2 0 of a soft circular

aperture formed with twisted GSM beam for different values of the initial transverse coherence width V g 0

Fig. 4. Normalized ghost image W  u 1

0, u 2 / W  u1 2

0, u 2

2 0 of a soft circular

aperture formed with twisted GSM beam for different values of the initial twist factor P0

#102691 - $15.00 USD

(C) 2009 OSA

Received 13 Oct 2008; revised 29 Nov 2008; accepted 2 Jan 2009; published 5 Feb 2009

16 February 2009 / Vol. 17, No. 4 / OPTICS EXPRESS 2461

Fig. 5. Normalized ghost image W  u 1

0, u 2 / W  u1 2

aperture formed with twisted GSM beam for different values of

0, u 2

2 0 of a soft circular

z

Figure 3 shows the normalized ghost image W  u1 0, u2 2 / W  u1 0, u2 0 2 of a soft circular aperture formed with twisted GSM beam for different values of the initial transverse coherence width V g 0 with O 632.8nm, V I 0 10mm, P0 0, z 200mm , z1 20mm , a1 0.1mm and N=30. One finds from Fig. 3 that the ghost image is closely related to transverse coherence width V g 0 . When the value of V g 0 is very small (see Fig. 3(a)), a

perfect ghost image of the object is formed. As the value of V g 0 increases, the image of the object gradually disappears (see Fig. 3 (b)-(d)). Our results are consistent with those reported for usual GSM beam in [40, 43, 52]. This phenomenon can be explained by the fact that as the initial transverse coherence increases, the intensity fluctuation correlation between two points or two detectors after propagation also decreases, which leads to the disappearance of the ghost image. Figure 4 shows the normalized ghost image W  u1 0, u 2 2 / W  u1 0, u 2 0 2 of a soft circular aperture formed with a twisted GSM beam for different values of the initial twist factor P0 with fixed V g 0 0.0001mm , at distance z 200mm from the source. It is clear from Fig. 4 that the twist phase has strong influence on the ghost image. A good quality of the image of the object is formed for small values of the twist factor (see Fig. 4(a) and (b)). As the absolute value of the twist factor increases gradually, the ghost image of the object disappears gradually, although the initial transverse coherence length is chosen to be very small. Figure 5 2 2 shows the normalized ghost image W  u1 0, u 2 / W  u1 0, u 2 0 of a soft circular aperture formed with the twisted GSM beam for different values of z with V g 0 0.0001mm . The other parameters for calculation of Figs. 4 and 5 are the same as in Fig. 2. We can find from Fig. 5 that the ghost image is also closely related to z. In our case, when the distance z  500mm , a ghost image with good quality can be obtained (see Fig. 5 (a) and (b)). As the distance z increases, the image of the object disappears gradually (see Fig. 5 (c) and (d)).

#102691 - $15.00 USD

(C) 2009 OSA

Received 13 Oct 2008; revised 29 Nov 2008; accepted 2 Jan 2009; published 5 Feb 2009

16 February 2009 / Vol. 17, No. 4 / OPTICS EXPRESS 2462

In order to fully understand the effect of the initial twist phase and of the distance z on ghost image we will now analyze their dependence on the state of coherence of the twisted GSM beam on propagation. Applying Eqs. (5) - (7), the cross-spectral density W of the twisted GSM beam after propagation at distance z in free space can be expressed in the form 1 § ik · exp ¨  u T M 1u ¸ , 1/ 2 1 2   © ¹ ª¬ Det(I  BM 0 ) º¼

W  u

(26)

where M 1 is the partially coherent complex curvature tensor at distance z

M 1

M

0

 B

1

§ m11 ¨ ¨ m21 ¨ m31 ¨ © m41

m12

m13

m22 m32 m42

m23 m33 m43

m14 · ¸ m24 ¸ , m34 ¸ ¸ m44 ¹

(27)

with B is given by Eq. (14). Comparing Eqs. (3), (4) and (26), we obtain, after some computations, the expression for the transverse coherence length V gz on propagation of the form 1/ 2

V gz

§ i · ¨ ¸ © km13 ¹

1/ 2

§ i · ¨ ¸ © km24 ¹

1/ 2

§ O 2 z2 O 2 z 2V g20 · 2 2 2 2 ¨¨ 2 2  V g 0  V g 0 P0 z  ¸ . 16S 2V I40 ¸¹ © 4S V I 0

(28)

It is now clearly seen that the transverse coherence width V gz of the twisted GSM beam on propagation is closely related to propagation distance z and twist factor P0 . Of course, parameter V gz reduces to its source value V g 0 when z

0 . Figure 6 shows the dependence

of the transverse coherence width V gz of the twisted GSM beam on propagation distance z for different values of the initial twist factor P0 with O from Fig. 6 that the transverse coherence length

632.8nm and V I 0

V gz

10mm . It is clear

of the twisted GSM beam increases

practically linearly with the increase of the propagation distance (note, however, that in Eq. (28) this dependence is not strictly linear) and it does so more rapidly as the absolute value of the initial twist factor P0 increases.

Fig. 6. Dependence of the transverse coherence length

V gz

of the twisted GSM beam on the

propagation distance z for different values of the initial twist factor P0

#102691 - $15.00 USD

(C) 2009 OSA

Received 13 Oct 2008; revised 29 Nov 2008; accepted 2 Jan 2009; published 5 Feb 2009

16 February 2009 / Vol. 17, No. 4 / OPTICS EXPRESS 2463

Now we can better understand the effect of the initial twist phase and of the distance z on the ghost image: although the initial spatial coherence width ( V g 0 ) of the twisted GSM beam is independent of the initial twist phase, the spatial coherence width ( V gz ) of the twisted GSM beam on propagation depends on the initial twist phase. Moreover, since the spatial coherence width of the twisted GSM beam always increases on propagation, doing so more rapidly with increasing values of the initial twist factor, the intensity fluctuation correlation between two detectors must decrease, forcing the ghost image to disappear. Such a conclusion is numerically demonstrated by Figs. 4 and 5. At this point we stress on the fact that in the preceding numerical calculations, we have not taken the term  I (u1 0) ! I (u 2 ) ! into consideration. Experimentally, one can measure both G (2)  u1 0, u 2 (i.e., I (u1 0) I (u 2 ) and  I (u1 0) ! I (u 2 ) ! [62]. Quantity W  u1

0, u 2 can be obtained by subtracting  I (u1 2

0) ! I (u 2 ) ! from G (2)  u1

0, u 2 .

If the product  I (u1 0) ! I (u 2 ) ! is not subtracted, then the twist phase will influence the visibility of the ghost image, which is defined by the formula V

W  u1

G (2)  u1

0, u 2 max 2

0, u 2 max

.

(29)

Thus we have shown that the visibility of the ghost image enhances with the increase of the initial transverse coherence width V gz of the beam, which is in the agreement with results obtained previously in Refs. [40] and [43]. Moreover, the visibility of the ghost image also enhances as the value of the initial twist factor or of the distance z increases, because the presence of the twist factor results in the increase in the transverse coherence width of a twisted GSM beam on propagation in free space.

5. Conclusion In conclusion, we have investigated ghost imaging with a partially coherent twisted GSM beam with the help of a tensor method. We have found that the ghost image formed with a twisted GSM beam is also closely dependent on the twist factor, in addition to the dependence on the imaging distance and on the initial coherence length. The image disappears gradually as either the absolute value of the twist factor or the imaging distance increases, but the visibility of the image enhances as the values of these parameters increase. This is caused by the fact that the transverse coherence width of a twisted GSM beam grows on propagation in free space, and that the non-zero initial twist phase stimulates its increase on propagation. Thus, in any practical situations involving the use of twisted GSM beams in ghost imaging it is necessary to take both the twist phase and the imaging distance into consideration, in addition to the initial transverse coherence width.

Acknowledgments Y. Cai gratefully acknowledges the support from the Alexander von Humboldt Foundation. Q. Lin gratefully acknowledges the supports from the Ministry of Science and Technology of China (grant no. 2006CB921403 & 2006AA06A204). O. Korotkova’s research is funded by the AFOSR (grant FA 95500810102).

#102691 - $15.00 USD

(C) 2009 OSA

Received 13 Oct 2008; revised 29 Nov 2008; accepted 2 Jan 2009; published 5 Feb 2009

16 February 2009 / Vol. 17, No. 4 / OPTICS EXPRESS 2464