## Abstract

Coherent vector beams with involved states of polarization (SOP) are widespread in the literature, having applications in laser processing, super-resolution imaging and particle trapping. We report novel vector beams obtained by transforming a Gaussian beam passing through a biaxial crystal, by means of the conical refraction phenomenon. We analyze both experimentally and theoretically the SOP of the different vector beams generated and demonstrate that the SOP of the input beam can be used to control both the shape and the SOP of the transformed beam. We also identify polarization singularities of such beams for the first time and demonstrate their control by the SOP of the input beam.

© 2015 Optical Society of America

## 1. Introduction

The state of polarization (SOP) is one of the fundamental signatures of light fields associated with their vectorial nature. In general, at each point in space the dynamics of the electric field vector of a wave can be described by an ellipse. This ellipse is known as polarization ellipse and it is characterized by the orientation of its major axis through the azimuth angle Φ ∈ [0, *π*] and by the ellipticity parameter *β* ∈ [−*π*/4, *π*/4] so that tan*β* is the ratio of the axes of the polarization ellipse. If *β* = 0 the light field is linearly polarized, while if *β* = ±*π*/4 the SOP will be circular (left handed for ‘−’ and right handed for ‘+’, if we consider an observer looking in the direction from which light is coming). Usually coherent light beams are homogeneously polarized, i.e. the SOP is identical for all points at any transverse plane along the beam propagation. However, there exist light beams possessing non-homogenous polarization, known as vector beams, such as the well known radial or azimuthal polarizations [1] or even beams with more involved polarization distributions [2–4]. The non-homogeneous polarization distribution of vector beams can lead to singular points where the SOP is exactly circular (C points), lines along which the SOP is linear (L lines) or disclinations where the instantaneous electric field is null [5–10]. Vector beams have been applied to laser material processing, optical imaging, atomic spectroscopy, and optical trapping (see [1] and the references therein), among many others.

