Following our earlier work (F. Flossmann et al., Phys. Rev. Lett. 95 253901 (2005)), we describe the fine polarization structure of a beam containing optical vortices propagating through a birefringent crystal, both experimentally and theoretically.We emphasize here the zero surfaces of the Stokes parameters in three-dimensional space, two transverse dimensions and the third corresponding to optical path length in the crystal. We find that the complicated network of polarization singularities reported earlier – lines of circular polarization (C lines) and surfaces of linear polarization (L surfaces) – can be understood naturally in terms of the zeros of the Stokes parameters.
©2006 Optical Society of America
Singular optics is an important branch of modern optics. In scalar optics, this is the study of optical vortices (phase singularities, or wave dislocations) [1, 2] and in polarization optics, polarization singularities (C points and lines, L lines and surfaces) [3, 4]. These singularities are geometric organizing centers in spatially extended optical fields which naturally arise in the propagation, interference and diffraction of optical fields , and the 2- and 3-dimensional singularity geometry reveals features of the underlying optical wave physics.
Optical vortices are generalized interference fringes, and generically occur whenever three or more waves are superposed , as well as in speckle patterns  and in beams carrying orbital angular momentum, such as Laguerre-Gauss (LG) beams [8, 9]. They are places (points in two dimensions, lines in three) where the scalar amplitude is zero, and in a loop around which the phase changes by ±2πn, where n is the integer vortex strength (charge).
Polarization singularities are places in coherent elliptic polarization fields where geometric polarization parameters cannot be mathematically defined. They are important in understanding the fine structure of beams whose polarization is spatially inhomogeneous, and their presence can give rise to significant physical effects . The polarization singularities which propagate stably in paraxial beams are C lines of circular polarization, where the major axis of the polarization ellipse is not defined, and L surfaces of linear polarization, where polarization handedness is not defined. These appear as C points and L lines respectively in the transverse plane. Polarization singularities occur in fields evolving from waveguides , polarization-varying microgratings , and light propagating through anisotropic media [13, 14, 15].
In Ref. , a physical connection between the optical vortices of scalar optics and polarization singularities of vector optics was investigated. A linearly polarized L beam with an on-axis vortex propagated through a birefringent crystal, and the outgoing polarization structure was interpreted in three dimensions (two transverse and a third, the optical thickness of the crystal, effectively a phase shift). It was found that the crystal birefringence unfolds the input optical vortex into to a complicated network of polarization singularities. A pair of C lines, of opposite index and handedness, forms a double helix structure, intertwined around an L surface, which is topologically identical to the periodic minimal surface known as Riemann’s examples or Riemann’s minimal surface . This pattern of singularities arises from the breaking of symmetries of the input beam (uniform linear polarization, intensity cylindrically symmetric) due to the anisotropy of the crystal.
The standard description of optical polarization employs the Stokes parameters, and their geometric representation on the Poincaré sphere. Here, we describe the behavior of Stokes parameters in the unfolding of a linearly polarized optical vortex beam through a birefringent crystal, and find that features of the polarization singularity geometry reported in Ref.  have a natural interpretation in terms of the zero contour surfaces of the Stokes parameters. We also present new experimental results on the structure of polarization singularities stemming from the unfolding of Gaussian beams with no vortices and higher order vortices on propagating through the crystal.
2. Interpretation using Stokes Parameters
Any state of paraxial, elliptic polarization can be described completely using the Stokes parameters [18, 19]. The complex electric field vector E=(Ex ,Ey ) describes the polarization ellipse, parameterized by its overall size (intensity), angle of major axis α, and a shape parameter ω, defined such that |tanω| is the ratio of minor and major ellipse axes, and sign ω is the ellipse handedness (see Figure 1(a)). Up to a phase factor, the polarization vector is therefore
where I is the overall intensity |E|2, α and ω are the ellipse parameters, and β and δ are a pair of equivalent parameters , which prove to be more useful in our analysis. The unnormalized Stokes parameters are defined [18, 19]
The terms in the third column denote the intensity of the scalar component of a particular polarization determined with the appropriate analyzer. The third set is in terms of the ellipse parameters α,ω, and the third, the complementary parameters β,δ. Since =++, any state of polarization can be represented on the surface of a sphere, with coordinates Si /S 0 (i=1,2,3)- this is the Poincaré sphere (Figure 1). From Eq. (2), angles 2α and 2ω are the azimuth and colatitude angles on the Poincaré sphere with S 3 as the axis. δ and β may be understood as the corresponding half-angles for azimuth and latitude with respect to S 1 as the axis. Physically, S 0 is simply the field intensity I, and S 3/S 0=sin 2ω is the ellipticity, signed by handedness: positive for right-handed elliptical polarization, negative for left-handed. When S 3=0, the polarization is linear, so the zero contour surface of S 3 is the L surface [1, 4, 20]. The angle α=arg(S 1+iS 2)/2 is the angle of the polarization ellipse major axis, and is undefined when S1=S2=0; this is the condition for circular polarization (where ω=±π/2, |S 3|=S 0). The polarization singularity positions are therefore determined by the zero surfaces of the three Stokes parameters: L surfaces by S 3=0, and C lines by the intersection of the surfaces S 1=0 and S 2=0.
The polarization phenomena described here are the result of a spatially inhomogeneous, linearly polarized monochromatic optical beam propagating through a birefringent crystal, and depend on three birefringence effects: the incident polarization is resolved into the two characteristic linear polarizations corresponding to that propagation direction; each of these two rays propagates with a particular speed, resulting in a net phase shift between the two; beam walkoff results in a transverse spatial separation of the two rays. The incident propagation direction is taken not to be close to the optic axes of the crystal, so the effect of conical diffraction  may be ignored, as well as whether the crystal is uniaxial or biaxial. The details of our experiment, involving uniaxial KDP, are described in Ref. ; they are summarized in Fig. 2. The Gaussian beam (with waist width w 0) is incident in the plane containing the surface normal and the optic axis, implying that the spatial shift between the two beam components is in the direction of one of the crystal eigenpolarizations, and the difference in geometrical path length is negligible compared to the difference in optical path length. The phase shift therefore arises from the difference in refractive index, i.e. Λ=2π(n 2-n 1)l/λ vac, where l is the average path length of the two rays. We will comment in the final section on what happens for crystals with optical activity and absorption, but in our experiment, these effects are negligible.
Following the notation of Ref. , the uniformly polarized beam incident on the crystal is E in(x,y)=ψ(x,y)d in, where ψ(x,y) is the complex scalar amplitude of the incident field, and d in is a unit vector representing the initial polarization, written in terms of β in and δ in by appropriate modification of Eq. (1). In the crystal, the two polarizations which propagate in the crystal are d 1,d 2 (defined to be in the x,y-directions), the path length-dependent phase shift is Λ and transverse shifts (in the x-direction) ±s result in a total outgoing field
The Stokes parameters of the outgoing field are easily computed, where ψ ±=ψ(x±s,y) :
By inspection of these equations, the following observations can be made:
- S 0 and S 1 are independent of Λ (since the only Λ-dependence in Eq. (3) is in the phase of the two Cartesian components);
- At x=±s,y=0, S 1≠0, but S 2=S 3=0;
- The only effect of shifting δ in is effectively to shift Λ;
- S 2,S 3 are identical, except for a quarter-period Λ-shift.
The final two points can be understood by interpreting the effect of a quarter-period phase difference as the effect of an appropriately oriented quarter-wave plate, which does not change S 1, but transforms S 2 to S 3, and S 3 to -S 2. In Stokes space, this is illustrated by a rotation of the Poincaré sphere by π/2 about the S 1-axis, corresponding to a change of π/4 in δ. This is represented in Fig. 3, in which experimental plots are compared for d in oriented at 45° (Fig. 3(a)) and right handed circular polarized (Fig. 3(b)), where δ in differs π/4 and Λ by π/2. On changing Λ by π, the pattern changes in a way analogous to a half-wave plate (as observed in Ref. ), namely S 1 is unchanged, but S 2 and S 3 change sign.
The comments above apply for any input field ψ(x,y). In Ref. , this was chosen to be a Gaussian beam carrying a strength 1 on-axis optical vortex (i.e. a Laguerre-Gauss Lmode), so ψ 1=(x+iy)exp(-(x 2+y 2)/). The polarization singularities are determined by the zero loci of the Stokes parameters S 1,S 2,S 3. Using Eq. (4), the zeros occur at the following x,y,Λ- coordinates (where, for simplicity here and hereafter, δ in=0) :
(The equations for S 2,S 3=0 assume that sin 2β≠0, i.e. d in≠d 1 or d 2.) As expected, the S 1=0 surface is independent of Λ. The solution for x and y is transcendental, depending only on the input β in and the splitting s. On the other hand, the surfaces S 2=0 and S 3=0 are independent of β in, and are mathematically easily determined: each constant-Λ section is a circle, with radius |(co) secΛ| and center x=0,y=(co) tan Λ. Such a circle for S 3=0 is plotted as the L line in Fig. 3(c); the experimental results, although topologically the same, do not give circles, probably due to lack of perfect rotational symmetry of the input beam. Various representations of the experimentally-determined ellipse parameters and zero contours of Stokes parameters, for a fixed choice of Λ, are plotted in Fig. 4.
The experimentally-determined zero surfaces of the Stokes parameters, in x,y,Λ-space, are plotted in Fig. 5. The surface for S 1=0 consists of a tube and a nearby sheet. The periodic surfaces S 2,S 3=0 are identical, and topologically the same as Riemann’s periodic minimal surface . The two surfaces intersect at (x,y)=(±s,0), i.e. at the position of the shifted vortex in each of the two waves. Although topologically the same as Riemann’s example, using the expressions in Eq. (5), it can be shown that the surfaces are not themselves minimal (i.e. the mean curvature does not vanish everywhere ). The intersection of the translation-invariant S 1=0 surface with the surface S 2=0 gives rise to the C line geometry described in Ref. : a double helix of C lines of opposite index and handedness (lying on the S 1=0 tube) and a transversally infinite C line (lying on the S 1=0 sheet). If β in changes (this is the angle of the polarization if d in is linearly polarized), the S 2,S 3=0 surfaces do not change, but the S 1 surface does; in particular, when β in=45°, and the C lines reconnect, the S 1 surface and tube touch along two self-intersection lines.
As described in Ref. , the complicated polarization singularity geometry in the previous example may be explained as the unfolding of the input optical vortex of the scalar field. When there is no vortex input, the polarization pattern is somewhat simpler. For instance, with a simple Gaussian beam as input, with ψ 0=exp(-(x 2 +y 2)/), the equations for the zeros of the Stokes parameters are
(Similar assumptions are made here as with Eq. (5): δ in is assumed to be zero, and sin 2β≠0.) The Λ-invariant S 1 surface is just a single plane with constant x, the S 2 and S 3 surfaces are planes occurring at constant Λ. Thus the output polarization field has Λ-layers alternating in handedness, with a single C line between two adjacent L surfaces. Since the C line lies in the transverse x,y-plane, its index is not defined , although its handedness is. There is no singularity structure threading through the pattern, reflecting the fact that there was no scalar vortex in the input field. The experimentally-determined zeros of S 1,S 2,S 3 for an input Gaussian beam are plotted in Fig. 6.
The input field can be modified to a Laguerre-Gauss beam L with any topological charge n (assumed positive here for simplicity), i.e. ψn (x,y)=(x+iy) n exp(-(x 2+y 2)/). In this case, the zeros of the Stokes parameters S 1 …S 3 are given by:
This expression generalizes Eqs. (5) (when n=1), and, with appropriate limits, (6) (when n=0). The zero contour for S 1 is topologically the same as the n=1 case described above, and the factor 1/n just scales the beam waist w 0. However, there are differences in the zero contours of S 2 and S 3. Although the sections of constant Λ are again circles, for each constant L there are now n distinct solutions, resulting in exactly n circles with different radii and centers. The corresponding surfaces are the same as those in the n=1 case, but occur with an n-fold multiplicity whilst being shifted in the Λ-direction and intersecting each other in a highly nongeneric and unstable way.
For n=2, the Stokes parameter zero surfaces are plotted in Fig. 7, both experimentally (7(a)) and theoretically (7(b)). There are now six C lines, with four coiling around each other, with an overall pitch of 4π instead of 2π.
An alternative representation of the polarization uses the ratio of Cartesian components of E, equivalent to stereographic projection of the Poincaré sphere in the positive S 1-direction, rather than the S 3 direction, as is more usual [18, 13],
All information about the polarization state is expressed by ρ; in particular, the Stokes zero contour surfaces are specified by Reρ=0 (S 2=0), Imρ=0 (S 3=0), and |ρ|=1 (S 1=0), and ρ has a zero or pole at the positions of the refracted vortices of the input field.
When the incident field has a regular Gaussian profile (possibly with n axial vortices embedded), the Gaussian factor gives rise to a factor of exp(-4xs/w 2); for the field ψ 01, this is multiplied by (x+s+iy)n/(x-s+iy) n .
Several general observations can now be made, working with the Poincaré sphere of Stokes parameters, with the non-usual convention of S1 as the axis of polar coordinates (i.e. with polar angles 2β,2δ). S 1 has this privilege as the crystal eigenpolarizations have S 1=±S 0,S 2,S 3=0. The behavior of all incoming polarizations with the same β(=arctan|d in,y/d in,x|) on the Poincaré sphere is the same. Also, the surface locus of polarizations lying on any great circle including the S 1-direction is the same as the surfaces S 2=0,S 3=0, up to a phase shift (corresponding to a surface Reρ exp(-iχ) for some phase angle χ). Finally, the loci of any pair of orthogonal polarizations corresponding to opposite points on the great circle perpendicular to the S 1-direction will be the same as the C lines, again up to a phase shift.
The unfolding of an optical vortex on propagation through a crystal, as explained here and in Ref. , demonstrates the interplay between optical vortices, the singularities of scalar optics, and C lines and L surfaces, the singularities of paraxial polarization optics. The effects described occur due to birefringence, and the resulting resolution of polarization, beam walkoff, and phase shift of the two rays; although we have emphasized the situation for uniaxial crystals, the only difference for biaxial crystals is that neither eigenpolarization corresponds to the direction of ray separation (walk-off). Of course, crystal optics has its own singularities [13, 14], particularly when the incoming ray direction coincides with an optic axis of the crystal; the question of how a vortex-carrying beam incident on the optic axis of a crystal is diffracted is outside the scope of the present work. Our experiment can be thought of as the reverse experiment of that described in Ref. , in which an optical vortex beam is synthesized from a polarization singularity in the pattern generated from propagation through a crystal (near the optic axis).
The periodicity of the phase shift Λ between the two rays is reminiscent of propagationperiodicity elsewhere in singular optics. As observed in Ref. , the configuration of polarization singularities here resembles the braided vortex lines of Ref. . It is natural to ask whether, with an appropriate choice of input field, birefringence alone can produce braided C lines. However, the following argument shows that a pigtail braid of three C lines, as constructed for vortices in Ref. , is not possible. C lines occur at the intersection of the surfaces S 2=0 and S 1=0. Since the latter must be independent of Λ, a (pigtail) braid requires a transverse section of this surface to be a figure-8 . However, such a braid consists of three filaments, requiring the S 2=0 surface to cross the Λ-invariant S 1=0 surface three times for each Λ. This is clearly impossible, as contour surfaces intersect generically an even number of times, proving that braided C lines are impossible using Eq. (3).
A common example of perturbing an optical vortex is the superposition of a screw dislocation with a copropagating plane wave of small amplitude [1, 6], in which the straight vortex line becomes helical. The situation we describe is a generalization of this in two ways . Firstly, the superposition of Eq. (3) is a superpositions of waves with orthogonal polarizations, rather than (implicitly) the same polarization, so the final pattern involves polarization singularities rather than phase singularities. Secondly, there is a vortex in each of the fields being superposed. If two orthogonally polarized beams were superposed, one with a vortex, the other a plane wave (with a phase shift Λ), then the zero contours of the Stokes parameters are a cylindrical tube invariant with Λ (S 1=0), and a pair of orthogonal helicoids along the same axis (S 2,S 3=0): this results in a double helix of C lines with opposite sign and handedness. Although illustrative, this example cannot be physically realized in an experiment involving crystal birefringence, as both fields must be the same up to a spatial and phase shift.
A further generalization of the effect described here would be to extend to crystals which possess optical activity or absorption (dichroism) . In these cases, we do not anticipate the rich polarization topology of pure birefringence to occur. The dielectric tensor of an optically active crystal is complex Hermitian, and the orthogonal eigenpolarizations for a non-singular propagation direction are in general elliptic. The behavior of the zero contours for generalized Stokes parameters, based on an orthogonal triad of polarizations on the Poincaré sphere (with one direction corresponding to the eigenpolarizations), is the same as that described in Section 2 for birefringent crystals without optical activity; however, the polarization singularities (C lines, L surfaces) do not possess the geometry described in Section 2 because the usual Stokes parameter triad describing them has no relationship to this generalized triad based on elliptic polarizations. Dichroic crystals do not, in general, have orthogonal eigenpolarizations, so the decomposition (3) fails.
The discussion at the end of Section 2, as well as the previous paragraph, implies that there is nothing special about the polarization state on C lines and L surfaces. Of course, from the viewpoint of the Poincaré sphere this may be true; however, the approach of singular optics is to deduce properties of the entire optical field from the behavior of the singularities, possible due to their geometrically singular nature (the polarization at a C line is rotationally invariant, the handedness on an L line is not defined). The remarkable properties of Stokes parameters and the Poincaré sphere remove this specialness, but only in the abstract, three-dimensional Stokes space: the geometric fine structure of the spatially inhomogeneous polarization field we describe is completely physical.
UTS acknowledges support from the Japan Society for the Promotion of Science (JSPS), and MRD acknowledges support from the Royal Society of London.
References and links
1. J. F. Nye, Natural focusing and fine structure of light: caustics and wave dislocations, IoP Publishing, 1999.
2. M. S. Soskin and M. Vasnetsov, “Singular optics,” Prog. Opt. 42219–276 (2001). [CrossRef]
3. J.F. Nye, “Lines of circular polarization in electromagnetic wave fields,” Proc. R. Soc. Lond. A 389279–290 (1983). [CrossRef]
4. M.R. Dennis, “Polarization singularities in paraxial vector fields: morphology and statistics,” Opt. Commun. 213201–221 (2002). [CrossRef]
5. M.V. Berry, “Singularities in waves and rays,” in R. Balian, M. Kléman, and J.-P. Poirier, editors, Les Houches Session XXV - Physics of Defects (North-Holland, 1981).
6. K. O’Holleran, M.J. Padgett, and M.R. Dennis, “Topology of optical vortex lines formed by the interference of three, four and five plane waves,” Opt. Express 143039–3044 (2006). [CrossRef] [PubMed]
7. M.V. Berry and M.R. Dennis, “Phase singularities in isotropic random waves,” Proc. R. Soc. Lond. A 4562059–2079 (2000). [CrossRef]
8. L. Allen, M. Beijersbergen, R.J.C. Spreeuw, and J.P. Woerdman, “Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes”, Phys. Rev. A 458185–8189 (1992). [CrossRef] [PubMed]
9. F. Flossmann, U.T. Schwarz, and M. Maier “Propagation dynamics of optical vortices in Laguerre-Gaussian beams,” Opt. Commun. 250218–230 (2005). [CrossRef]
11. J.V. Hajnal, “Observation on singularities in the electric and magnetic fields of freely propagating microwaves,” Proc. R. Soc. Lond. A 430413–421 (1990). [CrossRef]
13. M.V. Berry and M.R. Dennis, “The optical singularities of birefringent dichroic chiral crystals,” Proc. R. Soc. Lond. A 4591261–1292 (2003). [CrossRef]
14. M.V. Berry and M.R. Jeffrey, “Conical diffraction: Hamilton’s diabolic point at the heart of crystal optics,” Prog. Opt. in press.
15. Y. A. Egorov, T. A. Fadeyeva, and A. V. Volyar, “Fine structure of singular beams in crystals: colours and polarization,” J. Opt. A: Pure Appl. Opt. 6S217–S228 (2004). [CrossRef]
16. F. Flossmann, U. T. Schwarz, M. Maier, and M. R. Dennis, “Polarization singularities from unfolding an optical vortex through a birefringent crystal,” Phys. Rev. Lett. 95253901 (2005). [CrossRef] [PubMed]
17. B. Riemann, in Bernhard Riemann’s gesammelte mathematische Werke und wissenschaftlicher Nachlass, edited by H. Weber (Teubner, Leipzig, 1892) p. 301.
18. R.M.A. Azzam and N.M. Bashara, Ellipsometry and polarized light (North-Holland, 1977).
19. D.S. Kliger, J.W. Lewis, and C.E. Randall, Polarized light in optics and spectroscopy (Academic Press, San Diego, 1990).
20. I. Freund, A.I. Mokhun, M.S. Soskin, O.V. Angelsky, and I.I. Mokhun, “Stokes singularity relations,” Opt. Lett. 27545–547 (2002). [CrossRef]
21. D.J. Struik, Lectures on Classical Differential Geometry (Dover, 1988).
22. A. Volyar, V. Shvedov, T. Fadeyeva, A.S. Desyatnikov, D.N. Neshev, W. Krolikowski, and Y.S. Kivshar, “Generation of single-charge optical vortices with an uniaxial crystal” Opt. Express 143724–3729 (2006). [CrossRef]
23. M.R. Dennis, “Braided nodal lines in wave superpositions,” New J. Phys. 5134 (2003). [CrossRef]
24. F. Flossmann, Singularitäten von Phase und Polarisation des Lichts (PhD Thesis, University of Regensburg, 2006).