1. Introduction
In the last two decades, there has been a great interest in studying the natural singularities of scalar optics, namely, optical vortices. The vortices are rotational phase structures, with a singularity at center where the phase is undefined. Interest in such structures has become a separate topic in modern optics called Singular Optics [
1M. S. Soskin and M. Vasnetsov, “Singular optics,” Prog. Opt. 42, 219–276 (2001).
]. The concept has been extended to vector optics to study the polarization singularities, where the geometric polarization parameters of coherent fields cannot be defined [
2J. F. Nye, Natural Focusing and Fine Structure of Light: Caustics and Wave Dislocations (IoP Publishing, 1999).
]. The two types of polarization singularities are the C-points where the major axis of the polarization ellipse is not defined and the L-lines where the polarization handedness is not defined. These structures have been observed in the transverse plane of different fields [
3M. R. Dennis, “Polarization singularities in paraxial vector fields: morphology and statistics,” Opt. Commun. 213(4-6), 201–221 (2002).
–
5M. V. Berry, M. R. Dennis, and R. L. Lee, Jr., “Polarization singularities in the clear sky,” New J. Phys. 6, 162 (2004).
]. The polarization singularities have significant importance in understanding the fine details of beams with complex polarization structure.
Recently, there has been a tremendous interest in studying the connection between the singularities of scalar optics and vector optics using birefringent crystals [
6A. Desyatnikov, T. A. Fadeyeva, V. G. Shvedov, Y. V. Izdebskaya, A. V. Volyar, E. Brasselet, D. N. Neshev, W. Krolikowski, and Y. S. Kivshar, “Spatially engineered polarization states and optical vortices in uniaxial crystals,” Opt. Express 18(10), 10848–10863 (2010).
–
8F. Flossmann, U. T. Schwarz, M. Maier, and M. R. Dennis, “Stokes parameters in the unfolding of an optical vortex through a birefringent crystal,” Opt. Express 14(23), 11402–11411 (2006).
]. A detailed study of the spatial variation of polarization due to a circularly polarized input beam with Gaussian and Laguerre-Gaussian envelopes propagating through a birefringent crystal along and at an angle to its optic axis is presented in Ref [
6A. Desyatnikov, T. A. Fadeyeva, V. G. Shvedov, Y. V. Izdebskaya, A. V. Volyar, E. Brasselet, D. N. Neshev, W. Krolikowski, and Y. S. Kivshar, “Spatially engineered polarization states and optical vortices in uniaxial crystals,” Opt. Express 18(10), 10848–10863 (2010).
]. The formation of a complicated network of polarization singularities at the output of a birefringent crystal due to unfolding of a diagonally polarized optical vortex (OV) beam have been studied [
7F. Flossmann, U. T. Schwarz, M. Maier, and M. R. Dennis, “Polarization singularities from unfolding an optical vortex through a birefringent crystal,” Phys. Rev. Lett. 95(25), 253901 (2005).
,
8F. Flossmann, U. T. Schwarz, M. Maier, and M. R. Dennis, “Stokes parameters in the unfolding of an optical vortex through a birefringent crystal,” Opt. Express 14(23), 11402–11411 (2006).
] to show the topological interplay between scalar and vector optics. However, in these previous works, the polarization structures were studied at the output of a finite-length crystal where centers of the two beams are well separated. The spatial polarization structure near the unfolding point where the beam initially enters the crystal and splits into two eigen-beams has not been studied so far. In order to investigate the spatial polarization structure near the unfolding point one needs to attain the unfolding point at the output of the crystal.
In the present paper, we propose a novel setup to reconstruct the unfolding point at the crystal output and study the polarization structure near the unfolding point of a uniformly polarized OV beam. Using this setup, we obtain for the first time to our knowledge the polarization structure very near the unfolding/refolding point in a birefringent crystal. We also demonstrate that the separation between the two eigen-beams can be controlled, allowing observation of the beam structure further from the unfolding point.
At the input surface of a birefringent crystal, an initially diagonally polarized OV beam is split into two orthogonally polarized eigen-beams, namely, ordinary ray (O-ray) and extra-ordinary ray (E-ray). The point where this occurs is called the unfolding point. The resulting two eigen-beams are separated in space at the output of the crystal. In our setup a half wave plate (HWP) interchanges the polarization of the two beams, and propagation through a second identical crystal reconstructs the unfolding point at the output, folding back the two eigen-beams into a single beam. Such an arrangement will result in an interferometer of Mach-Zehnder type. The beam at the output of this optical vortex birefringent interferometer can be unfolded by rotating the second crystal with respect to the first crystal about its vertical axis. Small rotation of the second crystal produces a varying separation of the two eigen-beams at the output of the interferometer. The polarization of the beam at the output of the interferometer is characterized by measuring the Stokes parameters in the conventional way.
In previous works, the usefulness of the birefringent interferometer with Gaussian input beam is demonstrated to measure the surface deformations due to photo-thermal effects [
9X. D. Xu, P. K. Kuo, S. Y. Zhang, X. J. Shui, and Z. N. Zhang, “Application of an optical birefringence interferometer to photothermal detection,” Microw. Opt. Technol. Lett. 35(2), 140–143 (2002).
] and for signal processing techniques [
10C. Cheng, “The signal processing approach for the birefringent material based Mach-Zehnder interferometer design,” Proc. of IEEE, 48th Midwest Symposium on Circuits and Systems (Covington, Kentucky, 2005), 211–214, 10.1109/MWSCAS.2005.1594076. [CrossRef]
]. The birefringent interferometer has also been used as a stable interferometer for quantum information processing [
11J. L. O’Brien, G. J. Pryde, A. G. White, T. C. Ralph, and D. Branning, “Demonstration of an all-optical quantum controlled-NOT gate,” Nature 426(6964), 264–267 (2003).
]. In the present paper the birefringent Mach-Zehnder interferometer using OV beam is demonstrated and characterized for the first time to show its applicability in reconstruction of unfolding points and to study the spatially varying polarization structure near the unfolding point. Using the proposed configuration, the polarization structure of the OV beam very near and far from the unfolding point is created. It can also be observed from the results that it is easy to generate and control spatially varying polarization using the designed setup. The present study can also find potential applications in spectroscopy and microscopy dealing with spatially varying polarization structures.
The paper is organized as follows. Basic theoretical model regarding the present work is given in section 2. Section 3 describes the experimental setup and the results obtained using the stable optical vortex birefringent Mach-Zehnder interferometer (BR-MZI). The conclusions of the present work are given in section 4.
3. Experimental details and results
The experimental setup to demonstrate the proposed optical vortex birefringent interferometer is shown in
Fig. 1. The polarization of the laser beam from a He-Ne laser is set to 45 degrees using a polarizer (P
1). The diagonally polarized laser beam then enters a transmission type fork structured phase hologram (HG) made by electron beam lithography (EBL) [
14Y. Miyamoto, M. Masuda, A. Wada, and M. Takeda, “Electron-beam lithography fabrication of phase holograms to generate Laguerre-Gaussian beams,” Proc. SPIE 3740, 232–235 (1999).
]. The Laguerre-Gaussian (LG) beam with a topological charge of + 1 generated in the first order diffraction from the hologram is selected using an aperture (A). The size of the beam (2
w0) at the detector plane is estimated to be 1.37 mm. The singly charged and diagonally polarized LG beam is then allowed pass through a birefringent crystal (Cr
1) made of yttrium vanadate (YVO
4) with dimensions 30 mm (length) X 12 mm (width) X 10 mm (height) (from CASTECH, China) and its optic axis in the horizontal plane making an angle of 45 degrees to the propagation direction.
The ordinary and extraordinary refractive indices
no and
ne of the crystals are given by 1.9929 and 2.2154 respectively. The birefringence (Δ
n) of the crystals and the walk-off angle (
ρ) between the two eigen-beams are calculated to be 0.2225 and 6.04 degrees respectively at 630 nm [
15]. The two eigen-beams with orthogonal polarizations separated at the output of the first crystal due to birefringence are passed through a half wave plate (HWP) to interchange the polarizations of the two beams. The two beams after the HWP plate then enter a second crystal (Cr
2) identical to the first one with its optic axis parallel to the first crystal and mounted on a rotational stage. The two eigen-beams are folded back together to form a single beam at the output of Cr
2. Such an arrangement will result in a stable interferometer of Mach-Zehnder type.
The unfolding point occurring at the input surface of the first crystal is reconstructed at the output of the second crystal. Small rotation of the second crystal will result in a small separation of the two eigen-beams at the output of the interferometer. Similar rotation techniques have been previously used to control the relative phase between the two eigen-beams with a very small rotation angle [
7F. Flossmann, U. T. Schwarz, M. Maier, and M. R. Dennis, “Polarization singularities from unfolding an optical vortex through a birefringent crystal,” Phys. Rev. Lett. 95(25), 253901 (2005).
,
8F. Flossmann, U. T. Schwarz, M. Maier, and M. R. Dennis, “Stokes parameters in the unfolding of an optical vortex through a birefringent crystal,” Opt. Express 14(23), 11402–11411 (2006).
], and also for a different crystal orientation [
16T. Fadeyeva, Y. Egorov, A. Rubass, G. A. Swartzlander, Jr., and A. Volyar, “Indistinguishability limit for off-axis vortex beams in uniaxial crystals,” Opt. Lett. 32(21), 3116–3118 (2007).
]. The polarization structure of the resultant beam at the output of the interferometer is characterized by measuring the Stokes parameters using a quarter wave plate (QWP) and polarizer (P
2) as shown in
Fig. 1. The intensity distribution for different combinations of QWP and P
2 are recorded using a CCD (Model No. C5948, Hamamatsu) connected to a personal computer. As the BR-MZI is mounted on a translation stage it can be moved in and out of the input beam. Initially, the BR-MZI is moved out from the setup and we characterized the input beam. The polarization of the input OV beam is characterized by measuring the Stokes parameters using [
17T. Kihara, “Measurement method of Stokes parameters using a quarter-wave plate with phase difference errors,” Appl. Opt. 50(17), 2582–2587 (2011).
]
where
I (
β,γ) is the output intensity when the axes of the QWP and P
2 are at
β and
γ respectively as measured from horizontal direction.
Figure 2
shows the normalized stokes parameters of the input OV beam. S
0 is normalized to its maximum value and S
1, S
2, S
3 are normalized to S
0 such that
s1 = S
1/S
0,
s2 = S
2/S
0 and
s3 = S
3/S
0. The CCD images shown in the present paper are inverted in the vertical direction as indicated by the arrow showing the + x direction in the figures. It can be seen from
Fig. 2(a),
2(b),
2(c) and
2(d) that the polarization of the input beam is 45 degrees and the contribution of the horizontal/vertical and circular polarizations are minimal. The two angles 2α = arg (
s1 + is2) and δ = arg (
s2 + is3) calculated from the Stokes parameters are shown in
Fig. 2(e) and
2(f) respectively. The angle α is the azimuthal angle of the major axis of the polarization ellipse and δ is the phase difference between the horizontal and vertical components.
Figure 2(e) and
2(f) confirm that the major axis of the polarization ellipse is at 45 degrees and there is no phase difference between the two orthogonal components.
Fig. 2 Experimental results; (a) Normalized total intensity (S0); (b), (c) and (d) represents s1, s2 and s3 ; (e) and (f) represent 2α and δ respectively of the input beam.
Now the BR-MZI is introduced into the beam and the output beam is characterized. The crystal Cr
2 is rotated about its vertical axis (x-axis) to a maximum angle on one side where we see the two eigen-beams separated well in horizontal direction (y-axis) and then it is rotated to the other side in steps of one degree. At each rotation angle the output beam is captured through the CCD and we estimated the position of the O-ray and E-ray separately by calculating the center of gravity of the intensity distribution of each beam. The positions of the O-ray and E-ray at each rotation angle are plotted in
Fig. 3(a)
.
Fig. 3 (a) Position of O-ray and E-ray on CCD (b) Separation d between the two eigen-beams as a function of rotation angle of Cr2.
The relative separation
dr of the two eigen-beams with respect to the input beam size is defined as
Here
d is the separation between the two eigen-beams and 2
w0 is 1/e
2 beam width of the input beam at the detector plane.
Figure 3(b) shows the separation
d of the two eigen-beams as a function of rotation angle of the Cr
2. As the length and birefringence of the crystal (YVO
4) used in our experiment is large compared to that of reference [
7F. Flossmann, U. T. Schwarz, M. Maier, and M. R. Dennis, “Polarization singularities from unfolding an optical vortex through a birefringent crystal,” Phys. Rev. Lett. 95(25), 253901 (2005).
], we get relative separation of about 220% between the two eigen-beams at the output of a single crystal. This larger separation per crystal and larger angles of rotation of the crystal Cr
2 compared to reference [
7F. Flossmann, U. T. Schwarz, M. Maier, and M. R. Dennis, “Polarization singularities from unfolding an optical vortex through a birefringent crystal,” Phys. Rev. Lett. 95(25), 253901 (2005).
] allows us to control the separation over a wide range of values.
It can be seen from
Fig. 3(b) that the separation decreases linearly with rotation angle of Cr
2. The zero value of the separation indicates the unfolding point and it is observed from
Fig. 3(b) that this occurs between 346 and 347 degrees as read from meter of the rotational stage. The separation
d between the O-ray and E-ray is found to be around 30.2 μm and −29 μm at 346 degrees and 347 degrees respectively. The relative separation,
dr corresponding to 346 degrees and 347 degrees is calculated using
Eq. (5) to be 2.2% and −2.1% respectively at the output of the second crystal. The negative sign represents the interchange of the positions of the two-eigen beams after crossing the unfolding point. The dotted line in
Fig. 3(b) represents the linear fit of the data. From the linear fit it is found that the separation changes −32.28
0.82 μm per unit angle rotation of Cr
2. We note that the exact angle for the unfolding point depends on the initial separation at the input of the second crystal, which in turn depends on the orientation of the first crystal.
Now the Cr
2 is adjusted to the angle near the crossing point where the relative separation is measured to be 1.7% in the y-direction, which is very close to the unfolding point. The polarization of the output beam is characterized with the same method as for the input beam. The normalized Stokes parameters measured at the output of the BR-MZI are shown in
Fig. 4
. It can be seen from
Fig. 4(c) that in the upper half of the beam as displayed the polarization is 45 degrees, which is identical to the input beam polarization. However, the lower half of the beam is closer to right circular polarization (
Fig. 4(d)). This variation can also be seen in
Fig. 4(f) as a slope in δ, indicating relative tilt in the wave front of the two eigen-beams.
Figure 4(e) shows that 2α is more or less constant, indicating that the major axis of the polarization ellipse remains near 45 degrees throughout the beam.
Figure 4(b) shows that
s1 is relatively small and varies in the direction from top left to bottom right. Since
s1 is the balance between the intensities of the two eigen-beams, this variation is expected to coincide with the direction of the separation between the two eigen-beams, inferring that the two eigen-beams are separated in the x (vertical) direction as well as the y (horizontal) direction.
Fig. 4 Experimental results; (a) Normalized total intensity (S0); (b), (c) and (d) represents s1, s2 and s3; (e) and (f) represent 2α and δ respectively of the beam at the output of the interferometer for the relative separation of 1.7% in y-direction
In order to verify the effects of the tilt and the separation in x and y directions, we consider the following model which incorporate these details as well as the overall phase difference between the two eigen-beams. The tilt of the wave front will result in a phase variation in the direction of the relative wave vector of the two beams. To incorporate the wave vector changes due to the tilt we modify
Eq. (2) to get the general form as follows:
Here
d1 and
d2 represent the separation between the two eigen-beams in the x and y directions respectively.
θ and
ϕ represent the spherical coordinates corresponding to the tilt angle and Δ is the overall phase difference between the two-eigen beams. Using
Eq. (3) and
Eq. (6) the normalized Stokes parameters are calculated and matched with experimental results. From the calculations,
d1 and
d2 are found to give a relative separation of 2.5% and 1.7% respectively, and
θ and
ϕ are estimated to be 35x10
−6π and 0.1π respectively. The overall phase difference Δ between the two eigen-beams is found to be equal to 0.0012π. We note that the tilt can be a combined effect of a curved wave front and the separation between the two eigen-beams. The calculated results using
Eq. (3) and
Eq. (6) are shown in
Fig. 5
. It can be observed that the calculated results coincide with the experimental results.
Fig. 5 Simulation results using
Eq. (3) and
Eq. (6); (a) Normalized total intensity (S
0); (b), (c) and (d) represents
s1,
s2 and
s3; (e) and (f) represent 2α and δ respectively of the beam at the output of the interferometer for the relative separation of 1.7% in y-direction
The spatial polarization structure for a large rotation angle of Cr
2 which results in a large separation between the two eigen-beams is also investigated. The crystal Cr
2 is rotated such that the rotation results in the relative separation of −42% in the y-direction between the eigen-beams at the output of the interferometer. The spatial distributions of the measured normalized Stokes parameters are shown in
Fig. 6
. It can be observed that the large separation of the eigen-beams results in a non-uniform polarization across the beam cross section.
Fig. 6 Experimental results; (a) Normalized total intensity (S0); (b), (c) and (d) represents s1, s2 and s3; (e) and (f) represent 2α and δ respectively of the beam at the output of the interferometer for the −42% separation in y-direction
These polarization characteristics are identical to the results presented in the earlier research [
18U. T. Schwarz, F. Flossmann, and M. R. Dennis, “Topology of generic polarization singularities in birefringent crystals,” Topologica 2, 006 (2009), [CrossRef]
] for the large separation between the two eigen-beams using a single birefringent crystal. The experimental results are shown in
Fig. 6.
Figure 6(d) shows a horizontal line at center where S
3 = 0, indicating linear polarization (L-line). In
Fig. 6(e) two singular points can be observed which correspond to the positions where the major axis of the polarization ellipse is undefined. Such a situation will occur only for circular polarizations (C-points). The rotations of the major axis around the two C-points in
Fig. 6(e) have opposite handedness.
The experimental results are compared with the computer simulations for the relative separation of −2.5% in the x-direction and −42% in the y-direction and using
Eq. (3) and
Eq. (6). The computer calculations are shown in
Fig. 7
. It can be observed from
Fig. 6 and
Fig. 7 that the experimental results match well with the computer simulations, even though we have not included tilt and constant phase difference for these simulations. These results infer that for large relative separation the polarization structures are governed by the rotating phase of OVs and that the effect of the small wave front tilt and small overall phase difference is not as visible as for the case of small relative separation. The variation of
s1 is slightly tilted from the horizontal direction, showing that the effect of separation in the y-direction is dominant but there is also a relatively small contribution from the separation in the x-direction. The separation in the x-direction is similar in magnitude to the one observed in
Fig. 4 and
Fig. 5, but opposite in sign. These vertical separations may be caused by small imperfections and misalignment in the experimental setup.
Fig. 7 Computer simulations; (a) Normalized total intensity (S0); (b), (c) and (d) represents s1, s2 and s3; (e) and (f) represent 2α and δ respectively of the beam at the output of the interferometer for the −42% separation in y-direction
The present investigations show that the proposed birefringent interferometer in Mach-Zehnder configuration with an optical vortex input can be used to achieve the spatial separation of the eigen-beams very close to unfolding points and enable the study of the fine details of the spatial polarization structure very close to such points. The present setup can also be used to generate and control the spatial variation of polarization which results in polarization singularities in the output field. Further investigations on the sensitivity of the polarization structure to the relative separation between the eigen-beams and the evolution of the polarization singularities from uniform polarization are in progress.