In recent years, there has been a renewed interest in the conical refraction (CR) phenomenon ocurring in biaxial crystals, which can also provide different vector beams [11–23]. In CR, when a focused Gaussian input beam with waist radius *w*_{0} propagates along one of the optic axes of a biaxial crystal it appears transformed at the focal plane as a pair of concentric and well resolved bright rings split by a dark (Poggendorff) ring of geometrical radius *R*_{0}, as shown in Fig. 1(b). This intensity distribution is found as long as the condition *ρ*_{0} ≡ *R*_{0}/*w*_{0} ≫ 1 is satisfied [14,24, 25], while other light structures have been found for *ρ*_{0} ≤ 1.50 [20,21,23]. *R*_{0} is the product of the crystal length, *l*, and the CR semi-angle *α*, i.e. *R*_{0} = *lα* [15]. The CR semi-angle *α* depends on the principal refractive indices of the crystal as
$\alpha =\sqrt{\left({n}_{2}^{2}-{n}_{1}^{2}\right)\left({n}_{3}^{2}-{n}_{2}^{2}\right)/{n}_{2}^{2}}$, where we have assumed *n*_{1} < *n*_{2} < *n*_{3}. The non-uniform polarization distribution of a CR beam at the focal plane under Poggendorff splitting conditions is depicted by blue double arrows in Fig. 1(b). At any point of the CR rings, the electric field is linearly polarized with the azimuth varying so that every pair of diametrically opposite points have orthogonal polarizations. This polarization distribution depends only on the orientation of the plane of optic axes of the biaxial crystal. As a consequence, the intensity distribution along the ring depends on the input SOP, corresponding the one depicted in Fig. 1 to an incident circularly polarized beam. For *ρ*_{0} ≫ 1, the SOP of the CR beam at the focal plane is already well known [13,16,18,24].

Here, for the first time to our knowledge, we investigate the CR SOP out of the focal plane, including values of *ρ*_{0} ≤ 1.50, when Poggendorff fine splitting vanishes. By means of the Stokes vector formalism, we characterize the resulting novel vector beams of CR and show that the SOP of the input beam can be used to control both the shape and the SOP of the transformed CR beam.

The article is organized as follows. Section 2 is devoted to briefly review the basics of the Stokes vector and CR formalisms needed for the subsequent sections. In Sections 3 and 4, we report for *ρ*_{0} ≫ 1 and *ρ*_{0} ≤ 1.50, respectively, the generation of novel CR vector beams both theoretically and experimentally. We characterize these vector beams by measuring their Stokes parameters and propose methods to manipulate them as, for instance, making use of the SOP of the input beam. Finally, we summarize the main results of this work and discuss potential applications of these CR vector beams in different fields of optics.

## 2. Theoretical background

#### 2.1. Stokes vector formalism

The standard tool to analyze the SOP of a light beam is the Stokes vector: *S* = (*S*_{0}, *S*_{1}, *S*_{2}, *S*_{3}). For an electric field **E** = (*E _{x}*,

*E*) with intensity

_{y}*I*the Stokes parameters read [10]:

*I*

_{Φ}(Φ = 0°, 45°, 90°, 135°) indicates the intensity of linearly polarized light with azimuth Φ, and

*I*and

_{R}*I*indicate the intensity of right- and left-handed circularly polarized light, respectively. In what follows, we will use equations normalized to

_{L}*E*

^{2}, i.e. we will consider

*I*=

*E*

^{2}= 1. These definitions show that

*S*

_{0}account for the intensity of the light beam,

*S*

_{1}measures the amount of light which is linearly polarized (LP) in the vertical/horizontal basis,

*S*

_{2}does the same but with the diagonal basis and

*S*

_{3}relates the SOP in the right- and left- circularly polarized (CP) basis. The following equations show how the Stokes parameters relate to the azimuth and ellipticity of the polarization ellipse [10]:

#### 2.2. Conical refraction theory

The theoretical model predicting the transformation of an input beam once it passes along one of the optic axis of a biaxial crystal is based on the Belsky-Khapalyuk-Berry (BKB) integrals [11,13]. For a uniformly polarized input beam with electric field **E** and cylindrically symmetric 2D Fourier transform
$a(\kappa )=2\pi {\int}_{0}^{\infty}r{E}_{\text{in}}(r){J}_{0}(\kappa r)dr$, the electric field amplitude behind the crystal can be written as follows:

*C*=

*B*

_{1}cos (

*φ*+

*φ*) and

_{C}*S*=

*B*

_{1}sin (

*φ*+

*φ*).

_{C}**E**= (

*E*,

_{x}*E*) denotes the transverse Cartesian components of the output field,

_{y}**e**= (

_{0}*e*

_{x},

*e*

_{y}) is the unit vector of the input electric field, and

*φ*is the polar angle that defines the orientation of the plane of the crystal optic axes.

_{C}*Z*=

*z/z*and

_{R}*ρ*= (cos

*φ*, sin

*φ*)

*r/w*

_{0}define cylindrical coordinates whose origin is associated to the geometrical center of the light pattern at the focal plane (

*Z*= 0). ${z}_{R}=\pi {w}_{0}^{2}/\lambda $ denotes the Rayleigh length. Without loosing generality and for simplicity we will consider

*φ*= 0. For any other case, it would be enough to replace

_{C}*φ*by

*φ*+

*φ*in the following expressions.

_{C}*B*

_{0}and

*B*

_{1}are the main integrals of the Belsky–Khapalyuk–Berry (BKB) paraxial solution for CR [25]:

*η*=

*κw*

_{0}and

*J*is the

_{q}*q*-order Bessel function of the first type. Below we consider a fundamental Gaussian input beam with normalized transverse profile of the electric field amplitude $E(\rho )=\sqrt{2P/\pi w{(Z)}^{2}}\text{exp}(-{\rho}^{2})$ and 2D Fourier transform $a(\eta )=\sqrt{2P\pi /w{(Z)}^{2}}\text{exp}(-{\eta}^{2}/4)$, being $w(Z)={w}_{0}^{2}\sqrt{1+{Z}^{2}}$ and

^{th}*P*the power of the input beam. Note the central role

*ρ*

_{0}≡

*R*

_{0}/

*w*

_{0}in Eqs. (8) and (9) so it will determine the resulting CR intensity distribution. Accordingly, we will use

*ρ*

_{0}as the control parameter for our investigations. If

*ρ*

_{0}≫ 1 is fulfilled then the transverse intensity pattern at the focal plane of the system is formed by the two characteristic bright rings split by the dark Poggendorff one, as shown in Fig. 1(b). As the imaging plane is moved away from the focal plane, more involved structures are found including secondary rings and a bright central spot known as the Raman spot, see Fig. 1(a). The intensity maximum of the Raman spot is found at positions given by [26] Note that Eq. (10) is obtained under the assumption that

*ρ*

_{0}≫ 1, i.e. in the limit

*ρ*

_{0}→ ∞. Out of this approximation, the greater the value of

*ρ*

_{0}is, the better it is predicted by Eq. (10).

## 3. State of polarization for *ρ*_{0} ≫ 1

The first information that can be extracted from Eqs. (7)–(9) with respect to the SOP of the CR beam is that at *ρ* = 0 there is only contribution of *B*_{0}, since *B*_{1} ∝ *J*_{1}(*ηρ* = 0) = 0. Additionally, from Eqs. (8) and (9) note that the SOP of the *B*_{0} component is **e**_{0}. As a consequence, the center of the CR beam will possess always the same SOP as the input beam. This fact, that was already pointed out in [27–29], will be discussed with more detail below. To obtain the Stokes parameters of the CR beam, Eqs. (2)–(4) must be combined with Eqs. (7)–(9). For a CP input beam, the electric field and intensity beyond the crystal become:

**e**= (cos Φ, sinΦ).

_{0}For well resolved concentric rings with Poggendorff splitting, i.e. for *ρ*_{0} ≫ 1, Eqs. (8), (9), (13) and (16) show that a radially symmetric intensity pattern of CR is obtained only for a CP input beam. Instead, for a LP input beam, a crescent annular intensity pattern appears such that the zero intensity point is obtained for the ring position that possesses orthogonal polarization to the input beam. In both cases, the polarization distribution of the CR pattern is the same. Every point of the rings is linearly polarized and the azimuth rotates continuously along the ring so that every two diametrically opposite points have orthogonal polarizations.

Figure 2 presents the numerically obtained Stokes parameters at transverse sections of the CR beam (*ρ*_{0} = 10) at *Z* = 0 (first and second rows) and *Z* = 10.92 (third and fourth rows) obtained from a RHCP (first and third rows) and a LP (Φ = 45°) (second and fourth rows) Gaussian input beam. At the focal plane, see first two rows in Fig. 2, the SOP described by the Stokes parameters is the expected: symmetric pattern for the RHCP case and with a node at *φ* = 270° (Φ = 135°, since we have used *φ _{C}* = 0°). Last column in Fig. 2 demonstrates that at the focal plane the SOP of the CR beams, either RHCP or LP (Φ = 45°) (or any other) is linear, i.e.

*S*

_{3}= 0. In contrast the Raman spot,

*Z*= 10.92, all Stokes parameters are substantially different from zero, as shown in the last two rows in Fig. 2.

Stokes parameters in Fig. 2 clearly identify polarization singularities of CR beams. For RHCP input light, it is a C-point at the center of the CR beam, i.e. it is of circular polarization at any point of the beam center along propagation. The value of the *S*_{3} transverse pattern at its center is maximum, while the other two Stokes parameters *S*_{1} and *S*_{2} have zero values. For LP_{45°} input light, one can identify L-line singularity. In Fig. 2 this line can be identified as a vertical line of zero value of the Stokes parameters *S*_{3} and *S*_{1} and nonzero value of *S*_{2}. It belongs to the plane defined by the points of CR ring with linear SOP of LP_{Φ=45°} and LP_{Φ=135°}.

Media 1 and
Media 2 show, respectively, the spatial evolution of the Stokes parameters along the propagation direction for a RHCP and LP (Φ = 45°) Gaussian input beam. As we move away from the focal plane, the bright rings become wider and the intensity at the Poggendorff dark ring is no longer zero. The intensity of the outer ring decreases while the intensity of the inner ring increases. The outer ring expands and the inner ring becomes smaller in radius. As a result, at *Z* ≈ 6 a spot in the center of the pattern appears. Finally, the inner ring shrinks into a bright spot at *Z* = 10.92 for *ρ*_{0} = 10, corresponding to the Raman spot. At this point, there can be found an additional type of polarization singularity independently of the polarization state of an input beam. They are clearly distinguishable by expecting Stokes parameter *S*_{3} that defines degree of circular polarization of the field. Inspection of *S*_{3} far from the focal plane, for instance, at the Raman spot as demonstrated in the last two rows in Fig. 2, reveals alternating annular regions of RHCP and LHCP states. These annular regions of circular polarization are separated by circles of null intensity. At these circles the field is linearly polarized. Therefore, these polarization singularities can be called as L-circles. This behavior has been also observed in the focusing of radially polarized beams and explored in detail in Ref. [30]. For a LP (Φ = 45°) input beam, the central spot is broken by a line of null intensity (L-line) that connects the two points with LP and azimuth Φ = 45° and Φ = 135°, as commented before.

One feature that supplements demonstrate is that polarization singularities, i.e. C-point for CP input beam, L-lines for LP input beam and L-circles, are invariants of CR beam propagation behind the crystal. Another feature demonstrated by the supplement movies is the inversion of the center of the *S*_{1} and *S*_{2} parameters after the focal plane, which can be understood as a manifestation of the Gouy phase [31,32] of the CR beam. This effect will be more appreciable and discussed in the CR vector beams reported in the following section.

To test the validity of the obtained theoretical results, we have performed the corresponding experimental measurements. Figure 3 shows the experimental set-up. The input beam is obtained from a diode laser at 640nm coupled to a monomode fiber with a collimator (Monocrom C-6305C-4S), yielding a beam waist of *w*_{0} = 1.26mm. The *M*^{2} factor of the beam has been measured to be 1.0769, which is very close to an ideal Gaussian beam (*M*^{2} = 1). Therefore, in what follows we will refer to the beam emerging from the collimator as a Gaussian beam. After the collimator we place a linear polarizer to fix the state of polarization of the Gaussian beam. To generate the different input polarization states (LP with Φ = 45° and RHCP) we use a *λ*/2 and a *λ*/4 waveplates. We use lenses with different focal lengths to modify the waist radius *w*_{0} of the input beam and to reach different values of *ρ*_{0}. The beam passes along one of the optic axes of a biaxial crystal, placed always before the expected focal plane of the beam. An additional imaging lens is used to take different planes along the beam propagation direction and transfer them onto the CCD camera. Linear and circular polarizers are used to measure the Stokes parameters of the beam after being transformed by the CR phenomenon. We use a commercially available (CROptics) KGd(WO_{4})_{2} biaxial crystals with *α* = 16.9 mrad and length *l* = 10.5mm yielding CR ring radius of *R*_{0} = 180*μ*m. For more details about the material and the alignment procedure, see e.g. Ref. [23].

