## Abstract

Influence of the anisotropic tensorial electro-optic effect of *LiNbO*
_{3}:Fe photorefractive defocusing medium on propagation of a vortex beam is numerically and experimentally investigated. Characteristic behaviors are depicted by varying light polarization, sign of vortex angular momentum and propagation directions.

© 2008 Optical Society of America

## 1. Introduction

Vortices are entities that can be found in various physical systems with linear and nonlinear properties. In the field of optics, vortices are associated with phase singularities carried by optical beams [1]. The past decade has seen a renewed interest on optical vortices and various applications have been suggested, such as microlithography [2], optical trapping of particles and transfer of orbital angular momentum [3, 4] or soliton algebra [5]. Especially, propagation of optical vortex in nonlinear medium can lead to soliton vortex formation when beam diffraction is exactly compensated by a self-defocusing nonlinearity. For example soliton vortices have been demonstrated in medium with isotropic saturable thermal nonlinearity [6] or anisotropic photorefractive nonlinearity [7]. Vortex propagation in photorefractive medium was studied in detail [8, 9] and it was demonstrated, in particular, that intensity distribution due to the anisotropic nonlinear effect depends on the traveling vortex topological charge (angular momentum)[10, 11]. However only extraordinarily polarized light was considered in the literature which overlooks the influence of the directions of polarization and propagation of light in anisotropic materials.

In this paper we demonstrate significative influence of these parameters when propagation of light is performed in *LiNbO*
_{3}:Fe crystals. Influence of sign of angular momentum is also investigated. Experimental results are confirmed by a time-dependent three dimensional numerical model which primarily relies on the calculation of the charge-distribution in the photorefractive medium to deduce both transverse components of the space-charge field. The tensorial character of the electro-optic effect is then taken into account for computation of the refractive index perturbation. Numerical simulations and experiments reveal the contribution of the anisotropic properties of the material with respect to the vortex angular momentum and polarization of the light.

## 2. Theoretical background and numerical model

Vortex beam presents a ring-shaped intensity profile due to an helicoidal wavefront in the transverse plane that forms a singularity in the beam center. The optical phase varies from 0 to 2*mπ* around the central singularity, *m* represents the topological charge of the vortex. Electric field amplitude of a vortex beam is given by [12]:

where *r⃗*
_{⊥} represents the coordinate in the transverse plane and z the position along the propagation axis, *m* is the topological charge of the vortex, *θ* is the azimuth angle and
$\overrightarrow{A}=\overrightarrow{{A}_{0}}{\left(\frac{{r}_{\perp}}{{\omega}_{0}}\right)}^{\mid m\mid}\mathrm{exp}\left(\frac{-{r}_{\perp}^{2}}{2{\omega}_{0}^{2}}\right)$
is the optical field amplitude. To calculate the space charge field *E⃗* induced by the light intensity *I*, we solve the well known photorefractive system of equations [8, 13]:

where *I _{d}* is the equivalent dark intensity,

*N*,

_{A}*N*,

_{D}*N*

^{+}

_{D}and

*N*are respectively the densities of shallow acceptors, deep donors, deep ionized donors and free electrons.

_{e}*s*is the photoexcitation coefficient,

*γ*is the recombination constant, k

_{B}is the Boltzmann constant,

*µ*is the electron mobility and

*T*is the temperature. Photovoltaic current directed along

*LiNbO*

_{3}c-axis (i.e Z axis) is considered throughout

*β*which is a polarization dependent component of the photogalvanic tensor. [

_{ph}*ε*] is the static dielectric tensor,

*ρ*is the charge density and

*J⃗*is the current density. From this set of equations usual assumptions are made. Acceptors density is assumed to be greater than the free electrons density (

*N*≫

_{A}*N*). Moreover generation time of conducting electrons is neglected with respect to the characteristic evolution time of space charge, $\left(\frac{\partial {N}_{D}^{+}}{\partial t}=0\right)$ . Following these hypothesis and Eq.(2), Eq.(4) and Eq.(5) we obtain :

_{e}$$\phantom{\rule{.2em}{0ex}}\phantom{\rule{.2em}{0ex}}\phantom{\rule{.2em}{0ex}}\phantom{\rule{.2em}{0ex}}\phantom{\rule{.2em}{0ex}}\phantom{\rule{.2em}{0ex}}\phantom{\rule{.2em}{0ex}}\phantom{\rule{.2em}{0ex}}\phantom{\rule{.2em}{0ex}}\phantom{\rule{.2em}{0ex}}+\xi {E}_{\mathrm{ph}}\overrightarrow{\nabla}\left[\left({\tilde{N}}_{D}-{\tilde{N}}_{D}^{+}\right)I\right]\overrightarrow{c}\}$$

