We experimentally address the wave-vector and polarization dependence of the internal conical refraction phenomenon by demonstrating that an input light beam of elliptical transverse profile refracts into two beams after passing along one of the optic axes of a biaxial crystal, i.e. it exhibits double refraction instead of refracting conically. Such double refraction is investigated by the independent rotation of a linear polarizer and a cylindrical lens. Expressions to describe the position and the intensity pattern of the refracted beams are presented and applied to predict the intensity pattern for an axicon beam propagating along the optic axis of a biaxial crystal.
© 2013 Optical Society of America
Light beams propagating in optically uniaxial crystals, with two principal dielectric indexes, exhibit birefringence or double refraction, i.e. decomposition of an input beam into ordinary (o) and extraordinary (e) beams with orthogonal linear polarizations. The extraordinary beam is laterally shifted with respect to the ordinary one, being this shift equal to zero when the input beam propagates along the optic axis. The intensity distribution between the o- and e- beams, expressed by the Malus law, depends only on one parameter: the relative orientation of the crystal with respect to the polarization plane of the input beam.
Biaxial crystals (BCs), with three different principal dielectric indexes and two optic axes, also exhibit double refraction. Nevertheless, propagation of light beams along the optic axis is substantially different from the uniaxial case. Hamilton predicted that, in this case, the beam should propagate conically inside the crystal and emerge as a hollow cylinder yielding an annular transverse intensity pattern or light ring of internal conical refraction (CR) . Simultaneously to this prediction, Hamilton also demonstrated that an infinity of wave-vectors distributed on a hollow cone entering in the a biaxial crystal along the biradial propagate with common Poynting vector along the latter one, the so-called external conical refraction. Soon afterwards Hamilton’s prediction, Lloyd reported the first experimental observation of both internal and external CR  and since then few experimental works have further studied CR experimentally (see  and references therein). Recently, internal CR has been applied to the laser technique , for optical micromanipulation  and for free space optical communications . Firsts theoretical efforts to describe precisely the internal conical refraction phenomenon were started by Portigal and Burstein . The complete calculation according to Fourier analysis was performed by Lalor  and by Schell and Bloembergen  and later on rewritten by Belsky and Khapalyuk [9, 10], Belsky and Stepanov  and Berry and Jeffrey . From the experimental point of view, the full validation of the internal CR was given by Fève and coworkers . However, a simple analytical formulation capable of predicting the transverse intensity pattern at the Lloyd (focal) plane after a light beam propagates along or near the optic axis of a biaxial crystal (or a cascade of BCs [13, 14]) is still lacking.
In this article, we experimentally address this last need by analyzing the wave-vector and polarization dependence of the internal CR phenomenon. In the following we will refer to internal CR simply as CR. Our experimental scheme is based on the propagation of elliptical beams (EBs) along the optic axis of a BC. EBs do not possesses continuous cylindrical symmetry, at variance with Gaussian beams typically used in conical refraction experiments. In Section 2, we demonstrate that EBs do not produce the otherwise expected transverse annular intensity pattern when they propagate through the BC, but split into two refracted beams with orthogonal linear polarizations. Simple phenomenological expressions for the position and the intensity of the refracted beams are derived in Section 3. Worth noticing, we show that the position of the refracted beams is defined not only by the crystal orientation, but also by the input beam properties, in particular, by the orientation of the preselected plane of input wave-vectors. Additionally, we demonstrate that the intensity distribution between the two refracted beams can be well described with a single angular parameter. When the biaxial crystal is rotated, the intensity distribution between the two refracted beams follows the transformation law recently pointed out by Loiko et al. for CR filtered beams, that differs from the well known Malus law for double refraction in uniaxial crystals. To verify the usefulness of the phenomenological transformation rules, we apply them in Section 4 to calculate the resulting intensity pattern for an incident axicon beam and confirm our approach with experimental results. Section 5 summarizes the main results and conclusions of the work.
2. Conical refraction of spatially anisotropic beams
Conical refraction has been reported only for input beams with intensity pattern possessing continuous cylindrical (rotation) symmetry around the propagation axis, see Fig. 1(a) and Fig. 1(c). This has been motivated by Hamilton’s prediction of the CR within the ray optics description, where rays are always associated with cylindrically symmetric collimated beams. Real beams of finite size can be modeled as a bundle of rays with propagation directions isotropically distributed around the beam axis. On the one hand, according to the modified Hamilton theory , CR can be interpreted as a wave-vector and polarization dependent type of double refraction, when each ray of the input bundle refracts into two rays at the entrance crystal surface and all refracted rays form a cone. In this process, the two refracted rays go to opposite points on the CR ring. On the other hand, it has been also demonstrated that in CR the transverse wave-vector components are conserved [2, 13]. This means that all input rays (plane waves) with wave-vectors confined in a certain plane defined by an azimuthal angle ϕ refract to a plane with the same azimuthal angle at the Lloyd plane, i.e. into two opposite points on the CR ring.
Therefore, double refraction along the optic axis of a BC and, in particular, the wave-vector and polarization dependence of the CR, can be studied with linearly polarized beams composed from plane waves whose wave-vectors are confined in a certain plane. Such beams are known as cylindrical beams and elliptical beams (EBs) provide their finite size approximation. The latter ones can be obtained from collimated Gaussian beams focused by a cylindrical lens. Indeed, we have observed that an EB passing through the BC does not generate the CR ring, but splits into two beams (in case of circularly polarized input beams) oppositely placed along the otherwise expected CR ring and tangent to it, see Fig. 1(d). The form and orientation of the refracted beams resemble the transverse intensity pattern and orientation of the input EB. Below, using EBs we will obtain phenomenological laws that describe double refraction along the optic axis in biaxial crystals.
3. Transformation rules of conical refraction
In this section, we will investigate the position and the relative intensity of the two refracted beams produced when an input EB propagates along the optic axis of a BC. The experimental set-up shown in Fig. 2(a) is used. The initial circularly polarized Gaussian beam with waist radius of w = 1mm is obtained from a 640nm diode laser coupled to a monomode fiber. Then a linear polarizer is introduced to fix the polarization plane in a well defined direction. The resulting linearly polarized beam is focused by a cylindrical lens of 150mm focal length with its flat face oriented strictly perpendicular to the beam propagation direction. The cylindrical lens only focuses the Gaussian beam in one direction, so that it transforms its transverse circular shape to an elliptical one with a ratio 3/100 of the semi-axes of the ellipse. As a consequence, the divergence of the generated EB is different along the focused and unfocused directions (wf = 30μm, θf = 6.8mrad; wuf = 1000μm, θuf = 0.2mrad, where f and uf subscripts refer to focused and unfocused directions, respectively). The EB is characterized by its polarization plane, represented by the azimuthal angle ϕE, and by its plane of wave-vectors (or K-plane), represented by azimuthal angle ϕK; see Fig. 2(b). Different EBs are obtained by rotating either the cylindrical lens or the linear polarizer. The BC is 28mm long and it was cut from a monoclinic centrosymmetric KGd(WO4)2 crystal. Its polished entrance surfaces (cross-section 6×4mm2, parallelism 10 arc seconds) are perpendicular to one of the optic axes (misalignment angle < 1.5mrad) and positioned strictly perpendicular to the beam propagation direction, so that the incoming beam to the BC passes along its optic axis. Both the orientation of the optic axis of the BC and the cylindrical lens are well controlled in the θ and φ directions in 3D spherical coordinates by a micrometer positioning system. The resulting pattern is captured by a CCD camera at the Lloyd (focal) plane behind the BC .
3.1. Position of the refracted beams
First, we briefly outline the CR annular pattern observed at the Lloyd plane with an input Gaussian beam. While the complete CR ring is observed for circular polarization of the input Gaussian beam, as it is shown in Fig. 1(d), a crescent ring with one point of the ring being of zero intensity appears for linear polarization. In both cases, the center of the CR ring is laterally shifted with respect to the incident beam. This shift can be represented by a vector G = R0 (cosϕG, sinϕG) that belongs to the plane of the crystal optic axes and is perpendicular to them . Its modulus is equal to the ring radius |G| ≡ R0 (in our set-up R0 = 476μm), given by a product of the length and the conical refraction semiangle of the BC (17mrad for KGd(WO4)2) .
Now we will experimentally deduce the lateral shift of the refracted beams for EBs propagated through a BC. In Fig. 3 two series of images present the transverse intensity pattern at the Lloyd plane recorded varying either ϕE (ϕK = 0°) (a) or ϕK (ϕE = 0°) (b) from 0 to 157.5° in 22.5° intervals, when the crystal orientation remains fixed at ϕG = 0°. The geometric center of the refracted beams coincides with the center of the otherwise expected CR ring, since it is shifted from the position of the initial EB by the vector G. This shift is shown schematically in Fig. 3(c), and it is subtracted in Fig. 3(a) and Fig. 3(b).
From Fig. 3 it is clear that in a local frame with origin at the ring center, the position of the two refracted beams does not depend on the input beam polarization, ϕE, but rotates linearly with ϕK. Moreover, the azimuthal angles ϕ± of the refracted beams are defined by the wave-vector plane, ϕK, of the input EB, namely,Fig. 3(c). In other words, the two refracted beams are located at diagonal positions of the otherwise expected CR ring. EBs obtained from the same Gaussian beam have the same initial position S, i.e.S does not depend on ϕE,K,G in this case. As a final comment, Eqs. (2) generalize the geometrical approach .
3.2. Relative intensity distribution of the refracted beams
Now we report how the intensity of the input EB is distributed between the two refracted beams. In the experiments, three parameters can be varied independently: ϕG, ϕE and ϕK, associated to the crystal orientation represented by G and to the polarization and wave-vectors planes of the incident EBs. Below we show that only one combination of these angles governs the relative intensity distribution between the two refracted beams. With this purpose, we have repeated the experiments shown in Fig. 3(a) and Fig. 3(b) for different orientations of the BC. Symbols (black crosses and red circles) in Fig. 4(a) and Fig. 4(b) show the corresponding experimental results for the intensities I± of the two refracted beams normalized with respect to the intensity of the incident beam. Black solid and red dashed curves represent their analytical fittings given by the following expressions:Eq. (2).Eq. (5) it follows that, with respect to the relative energy distribution the only significant parameters of the input EB is ϕχ. Therefore, the intensity splitting under CR can be expressed in terms of the difference between the parameter ϕχ and the orientation of the BC, given by ϕG. We would like to highlight here that the expression ϕχ = 2ϕE − ϕK that we have obtained experimentally, corrects the theoretical result derived in , where a plus sign was considered instead of a minus one. Eqs. (2) and (5) describe also recent experimental results on the propagation of CR filtered beams  along the optic axis of a BC. In that case, ϕχ = ϕG, where ϕG denotes the orientation of an initial biaxial crystal that produces a CR ring with lateral shift given by vector G0. Then, a pinhole placed at angle ϕ on the CR ring produces a CR filtered beam with polarization and plane of wave-vectors ϕK = ϕ that also splits into two beams after passing along the optic axis of a biaxial crystal.
Eqs. (2) and (5) constitute the Transformation Rules of Conical Refraction. For different polarization of the input beams they allow explaining the ratio of intensities for any pair of diagonally opposite points of the CR ring. These results corroborate with experimental observations that the complete CR ring appears only for beams with azimuthally continuous symmetric distribution of wave-vectors. As a proof of usefulness of the formalism introduced in this work, in the next section we apply the derived transformation rules to an axicon input beam.
4. Application of the transformation rules of CR to an axicon beam
The experimental set-up is sketched in Fig. 5. A linearly polarized conical beam is prepared when an initial circularly polarized Gaussian beam passes through a linear polarizer and an axicon lens. The axicon beam is focused by a spherical lens into the BC. The resulting pattern is captured with a CCD camera at the focal plane. As in previous experiments, the orientation of the optic axis of the BC, the spherical lens and the axicon lens are well controlled in the θ and φ directions in 3D spherical coordinates by a micrometer positioning system. Each infinitesimally thin azimuthal sector of the axicon lens characterized by azimuthal angle ϕ forms a thin prism that produces a wave with particular wave-vector whose transverse projection comprise an angle ϕk = ϕ. The axicon lens generates, therefore, a continuous collection of beams with ϕ ∈ [0, 2π). After the axicon, the refracted beams, following Snell law, have the same inclination angle θ0 with respect to the z-axis. At the focal plane of the lens they form a ring such that each point can be characterized by wave-vector plane, ϕk, and polarization plane, ϕE. In other words, each point of the axicon ring is an EB. In this case all these EB have their polarization plane fixed at ϕE and their wave-vector plane ϕ is varying continuously along the ring as shown in Fig. 6(a). Behind the BC, the refraction pattern can be calculated by applying Eq. (2) to every point of the input axicon annular beam taking into account the initial positions as given by S(ϕ) = Rax (cosϕ, sinϕ) (where Rax is the radius of the axicon light ring). Therefore, from Eq. (2) one can obtain the refracted pattern for an axicon beam:Eq. (5). All points of the incident axicon beam have the same intensity, which is distributed between the two refracted rings as follows: Fig. 6(a) shows the intensity pattern of the axicon beam with polarization and K-plane distribution. The experimental refraction pattern behind the BC is shown in Fig. 6(b). The pattern is formed by two concentric rings oppositely polarized, with polarization distribution analogous to that one obtained in a cascaded CR configuration . Polarization of the ring points with minimum (zero) and maximum intensities are orthogonal and parallel respectively to that of the input axicon beam. Fig. 6(c) presents the theoretical prediction of the light refraction pattern from Eqs.(6) and (7). Comparison between Fig. 6(b) and Fig. 6(c) shows that the theoretical prediction and the experiment (see Fig. 6(b)) are in complete agreement.
It is well known that an input Gaussian beam produces a light ring after propagating along the optic axis of a biaxial crystal. In this paper we have experimentally addressed the wave-vector and polarization dependence of the internal conical refraction phenomenon in biaxial crystals. We have experimentally proved that this phenomenon is only one particular case of refraction associated to beams possessing continuous axial cylindrical symmetry. In fact, we have reported that elliptical beams – obtained by focusing a Gaussian beam with a cylindrical lens – do not generate the complete characteristic light ring of conical refraction when they propagate along an optic axis of a biaxial crystal but split, in general, into two beams. We have phenomenologically obtained expressions for the positions and relative intensities of the resulting refracted beams. These expressions play an analogous role as the Malus law but for biaxial crystals and can be used to predict the transverse intensity pattern of much more involved incident beams as it has been demonstrated here for an input axicon beam.
The authors gratefully acknowledge financial support through Spanish MICINN contracts FIS2010-10004-E and FIS2011-23719, and the Catalan Government contract SGR2009-00347. A. T. acknowledges financial support through grant AP2010-2310 from the MICINN. Dr. José Luis Martínez is also acknowledged for his helpful comments.
References and links
1. J. G. O’Hara, “The prediction and discovery of conical refraction by William Rowan Hamilton and Humphrey Lloyd,” Proc. Roy. Ir. Acad. 82 (2), 231–257 (1982).
2. M. V. Berry and M. R. Jeffrey, “Conical diffraction: Hamiltons diabolical point at the heart of crystal optics,” Prog. Opt. 50, 13–50 (2007). [CrossRef]
4. D. P. O’Dwyer, K. E. Ballantine, C. F. Phelan, J. G. Lunney, and J. F. Donegan, “Optical trapping using cascade conical refraction of light,” Opt. Express 20, 21119–21125 (2012). [CrossRef]
5. A. Turpin, Yu. V. Loiko, T. K. Kalkandjiev, and J. Mompart, “Free-space optical polarization demultiplexing and multiplexing by means of conical refraction,” Opt. Lett. 37, 4197–4199 (2012). [CrossRef] [PubMed]
6. D. L. Portigal and E. Burstein, “Internal Conical Refraction,” J. Opt. Soc. Am. 59, 1567–1573 (1969). [CrossRef]
7. E. Lalor, “An Analytical Approach to the Theory of Internal Conical Refraction,” J. Math. Phys. 13, 449–454 (1972). [CrossRef]
8. A. J. Schell and N. Bloembergen, “Laser studies of internal conical diffraction. I. Quantitative comparison of experimental and theoretical conical intensity distribution in aragonite,” J. Opt. Soc. Am. 68, 1093–1098 (1978). [CrossRef]
9. A. M. Belskii and A. P. Khapalyuk, “Internal conical refraction of bounded light beams in biaxial crystals,” Opt. Spectrosc. 44, 436–439 (1978).
10. A. M. Belskii and A. P. Khapalyuk, “Propagation of confined light beams along the beam axes (axes of single ray velocity) of biaxial crystals,” Opt. Spectrosc. 44, 312–315 (1978).
11. A. M. Belsky and M. A. Stepanov, “Internal conical refraction of light beams in biaxial gyrotropic crystals,” Opt. Commun. 204, 1–6 (2002). [CrossRef]
12. J. P. Fève, B. Boulanger, and G. Marnier, “Experimental study of internal and external conical refractions in KTP,” Optics Commun. 105, 243–252 (1994). [CrossRef]
13. M.V. Berry, “Conical diffraction from an N-crystal cascade,” J. Opt. 12, 075704 (2010). [CrossRef]
14. A. Abdolvand, “Conical diffraction from a multi-crystal cascade: experimental observations,” Appl. Phys B 103, 281–283 (2011). [CrossRef]
15. Y. Loiko, M. A. Bursukova, T. K. Kalkanjiev, E. U. Rafailov, and J. Mompart, “Fermionic transformation rules for spatially filtered light beams in conical refraction,” in Complex Light and Optical Forces V, D. L. Andrews, E. J. Galvez, and J. Glückstad, eds., Proc. SPIE7950, 79500D (2011). [CrossRef]
16. T. K. Kalkandjiev and M. A. Bursukova, “Conical refraction: an experimental introduction,” in Photon Management III, J. T. Sheridan and F. Wyrowski, eds., Proc. SPIE6994, 69940B (2008). [CrossRef]