Figure 2(b) shows the obtained experimental Stokes parameters for *ρ*_{0} = 10.81. The experimental results agree with the theoretical results presented in Fig. 2. Discrepancy has been observed only for the *S*_{3} parameter for the case of a LP input beam (see last image of second rows). This can be explained in terms of the experimental error introduced by the polarization state detector elements (PSA) used, that disturb the beam shape and its position, which is central for the quality of the experimental results. Additionally it must be taken into account that the PSA and PSG used are designed for a working wavelength of 633nm, while our laser system has a wavelength of 640nm.

## 4. State of polarization for *ρ*_{0} ≤ 1.50

For *ρ*_{0} ≤ 1.50 CR patterns are significantly different from the double bright concentric rings with clear Poggendorff splitting (occurring for *ρ*_{0} ≫ 1 as shown in the previous section). The region *ρ*_{0} ≤ 1.50 has been explored recently [20, 21, 23] showing that CR can be used to new CR lasers [22], to increase the directivity of laser beams [20], to generate a super-Gaussian beam [23], to create a three dimensional dark focus [21] and even to develop a novel scheme for super-resolution microscopy [35, 36]. However, in all these works, no deep insight about the SOP of the generated CR beams has been provided. In what follows by considering the Stokes parameters we uncover the evolution of the SOP and polarization singularities of the CR beams and demonstrate how they depend on the SOP of the input beam. Figure 4 shows the main features and general view of CR beams with *ρ*_{0} = [1.50, 0.92, 0.45]. The cross-section of the CR transverse intensity pattern at the focal plane and along the axial direction *Z* at the beam center are shown in Figs. 4(a) and (b), respectively. Figs. 4(c)–(e) are 2D density plots of the intensity of the CR beams in the *Z*–*ρ* plane.