*Ñ _{e}*,

*Ñ*and

_{D}*Ñ*

^{+}

_{D}are free electrons, donors and ionized donors densities normalized to

*N*,

_{A}*Ñ*

^{+}

_{D}=1+

*$\tilde{\rho}$*where

*$\tilde{\rho}$*is the space charge density normalized to e

*N*, $\xi =\frac{s}{\gamma {N}_{A}}$ and ${E}_{\mathrm{ph}}=\frac{{\beta}_{\mathrm{ph}}\gamma {N}_{A}}{e\mu s}$ is the photovoltaic field.

_{A}To deduce the space charge distribution evolution induced by the vortex intensity the above system of equations is solved using a numerical iterative method. In conjunction, space charge field distribution is deduced by integration of Eq.(4) that gives the electrical field produced by charges distribution *ρ*(*r⃗*)d*V* in the medium volume *V* (d*V* being an elementary volume) which is easier to solved than Eq.(3).

In this model discrete Fourier Transform is used to resolve the space charge field in 3D. Two configurations are considered. In the first case vortex beam propagates along X-crystallographic *LiNbO*
_{3} axis while in the second case propagation along Y axis is analyzed. In both case propagation is perpendicular to the c-axis to benefit from a large photovoltaic effect but the refractive index perturbation Δ*n* induced by the space charge field differs. Furthermore, the tensorial nature of the electro-optic effect gives rise to an anisotropic space charge field. For 3*m* group, like *LiNbO*
_{3}, the refractive index modulation for ordinary (X or Y-polarized) and extraordinary (Z-polarized) beams are given by [14]. For propagation along X axis:

For propagation along Y axis, Δ*n _{Z}* for extraordinary polarization is identical to Eq.(11) while for X-polarized beam it yields to:

*r*
_{22}, *r*
_{13} and *r*
_{33} are the electro-optical coefficients and n_{o} and n_{e} are respectively the ordinary and extraordinary indices. E_{Y} and E_{Z} are the components of the space charge field *E⃗* (*r⃗*). Note that E_{X} does not play a significative role for the chosen configurations. Vortex propagation can now be calculated for the different situations of interest by classical split-step Fourier method taking into account the time dependent photorefractive effect. For the numerical study we consider a CW 200*µ*W vortex beam at 473 nm focused to a 26*µ*m FWHM at the entrance of a 9 mm long *LiNbO*
_{3}:Fe crystal. The input and output beams are presented on Fig. 1(a) and Fig. 1(b) in linear regime. The profiles are identical whatever the sign of angular momentum and the directions of propagation and polarization. The maximum intensity to dark irradiance ratio is set to
$\frac{{I}_{max}}{{I}_{d}}=260$
. Photovoltaic field is set to *E _{ph}*=-7,7.10

^{4}

*V*.cm

^{-1},

*N*=5.10

_{D}^{16}cm

^{-3}, and $\frac{{N}_{D}}{{N}_{A}}=1.1$ and other parameters are extracted from [15]. Vortex topological charges

*m*=1 and

*m*=-1 for both polarizations are considered.

The first studied case corresponds to the familiar configuration of an extraordinary polarized vortex [10]. Results, independent of propagation direction, are presented in Fig. 2 at a typical instant of the formation process. The intensity distribution of the vortex beam and the corresponding refractive index modulation at the output of the crystal are depicted for a topological charge *m*=+1 (Figs. 2(a,b)) and for a topological charge *m*=-1 (Figs. 2(c,d)). We specify that the phase of a positively charged vortex rotates anti-clockwise as the beam travels. Because of the intrinsic anisotropy of the photorefractive effect the beam spreads along the c-axis thanks to the deeper index modulation present in the upper and lower part of the beam as depicted in Figs. 2(b,d). Elliptical shapes of vortex and dark core are observed for both topological charges. However the sign of the topological charge impacts the light distribution. Indeed two bright areas whose location are dependent on vortex charge are present on each side of the beam. Symmetrical light distribution relative to the c-axis is observed for opposite topological charge (Fig. 2(c)).

