Nonlinear dynamics of two-color optical vortices in lithium niobate crystals

We study experimentally the nonlinear dynamics of two-color optical vortex beams in the presence of second-harmonic generation combined with the effects of photoand thermal refraction, as well as selfand induced-phase modulation. We use an iron-doped lithium niobate crystal as a nonlinear medium for the vortex propagation and observe experimentally, depending on the laser wavelength, a decay of a double-charge vortex, splitting and reshaping of background beam, pattern formation, and controllable nonlinear rotation of a vortex pair. © 2008 Optical Society of America OCIS codes: (190.4420), Nonlinear optics, transverse effects in; (190.0190), Nonlinear optics; (190.5330), Photorefractive nonlinear optics; (190.2620), Frequency conversion; (190.5940) Self-action effects. References and links 1. D. G. Grier, “A revolution in optical manipulation,” Nature 424, 810–816 (2003). 2. G. Foo, D. M. Palacios, and G. A. Swartzlander, “Optical vortex coronagraph,” Opt. Lett. 30, 3308–3310 (2005). 3. S. Furhapter, A. Jesacher, S. Bernet, and M. Ritsch Marte, “Spiral interferometry,” Opt. Lett. 30, 1953–1955 (2005). 4. G. Molina Terriza, J. P. Torres, and L. Torner, “Twisted photons,” Nature Physics 3, 305–310 (2007). 5. L. Allen, M. J. Padgett, and M. Babiker, “The orbital angular momentum of light,” Prog. Opt. 39, 291–372 (1999). 6. M. S. Soskin and M. V. Vasnetsov, “Singular optics,” Prog. Opt. 42, 219–276 (2001). 7. A. S. Desyatnikov, Yu. S. Kivshar, and L. Torner, “Optical vortices and vortex solitons,” Prog. Opt. 47, 291–391 (2005). 8. W. J. Firth and D. V. Skryabin, “Optical solitons carrying orbital angular momentum,” Phys. Rev. Lett. 79, 2450– 2453 (1997). 9. D. V. Skryabin and W. J. Firth, “Dynamics of self-trapped beams with phase dislocation in saturable Kerr and quadratic nonlinear media,” Phys. Rev. E 58, 3916–3930 (1998). 10. L. Torner and D. V. Petrov, “Azimuthal instabilities and self-breaking of beams into sets of solitons in bulk second-harmonic generation,” Electron. Lett. 33, 608–610 (1997). 11. L. T. Vuong, T. D. Grow, A. Ishaaya, A. L. Gaeta, G. W. Hooft, ’t, E. R. Eliel, and G. Fibich, “Collapse of optical vortices,” Phys. Rev. Lett. 96, 133901–4 (2006). 12. V. Tikhonenko, J. Christou, and B. Luther-Davies, “Spiraling bright spatial solitons formed by the breakup of an optical vortex in a saturable self-focusing medium,” J. Opt. Soc. Am. B 12, 2046–2052 (1995). 13. Yu. S. Kivshar and B. Luther-Davies, “Dark optical solitons: physics and applications,” Phys. Rep. 298, 81–197 (1998). 14. Yu. S. Kivshar and G. P. Agrawal, Optical Solitons: From Fibers to Photonic Crystals (Academic Press, San Diego, 2003). 15. D. Rozas, Z. S. Sacks, and G. A. Swartzlander, “Experimental observation of fluidlike motion of optical vortices,” Phys. Rev. Lett. 79, 3399–3402 (1997). #91939 $15.00 USD Received 22 Jan 2008; revised 31 Mar 2008; accepted 31 Mar 2008; published 3 Apr 2008 (C) 2008 OSA 14 April 2008 / Vol. 16, No. 8 / OPTICS EXPRESS 5406 16. E. B. Sonin, “Vortex oscillations and hydrodynamics of rotating superfluids,” Rev. Mod. Phys. 59, 87–155 (1987). 17. D. Rozas, C. T. Law, and G. A. Swartzlander, “Propagation dynamics of optical vortices,” J. Opt. Soc. Am. B 14, 3054–3065 (1997). 18. D. Neshev, A. Dreischuh, M. Assa, and S. Dinev, “Motion control of ensembles of ordered optical vortices generated on finite extent background,” Opt. Commun. 151, 413–421 (1998). 19. G. Indebetouw, “Optical vortices and their propagation,” J. Mod. Opt. 40, 73–87 (1993). 20. F. S. Roux, “Dynamical behavior of optical vortices,” J. Opt. Soc. Am. B 12, 1215–1221 (1995). 21. Yu. S. Kivshar, J. Christou, V. Tikhonenko, B. Luther-Davies, and L. M. Pismen, “Dynamics of optical vortex solitons,” Opt. Commun. 152, 198–206 (1998). 22. G. H. Kim, H. J. Lee, L. U. Kim, and H. Suk, “Propagation dynamics of optical vortices with anisotropic phase profiles,” J. Opt. Soc. Am. B 20, 351–359 (2003). 23. J. Christou, V. Tikhonenko, Yu. S. Kivshar, and B. LutherDavies, “Vortex soliton motion and steering,” Opt. Lett. 21, 1649–1651 (1996). 24. B. Luther-Davies, R. Powles, and V. Tikhonenko, “Nonlinear rotation of 3-dimensional dark spatial solitons in a gaussian laser-beam,” Opt. Lett. 19, 1816–1818 (1994). 25. A. Berzanskis, A. Matijosius, A. Piskarskas, V. Smilgevicius, and A. Stabinis, “Conversion of topological charge of optical vortices in a parametric frequency converter,” Opt. Commun. 140, 273–276 (1997). 26. I. V. Basistiy, V. Y. Bazhenov, M. S. Soskin, and M. V. Vasnetsov, “Optics of light-beams with screw dislocations,” Opt. Commun. 103, 422–428 (1993). 27. K. Dholakia, N. B. Simpson, M. J. Padgett, and L. Allen, “Second-harmonic generation and the orbital angular momentum of light,” Phys. Rev. A 54, R3742–R3745 (1996). 28. M. S. Soskin and M. V. Vasnetsov, “Nonlinear singular optics,” Pure Appl. Opt. 7, 301–311 (1998). 29. A. Dreischuh, G. G. Paulus, F. Zacher, F. Grasbon, D. Neshev, and H. Walther, “Modulational instability of multiple-charged optical vortex solitons under saturation of the nonlinearity,” Phys. Rev. E 60, 7518–7524 (1999). 30. S. M. Saltiel, A. A. Sukhorukov, and Yu. S. Kivshar, “Multistep parametric processes in nonlinear optics” Prog. Opt. 47, 1–73 (2005). 31. K. H. Yang, J. R. Morris, P. L. Richards, and Y. R. Shen, “Phase-matched far-infrared generation by optical mixing of dye laser beams,” Appl. Phys. Lett. 23, 669–671 (1973). 32. G. D. Miller, R. G. Batchko, W. M. Tulloch, D. R. Weise, M. M. Fejer, and R. L. Byer, “42%-efficient single-pass cw second-harmonic generation in periodically poled lithium niobate,” Opt. Lett. 22, 1834–1836 (1997). 33. K. Noguchi, O. Mitomi, and H. Miyazawa, “Millimeter-wave Ti : LiNbO3 optical modulators,” J. Lightwave Technol. 16, 615–619 (1998). 34. G. C. Valley, M. Segev, B. Crosignani, A. Yariv, M. M. Fejer, and M. C. Bashaw, “Dark and bright photovoltaic spatial solitons,” Phys. Rev. A 50, R4457–R4460 (1994). 35. D. Runde, S. Brunken, C. E. Rutter, and D. Kip, “Integrated optical electric field sensor based on a Bragg grating in lithium niobate,” Appl. Phys. B 86, 91–95 (2007). 36. F. S. Chen, J. T. LaMacchia, and D. B. Fraser, “Holographic storage in lithium niobate,” Appl. Phys. Lett. 13, 223 (1968). 37. X. An, D. Psaltis, and G. W. Burr, “Thermal fixing of 10,000 holograms in LiNbO3 : Fe,” Appl. Opt. 38, 386–393 (1999). 38. D. N. Nikogosyan, Properties of Optical and Laser-Related Materials: A Handbook (Wiley, Chichester, UK, 1997). 39. L. Arizmendi, “Photonic applications of lithium niobate crystals,” Phys. Status Solidi A 201, 253–283 (2004). 40. F. Pettazzi, V. Coda, M. Chauvet, and E. Fazio, “Frequency-doubling in self-induced waveguides in lithium niobate,” Opt. Commun. 272, 238–241 (2007). 41. R. G. Batchko, M. M. Fejer, R. L. Byer, D. Woll, R. Wallenstein, V. Y. Shur, and L. Erman, “Continuous-wave quasi-phase-matched generation of 60 mW at 465 nm by single-pass frequency doubling of a laser diode in backswitch-poled lithium niobate,” Opt. Lett. 24, 1293–1295 (1999). 42. B. Crosignani, M. Segev, D. Engin, P. Diporto, A. Yariv, and G. Salamo, “Self-trapping of optical beams in photorefractive media,” J. Opt. Soc. Am. B 10, 446–453 (1993). 43. G. C. Valley, M. Segev, B. Crosignani, A. Yariv, M. M. Fejer, and M. C. Bashaw, “Dark and bright photovoltaic spatial solitons,” Phys. Rev. A 50, R4457–R4460 (1994). 44. M. Taya, M. C. Bashaw, M. M. Fejer, M. Segev, and G. C. Valley, “Observation of dark photovoltaic spatial solitons,” Phys. Rev. A 52, 3095–3100 (1995). 45. D. N. Neshev et al., “Nonlinear spectral-spatial control and localization of supercontinuum radiation,” Phys. Rev. Lett. 99, 123901–4 (2007). 46. Z. Chen, M. Segev, D. W. Wilson, R. E. Muller, and P. D. Maker, “Self-trapping of an optical vortex by use of the bulk photovoltaic effect,” Phys. Rev. Lett. 78, 2948–2951 (1997). 47. R. DeSalvo, A. A. Said, D. J. Hagan, E. W. Van Stryland, and M. SheikBahae, “Infrared to ultraviolet measurements of two-photon absorption and n(2) in wide bandgap solids,” IEEE J. Quantum Electron. 32, 1324–1333 (1996). 48. D. H. Jundt, M. M. Fejer, and R. L. Byer, “Optical-properties of lithium-rich lithium-niobate fabricated by vapor #91939 $15.00 USD Received 22 Jan 2008; revised 31 Mar 2008; accepted 31 Mar 2008; published 3 Apr 2008 (C) 2008 OSA 14 April 2008 / Vol. 16, No. 8 / OPTICS EXPRESS 5407 transport equilibration,” IEEE J. Quantum Electron. 26, 135–138 (1990). 49. S. Lan, M. F. Shih, G. Mizell, J. A. Giordmaine, Z. G. Chen, C. Anastassiou, J. Martin, and M. Segev, “Secondharmonic generation in waveguides induced by photorefractive spatial solitons,” Opt. Lett. 24, 1145–1147 (1999). 50. A. D. Boardman, W. Ilecki, and Y. Liu, “Spatial solitons in a photorefractive medium sustaining second-harmonic generation,” J. Opt. Soc. Am. B 19, 832–838 (2002). 51. J. F. Reintjes, Nonlinear Optical Parametric Processes in Liquids and Gasses (Academic Press, New York, 1984). 52. S. Dinev and A. Dreischuh, “Sum-frequency generation in the xuv,” Opt. Quantum Electron. 23, 91–97 (1991). 53. A. Dreischuh, U. Reiter Domiaty, D. Gruber, L. Windholz, and S. Dinev, “Nonlinear alignment between conical emissions generated in a four-wave parametric mixing process,” Appl. Phys. B 66, 175–180 (1998). 54. M. F. Shih, Z. G. Chen, M. Mitchell, M. Segev, H. Lee, R. S. Feigelson, and J. P. Wilde, “Waveguides induced by photorefractive screening solitons,” J. Opt. Soc. Am. B 14, 3091–3101 (1997). 55. S. Lan, C. Anastassiou, M. Segev, M. Shih, J. A. Giordmaine, and G. Mizell, “Tuning of second-harmonic generation in waveguides induced by photorefractive spatial solitons,” Appl. Phys. Lett. 77, 2101–2103 (2000). 56. S. Orlov, A. Yariv, and M. Segev, “Nonlinear self-phase matching of optical second harmonic generation in lithium niobate,” Appl. Phys. Lett. 68, 1610–1612 (1996). 57. C. Lou, J. Xu, H. Qiao, X. Zhang, Y. Chen, and Z. Chen, “Enhanced second-harmonic generation by means of high-power confinement in a ph