Figure 5 presents density plots of the numerically calculated Stokes parameters for (a) *ρ*_{0} = 1.50, (b) *ρ*_{0} = 0.92 and (c) *ρ*_{0} = 0.45 at *Z* = 0 (first and second rows) and out of the focal plane (third and fourth rows). For *ρ*_{0} = 1.50 and *ρ*_{0} = 0.92 these planes give the axial intensity maximum, while for *ρ*_{0} = 0.45 we have considered the plane where the cross-section area of the beam is doubled [23]. First and third rows refer to a RHCP Gaussian input beam while second and fourth rows present the case of a LP (Φ = 45°) Gaussian input beam. For *ρ*_{0} ≤ 1.50 the position of the Raman spots are not well determined by Eq. (10) and the position of the axial intensity maxima must be determined for each particular case of *ρ*_{0}. Its value (*Z*) is indicated at each image.

While the transverse distribution for intensity and polarization of CR beam look different for the *ρ*_{0} ≤ 1.50 case with respect to *ρ*_{0} ≫ 1, it should be noted that singularities, i.e. topological structure of the CR beam, remains similar. It means that (i) the CR beam center constitutes a C-point for the case of RHCP input beam (see the first and third rows in Fig. 5), (ii) there is L-line polarization singularity for the LP_{Φ=45°} input beam (see the second and fourth rows in Fig. 5) and (iii) there are L-circle singularity for input beam with arbitrary SOP (see the third and fourth rows in Fig. 5).

