Photonics Research, 2019, 7 (12): 12001485, Published Online: Nov. 21, 2019   

Generation of coherence vortex by modulating the correlation structure of random lights Download: 612次

Author Affiliations
1 College of Optical and Electronic Technology, China Jiliang University, Hangzhou 310018, China
2 National Demonstration Base for International Science and Technology Cooperation on Micro/Nano-Fabrication & Optoelectronic Detection, China Jiliang University, Hangzhou 310018, China
Abstract
A coherence vortex (CV) carrying topological-charge information in its correlation dimension is a new option for optical manipulation and communication. CV generation by directly modulating the correlation function enables a way to control the light field in this dimension. However, few experimental realizations on this issue have been reported because of the difficulty in phase modulation when the light arrays are of low coherence. In this paper, we propose a method for generating a CV by utilizing partially coherent light arrays. A proper design of random arrays at the input plane leads to a complex CV field at the output plane after free-space propagation. This generation mechanism works well for beamlets of low coherence.

1 INTRODUCTION

A coherence vortex (CV) [13" target="_self" style="display: inline;">–3] is a common phenomenon in various laser speckles [48" target="_self" style="display: inline;">8]. The study of CVs began in 1998, when Gori et al. proposed that a partially coherent light source would possess helicoidal modes in its spectral degree of coherence [9]. For speckle fields that are partially coherent lights, the random fluctuations of the light fields tend to hide the singular points of the averaged intensity [10]. However, in the picture of two-point correlation, the singularities survive even when the fields are of low coherence [11]. A CV field, which is a special optical vortex, possesses a spiral phase structure in a two-point correlation function formed as cross-spectral density (CSD) [12]. The accumulated phase change around the singularity of the helical correlation structure is 2mπ, where m is the topological charge of the CV field and is linked to the orbital angular momentum (OAM) of the photon according to m (where is the Planck constant). The two-point correlation is responsible for the coherence of light, hence the term CV. A CV is a point singularity in the coherence of the optical field where the intensity is unnecessarily zero. The benefit from the low coherence is that a CV is less influenced than an optical vortex when it propagates in turbulence [13] and through obstacles [14]; thus, it has been applied in areas such as object identification and imaging [7,8,15], reduction in turbulence-induced degradation during propagation [16], optical manipulation [17,18], and self-reconstruction in obstacle imaging [14,19]. In quantum optics, the correlated photons can be used to transmit and process information with respect to the physical dimension of correlation [20]. Similar to its quantum counterpart, a partially coherent laser beam also possesses a spatial correlation of freedom, which is called the spectral degree of coherence and is embedded in CSD functions [21]. This property enables the use of the correlation of freedom of the light field to perform imaging tasks [22,23]. Information on correlation structures can also be transferred by CV beams when encoded with topological charges [24,25]. In communication, the generation and multiplexing strategy of multiple information channels, i.e., propagating several transmission modes within a common medium and without interacting with one another, is important to increase the transmission capacity of light [26,27]. Accordingly, correlation structures provide a routine to increase the transmission capacity with consideration of coherence as an inherent property of light beams.

Recently, the combination of mixed CV fields in the low coherence limit has been distinguished by designing a ghost imaging system [28]. Given that CV is a special mode in the coherence dimension, CV communication is possible as long as the generation of various coherence structures has been conducted. As a result, methods of CV generation and correlation modulation warrant investigation.