Introduction
The study of optical vortices (OVs) has received a special attention in recent years due to a variety of potential applications of the vortex beams including particle micro-manipulation [1], imaging [2], interferometry [3], and quantum information [4].An optical vortex is associated with a light beam whose phase varies in a screw-type way around an isolated singular point.This characteristic helical phase of an OV is described by the multiplier exp(imθ ), where θ is the azimuthal coordinate, and the integer number m stands for the vortex topological charge.The helical phase results in a non-zero orbital angular momentum of the light beam [5,6] which strength depends on the vortex charge.
The optical vortices are associated with doughnut beams or Laguerre-Gaussian (LG) modes.As such they are eigenmodes of the free-space equation and propagate without any change of their shape.However, in nonlinear media the vortex propagation is affected by the beam induced change of the refractive index as well as by the phase and index gradients induced by other beams [7].While in nonlinear media with the self-focusing nonlinearity the vortex beam becomes unstable due to the azimuthal modulational instability [8,9,10,11,12], in selfdefocusing media a doughnut beam can form a dark vortex soliton with a strongly localized self-trapped vortex core [7,13,14].Analogously to point vortices in fluid dynamics [15,16], optical vortices of the same charge with localized cores orbit each other at the rate which is inversely proportional to the squared vortex separation distance [17,18].The nonlinearity of the material, being the key for manipulation of the vortex dynamics [17,19,20,21,22], can alter dramatically this type of the vortex fluid-like motion, including the vortex steering [23] and the rotation of a vortex pair [18,24].One of the important nonlinear processes which affects strongly the vortex propagation in a nonlinear medium is the parametric frequency conversion.In this process the orbital angular momentum of the vortex beam is conserved.Thus, when the fundamental frequency (FF) beam carries an orbital angular momentum, the generated harmonic beams have to carry an orbital angular momentum as well, resulting in the generation of a two-color optical vortex.In the general case of the frequency mixing of two LG modes with the topological charges m 1 and m 2 , the generated vortex has a topological charge of m 3 = m 1 + m 2 [25].Therefore, in the parametric process such as the second-harmonic generation (SHG) the topological charge of an optical vortex should double [5,7,26,27,28].Higher-charge vortices however, are structurally unstable and weak perturbation can initiate a decay of the second-harmonic (SH) vortex into a pair of single-charge vortices [29].Such perturbations can arise from the spatial overlap of the first-and second-harmonic beams and the complex field gradients due to back conversion from the SH beam to the FF beam.
In typical crystals used for SHG, a variety of additional nonlinear processes, such as photorefraction, thermal focusing/defocusing, and Kerr-type self-phase modulation, can influence the nonlinear propagation of a vortex beam.When ultrashort pulses of high-peak intensities are used, the frequency conversion processes can also be accompanied by cascaded parametric [30] and higher-order nonlinear processes (e.g.self-and induced-phase modulation).The rich dynamics of such combined nonlinear interactions calls for a systematic study of the vortex dynamics in the process of SHG accompanied by the action of additional effects.
In this paper we study experimentally the nonlinear dynamics of two-color optical vortices generated in the process of non-phase matched SHG and affected by the combined effect of several different nonlinear processes such as thermal focusing, photorefractive effect, and selfand induced phase modulation.We use an iron-doped lithium niobate crystal as a nonlinear medium for the vortex propagation and observe a variety of interesting nonlinear effects including a decay of a double-charge vortex, splitting and reshaping of the background beam, pattern formation, and nonlinear rotation of the a vortex pair.