If we consider the special case of *ρ*_{0} = 0.92, the central point at the focal plane is a null-intensity point and therefore all Stokes parameters are also 0 at the beam center. Vanishing intensity at the beam center leads to another feature: at the focal plane the sign of the circular polarization state associated with the Stokes parameter *S*_{3} is the opposite with respect to the input beam. In other words, at the focal plane the CR beam is predominantly LHCP for the RHCP input beam that we consider in this article (see first row in Fig. 5 for the case of *ρ*_{0} = 0.92). This feature is easily deduced from the mathematical formulation, since for this value of *ρ*_{0} the value of the integral *B*_{0} at any radial point is null, which leads to dominance of the contribution of integral *B*_{1} associated with a SOP orthogonal to the input one. For the case of a LP (Φ = 45°) input beam and for all values of *ρ*_{0} investigated, at the focal plane *Z* = 0 the transverse patterns are crescent-like, with the point of intensity minimum at a point, which is diagonally opposite to the point with maximum intensity. In the case of LP_{Φ=45°} input beam that we consider, the intensity minimum is observed at the bottom, corresponding to the point of LP with Φ = 135°. In contrast to the RCHP case, for all three values of *ρ*_{0} the *S*_{3} parameter is null, which means that the patterns are completely LP. Out of the focal plane *S*_{3} ≠ 0 except at the beam center, where there is an L-line connecting the points with LP Φ = 45° and Φ = 135°.

Additional features of CR beams are revealed by considering the evolution of the transverse pattern of the Stokes parameters along CR beam propagation, as shown in
Media 3,
Media 4,
Media 5,
Media 6,
Media 7,
Media 8 for a RHCP and LP (Φ = 45°) input beam for *ρ*_{0} = 1.50, *ρ*_{0} = 0.92 and *ρ*_{0} = 0.45. For RHCP input light (
Media 3,
Media 5,
Media 7) the focal plane (*Z* = 0) is a symmetry plane for *S*_{3}. In contrast, *S*_{1} and *S*_{2} are rotated roughly 180° before and after the focal plane, which is associated to the Gouy phase in vector beams [31, 32]. For a LP (Φ = 45°) input beam (
Media 4,
Media 6,
Media 8) the Stokes parameters *S*_{1} and *S*_{2} are symmetric with respect to the focal plane and now *S*_{3} suffers from a phase shift of 180°. In general, analogous results are obtained with focused Gaussian beams: due to the Gouy phase shift, that induces an additional phase in the electric field of 180° at the focal plane, the *S*_{1} and *S*_{2} of a Gaussian beam rotate around the *Z* axis by 180°. However, note that at variance with the results expected for a focused Gaussian beam, the *S*_{1} and *S*_{2} parameters of CR beams obtained from CP input beams not only rotate along the *Z* axis but also twist around it at the focal plane. CR beams possess non-integer orbital angular momentum (OAM) that is generated due to the spin-orbit coupling provided by the biaxial crystal [37]. The value and handiness of the OAM depend both on the ellipticity of the input beam and on the *ρ*_{0} parameter. For *ρ*_{0} ≫ 1 and a CP beam, i.e. with spin angular momentum (SAM) *σ _{s}* = ±1 (+ for LHCP and − for RHCP), the CR beam carries total optical angular momentum ±1/2 and zero SAM. For values of

*ρ*

_{0}out of that condition, the spin-orbit coupling differs. For instance, for

*ρ*

_{0}= 0.924 the SAM of the CR beam has opposite handiness with respect to the input beam, which explains the results discussed above. As reported in Ref. [37], for

*ρ*

_{0}= 1.50 the OAM of the CR beam is maximum compared with any other value of

*ρ*

_{0}. At the focal plane the twist is more appreciable at the central spot. Additionally, the velocity of rotation of the

*S*

_{1}and

*S*

_{2}parameters of that region as the beam approaches to the focal plane (

*Z*= 0) is much faster than the rotation velocity of the outer ring. For the value

*ρ*

_{0}= 0.445 the beam carries almost null OAM, which explains the small twist of the

*S*

_{1}and

*S*

_{2}parameters at the focal plane and its similarities with a homogeneously CP Gaussian beam. As a summary, the combination of both SAM and OAM of the CR vector beams leads to the rotation and twisting of the