Previous works have revealed the generation of various structured lights that utilize the combination of special optical beams [2934" target="_self" style="display: inline;">34], and sometimes ultrafast pulses [35]. Applications have been demonstrated in optical manipulation of microparticles [36,37]. Among them [29,30], it is pointed out that an optical vortex with a net OAM can be produced by completely coherent light arrays during free-space propagation. Whether a CV can be produced in the same way is unclear. In fact, these vortices are produced differently. An optical vortex exists in the phase of the wave function, whereas a CV is embedded in spiral structures in the phase of the correlation function [38]. In the former case, the decrease in coherence or the addition of turbulence results in OAM spectrum dispersion [5,39]. Even in free-space propagation, some special slits that cause wavefront distortion by diffraction also result in the distortion of OAM spectra [40,41]. For the latter case in which spiral structures exist in the dimension of coherence, the variation in the CV spectrum during propagation has not yet been studied. To the best of our knowledge, few works have focused on CV generation by directly configuring the phase of the correlation function.

2 THEORETICAL ANALYSIS

We consider a partially coherent light (PCL) array with N fundamental beams. At the source plane (z=0), the mixed field is E(x,y;0)=eiδj=1NEj(x,y;0),where δ is a random phase and Ej is the light field of the j-indexed beamlet, Ej(x,y;0)=exp[(xaj)2+(ybj)2w2]exp(iφ0j),and w is the waist width of the beamlet. φ0j=L0αj is an initial phase, αj=π(2j1)/N is the azimuthal angle of the j-indexed beamlet, and L0 is a constant. We restrict that 0L0<N due to periodicity. Parameters aj and bj illustrate the displacement of the center of the j-indexed beamlet and are related to the radius r0 of the radial array according to aj=r0cosαj and bj=r0sinαj. The positions of the fundamental beamlets are shown in Fig. 1(a), where the number of the fundamentals is chosen as N=3 for concreteness.

Fig. 1. (a) Schematic of a PCL array at the source with N=3. Relative CV spectrum of the array at (b) z=10zr and (c) z=20zr. The beam parameters are λ=632.8  nm, w=0.157  mm, σ=0.1w, r0=1.5w, and L0=1. Here, zr=πw2/λ=0.122  m is the Rayleigh distance.

下载图片 查看所有图片

The statistical distribution of the random phase δ corresponds to a correlation within the combined field, E*(ρ1)E(ρ2)=E*(ρ1)E(ρ2)g(ρ1,ρ2). Specifically, we use a Schell model [42] to mathematically represent the correlation function, i.e., g(ρ1,ρ2)=exp(|ρ1ρ2|2/2σ2), where σ is the transverse coherence width of the beam. Hence, the CSD of the PCL array reads W(ρ1,ρ2;0)=E*(ρ1;0)E(ρ2;0)=s=1Nq=1NWsq(ρ1,ρ2;0),with Wsq(ρ1,ρ2;0)=exp[(x1as)2+(y1bs)2w2]exp[(x2aq)2+(y2bq)2w2]×exp(iφsq)g(ρ1,ρ2;0),where ρ1(x1,y1) and ρ2(x2,y2) are two arbitrary spatial points expressed in Cartesian coordinates at the source plane; * indicates the complex conjugate; and φsq=L0(αqαs) is the phase difference between the numbering index s and q.

According to coherence theory and the paraxial propagation law, the propagated CSD at z>0 can be obtained by [43] Wsq(r1,r2;z)Wsq(ρ1,ρ2;0)h1*(r1,ρ1)h2(r2,ρ2)d2ρ1d2ρ2,where r1(u1,v1) and r2(u2,v2) are two arbitrary spatial points at the output plane, and hp with p=1,2 is the response function: hp(rp,ρp)=exp[ik2B(Axp22xpup+Dup2)+ik2B(Ayp22ypvp+Dvp2)],where k=2π/λ is the wavenumber and λ is the wavelength, and A, B, and D are transfer-matrix elements of the optical system. For the case of free-space propagation, A=D=1 and B=z. Substituting Eqs. (3) and (4) into Eq. (5) yields the CSD function of the PCL array [Fig. 1(a)] at the output plane as follows: W(r1,r2;z)=s=1Nq=1NWsq(r1,r2;z),where Wsq(r1,r2;z) reads Wsq(r1,r2;z)=π2M1M2(1λ|B|)2exp(iφsq)exp(2r02w2)exp(N22+N124M2)×exp[14M1(2bsw2+ikB*v1)2+14M1(2asw2+ikB*u1)2]×exp[ikD*2B*(u12+v12)+ikD2B(u22+v22)],with M1=1w2+12σ2+ikA*2B*,M2=1w2+12σ214σ4M1ikA2B,N1=1w2(2aq+asσ2M1)+ik(u12σ2M1B*u2B),N2=1w2(2bq+bsσ2M1)+ik(v12σ2M1B*v2B).

To examine the existence of CVs in the propagated radial array, the output CSD can be projected to OAM modes [44,45]: P(L)|0002π02πW(r1,r2;z)exp(iL·Δθ)r1r2dr1dr2dθ1dθ2|,where r1(r1,θ1) and r2(r2,θ2) are re-expressed in polar coordinates, and Δθ=θ1θ2. P(L) is a weight coefficient that represents the probability for an OAM mode with the topological charge as L. The center correlation function of the field is the CSD taken at points of (0, 0) and (r2,θ2), i.e., W(0,0,r2,θ2). On this basis and by substituting Eqs. (6)–(8) into Eq. (9), under the assumption of the paraxial propagation w2λz and low coherence source σw, we obtain the mode coefficient as P(L)(12π)|s=1Nq=1NPsq(L)|,where Psq(L)=Q1×exp[i(φsqL·αsq(+))]exp(4r02cos2αsq()w4M2)×n=0[4k2r02cos2αsq()w4M2k(k4iM2z)]n+|L|/2Γ(n+|L|/2+1)Q2×n!Γ(n+|L|+1),with the parameters Q1 and Q2 as Q1=2π3M1M2λ2z2exp(2r02w02+r02M1w04),Q2=π2M2λ2z2iπλz,and αsq(+)=(αs+αq)/2, αsq()=(αqαs)/2, Γ(·) as the Euler gamma function. During the derivation of Eq. (11), the following formulas are used [45,46]: 02πexp[inϕ1+ξcos(ϕ1ϕ2)]dϕ1=2πexp(inϕ2)In(ξ),In(ξ)=k=01k!Γ(k+n+1)(ξ2)2k+n,with  Re(n)>1,0x2m+1exp[qx2]dx=12qm+1Γ(m+1),with  Re(q)>0  and  m>1.

In Eq. (11), the magnitude of the series summation converges over n when σw and w2λz are satisfied. Therefore, a finite summation can be used to replace the infinite one during the calculation of Psq(L) under low-coherence conditions. According to the formula of P(L), one can simulate the OAM spectrum of the PCL array at different propagation distances. Numerical simulations reveal that P(L) comprises items at L=mN with m integers. As shown in Figs. 1(b) and 1(c), the largest contribution to the CV spectrum originates from the item with L=0. Side modes are distributed symmetrically at L=mN with m integers. This symmetry does not vary with the change in propagation distance. Hence, the total OAM of the combined field is always zero.

Figure 2 simulates the averaged intensity distribution I(r,z) and the correlation phase distribution ψ(x,y) of a radial array of three fundamentals (N=3). ψ(x,y) is the phase of the center correlation function W(0,0,x,y). The simulation process is as follows. To obtain the averaged intensity in Figs. 2(a) and 2(b), one makes use of Eq. (6) by setting r1=r2=r, i.e., W(r,r;z)=I(r,z). To construct the correlation phase in Figs. 2(c) and 2(d), one extracts the phase of W(0,0,x,y) by setting r1=(0,0) and r2=r=(x,y) in Eq. (6). The numerical results of Figs. 2(a) and 2(b) show that the combined field cannot maintain its initial irradiance distribution after far-field propagation. Figures 2(c) and 2(d) show that phase modulation among different beamlets at the source does not lead to a spiral phase structure at the far field.

Fig. 2. Numerical simulation of (a) spatial distribution of the averaged intensity for N=3 and z=0; (b) spatial distribution of the averaged intensity for N=3 and z=20zr; (c) correlation phase of a radial array for N=3 and z=0; and (d) correlation phase of a radial array for N=3 and z=20zr. The other beam parameters are the same as those in Fig. 1.

下载图片 查看所有图片

The numerical simulations correlating the PCL array in a Schell model of Gaussian function cannot produce a CV. The reason is that in free-space propagation, the irradiance distributions of PCL beams at the far field are determined by their correlation functions [47]. The correlation among different spatial points ruins the initial phase configuration among fundamental beamlets. Hence after far-field propagation, a vortex phase structure cannot be formed at the output plane. Although many coherence models have been proposed in the literature, experimental demonstrations are limited. Methods for the direct modulation of the complex field in the physical dimension of correlation are required.

3 EXPERIMENTAL METHOD

We present an experiment to demonstrate a strategy of generating a CV utilizing the previous radial arrays by modulating the correlation structure itself. Our approach is based on the work of Ref. [48] where a spatial light modulator (SLM) is utilized to produce a partially coherent light. Different from the work of Ref. [48], in this paper a PCL array is produced which provides a way to modulate the correlation functions azimuthally.

The experimental setup is shown in Fig. 3(a) [microscope objective (MO) is excluded], which digitally generates PCL fields by utilizing a spatial light modulator [48]. The PCL field of interest is a radial array with N fundamental beamlets [Eqs. (1) and (2)]. Here, we set N=3. A linearly polarized light from a He–Ne laser source (λ=632.8  nm, 2 mW; Newport R-30989) is expanded with a beam expander and then goes through a beam splitter before shining on the SLM (phase modulation type, 1280×720, 6.3 μm, 60 Hz; UPOLabs, HDSLMM63R). The polarization of the incident light matches the SLM requirement. The first order of the reflected light of the SLM is selected by a pinhole (P1) and then propagates through a neutral density filter. The light is reflected by a reflecting mirror to increase the propagation distance on the optical table and is then detected by a camera. The free-space propagation distance from the SLM to the camera is z. The SLM displays a movie at a certain rate. The movie consists of a series of computer-generated holographs (CGHs). The kth frame of the movie is taken as an example. It is a CGH that is designed to produce a superposed field E(k) formulated as [48] E(k)(x,y)=p=1M1[E1(p)(x,y)+E2(p)(x,y)+E3(p)(x,y)]exp(iδp),where δp is a random phase that ranges from 0 to 2π, and Ej(p) is the jth beamlet (j=1,2,3) of the pth array (p=1,2,,M1) in the superposed field E(k) [Eq. (12)]. We have Ej(p)=uj(p)exp(iφ0j) with the amplitude defined as uj(p)(x,y)=exp{[xajxj(p)]2+[ybjyj(p)]2w2}.Parameters aj, bj, φ0j, and w are the same as those in Eq. (2). The shifted coordinates (xj(p),yj(p)) representing the positions of the perturbation centers are random numbers that are uniformly distributed within a circle of radius c0. Notably, the fundamental beamlets are perturbed separately. In other words, the perturbation centers (xj(p),yj(p)) have no connection with each other. We denote this kind of randomness modulation as independent perturbation. Figures 3(b) and 3(c) show examples of the CGHs displayed on the SLM to produce the PCL array under different coherence conditions. Notably, the CGH modulation for the three beamlets merges together with increased c0. The parameter c0 is inversely proportional to the coherence width of the fundamental PCL beamlet [48]. If the array is of total coherence, then c0=0 and δp=0; if the array is of poor coherence, then c0>w and δp is a random phase. As shown in Fig. 3(d), the instantaneous intensity of the array with c0=3w, δp(0,2π) appears as speckles. As shown in Fig. 3(e), which is calculated according to Eq. (12), the phase of the field E(k)(x,y) has a random distribution. The speckles and phases change randomly with time, which guarantees that the light field is not totally coherent in the time domain.

Fig. 3. (a) Experimental setup to generate a CV with a radial array of N fundamental PCL beams. Examples of the CGH displayed on the SLM for the production of the complex field (N=3) under the following conditions: (b) total coherence with c0=0, δp=0; and (c) poor coherence with c0=3w, δp(0,2π). (d) Instantaneous intensity I(k)(x,y) of the PCL array produced with c0=3w, δp(0,2π). (e) Phase distribution of the complex light field E(k)(x,y) produced with the CGH of (c). The array is in the random state of independent perturbation. The beam parameters are w=0.157  mm, r0=1.5w, and L0=1. BE, beam expander (Newport T81-3X); BS, beam splitter; SLM, spatial light modulator; NF, neutral density filter (Absorptive, optical density 1.0); RM, reflecting mirror; P1, pinhole; MO, microscope objective; CMOS, complementary metal oxide semiconductor camera; PC1, PC2, personal computers.

下载图片 查看所有图片

By incoherent superposition of a large number of realizations of the field E(k), where k=1,2,,M2, the averaged CSD at the detection plane can be obtained according to W(x1,y1,x2,y2)=E(k)*(x1,y1)E(k)(x2,y2)=1M2k=1M2E(k)*(x1,y1)E(k)(x2,y2).The magnitude of the CSD can be obtained according to the Hanbury Brown and Twiss (HBT) effect [25]. A camera records the intensity pattern from each CGH frame, i.e., I(k)(x,y). With the data of the instantaneous intensities, the second-order correlation can be obtained by computation, G(2)(x1,y1,x2,y2)=I(k)(x1,y1)I(k)(x2,y2). Accordingly, the magnitude of the CSD can also be computed, |W(x1,y1,x2,y2)|2=G(2)(x1,y1,x2,y2)I(k)(x1,y1)·I(k)(x2,y2). To reconstruct the cross-correlation function (CCF), i.e., W(x,y,x,y), the transformation (x1,y1)(x,y), (x2,y2)(x,y) is necessary in Eq. (14).

In the experiment, the number of coherently superposed fields is M1=500. The frame rate of the SLM is 60 fps. The shooting frequency of the camera should be below 60 Hz; for example, we have utilized 19.8 and 49.6 Hz for the setting of the cameras. The number of CGH frames playing in the SLM movie is 1000. The number for ensemble average in Eq. (14) is M2=200. In the experiment of the production of the PCL array, the power ratio of the first-order diffracted beam from the SLM is around 10%.

Figure 4 presents the experimental results of the averaged intensity of the independently perturbed random array at different propagation distances. The parameter z presents the distance from the SLM to the camera. As shown in Fig. 4(a), the random array has a single-beam shape at a near field (z=0.2  m1.6zr). At a longer propagation distance z=0.57  m4.7zr, a dark center appears inside the light field, as shown in Fig. 4(b). As the generated beam propagates to z=1.14  m9.3zr [Fig. 4(c)], the dark hollow structure maintains its shape and the center of the beam darkens. The structure appearing at a far field reveals the shape of the correlation function at the source. The far-field property has been influenced by the randomness at the source [Eqs. (12) and (13)]. In other words, the perturbation method and the initial phase configuration among different beamlets modulate the correlation structure of the field. In view of applications in optical trapping or manipulation, one may use a microscope objective to compress the vortex structure. As shown in Fig. 3(a) (MO is included), after a far-field propagation (z=1.14  m) the CV is focused with a microscope objective (10×, NA=0.25) with a focal length of 17.1 mm and a working distance of 7.3 mm. Figure 4(d) presents the averaged intensity of the CV at the focal plane, where the diameter of the vortex ring is compressed to approximately 250 μm and the maximum intensity at the ring is 5.17  mW/mm2. Comparatively, the sizes of the central hollow section in Figs. 4(b) and 4(c) are 650 μm and 1.46 mm, respectively. The intensities at the ring of Figs. 4(b) and 4(c) are 1.65  mW/mm2 and 0.34  mW/mm2, respectively. The data of the peak intensities in Figs. 4(b)4(d) is obtained when the neutral density filter is not used, and the total power of the incident beam on the detection plane is 2 mW.

Fig. 4. Experimental results for the averaged intensity of the independently perturbed PCL array at different propagation distances: (a) z=0.2  m, (b) z=0.57  m, and (c) z=1.14  m. (d) The averaged intensity of the PCL array focused by the MO with a propagation distance z=1.14  m. The beam parameters are the same as those in Fig. 3(c). In the 2D graphs, the x and y axes are in the unit of pixels, and 1 pixel is 5.2 μm. In the 1D graphs, the vertical axes representing the averaged intensity are normalized. The diameters of the dark centers are marked in the 1D graphs. The camera is a CinCam CMOS-1201 (1280×1024 at 19.8 fps, 5.2 μm).

下载图片 查看所有图片

Figure 5 presents the experimental results for the propagation of the array with different perturbation models. Figures 5(a) and 5(b) show the case of a phase modulated radial array with total coherence, where the array holds the shape of fundamental beamlets from z=0.57  m to 1.14 m. Figures 5(c) and 5(d) are the case of the independent perturbation presented in Fig. 4. A comparison of Figs. 5(a) and 5(b) with Figs. 5(c) and 5(d) demonstrates that when the number of the fundamental beamlets is small, i.e., N=3, the totally coherent array does not form a dark hollow structure according to propagation, but the specially designed random array does. In Figs. 5(e) and 5(f), we use an uniform perturbation to replace the independent one, which means that the perturbation centers in Eq. (13) transform as (xj(p),yj(p))(x(p),y(p)). This finding means that the fundamental beamlets of the array are perturbed uniformly. In the uniform-perturbation model, a dark hollow beam is not produced. A comparison of Figs. 5(e) and 5(f) with Figs. 5(c) and 5(d) shows that the independent perturbation model leads to the hollow structure in the far field. To verify the effect of the initial phase on the propagation of the independently perturbed array, Figs. 5(g) and 5(h) present a case without phase modulation L0=0, which reveals that the propagated random array is in a single-beam shape and no vortex structure is formed. A comparison of Figs. 5(g) and 5(h) with Figs. 5(c) and 5(d) shows that the dark hollow structure at the far field is a direct result of the phase modulation at the source.

Fig. 5. Experimental results of the PCL array under different perturbation conditions. (a) and (b) are for c0=0, δp=0 and L0=1 with (a) z=0.57  m, (b) z=1.14  m; (c) and (d) are for the case of independent perturbation with c0=3w and L0=1 at (c) z=0.57  m, (d) z=1.14  m; (e) and (f) are for the case of uniform perturbation with c0=3w and L0=1 at (e) z=0.57  m, (f) z=1.14  m; and (g) and (h) are for the independent perturbation with c0=3w and L0=0 at (g) z=0.57  m, (h) z=1.14  m. Other parameters are w=0.157  mm and r0=1.5w. The horizontal and vertical axes representing the x and y axes are scaled in pixels with 1  pixel=5.2  μm. The camera is a CinCam CMOS-1201 (1280×1024 at 19.8 fps).

下载图片 查看所有图片

One may notice that the averaged intensities in Figs. 5(c) and 5(d) are not symmetric. The reason is that at the far field the irradiance distribution of a partially coherent light beam is a result of its correlation function at the source plane. The total correlation function of the radial array is a superposition of the correlations of the beamlets. High-ordered interferences among the beamlets of the array lead to the irradiance distribution of the combined field at the far field. When the number of the beamlets is as small as three, the interfering patterns are not symmetrically distributed and as a result Figs. 5(c) and 5(d) present asymmetric distributions of the averaged intensity. Increasing the number of beamlets at the source, the symmetry of the generated dark hollow structure would be improved.

The findings based on Fig. 5 demonstrate that the PCL array designed according to the strategy of Eqs. (12) and (13) can produce a dark hollow beam in the far field despite the poor coherence and the small number of beamlets at the source plane. We then further confirm that the dark hollow structure generated by this random array possesses a CV.

The dark rings inside the contour graph of the CCF are known to indicate the existence of an OAM in a partially coherent vortex beam [25]. For a beam with good coherence, the existence of the vortex can also be revealed by the ring dislocations in the Fourier transform of the intensity [49]. Figures 6(a)6(c) present the computational realizations of the second-order correlation G(2)(x,y,x,y), the magnitude of the CCF |W(x,y,x,y)|, and the Fourier transform of the averaged intensity of the random array in Fig. 5(d). A dark ring is observed in the distributions of the second-order correlation, the CCF, and the Fourier transformed intensity. Thus, the results of Figs. 6(a)6(c) indicate a CV embedded in the dark hollow beam produced by the special random array designed by Eqs. (12) and (13). The topological charge is 1. Notably, the computational quantities of G(2)(x,y,x,y) and |W(x,y,x,y)| are obtained based on the HBT relation and the instantaneous intensities obtained experimentally [Eq. (14)]. To demonstrate the correctness of the computational realizations, we perform wavefront folding interferometry [setup not shown in Fig. 3(a)] in which two Dove prisms are used to produce the correlation between the fields that are diametrically opposed [48]. A Mach–Zehnder interferometer is inserted into the path between P1 and the detector. The experimental result of the interferometry is shown in Fig. 6(d). A ring dislocation is embedded in the interference fringes. This result agrees with that of the computational realizations in Fig. 6(a).

Fig. 6. Computational realizations of (a) the second-order correlation G(2)(x,y,x,y) of the random array; (b) the magnitude of the CCF, |W(x,y,x,y)|; and (c) the Fourier transform of the averaged intensity of the array. (d) Experimental realization of the second-order correlation obtained with a Mach–Zehnder interferometer with two Dove prisms. The parameters of the array are the same as those in Fig. 5(d). In (a) and (b), the shooting camera is a TUCSEN ISH500 (600×800 at 49.6 fps), and 1 pixel is 2.2 μm; while in (d), the shooting camera is a CinCam CMOS-1201 (1280×1024 at 19.8 fps), and 1 pixel is 5.2 μm. (c) is in spatial frequency domain where 1 pixel is 1/(5.2×1280)  μm1.

下载图片 查看所有图片

The finding in Fig. 6 confirms the existence of a CV in the light field generated by the random array. Two approaches are adopted for the justification. One approach relies on measuring the instantaneous intensity of the structured light and then reconstructing the second-order correlation function. The magnitude of the CCF can be computed according to the HBT relation and Eq. (14). The result for the second-order correlation is shown in Fig. 6(a). The second approach is to measure the interference pattern in a Mach–Zehnder interferometer with two Dove prisms inverting the beams with respect to the horizontal and vertical axes, respectively. The output intensity of the interferometer is the second-order correlation function between (x, y) and (x,y), as shown in Fig. 6(d). From the results of both approaches, i.e., Figs. 6(a) and 6(d), one can observe a dark ring embedded in the interference fringes. It means that a CV exists inside the structured light field.

It is worth noting that there are other situations that spatially incoherent sources are used to create vortices [50,51] where people utilized rotating ground-glass disk (RGGD) to destroy the coherence of the light source and produce special vortex speckles. The difference between the RGGD method and the one proposed here is that in this paper the coherence of the complex light field is designed in a digital way. Correlation modulations can be carried out in the azimuthal direction. Moreover, fewer optical elements are needed in the experiment setup.

Considering the limitations of this method, the following issues are listed. (1) The quality of the structured light produced with the method here is limited by the spatial resolution of the digital modulator, i.e., SLM; and (2) the temporal evolution of the light source is controlled by the CGHs, and hence the duration of the light source relies on the amount of the calculation carried out in advance.

4 CONCLUSION

We proposed an experimental strategy to modulate the correlation structure of light fields in a digital and controllable way. By utilizing PCL arrays, we provided a routine to control the correlation structure according to phase configuration. A vortex-beam structure that cannot be produced using conventional laser arrays nor uniformly perturbed PCL arrays was evidenced through experiments by digitally modulating the randomness of the array in a special way. The CCF and the Fourier transformed contours of the produced beam revealed the existence of a CV embedded therein. Although we demonstrated this method with the use of an SLM in the experiment, it is not limited to this device. Considering the emergence of digital coding metamaterials [52] that also manipulate light fields in a programmable way, more optional routines for experimental realization may be available.

The method in this study is a general one and can be useful in other studies on optical manipulation and communication where the correlation structure plays a key role.

References

[1] D. M. Palacios, I. D. Maleev, A. S. Marathay, G. A. Swartzlander. Spatial correlation singularity of a vortex field. Phys. Rev. Lett., 2004, 92: 143905.

[2] W. Wang, M. Takeda. Coherence current, coherence vortex, and the conservation law of coherence. Phys. Rev. Lett., 2006, 96: 223904.

[3] W. Wang, Z. Duan, S. G. Hanson, Y. Miyamoto, M. Takeda. Experimental study of coherence vortices: local properties of phase singularities in a spatial coherence function. Phys. Rev. Lett., 2006, 96: 073902.

[4] O. S. Magaña-Loaiza, M. Mirhosseini, R. M. Cross, S. M. H. Rafsanjani, R. W. Boyd. Hanbury Brown and Twiss interferometry with twisted light. Sci. Adv., 2016, 2: e1501143.

[5] M. P. J. Lavery, C. Peuntinger, K. Günthner, P. Banzer, D. Elser, R. W. Boyd, M. J. Padgett, C. Marquardt, G. Leuchs. Free-space propagation of high-dimensional structured optical fields in an urban environment. Sci. Adv., 2017, 3: e1700552.

[6] V. Kumar, B. Piccirillo, S. G. Reddy, R. P. Singh. Topological structures in vector speckle fields. Opt. Lett., 2017, 42: 466-469.

[7] Z. Yang, O. S. Magaña-loaiza, M. Mirhosseini, Y. Zhou, B. Gao, L. Gao, S. Mohammad, H. Rafsanjani, G.-l. Long, R. W. Boyd. Digital spiral object identification using random light. Light Sci. Appl., 2017, 6: e17013.

[8] A. B. Ortega, S. Bucio-Pacheco, S. Lopez-Huidobro, L. Perez-Garcia, F. J. Poveda-Cuevas, J. A. Seman, A. V. Arzola, K. Volke-Sepúlveda. Creation of optical speckle by randomizing a vortex-lattice. Opt. Express, 2019, 27: 4105-4115.

[9] F. Gori, M. Santarsiero, R. Borghi, S. Vicalvi. Partially coherent sources with helicoidal modes. J. Mod. Opt., 1997, 45: 539-554.

[10] T. D. Visser, G. Gbur, E. Wolf. Effect of the state of coherence on the three-dimensional spectral intensity distribution near focus. Opt. Commun., 2002, 213: 13-19.

[11] G. Gbur, T. D. Visser, E. Wolf. ‘Hidden’ singularities in partially coherent wavefields. J. Opt. A, 2004, 6: S239-S242.

[12] G. Gbur, T. D. Visser. Coherence vortices in partially coherent beams. Opt. Commun., 2003, 222: 117-125.

[13] J. Yu, X. Zhu, F. Wang, D. Wei, G. Gbur, Y. Cai. Experimental study of reducing beam wander by modulating the coherence structure of structured light beams. Opt. Lett., 2019, 44: 4371-4374.

[14] X. Liu, X. Peng, L. Liu, G. Wu, C. Zhao, F. Wang, Y. Cai. Self-reconstruction of the degree of coherence of a partially coherent vortex beam obstructed by an opaque obstacle. Appl. Phys. Lett., 2017, 110: 181104.

[15] G. Kulkarni, R. Sahu, O. S. Magaña-loaiza, R. W. Boyd, A. K. Jha. Single-shot measurement of the orbital-angular-momentum spectrum of light. Nat. Commun., 2017, 8: 1.

[16] F. Wang, X. Liu, Y. Cai. Propagation of partially coherent beam in turbulent atmosphere: a review. Prog. Electromagn. Res., 2015, 150: 123-143.

[17] C. Zhao, Y. Cai, X. Lu, H. T. Eyyuboglu. Radiation force of coherent and partially coherent flat-topped beams on a Rayleigh particle. Opt. Express, 2009, 17: 1753-1765.

[18] J. M. Aunon, M. Nieto-Vesperinas. Partially coherent fluctuating sources that produce the same optical force as a laser beam. Opt. Lett., 2013, 38: 2869-2872.

[19] C. R. Alves, A. J. Jesus-Silva, E. J. S. Fonseca. Robustness of a coherence vortex. Appl. Opt., 2016, 55: 7544-7549.

[20] N. Montaut, O. S. Magaña-Loaiza, T. J. Bartley, V. B. Verma, S. W. Nam, R. P. Mirin, C. Silberhorn, T. Gerrits. Compressive characterization of telecom photon pairs in the spatial and spectral degrees of freedom. Optica, 2018, 5: 1418-1423.

[21] GburG.VisserT. D., “The structure of partially coherent fields,” in Progress in Optics (Elsevier B.V., 2010), Vol. 55, pp. 285341.

[22] Y. Altmann, S. McLaughlin, M. J. Padgett, V. K. Goyal, A. O. Hero, D. Faccio. Quantum-inspired computational imaging. Science, 2018, 361: eaat2298.

[23] Y. Shao, X. Lu, S. Konijnenberg, C. Zhao, Y. Cai, H. P. Urbach. Spatial coherence measurement and partially coherent diffractive imaging using self-referencing holography. Opt. Express, 2018, 26: 4479-4490.

[24] X. Lu, C. Zhao, Y. Shao, J. Zeng, S. Konijnenberg, X. Zhu, S. Popov, H. P. Urbach, Y. Cai. Phase detection of coherence singularities and determination of the topological charge of a partially coherent vortex beam. Appl. Phys. Lett., 2019, 114: 201106.

[25] R. Liu, F. Wang, D. Chen, Y. Wang, Y. Zhou, H. Gao, P. Zhang, F. Li. Measuring mode indices of a partially coherent vortex beam with Hanbury Brown and Twiss type experiment. Appl. Phys. Lett., 2016, 108: 051107.

[26] N. Zhao, X. Li, G. Li, J. M. Kahn. Capacity limits of spatially multiplexed free-space communication. Nat. Photonics, 2015, 9: 822-826.

[27] J. Wang. Advances in communications using optical vortices. Photon. Res., 2016, 4: B14-B28.

[28] J. Chen, Y. Li. Discrimination of incoherent vortex states of light. Opt. Lett., 2018, 43: 5595-5598.

[29] A. S. Desyatnikov, Y. S. Kivshar. Rotating optical soliton clusters. Phys. Rev. Lett., 2002, 88: 053901.

[30] L.-G. Wang, L.-Q. Wang, S.-Y. Zhu. Formation of optical vortices using coherent laser beam arrays. Opt. Commun., 2009, 282: 1088-1094.

[31] R. K. Singh, A. M. Sharma, P. Senthilkumaran. Vortex array embedded in a partially coherent beam. Opt. Lett., 2015, 40: 2751-2754.

[32] Z. Mei, O. Korotkova. Sources for random arrays with structured complex degree of coherence. Opt. Lett., 2018, 43: 2676-2679.

[33] J. Zeng, X. Liu, F. Wang, C. Zhao, Y. Cai. Partially coherent fractional vortex beam. Opt. Express, 2018, 26: 26830-26844.

[34] G. Piquero, M. Santarsiero, R. Martínez-Herrero, J. C. G. de Sande, M. Alonzo, F. Gori. Partially coherent sources with radial coherence. Opt. Lett., 2018, 43: 2376-2379.

[35] L. Rego, K. M. Dorney, N. J. Brooks, Q. L. Nguyen, C.-T. Liao, J. San Román, D. E. Couch, A. Liu, E. Pisanty, M. Lewenstein, L. Plaja, H. C. Kapteyn, M. M. Murnane, C. Hernández-García. Generation of extreme-ultraviolet beams with time-varying orbital angular momentum. Science, 2019, 364: eaaw9486.

[36] A. Lizana, H. Zhang, A. Turpin, A. V. Eeckhout, F. A. Torres-Ruiz, A. Vargas, C. Ramirez, F. Pi, J. Campos. Generation of reconfigurable optical traps for microparticles spatial manipulation through dynamic split lens inspired light structures. Sci. Rep., 2018, 8: 11263.

[37] J. Arlt, M. J. Padgett. Generation of a beam with a dark focus surrounded by regions of higher intensity: the optical bottle beam. Opt. Lett., 2000, 25: 191-193.

[38] Y. D. Liu, C. Q. Gao, M. W. Gao, F. Li. Coherent-mode representation and orbital angular momentum spectrum of partially coherent beam. Opt. Commun., 2008, 281: 1968-1975.

[39] S. Fu, C. Gao. Influences of atmospheric turbulence effects on the orbital angular momentum spectra of vortex beams. Photon. Res., 2016, 4: B1-B4.

[40] H. Wang, L. Liu, C. Zhou, J. Xu, M. Zhang, S. Teng, Y. Cai. Vortex beam generation with variable topological charge based on a spiral slit. Nanophotonics, 2019, 8: 317-324.

[41] Y. Yang, X. Zhu, J. Zeng, X. Lu, C. Zhao, Y. Cai. Anomalous Bessel vortex beam: modulating orbital angular momentum with propagation. Nanophotonics, 2018, 7: 677-682.

[42] MandelL.WolfE., Optical Coherence and Quantum Optics (Cambridge University, 1995).

[43] Q. Lin, Y. Cai. Tensor ABCD law for partially coherent twisted anisotropic Gaussian-Schell model beams. Opt. Lett., 2002, 27: 216-218.

[44] L. Chen, J. Lei, J. Romero. Quantum digital spiral imaging. Light Sci. Appl., 2014, 3: e153.

[45] Y. Yuan, D. Liu, Z. Zhou, H. Xu, J. Qu, Y. Cai. Optimization of the probability of orbital angular momentum for Laguerre-Gaussian beam in Kolmogorov and non-Kolmogorov turbulence. Opt. Express, 2018, 26: 21861-21871.

[46] GradshteynI. S.RyzhikI. M., Table of Integrals, Series, and Products (Elsevier/Academic, 2007).

[47] Y. T. Zhang, C. L. Ding, L. Z. Pan, Y. J. Cai. Laser arrays of partially coherent beams with multi-Gaussian correlation function. J. Quant. Spectrosc. Radiat. Transfer, 2018, 218: 1-11.

[48] B. Perez-Garcia, A. Yepiz, R. I. Hernandez-Aranda, A. Forbes, G. A. Swartzlander. Digital generation of partially coherent vortex beams. Opt. Lett., 2016, 41: 3471-3474.

[49] S. Prabhakar, A. Kumar, J. Banerji, R. P. Singh. Revealing the order of a vortex through its intensity record. Opt. Lett., 2011, 36: 4398-4400.

[50] F. Wang, S. Zhu, Y. Cai. Experimental study of the focusing properties of a Gaussian Schell-model vortex beam. Opt. Lett., 2011, 36: 3281-3283.

[51] Y. Cai, Y. Chen, F. Wang. Generation and propagation of partially coherent beams with nonconventional correlation functions: a review invited. J. Opt. Soc. Am. A, 2014, 31: 2083-2096.

[52] X. Wan, Q. Zhang, T. Yi Chen, L. Zhang, W. Xu, H. Huang, C. Kun Xiao, Q. Xiao, T. J. Cui. Multichannel direct transmissions of near-field information. Light Sci. Appl., 2019, 8: 60.

Min-Jie Liu, Jun Chen, Yang Zhang, Yan Shi, Chun-Liu Zhao, Shang-Zhong Jin. Generation of coherence vortex by modulating the correlation structure of random lights[J]. Photonics Research, 2019, 7(12): 12001485.

本文已被 1 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!