Experimental arrangements and parameters
In our experiments we generate OVs by passing the laser beam through a computer generated hologram (CGH) and selecting the '+1' or '-1' diffraction order after the hologram as shown in Fig. 1(top).The hologram is fabricated by etching of the corresponding vortex interference pattern onto a glass substrate.The holographic grating has a period of 80 μm and a diffraction efficiency of approximately 30% in the first order.The selected diffraction order is then focused by a lens onto an iron-doped Lithium Niobate (LiNbO 3 ) crystal with dimensions 3 × 15 × 20 mm.The crystal is cut with c-axis parallel to its longest side [Fig.1(top)].The output of the crystal is then monitored with a CCD camera and digitized on a computer.
The LiNbO 3 is a negative uniaxial crystal that exhibits a combination of electro-optical, piezoelectric, pyroelectric, and photorefractive (PR) properties which makes it attractive for applications in frequency conversion [31,32], optical modulators [33], optically-written waveguides [34], sensors [35], and holographic recording [36,37].The LiNbO 3 is a well established optical material [38,39] with good nonlinear properties.It is most commonly used as a medium with quadratic nonlinearity, but other effective third-order nonlinearities, including photorefractive, thermal, and pure electronic Kerr-type nonlinearity are well pronounced.
The crystalline structure of the LiNbO 3 is 3m thus it lacks inversion symmetry allowing the process of SHG.Obtaining proper phase matching conditions for wavelengths in the visible and near infra-red however, is practically impossible in LiNbO 3 .For example, for degenerate phase matching of SHG at 1064 nm, temperatures of the order of -5 • C are required [40].At other temperatures, the birefringent phase matching leads to spatial walk-off between the fundamental frequency and the second harmonic.To satisfy the phase-matching conditions in the degenerate case of propagation along the optical axis, the technique of periodic poling is exploited to ensure quasi-phase-matching.The LiNbO 3 appears well suitable for periodic poling and various techniques for such poling have been demonstrated [39,41].
LiNbO 3 crystals doped with certain metal ions (e.g., iron or copper) exhibit enhanced photovoltaic effect at photosensitive wavelengths.The photovoltaic effect is related to photoexcitation of free electrons and their drift along the optical c-axis.The resulting space-charge field, E pv , creates via electro-optic effect a nonlinear refractive index change, Δn, which is negative, saturable (as a function of light intensity), and anisotropic.When the polarization of the input light is extraordinary (along the c-axis) the largest electro-optic coefficient r 33 contributes to the decrease of n.In this case the photorefractive index change is Δn = −(1/2)n 3 r 33 E pv I/(I + I d ), where I, and I d are the light intensity and dark irradiance, respectively [42,43,44].The spectral sensitivity of the LiNbO 3 is mainly determined by the dependence of absorption on wavelength and can be approximated by the Gaussian function [45] with λ > λ b = 400 nm and λ w = 150 nm.The defocusing photovoltaic nonlinearity in LiNbO 3 can support vortex solitons.The first observation of vortex self-trapping was reported in Fe-doped LiNbO 3 in Ref. [46].LiNbO 3 also exhibits third order Kerr-type nonlinearity, which is positive and approximately three times larger than the nonlinearity of BK7 glass.For a wavelength of 1.06 μm the value of the nonlinear coefficient is n 2 = 9.1 × 10 −20 [m 2 /W] [47].Even though this is a small value, the contribution of the Kerr nonlinearity can become appreciable when high power short laser pulses are used.The Kerr effect is mainly expressed in self-phase modulation of the pulse but, as can be seen below, can dramatically affect vortex propagation.
Due to the non-negligible absorption of our crystal in the visible and near infrared region (400 − 800 nm) [Fig.1(bottom)], part of the laser beam power can be absorbed during the process of linear or two-photon absorption.The absorption of a portion of the laser power leads to heating of the sample and change of its refractive index via thermooptic effect.In the case of LiNbO 3 the thermooptic coefficients are 1.4× 10 −6 K −1 for the ordinary and 39× 10 −6 K −1 for the extraordinary polarization components, respectively at a wavelength of 1.06 μm [48].This thermooptic effect results in an effective increase of the material refractive index and therefore self-focusing effect of the beam.The increased light intensity due to the process of thermal self-focusing on the other hand facilitates the increase of the other nonlinear effects like SHG and Kerr-nonlinearity.
The combination of different nonlinearities in LiNbO 3 has been studied in the terms of SHG enhancement via photorefractive soliton induced waveguides [49].When the SH wave is generated in a crystal with non-negligible photorefraction the beam dynamics is far more rich [50].As seen from Eq. ( 1), depending on the FF wavelength, λ FF , the photorefraction could be initiated by the second harmonic (wavelength λ SH ), or to a different extent, by both the first and second harmonics.As a result, the phase matching for SHG changes in time and space.This effect appears similar to the influence of the Kerr nonlinearity on the four-wave mixing process [51,52,53].
In optical soliton-induced waveguides the nonlinear frequency conversion relies on phase matching between the guided modes, which have different propagation constants than the free propagating waves.This difference in propagation constants enables photorefractive tuning of the SHG process by modifying the properties of the guided modes.The first proposal of this idea (when the waveguides are induced by bright screening solitons) [54], followed by successful SHG [49] and SH wavelength tuning [55] in KNbO 3 was recently confirmed in LiNbO 3 [40].It is also shown [56] that in Fe-doped LiNbO 3 crystals used in photovoltaic regime, nonlinear refractive index change originating from the second-harmonic-induced space-charge field can give rise to SH self-phase matching and effective increase of the SHG efficiency.
The experimental observation of enhanced SH in photovoltaic-dark-soliton-induced waveguide [57] is qualitatively similar to the observations with bright solitons [49,55,56].In the experiment of Ref. [57] the dark soliton is created with a weak continuous-wave background laser beam with an one-dimensional odd dark stripe.Femtosecond pulse train is used as fundamental wave, while both the FF and SH radiation are able to cause photorefraction and to destroy the optically-written waveguide when the dark soliton is absent.
Here we explore a somewhat different geometry with a two-dimensional dark dislocation which is embedded already in the FF beam.In addition, we also employ three different laser sources enabling us to distinguish between different nonlinear phenomena pronounced in the LiNbO 3 crystal.
The parameters of the three different laser sources used in our experiments are summarized in Table 1.The lasers differ in wavelength and power, allowing us to observe independently the influence of the different nonlinear effects on vortex dynamics.In the first experiments, refereed below as Case 1, we use a Ti:Sapphire femtosecond oscillator (Mira 900, Coherent) emitting at a wavelength of 845 nm.In the second case, Case 2, a picosecond mode-locked When the LiNbO 3 crystal is oriented with its c-axis perpendicular to the FF beam polarization (ordinary) the propagation lengths are either l cryst = 3 mm or 15 mm [Fig.1(top)].The estimated pulse walk-off lengths due to group velocity mismatch (GVM) between the FF and the SH are [38] for the three lasers we have three different propagation regimes: L W l cryst for Case 1 (λ FF = 845 nm), L W > l cryst for Case 2 (λ FF = 1064 nm), and L W ∼ l cryst for Case 3 (λ FF = 1464 nm).The long L W in Case 2 is mainly due to the longer duration of the pulse (18 ps).
The second-order dispersion (SOD) and the FF output pulse durations at the exit of the sample are estimated using refractive index data from [38] according to the relations SOD = λ 3 0 l cryst /(2πc 2 )(d 2 n/dλ 2 ).Then the pulse broadening due to SOD can be estimated as valid in a Gaussian pulse shape approximation with τ in being the input full-width pulse duration at half maximum) [58].Our estimates for the pulse broadening (see Table 1) show that in all three cases this effect is weak and can be largely neglected.
For the process of SHG an important parameter is the coherence length, l coh , defined as l coh = π/|Δk|, where Δk = 4πn(ω)/λ FF − 2πn(2ω)/λ SH .As shown in Table 1, in all cases l coh << l cryst , therefore the process of SHG is realised under the condition of strong phase mismatch.The relatively long femtosecond pulses allow us to employ a CGH for the generation of the OV beams without noticeable OV beam smearing [59,60], pulse-front tilt [61] and associated pulse lengthening [62].The problem of topological dispersion [63] (frequency-dependent OV topological charge) was absent as seen in Figs.4(a,d) and the length needed to separate the desired first-order diffracted OV beam was only adjusted depending on the FF laser wavelength.