*S*

_{1}and

*S*

_{2}Stokes parameters when CP input beams into the biaxial crystal are considered.

We have performed experiments in order to prove theoretical findings reported above. Figure 6 shows the experimentally measured Stokes parameters for *ρ*_{0} = 1.48 (box (a)), *ρ*_{0} = 0.95 (box (b)) and *ρ*_{0} = 0.44 (box (c)). For these experiments, the same set-up shown in Fig. 3 was used but in this case taking a 2.3mm long KGd(WO_{4})_{2} biaxial crystal yielding CR ring radius of *R*_{0} = 39 *μ*m. In general, the theoretical predictions agree with the experimental results. Again, the *S*_{3} parameter for a LP input beam (last image of second and fourth rows) is the measurement that mostly differs with respect to the numerical predictions. In addition to the experimental difficulties commented above, i.e. the experimental error introduced by the polarization state detector elements used that disturb the beam shape and its position (being these ones particularly significant for large values of *w*_{0}), here it must be also taken into account the fact that small changes in the *ρ*_{0} can modify quantitatively the CR pattern.

Finally, Fig. 7 presents the calculated transverse patterns and SOP (blue curves) of the vector beams analyzed in this work for a RHCP (first row) and a LP (Φ = 45°) (second row) input Gaussian beam at the focal plane. As it can be observed, for low values of *ρ*_{0} the SOP of the generated vector beams is rich and include regions with different ellipticity and azimuth, leading to Poincare beams [38].

## 5. Conclusions

In summary, we have studied in detail the SOP of conical refraction based vector beams for a wide range of *ρ*_{0} and for different SOP of the input beam. We have determined the Stokes parameters of the CR beam at different transverse planes along the beam propagation direction and we have shown that both the shape and the SOP of the transformed beams depend on the SOP of the input beam. We have shown that the polarization distribution formed by orthogonal polarizations at any two radially opposite points of the pattern and usually associated to the CR phenomenon remains relevant under the condition *ρ*_{0} ≫ 1 and at the focal plane only. For *ρ*_{0} ≤ 1.50 and CP input beams we have found that CR beams demonstrate more involved structure of non-homogeneously elliptically polarized states not only with different azimuth but also with different ellipticity. In contrast, for LP input beams, the SOP of the CR has been reported to be completely linear at the focal plane and with variable ellipticity and azimuth out of it. Agreement between the theoretical predictions and the experimental results has been obtained.

Additionally, we have demonstrated for the first time to our knowledge, experimental results on polarization singularities of CR beams and demonstrate how types of polarization singularities can be changed by and controlled with varying SOP of the input light beams. Such polarization singularities as C-points, L-lines and L-circles have been identified for CR beams.

The reported results can be particularly interesting for experiments with tightly focused beams [1], for the generation of novel polarizations in CR [19], in optical micromanipulation [26, 33, 34], mode conversion between Heremite–Gauss-like beams and Laguerre–Gauss-like beams [39] and in super-resolution imaging [35, 36]. It is also promising the generation of polarization-tunable potentials to inject, extract and manipulate ultra-cold atoms [25, 40]. Additionally, by taking into account that the CR beams posses orbital angular momentum [37] and that the biaxial crystals used are transparent to a wide frequency range, the presented technique could be an alternative for the production of high-frequency vortex beams for molecular-scale super-resolution microscopy [41]. Finally, by using quantum sources, CR vector beams can be a tool for testing new Bell-like inequalities with hybrid polarization-momentum states useful in quantum information technologies [42].

## Acknowledgments

We acknowledge financial support from the Spanish Ministry of Science and Innovation (MICINN) (contracts FIS2011-23719 and FIS2012-39158-C02-01). A.T. acknowledges financial support through the grant AP2010-2310.

## References and links

**1. **Q. Zhan, “Cylindrical vector beams: from mathematical concepts to applications,” Adv. Opt. Photon. **1**, 1–57 (2009). [CrossRef]

**2. **D. Maluenda, I. Juvells, R. Martínez-Herrero, and A. Carnicer, “Reconfigurable beams with arbitrary polarization and shape distributions at a given plane,” Opt. Express **21**, 5424–5431 (2013)., [CrossRef]

**3. **X.-L. Wang, J. Ding, W.-J. Ni, C.-S. Guo, and H-T. Wang, “Generation of arbitrary vector beams with a spatial light modulator and a common path interferometric arrangement,” Opt. Lett. **32**, 3549–3551 (2007). [CrossRef] [PubMed]

