Sorting of spatially incoherent optical vortex modes.

Coherent optical vortices have promising applications in quantum and classical optical communication. They add new degrees of freedom to code information. In this context, to implement a tool enabling sorting of spatially multiplexed vortex states is fundamental. By other hand, spatially incoherent vortices can be more robust in propagation through noise media, such as turbulent atmosphere or obstacles that block part of the light. Therefore, in this work we propose directly applying a high-resolution sorting scheme to spatially incoherent vortex states.

Introduction

Since the discovery of the orbital angular momentum (OAM) of light or optical vortices[1], one application that have attracted more attention recently is its use as a new degree of freedom to burst optical communications in the classical[2,3] and quantum regimes[4,5]. Several orthogonal modes can be exploited to spatial multiplexing and increase the transmission capacity in optical communication systems. However, the most used optical beam is the OAM beam called Laguerre-Gauss (LG) beam[6] which can be decomposed in terms of orthogonal components allowing being space-divided with low inter-modal crosstalk among multiple modes[2,3,7]. This beam has a finite spatial extend and due to its azimuthal phase dependence can be sorted[8-10], being ideal to optical communications.

There are several proposals for OAM sorting. For example, the interferometry methods[11,12], which require increasing complexity of the scheme with increasing number of modes to be sorted. However, the best approach to OAM sorting is based on a ray-optics coordinate transformation. To date, there are two types of such mode sorting. The first is based on log-polar transformation[13] which can map OAM modes with an azimuthal phase to plane waves with tilts proportional to the vortex topological charge (TC)[8]. This approach is simple and efficient, but there is a limited separation of adjacent modes. One work improved this result at the cost of adding an additional phase-only hologram known as fan-out element that makes copies of the input optical fields and the interference of these copies can produce a better resolution[9], but it also increase the complexity of the scheme. The second type of ray-optics coordinate transformation is the spiral transformation that maps spirals to parallel lines. In this case, the tilt of the plane waves is theoretically unlimited and, in practice, the resolution of the mode separation is considerably improved[10].

Spatially incoherent optical vortices, are vortices originated by passing an incoherent light through a hologram designed to generate optical vortices[14] or when a coherent optical vortex is scattered by a rough surface[15-18]. Those works have shown that the generated light presents a vortex in the coherence function, the so-called coherence vortex[19]. Remarkable, the lower coherence vortices are less influenced than those of higher coherence during propagation through turbulent atmosphere[20] and through obstacles[21]. Therefore, discriminating incoherent optical vortices with different superimposed TCs is a topic of current interest for optical communication through turbulent atmosphere[22,23].

Results

Concept and principle

The spiral transformation method has been successfully applied for high resolution sorting of coherent OAM modes of light[10]. In this work we have applied it to OAM mode sorting of spatially incoherent OAM modes. Figure 1 illustrates the main idea of the present work. A spatially incoherent light, which actually is a speckle field, is used with a hologram designed to generate LG beams. Therefore, incoherent LG beams are generated and their optical Fourier transform[24] is projected over the transformation phase mask . The resulting field has its optical Fourier transform projected over the correction phase mask . A final optical Fourier transform generates a final speckle field. The intensity of the final speckle field fluctuates randomly, but it still has the information of the sorted modes. This information is recovered from the correlations of the spatial intensity fluctuations[25].

Figure 1

The general idea of the proposal: a spatially incoherent light represented by a speckle field, instead of coherent light, is used to generate OAM modes. The incoherent modes are passed through a usual system of coherent mode sorting and a final speckle field is produced from which we have to extract the information of the sorted modes.

The spiral transformation is generated by a phase mask, which is described by the following equation[10],where and is the light wavelength, is the focal length of a lens put at a distance behind the phase element to perform its Fourier transform[24], and , and are constants. The and are the polar coordinates but the variable is not the conventional one since it is unlimited,whereand is the floor function. A second phase element is used as a phase-correction mask such that its local gradient at a point (u; v) compensate the slope of the ray arriving at this point from the corresponding point (x; y) in the previous phase element[10]. This phase-correction element is given by,where the constants are the same defined for the function and are the spatial coordinates at the plane of the phase-correction element.

Experimental setup

The Fig. 2 displays the used experimental setup. A laser beam is expanded and collimated by a telescope formed through lenses specified by its focal length and , where a rotating ground glass disc (RGGD) is placed between these lenses near the focus of the lens of focal length . In order to get the best mode distinguishability we make use of the control of the source size[26,27] slightly dislocating the RGGD from the focus. The obtained light illuminates the right half of which was half divided in order to contain, in the right, the hologram to generate LG beams[28] and, in the left, the phase. The generated LG beam is Fourier transformed and projected over by a lens of focal length .

Figure 2

Experimental setup: and are spatial light modulators, CCD stands for charge coupled device camera, to are mirrors, and are beam splitters, the focal length of the lenses are , , , , .

After the second time reflecting in the the LG beams acquire a phase and a lens of focal length projects the Fourier transform of this new field over the which contain the correction phase . A lens of focal length projects the Fourier transform of the field reflected from the SLM2 on the CCD camera. To calculate the phase elements and we have used , and .

Theory

We present the theoretical formulation of the conceptual idea of the experiment with spatially incoherent light. The speckle field generated in the RGGG hits the first half of the SLM1 containing a hologram that encode the signal field , i.e., LG modes. The resulting field is Fourier transformed by a lens of focal length ,