We now consider the case of an ordinary polarized vortex. Figures 3(a,b) shows the key configuration of a beam propagating in the X-axis direction. The beam ellipticity is less pronounced compare to the case of an extraordinary polarized beam due to a weaker electro-optical coefficient *r*
_{13} and consequently a smaller photorefractive index change. But most importantly the two bright spots arrangement totally differs. While they are located on a quasi horizontal diameter for a positively charged vortex (Fig. 3(a)) they instead are titled approximatively at 30 degrees for a negatively charged vortex (Fig. 3(c)). Consequently, the intensity distributions for both charges are no more symmetric relative to the c-axis. The reason of this revealed behavior holds in the index distribution presented in Figs. 3(b,d). Indeed index modulation now depends on both transverse components *E _{Y}* and

*E*of the space charge field (see Eq.(10)).

_{Z}*E*component has a symmetric distribution relative to the Y-axis while the

_{Z}*E*component is anti-symmetric relative to both Y and Z-axis. Moreover since electro-optical coefficients

_{Y}*r*

_{13}and

*r*

_{22}have similar values it yields that conjunction of both transverse space charge field components induces a complex refractive index distribution. To highlight the major role of the space charge field transverse component perpendicular to c-axis, a Y-propagating ordinary polarized vortex is also analyzed. Note that Δ

*n*only depends on

_{X}*E*(Eq.(12)). Calculated vortex intensity distribution is represented in Fig. 4 with corresponding index modulation. Comparison of light distributions for

_{Z}*m*=+1 (Fig. 4(a)) and

*m*=-1 (Fig. 4(c)) shows that they are symmetric with respect to the c-axis like for extraordinary polarized vortices depicted in Fig. 2.

## 3. Experimental results

Figure 5 shows the experimental setup. A CWlaser beam at 473 nm is enlarged and encoded by reflection on a phase hologram. Vortex topological charge sign is chosen by spatially selecting the +1 or -1 diffraction order. Polarization and intensity of the vortex are controlled by means of a polarizer and a half wave plate. The vortex beam is focused to form a 26*µ* m (FWHM) spot at the input face of a *LiNbO*
_{3}:Fe sample (Fig. 6(a)). The 9 mm square sample is cut from a 1 mm thick c-oriented iron doped *LiNbO*
_{3} wafer allowing vortex propagation either along the X or Y axis. Exit face of the crystal is imaged onto a CCD camera. Background illumination of the sample with a white light source is used to adjust the intensity to a dark irradiance ratio
$\frac{I}{{I}_{d}}\simeq 260$
like in numerical simulations. Beam power is set to 200*µ*W.

In linear regime a circular 150*µ*m FWHM vortex beam (Fig. 6(b)) is present at the exit face of the crystal. To start with, extraordinary polarized vortices are used for both sign of the topological charge for X-propagation. Results are presented in Fig. 6(c) and 6(d) for respectively *m*=+1 and *m*=-1 for an induction time of 1000 seconds sufficient to observe beam reshaping. Vortex stretching in the c-axis direction, peripheral brighter areas and strong ellipticity of the dark core are observed in accordance with numerical prediction. In addition, results for positive topological charge can be retrieved from negative topological charge by the symmetry relative to the c-axis of light distribution and conversely. Measurements for Y-propagating shows identical behavior.

Figure 7 depicts experimental results for ordinary polarized vortices. Anisotropic diffraction of the beam forming two bright lobes and slight elliptical shape of the dark core can be observed. More interestingly, for propagation along X-axis, light distribution for a positive topological charge (Fig. 7(b)) presents a vertical symmetry while for a negative topological charge light distribution is clearly rotated (Fig. 7(d)). Such an asymmetric dependence versus topological charge is predicted from our numerical model (see Figs. 3(a,c)) and along with results in Fig. 6 it demonstrates the influence of the electro-optic effect anisotropy on propagation of an ordinary polarized beam with respect to the sign of the angular momentum.

Finally, for propagation along Y-axis, influence of the topological charge on the beam distribution tilt is exhibited in Fig. 7(a,c). For such configuration vortex is slightly tilted throughout the nonlinear photorefractive effect. As predicted by theory and contrary to X-propagation, vertical symmetry is observed between *m*=+1 and *m*=-1 charged vortex distributions. Results depicted in Fig. 7 clearly reveal that each pair angular momentum-propagation direction gives rise to a different intensity distribution.

## 4. Conclusion