**4. **A. 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**, 10848–10863 (2010) [CrossRef] [PubMed]

**5. **J. F. Nye and M. V. Berry, “Dislocations in wave trains,” Proc. R. Soc. London, Ser. A **336**, 165–190 (1974). [CrossRef]

**6. **M. R. Dennis, K. O Holleran, and M. J. Padgett, “Singular optics: Optical vortices and polarization singularities,” Prog. Opt. **53**, 293–363 (2009). [CrossRef]

**7. **M. R. Dennis, R. P. King, B. Jack, K. O Holleran, and M. J. Padgett, “Isolated optical vortex knots,” Nature Phys. **6**, 118–121 (2010). [CrossRef]

**8. **D. Kleckner and W. T. M. Irvine, “Creation and dynamics of knotted vortices,” Nature Phys. **9**, 253–258 (2013). [CrossRef]

**9. **I. Freund, “Ordinary polarization singularities in three-dimensional optical fields,” Opt. Lett. **37**, 2223–2225 (2012). [CrossRef] [PubMed]

**10. **M. Born and E. Wolf, *Principles of Optics: Electromagnetic Theory of Propagation, Interference and Diffraction of Light*, 7 (expanded) ed. (Cambridge University Press, Cambridge, 1999). [CrossRef]

**11. **A. M. Belskii and A. P. Khapalyuk, “Internal conical refraction of bounded light beams in biaxial crystals,” Opt. Spectrosc. (USSR) **44**, 436–439 (1978).

**12. **A. M. Belsky and M. A. Stepanov, “Internal conical refraction of coherent light beams,” Opt. Communs. **167**, 1–5 (1999). [CrossRef]

**13. **M. V. Berry and M. R. Jeffrey, “Conical diffraction: Hamilton’s diabolical point at the heart of crystal optics,” Prog. Opt. **50**, 13–50 (2007). [CrossRef]

**14. **M. V. Berry, “Conical diffraction asymptotics: fine structure of Poggendorff rings and axial spike,” J. Opt. A: Pure Appl. Opt. **6**289–300 (2004). [CrossRef]

**15. **T. K. Kalkandjiev and M. Bursukova, “The conical refraction: an experimental introduction,” Proc. SPIE **6994**, 69940B (2008). [CrossRef]

**16. **A. Turpin, Yu. V. Loiko, T. K. Kalkandjiev, and J. Mompart, “Wave-vector and polarization dependence of conical refraction,” Opt. Express **21**4503–4511 (2013). [CrossRef] [PubMed]

**17. **A. Turpin, Y. V. Loiko, T. K. Kalkandjiev, and J. Mompart, “Multiple rings formation in cascaded conical refraction,” Opt. Lett. **38**, 1455–1457 (2013). [CrossRef] [PubMed]

**18. **G. S. Sokolovskii, D. J. Carnegie, T. K. Kalkandjiev, and E. U. Rafailov, “Conical Refraction: New observations and a dual cone model,” Opt. Express **21**, 11125–11131 (2013). [CrossRef] [PubMed]

**19. **C. F. Phelan, J. F. Donegan, and J. G. Lunney, “Generation of a radially polarized light beam using internal conical diffraction,” Opt. Express **19**, 21793–21802 (2011). [CrossRef] [PubMed]

**20. **V. Peet, “Improving directivity of laser beams by employing the effect of conical refraction in biaxial crystals,” Opt. Express **18**, 19566–19573 (2010). [CrossRef] [PubMed]

**21. **Yu. V. Loiko, A. Turpin, T. K. Kalkandjiev, E. U. Rafailov, and J. Mompart, “Generating a three-dimensional dark focus from a single conically refracted light beam,” Opt. Lett. **38**, 4648–4651 (2013). [CrossRef] [PubMed]

**22. **Yu. V. Loiko, G. S. Sokolovskii, D. Carnegie, A. Turpin, J. Mompart, and Edik U. Rafailov, “Laser beams with conical refraction patterns,” Proc. SPIE **8960**89601Q (2014). [CrossRef]

**23. **A. Turpin, Yu. V. Loiko, T. K. Kalkandkiev, H. Tomizawa, and J. Mompart, “Super-Gaussian conical refraction beam,” Opt. Lett. **39**, 4349–4352 (2014). [CrossRef] [PubMed]

**24. **A. Peinado, A. Turpin, A. Lizana, E. Fernndez, J. Mompart, and J. Campos, “Conical refraction as a tool for polarization metrology,” Opt. Lett. **38**, 4100–4103 (2013). [CrossRef] [PubMed]