The field is projected in the second half of the SLM1 containing the phase . The resulting field is Fourier transformed by a lens of focal length ,

The field in Eq. (6) hits the SLM2 which contains the phase . The resulting field is Fourier transformed by a lens of focal length producing the final field at the CCD plane,

The reference field is just a Gaussian beam codded in the hologram and it propagates through the same optical components. In order to allow an analytical calculation of the correlation function we suppose, without loss of generality, that the for this field the transformations phases are . Therefore, the correlation function is written as[25],where the symbol “*” stands for the complex conjugated field and means an ensemble average.

We consider that the speckle fields are delta-correlated, i.e., its cross-correlation at the plane of the hologram that generates the signal and reference fields is . Therefore, the integral in and can be evaluated resulting in

Equation (9) represents a convolution between the Fourier transformed sign and reference fields . Since the waist size of the reference beam is bigger than that of the signal beam, we can make the approximation implying that . We also note that the integral in in Eq. (8) results . Using all of these results Eq. (8) simplifies to,

Using and we realize that the correlation function is identical with the result that could be obtained for the coherent fields that arrives at the final plane at the CCD[10]. There is no analytical result for the integral in Eq. (10) therefore, in the next section, we present a numerical simulation which is equivalent to the effect described by this integral.

Numerical simulation and experimental results

For the numerical simulation we have used the band-limited angular spectrum method[29] to propagate in the free space between the lenses and the SLMs, where each lens of focal length is represented by a quadratic phase . Therefore, the field is multiplied by this quadratic phase to simulate the crossing of a lens. Similarly, the beam is multiplied by a phase or to simulate the reflection of the field in the SLMs. We have first to generate a speckle field which is multiplied by the LG beams and then propagated though the lens of focus until the transformation phase , and after that it follows the path through the lens of focus until the correction phase , and then through the lens until the CCD camera. To generate a speckle field we have multiplied a Gaussian function by a random phase and performed a Fourier transform of it. The size of the speckles where very small such that we can consider the field as a delta-correlated field. The reference beam is generated just passing a zero order LG beam through the system.

In the experiment, we measured the signal intensity followed by the reference intensity with RGGD stopped and then rotate the RGGD a small angle and stop it again to measure the next pair of signal and reference beams. For each pair of measurement we need to calculate the cross-correlation between the measured signal and reference intensities and we average over 100 measurements. Thanks to the Reed’s momentum theorem[30] the intensity cross-correlation can be written as[31,32],where is a background and is given by Eq. 10.

Figure 3 shows samples of the numerically calculated (first line) and experimentally measured speckle patters (second line), for the signal (first column) and reference beams (second column) at the CCD plane.

Figure 3

Random intensities patterns observed for the signal (, first column) and reference beams (, second column) at the detection plane. First line simulation and second line experiment. The windows are squares of .

Figure 4, in the first and second columns, presents the results for the intensity cross-correlation patterns for the simulation and experiment, respectively. Although the signal and reference beams fluctuates randomly, the correlation function between them tends to a well-defined pattern in agreement with the prediction of the theory. We can observe that another prediction of the theory is also confirmed, i. e., despite of a background, the patterns for the cross-correlation was the same that the intensities of the coherent fields in ref. [10].

Figure 4

Simulation (first column) and experimental results (second column) for the intensity cross-correlation functions . The windows are squares of .

Figure 5 shows transversal profiles taken at the center of the simulation (first column of Fig. 4) of the cross-correlation functions of the spatially incoherent field, and the same for the experimental results (second column of Fig. 4). In the experiment and simulation of the cross-correlation functions we have averaged over 100 calculations of the cross-correlation through 100 realization of the pair of random pattern (see Fig. 2) and have performed a background subtraction[33], resulting . The distance between two adjacent peaks is [10]. The distance between adjacent peaks obtained directly by the image in the CCD camera is , close to that value. Using this distance we can extract the topological charge from the distance of each mode to the central zero order mode. Therefore, the only previous calibration that we have to do is determining the position the zero order mode. We observe a good agreement between simulation and experiment indicating that the sorting scheme is viable with spatially incoherent light.

Figure 5

Profiles for the simulation (fist column) and experiment (second column) for the field cross-correlation function .

Conclusion

We have observed that it is possible the use the spatial degrees of freedom of the coherence function to spatially multiplex coherence vortex to encode information for optical communications. Besides we note that it is possible to implement the sorting of spatially multiplexed spatially incoherent vortex states. Considering that spatially incoherent optical vortices beams are robust against scattering by the atmospheric turbulence and opaque obstacles, it may constitute a promising alternative for codding and transmitting information.

Methods

Devices in the experimental setup

A have used a laser model Ultralasers MSL-FN-532-200mW operating at 532 nm. The SLM1 used in the experiments is a Holoeye LETO and the SLM2 is a Hamamatsu X10468-01, both are phase-only SLMs based on reflective LCOS microdisplays. The Holoeye LETO SLM have a spatial resolution of 1920 × 1080 pixels and a pixel size of and the Hamamatsu SLM have a spatial resolution of 600 × 792 pixels and a pixel size of . The resolution of the phase patterns displayed on the Holoeye SLM are 960 × 960 pixels, square shape and displayed side by side at the center. The resolution of the phase pattern displayed on the Hamamatsu SLM is 600 × 600 pixels displayed at the center. The camera used in the experiment is a PixeLink PL-B781F.