Characterization of the effective nonlinearity
The input 140 fs pulses broaden negligibly (by 2 − 3 fs) along the LiNbO 3 crystal for both FF wave polarizations.The GVM length is much shorter than the sample length (l cryst = 3 mm), hence the SHG at normal incidence is strongly influenced by the group-velocity mismatch.This is well expressed in Fig. 2(a) which shows the calculated [64] output pulse profiles of the FF and SH fields.The intensities are normalized to their maximal value.From the plot we see that the SH is emitted in two separate pulses [65] one of which is locked to the FF pulse and the other is delayed due to GVM.
In order to check the phase-matching conditions experimentally, we vary the FF beam polarization at the input facet of the crystal.After a polarisation filtering, the ordinary and extraordinary components [SH(o) and SH(e)] of the generated SH are recorded.The strong, nearly symmetric oscillations and the low conversion efficiency [Fig.2(b)] are indicative for the nonphase-matched regime of the SHG.Each individual measurement in Fig. 2(b) is done at a fresh location of the crystal and at a reproducible waiting time of 2 s.According to Eq. ( 1) both the pump and the generated SH are able to initiate photorefraction, since the lower (three orders of magnitude) photosensitivity at the FF wavelength is compensated by the (five orders of magnitude) stronger FF intensity.In Fig. 3 we present the time evolution of the far-field FF (upper rows) and SH beam patterns (lower rows) when excited by a Gaussian FF beam.The crystal c-axis direction is vertical.Fig. 3 (oo-e)-regime.The infrared FF beam (λ FF = 845 nm) becomes elongated along the c-axis and this is stronger pronounced for extraordinary FF beam polarization.The formation of side lobes in the FF beam in this case [see Fig. 3(a), lower row, t=16s] we explain by two mechanisms of the SH beam pattern development [56].
At the initial stage (see Fig. 3, virtual time t = 0 s) the SH beam appears ring-shaped.The relatively fast (of the order of milliseconds), initial isotropic and polarization-independent response in the material indicates refractive index changes of a thermal origin [66] and results in a formation of a bright symmetric ring [not clearly visible in Fig. 3(a,   modulation Δn ∼ 3 × 10 −4 .Later, the photorefractive nonlinearity clearly dominates, resulting in a strong asymmetry in the transverse plane and development of a complex beam structure simultaneously in the FF and SH beams.The dominating photorefractive effect also leads to a total negative index change. When both the FF and SH waves are extraordinary [Fig.3(a)], the strong PR nonlinearity reshapes the ring and a complex pattern is formed in the far field.The pronounced asymmetry in the ±c-direction [see Fig. 3(a), lower row, t = 16 s] is due to the diffusion field [56].
With ordinary polarized FF waves [Fig.3(b)], in the SH wave we observe side-lobes (t = 2 s) developing in a broad (of the order of the frame window, unfortunately not well visible) halolike structure (for t ∼ 4 − 8 s) surrounding the elongated (along c-axis) main part of the beam.
In summary to this subsection, in the Case 1, the nonlinear processes which accompany the SHG are the thermal positive index change at a short millisecond time scale, followed by much stronger photorefraction from both FF and SH fields on a time scale of a few seconds.

Vortex dynamics
As a next step we encode a single-charged OV in the FF beam and focus it in the middle of the LiNbO 3 crystal by a lens with f = 50 mm.As such, we obtained average SH power of 0.42 μW for extraordinary FF beam polarization at average FF power of 120 mW.In Fig. 4(upper row) we show far-field energy density distributions of the SH OV beam at an initial (a), intermediate (b), and final evolution stages (c) for a total time of approximately 30 minutes.The data are recorded with a CCD-camera build-in in a GRENOUILLE frequency-resolved optical gating device (Swamp Optics) used in a spatial mode.The same device and alignment was used to estimate that the material dispersion of the hologram substrate does not increase noticeably the FF pulse width.As seen in Fig. 4

Characterization of the effective nonlinearity
For this experimental case, first we carried out benchmark measurements of the origin of the photorefraction in the crystal and its influence on the SHG phase-matching.Nonlinear selfphase matching of the SHG in LiNbO 3 at 1064 nm has been demonstrated with a Q-switched mode-locked laser emitting 100 ps pulses [56].In this process, the presence of SH changes spatially the refractive index of the material, thus causing a change in the SH phase matching conditions.This change leads to improved conversion efficiency with stronger SH component resulting in a feedback self-matching mechanism.In contrast to Ref. [56], here we use a crystal oriented strictly at a normal incidence (with c-axis perpendicular to the incoming beam) and a mode-locked Nd:YVO 4 laser emitting shorter (18 ps) pulses at a repetition rate of 1.5 MHz [67].Due to this pulse duration the generated SH pulse always overlaps with the FF pulse inside the crystal as shown in Fig. 5(a).The average power of the generated extraordinary SH at the output facet of the crystal (l cryst = 3 mm) is shown in Fig. 5(b) as a function of the FF beam input polarization.With an input average power P FF = 30 mW the ordinary SH always remained below 7 nW, whereas the polarization-dependent extraordinary SH showed pronounced oscillations with input polarisation and reveal that the only possible type of interaction is oo-e.Each measurement is done at a fresh location of the crystal after a waiting time of 5 s.
In Fig. 6 we show the increase of the SH(e) power with time when the FF beam is focused in the center (a) of the sample, on its input (b), and output facets (c).The crystal itself was rotated to use longer beam propagation (l cryst = 15 mm).Since only the focusing conditions are changed from curve to curve, the enhancement of the SH(e) average power with time is a clear signature of photovoltaic self-phase matching [56].In order to determine wheatear the origin of photorefraction is coming from the FF or the SH wave, we performed a test experiment where the beam of a cw Nd:YAG laser (Coherent 1064-100) with three times higher output power (90 mW) was focused on the crystal with the same polarization.Since no signature of cw beam reshaping or fanning was observed, we are convinced that the self-phase-matching effect (Fig. 6) is due to the photovoltaic field induced by the extraordinary SH light which changes (via photorefraction) the refractive indices and phase-matching conditions.Thus for this case, we can conclude that along with the SHG process we have a strong contribution of photorefractive nonlinearity induced entirely by the SH field at 532 nm.

Vortex dynamics
As the next step, we generate an optical vortex in the FF beam and focus it onto the LiNbO 3 crystal with a short focal length lens f = 38 mm.As expected, the generated double-charge vortex in the SH field splits into a pair of two single-charge vortices, which start rotating as the photorefraction influences the refractive index.For the (oo-e) type SHG process, the SH is relatively strong, inducing strong photorefraction.As a result, the SH beam reshapes quickly breaking its symmetry, thus preventing clear identification of the angular positions of the vortices.When the input polarization is changed to 45 • with respect to the c-axis [see Fig. 5(b)] the photovoltaic self-phase matching slows down keeping its monotonic and saturable character.The rotation angle of the OVs vs. time is shown in Fig. 7.The solid curve is a log-normal fit intended to guide the eye only.At the initial stage of evolution the vortex pair rotates quickly (in a few seconds) to approximately 35 • , after which it continues rotating monotonically at a slower rate.The inset in Fig. 7 shows the initial (t = 0; left frame) and final (right frame; t = 30 min) orientation of the OV pair.Similar to the behavior described in Case 1, in the steady state the OVs align along the c-axis of the LiNbO 3 crystal.5. Case 3: λ FF = 1464nm

Characterization of the effective nonlinearity
In the measurements described below we use high-power pulses with the duration of 150 fs at 1464 nm and repetition rate 250 Hz.Even though in this case the average laser power is much lower, the pulse peak power is four orders of magnitude higher than in the previous two cases.
The low average power and the long FF wavelength ensured that neither the FF beam nor the SH beam trigger photorefraction.Thus, in addition to the SHG process we encounter only the third-order Kerr-type nonlinearity.In this experiment we keep the LiNbO 3 crystal at normal incidence with respect to the laser beam (l cryst = 3 mm) and orient it with its c-axis parallel to the FF beam polarization (ee-e type of interaction).Unlike the previous cases, here we use an OV embedded in the FF beam even for characterisation of the nonlinearity.The input OV beam is focused by a f = 150 mm lens inside the LiNbO 3 crystal.After filtering out the FF beam, a second lens of a focal length 200 mm is used to re-collimate the SH onto a CCD camera.The lens, camera, and LiNbO 3 crystal are then translated together (similar to the experiments in Ref. [15]) thus varying the position of the input FF beam waist with respect to the crystal.Additionally, we can record the spectra of the SH beam exiting the crystal by a fibre spectrometer (Ocean Optics, HR2000).
In the SH spectra recorded at different positions of the LiNbO 3 crystal, we observe strong dependence of the SH spectral width, Δλ , ranging from 10.5 nm (at large waist offset, but still detectable SH) to Δλ MAX = 24.1 nm (at the longitudinal position z = 0 inside the crystal).In Fig. 8(a) (solid diamonds) we show the normalized SH spectral width Δλ /Δλ MAX for the half of the recorded data.The reason for this spectral broadening is the third-order nonlinear process of induced-phase modulation (IPM) from the FF pulse to the SH.Due to the mild focusing of the FF beam, we can neglect the linear diffraction inside the 3 mm long crystal.Then by using the fit of the experimental data [Fig.8 curve in Fig. 8(a) visualizes the beam broadening with increasing offset of the crystal from the beam waist which decreases the FF beam intensity.The reduced intensity decreased the effect of the IPM on the SH and consequently the normalized SH spectral width Δλ /Δλ MAX .Since the estimated GVM length L W (e-e) is approximately equal to the sample length and n gr e (ω) < n gr o,e (2ω), the generated SH will be influenced mainly by the trailing wing of the FF pulse.Therefore, in the case of a positive third-order nonlinearity, the IPM-assisted spectral broadening has to be asymmetric.In this case one can expect an induced shift of the SH centerof-gravity wavelength (λ c ) and stronger oscillations in one of the wings of the SH spectra [68,69].In Fig. 8(b,c) such signatures are clearly visible when comparing the SH spectra (solid curves) recorded near [Fig.8(b)] or at the FF beam waist [Fig.8(c)] with the spectrum of the weak SH signal at a large offset (x/L D = 3) from the FF beam waist (dashed curves).Because of the short propagation path-length in the crystal, it is unlikely that the observed red shift Δλ c ≈ 4 nm (Δν ≈ 74 cm −1 ) of the spectrum is due to the intrapulse Raman scattering [70,71].

Vortex dynamics
In the spatial domain, the presence of OV in the FF beam leads to the generation of a doublecharge OV in the SH beam, which splits into a pair of two OVs.The fast response of the electronic Kerr-type nonlinearity allows us to study the process of the vortex rotation (similar to Case 2), but for a positive nonlinearity and with respect to the position of the beam focusing.In Fig. 9(a) we show the retrieved angle of the rotation of an OV pair in the SH field as a function of the longitudinal position z of the crystal.The dashed curve is a sigmoidal fit intended to guide the eye only.The sequence of frames shows the far-field images recorded experimentally at the crystal positions as marked by dashed-blue circles.Despite the focusing Kerr nonlinearity, the azimuthal instability of the vortex beam [11] is not pronounced in our experiments due to a short propagation length and mild focusing inside the crystal.
In order to compare the rotation angle in the nonlinear regime with the rotation due to linear Guoy phase shift, we normalize the longitudinal coordinate to the effective Rayleigh diffraction length L D .In this way, we can add in Fig. 9(a) the solid curve for the calculated linear Guoy phase shift at the distances of interest.For the sake of clarity one has to mention that both the GVM length L W and the Rayleigh diffraction length L D are comparable to the sample length  l cryst .Since a significant OV rotation of ≈ 90 • is achieved over a distance z/L D ∼ 1.5, whereas the Gouy shift over this distance is only 40 • (in agreement with Ref. [17]) the results strongly indicate a nonlinear origin of the vortex pair rotation.As seen in Fig. 9(b), the SH OVs rotation with increasing the offset of the crystal from the FF beam waist is accompanied by changes in their radial positions.To our knowledge, this is the first experimental observation of the rotation of a vortex pair due to the effect of pure Kerr-type nonlinearity.

Conclusions
We have presented a systematic experimental study of the nonlinear dynamics of two-color optical vortices in the presence of non-phase-matched second-harmonic generation and additional nonlinear effects of photo-and thermal refraction, self-and induced-phase modulation.
We have revealed several different scenarios of the vortex evolution.When the medium is photosensitive to both the fundamental and second-harmonic wavelengths (Case 1), the different time-scales of the thermal and photovoltaic nonlinearity lead to an alignment of the two colors of the vortex beam in near field and to the generation of complex spatial patterns in far field.However, when photorefraction is initiated by the generated second-harmonic beam only (Case 2), the nonlinear self-phase matching is accompanied by rotation of a two-color vortex by 90 • .Finally, we have achieved controllable vortex rotation with two-fold nonlinear increase, as well as spectral broadening of the femtosecond vortices due to the pump-induced phase modulation with possible successive femtosecond vortex pulse compression (Case 3).

Fig. 1 .
Fig. 1.Experimental arrangement (top) and absorption of the Fe:LiNbO 3 crystal (bottom).Solid and dashed lines indicate the FF and SH wavelengths corresponding to the three different cases in the experiment.CHG: computer generated hologram, L: focusing lens.

Fig. 2 .
Fig. 2. (a) Calculated [64] pulse profiles of the FF (solid line) and SH (dashed line) normalized to their maximum intensity at the crystal output.(b) Average power of the ordinary and extraordinary second harmonic [SH(o) and SH(e), respectively] vs. angle of rotation of the FF beam polarization.Focusing by f = 50 mm lens, l cryst = 3 mm.Specific polarization components of interacting waves are marked with vertical lines.Solid lines: spline to the experimental points.
0 s)].At this initial stage the Gaussian input FF beam also creates photorefractive carriers with symmetric distribution, however it takes longer time for these charges to separate and create space-charge field.The formation of a single bright ring at this initial stage indicates a positive refractive index Far field Near field

Fig. 4 .
Fig. 4. SH far-field (upper row) and near field (lower row) energy density distributions of the OV beam at an initial (a,d), intermediate (b,e), and final evolution stages (c,f).The arrows indicate the position of the two vortices.

Fig. 5 . 4 .
Fig. 5. (a) Calculated [64] pulse profiles of the FF (solid line) and SH (dashed line) normalized to their maximum intensity at the crystal output.(b) Average power of the extraordinary second harmonic vs. angle of rotation of the FF beam polarization.Focusing by f = 38 mm lens; l cryst = 3 mm).Specific polarization components of interacting waves are marked with vertical lines.Solid line: spline to experimental points.

Fig. 6 .
Fig.6.Growth of the SH(e) mean power with time.Type of interaction: (oo-e).Gaussian FF beam is focused in (a) the center of the sample (l cryst = 15 mm), on its input (b), and output facets (c) with a f = 38 mm AR-coated lens.

Fig. 9 .
Fig. 9. (a) Rotation angle of a OV pair vs. crystal position.Squares -experimental data, dashed curve -sigmoidal fit.Solid curve -calculated rotation angle due to linear Guoy phase shift only.(b) Projection of the vortex trajectories onto the CCD camera.Lower row: experimental far-field images corresponding to the points on the graph marked with circles.

Table 1 .
Parameters of the lasers and the Fe 3+ :LiNbO 3 crystal.GVM -group velocity mismatch, SOD -second order dispersion.
coh , μm 1.57 (ee-e) 131.0 (oo-e) 8.3 (ee-e) laser at 1064 nm was employed, and in the last case, Case 3, we used a Ti:Sapphire amplified system (Clark-MXR 2001) and a parametric amplifier (Light Conversion TOPAS) emitting at a wavelength of 1464 nm.