**25. **A. Turpin, J. Polo, Yu. V. Loiko, J. Küber, F. Schmaltz, T. K. Kalkandjiev, V. Ahufinger, G. Birkl, and J. Mompart, “Blue-detuned optical ring trap for Bose-Einstein condensates based on conical refraction,” in Press.

**26. **A. Turpin, V. Shvedov, C. Hnatovsky, Yu. V. Loiko, J. Mompart, and W. Krolikowski, “Optical vault: A reconfigurable bottle beam based on conical refraction of light,” Opt. Express **21**, 26335–26340 (2013). [CrossRef] [PubMed]

**27. **C. F. Phelan, D. P. ODwyer, Y. P. Rakovich, J. F. Donegan, and J. G. Lunney, “Conical diffraction and Bessel beam formation with a high optical quality biaxial crystal,” Opt. Express **17**, 12891–12899 (2009). [CrossRef] [PubMed]

**28. **V. Peet and D. Zolotukhin, “Free-space evolution of focused Gaussian beams transformed by conical diffraction in a biaxial crystal,” Opt. Commun. **283**3011–3016 (2010). [CrossRef]

**29. **V. Peet, “The far-field structure of Gaussian light beams transformed by internal conical refraction in a biaxial crystal,” Opt. Commun. **311**150–155 (2010). [CrossRef]

**30. **R. W. Schoonover and T. D. Visser, “Polarization singularities of focused, radially polarized fields,” Opt. Express **14**, 5733–5745 (2006). [CrossRef] [PubMed]

**31. **R. W. Boyd, “Intuitive explanation of the phase anomaly of focused light beams,” J. Opt. Soc. Am. **70**, 877–880 (1980). [CrossRef]

**32. **G. M. Philip, V. Kumar, G. Milione, and N. K. Viswanathan, “Manifestation of the Gouy phase in vector-vortex beams,” Opt. Lett. **37**, 2667–2669 (2012). [CrossRef] [PubMed]

**33. **D.P. ODwyer, C.F. Phelan, K.E. Ballantine, Y.P. Rakovich, J. G. Lunney, and J.F. Donegan, “Conical diffraction of linearly polarised light controls the angular position of a microscopic object,” Opt. Express **18**, 27319–27326 (2009). [CrossRef]

**34. **V. Shvedov, A. R. Davoyan, C. Hnatovsky, N. Engheta, and W. Krolikowski, “A long-range polarization-controlled optical tractor beam,” Nature Photon. **21**, 26335–26340 (2014).

**35. **S. Rosen, G. Y. Sirat, H. Ilan, and A. J. Agranat, “A sub wavelength localization scheme in optical imaging using conical diffraction,” Opt. Express **21**, 10133–10138 (2013). [CrossRef] [PubMed]

**36. **G.Y. Sirat, S. Shorte, L.P.O. Braitbart, L. Moisan, J.Y. Tinevez, J. Caron, and C. Fallet, “Conical diffraction based super-resolution system for fluorescence microscopy: system Description and demonstration visualizing biological objects,” Microscopy and Microanalysis **13**, 1178–1179 (2013).

**37. **M. V. Berry, M. R. Jeffrey, and M. Mansuripur, “Orbital and spin angular momentum in conical diffraction,” J. Opt. A: Pure Appl. Opt. **7**, 685 (2005). [CrossRef]

**38. **A. M. Beckley, T. G. Brown, and M. A. Alonso, “Full Poincare beams,” Opt. Express **18**, 10777–10785 (2010). [CrossRef] [PubMed]

**39. **V. Peet, “Biaxial crystal as a versatile mode converter,” J. Opt. **12**, 095706 (2010). [CrossRef]

**40. **Yu. V. Loiko, V. Ahufinger, R. Menchon-Enrich, G. Birk, and J. Mompart, “Coherent injecting, extracting, and velocity-filtering of neutral atoms in a ring trap via spatial adiabatic passage,” Eur. Phys. J. D **68**, 147 (2014). [CrossRef]

**41. **G. Gariepy, J. Leach, K. T. Kim, T. J. Hammond, E. Frumker, R. W. Boyd, and P. B. Corkum, “Creating High-Harmonic Beams with Controlled Orbital Angular Momentum,” Phys. Rev. Lett. **113**, 153901 (2014). [CrossRef] [PubMed]

**42. **R. Fickler, R. LApkiewicz, S. Ramelow, and A. Zeilinger, “Quantum entanglement of complex photon polarization patterns in vector beams,” Phys. Rev. A **89**, 060301(R) (2014). [CrossRef]