Impact of the electro-optic effect on propagation of vortex beam in photorefractive media is established. Demonstration is realized in self-defocusing *LiNbO*
_{3}:Fe photorefractive crystal. The influence of tensorial nature of electro-optic effect gives rise to a beam pattern characteristic of the polarization-momentum pair. Propagation along *LiNbO*
_{3} X-crystallographic axis is found to be adequate for the demonstration. A full 3D photorefractive numerical model is used to successfully model the experimental observations. It shows that the observed phenomena is due to both the anisotropy of the space charge field and the rotating phase carried by the vortex. Influence of the two transverse components of the space charge field have to be taken into account to understand the vortex perturbation.

Note that this fundamental interaction effect can be viewed as an analysis technique. Indeed for given crystal configuration, polarization and momentum of vortex can be deduced from observed intensity pattern or otherwise, for given vortex properties, crystal orientation could be assessed.

## Acknowledgments

Authors would like to thanks Bruno Wacogne for preparing the holographic mask.

## References and links

**1. **A. S. Desyatnikov, Yu S. Kivshar, and L. Torner, “Optical vortices and Vortex Solitons” in *Progress in Optics*, Vol. 47, Ed. E. Wolf (Elsevier, Amsterdam, 2005)

**2. **Z. Jaroszewicz and A. Kolodziejczyk, “Zone plates performing generalized Hankel transforms and their metrological applications,” Opt. Commun. **102**, 391–396 (1993). [CrossRef]

**3. **H. He, M. E. J. Friese, N. R. Heckenberg, and H. Rubinsztein-Dunlop, “Direct observation of transfer of angular momentum to absorptive particles from a laser beam with a phase singularity,” Phys. Rev. Lett. **75**, 826–829 (1995). [CrossRef]

**4. **H. He, N. R. Heckenberg, and H. Rubinsztein-Dunlop, “Optical particle trapping with higher-order doughnut beams produced using high efficiency computer generated holograms,” J. Mod. Opt. **42**, 217–223 (1995). [CrossRef]

**5. **S. Minardi, G. Molina-Terriza, P. Di Trapani, J. P. Torres, and L. Torner, “Soliton algebra by vortex-beam splitting,” Opt. Lett. **26**, 1004–1006 (2001). [CrossRef]

**6. **A. Dreischuh, G. G. Paulus, F. Zacher, F. Grasbon, and H. Walther, “Generation of multiple-charged optical vortex solitons in a saturable nonlinear medium,” Phys. Rev. E **60**, 6111–6117 (1999). [CrossRef]

**7. **Z. Chen, M. Feng Shih, M. Segev, D. W. Wilson, R. Muller, and P. D. Maker, “Steady-state vortex-screening solitons formed in biased photorefractive media,” Opt. Lett. **22**, 1751–1753 (1997). [CrossRef]

**8. **M. Segev, B. Crosignani, A. Yariv, and B. Fischer, “Spatial solitons in photorefractive media,” Phys. Rev. Lett. **68**, 923–926 (1992). [CrossRef]

**9. **M. Segev, G. C. Valley, B. Crosignani, P. DiPorto, and A. Yariv, “Steady-State Spatial Screening Solitons in Photorefractive Materials with External Applied Field,” Phys. Rev. Lett. **73**, 3211–3214 (1994). [CrossRef]

**10. **A. V. Mamaev, M. Saffman, and A. A. Zozulya, “Time-dependent evolution of an optical vortex in photorefractive media,” Phys. Rev. A **56**, 1713–1716 (1997). [CrossRef]

**11. **A. V. Mamaev, M. Saffman, and A. A. Zozulya, “Vortex evolution and bound pair formation in anisotropic nonlinear optical media,” Phys. Rev. Lett. **77**, 4544–4547 (1996). [CrossRef]

**12. **I. V. Basistiy, M. S. Soskin, and M. V. Vasnetsov, “Optical wavefront dislocations and their properties,” Opt. Commun. **119**, 604–612 (1995). [CrossRef]

**13. **V. L. Vinetskii and N. V. Kukhtarev, “Wave front convolution in 4-wave interaction inmedia with nonlocal nonlinearity,” Sovi. Phys. JETP Letters **30**, 6 (1979).

**14. **P. Yeh, *Introduction to photorefractive nonlinear optics* (Wiley-Interscience, New York, 1993).

**15. **M. Simon, S. Wevering, K. Buse, and E. Krätzig,“The bulk photovoltaic effect of photorefractive *LiNbO*_{3}:Fecrystals at high light intensities,”J. Phys. D **30**, 144–149 (1997). [CrossRef]