Thermal lensing-induced bifocusing of spatial solitons in Kerr-type optical media
Thermo-optical effects cause a bifocusing of incoming beams in optical media, due to the birefringence created by a thermal lens that can resolve the incoming beams into two-component signals of different polarizations. We propose a non-perturbative theoretical description of the process of formation of double-pulse solitons in Kerr optical media with a thermally-induced birefringence, based on solving simultaneously the heat equation and the propagation equation for a beam in a one-dimensional medium with uniform heat flux load. By means of a non-isospectral Inverse Scattering Transform assuming an initial solution with a pulse shape, a one-soliton solution to the wave equation is obtained that represents a double-pulse beam which characteristic properties depend strongly on the profile of heat spatial distribution.
pacs:42.65.Tg, 42.70.Gi, 44.10.+i
In dielectric media with thermo-optical
effects, the modulations of incoming beams can lead to a wavefront
distortion dist1 (); dist3 (); dist4 ()
reflecting their instability. Generally
this instability gives rise to a depolarization of a high-power
field hans (); menzel1 (); moshe1 (); moshe2 (); wang (); resolv2 ()
due to a thermally-induced birefringence which is
attributed biref1 (); biref2 (); biref3 () to a
change in the refractive index of the medium. It has been
established dist1 () that
this change in the refractive index originates from heat deposition in the
propagation medium, resulting in a space-dependent
temperature gradient (so-called thermal lensing) dist1 (); hans (); koechner (). For
media with linear indices,
the thermal lens leads to a drastic change in the
the beam axis so that the resulting depolarization can strongly degrade
the beam quality requiring thermal lensing compensation. However in nonlinear
media such as Kerr media, the nonlinearity can be a relevant
self-compensaton factor alex () stabilizing
generated by the thermal birefringence.
While several materials exhibiting thermo-optical effects are known in the literature dist1 (), a most investigated one is the solid-state laser Nd:YAG dist1 (); hans (); koechner (); eich (). This material is represented as a rod crystal with a cylindrical geometry, where the change in temperature induces thermal distortion of incoming laser beams resolv2 (). For this particular material, several theoretical attempts have been made to formulate the spatial profile of the temperature gradient along the rod exploiting available experimental data. In particular, in refs. dist1 (); koechner (); eich () it was found that in the cylindrical rod configuration where the heat is generated at a constant rate koechner (); resolv2 (), a quadratic spatial distribution gives a very good description of the experimentally observed birefringence and the resulting beam bifocusing resolv2 (); maik ().
But thermo-optical processes are actually common to a broad class of materials, not just solid-state lasers. Indeed, photonic crystals and optical fibers (including laser fibers) displaying nonlocal thermal and photothermal properties have been considered in the recent past, from both experimental and theoretical points of view chen (); martha (). These materials share in common the fact that heat resulting from the input pump source causes physical variation of the material. Namely, the material expands with the heat load due to a stress gradient formed which produces space-dependent birefringence in the material. Quite remarkably when the thermal gradient is strong enough the thermally-induced birefringence can resolve a polarized high-power input beam into two-component beams hans (); resolv2 (); resolv1 (), moreover upon recombination after traversing the bulk the two beams with two polarization directions are no longer in phase with one another such that the polarization state of the input beam cannot be recovered dist4 (); resolv1 ().
The thermally-induced depolarization phenomenon together with the resulting beam bifocusing have been investigated experimentally and theoretically by several authors hans (); menzel1 (); moshe1 (); moshe2 (); wang (); resolv2 (). In photonic crystals with a nonlinear Kerr-type dielectric susceptibility competing with diffraction, nonlocal thermal properties of the propagation medium have very recently been shown to induce twin-mode spatial solitons so-called dipole solitons yak (); cohen (); skupa (); bambi (). The robustness of these double-pulse structures provides strong indication that the competition between nonlinearity and thermal lensing might be a stabilizing factor for the double-polarized modes, potentially observable in several optical materials including solid-state and optical fibers with thermally-induced birefringence.
The present work aims at proposing a non-perturbative description of the generation of a double-polarized pulse beam in optical materials with thermally-induced linear birefringence, paying particular attention to nonlinear optical media with a quadratic spatial profile of temperature distribution. Our objective is to provide a consistent theoretical framework for understanding better the effects of thermal lensing on beam propagation in nonlinear optical materials, of which the YAG fiber laser which temperature profile was exprimentally established dist1 (); koechner (); eich (). We will first demonstrate, by solving the heat equation, that the topography of temperature distribution is quadratic in the case of a one-dimensional () anisotropic material uniformly loaded at a constant heat flux rate. Next, using the relationship between the temperature gradient and the induced inhomogeneous optical index, we shall derive an effective refractive index for the material by combining the thermally-induced optical index and a Kerr-type refractive index reflecting the intrinsic nonlinearity of optical properties of the anisotropic material. With the help of this effective refractive index we shall formulate the propagation equation for beams in the thermal nonlinear medium. Note that this equation is a Nonlinear Schrödinger equation (NLSE) with a repulsive external quadratic potential, that has already been derived in ref. cari () but solved perturbatively. Here instead, a non-perturbative treatment will be proposed based on a non-isospectral Inverse Scattering Transform (IST) method rada (); dika1 () with emphasis on IST’s initial solution being single-pulse shaped.
Thermal lensing can be described in simple words as a thermo-optical process associated with a weak absorption of an input beam that induces a nonzero temperature gradient across the material, leading to a spatial variation of its refractive index dist1 (); koechner (); cohen (). Recent experimental as well as theoretical developments on this process suggest that the underlying mechanism involves a local refractive index change which increases linearly cohen (); eich () with the temperature change i.e. , where refers to the thermo-optic coefficient cohen (); sowade (). Thus when the optical beam of a uniform thermal load gets slightly absorbed and heats the material, this produces heat that is conveyed by the electromagnetic wave. If denotes the uniform heat flux density and the heat conductivity coefficient, the heat diffusion in the material along a preferred direction (for anisotropic materials of current interest), driven by the uniform heat flux load, is determined by the heat equation:
Since eq. (1) is key to the current analysis we consider its most general solution given by:
where , and being two arbitrary real constants.
Formula (2) is
consistent with the quadratic law of temperature variation found for most laser
fibers in the presence of a uniform
head load. More specifically in YAG fiber lasers this law
is common dist1 (); koechner (); eich () and is consistent with the optical
bifocalility associated with a thermally-induced
birefringence, that promotes double-polarized laser beams from an input laser
For the sake of simplicity we require the temperature gradient to be zero and temperature to take a bare value at (i.e. the ambient temperature). The change of temperature in the material along the axis, hereafter assumed to be the axis of beam propagation, then reads:
With the last formula we derive the following expression for the local refractive index change:
Now if the intrinsic optical properties of the material are dominated by Kerr-type phenomena, the homogeneous part of the refractive index can be expressed as where is the beam intensity. With the help of (4) we can readily define an effective refractive index for the thermal nonlinear material viz:
Assuming that the wave motion is fast along the axis of anisotropy (i.e. ) but very slow along cohen (), the paraxial approximation on the wave equation for an electromagnetic field leads to:
As already stressed eq. (6) is actually not new, indeed the same equation
was obtained cari () for the same problem but solved
following the collective-coordinate method. To this last point, eq. (6) is an inhomogeneous NLSE and so can in principle be solved using the
collective-coordinate method. However this is a perturbative
method and consequently requires that the thermal lensing
is sufficiently weak, so that the Kerr nonlinearity remains the main
governing factor in the modulation and stability of signals in the
thermal nonlinear medium. So to say any
input beam sent in the medium must be modulated into a signal of
permanent single-pulse shape, with eventually an acceleration or slow down of
pulse due to the
thermal lensing. In fact this
consideration is very far from any acceptable consistency with the physics of
process under study, which the double-polarization of the incident beam is a
most salient aspect.
Being interested in a solution to eq. (6) which is more consistent with experiments we try for a non-perturbative approach. Remark to start that this equation can be rewritten in the following form:
corresponding to a NLSE with an external potential
which is quadratic in with a maximum at . The physics behind this
quadratic potential is contained both in its profile and the parameter
defined in (4), which determines the strength of the
thermal birefringence on the beam shape. One remarkable side of this
physics emerges from the assumption of the heat flux density and the heat
conductivity as being fixed, such that appears to be
increasing with an increase of the thermo-optic coefficient . Thus when
increased the curvature of the quadratic profile of heat distribution in the
material is more and more pronounced so that
the effect of thermal lensing on beam modulation is more and more strong. In
fact the external potential is expected to be much effective
on the beam position and according to the profile of this potential given
the centre of mass of the input beam should experience a trapping
the manifestly expulsive quadratic potential.
When is large enough such that the potential field erected by the thermal lensing process on the beam path is strongly localized, its contribution must be fully taken into consideration. In this last respect we follow the non-isospectral IST technique proper rada (); dika1 () to equations of this specific kind, considering an initial signal of a permanent single-pulse shape along at . For eq. (6) this technique leads to the following one-soliton solution:
If the ”sech” function in (9) reminds a pulse signal, the complicated dependence of its prefactor in and clearly suggests an actually complex pulse structure for the beam on the plane. To gain insight about what this dependence implies for the signal profile, in fig. 1 we plotted the beam amplitude as a function of at two different positions in the direction transverse to the beam propagation, and for four distinct values of .
More explicitely the four left graphs in fig. 1 represent versus at
for , , and , while the left graphs
for same values of but at . As one sees, the signal intensity is
a two-component pulse which intensities are strongly
dependent on the strength of the thermal lens potential.
It is quite noticeable on exploring the six figures
that when increases, the intensities of the two-component pulse
their width at half tails diminish. The last behaviour is
with the dependence of parameter defined in (15), which
indeed represents the average spatial extension of the pulse along the
Another relevant feature emerging from graphs of fig. 1 is the fact that when is decreased for a fixed value of , the two constituents pulses in the double-polarized beam preserve their shapes but their peak positions (i.e. centres) are gradually shifted. The last feature is more transparent in the left graphs corresponding to a relatively large value of . In fact, the last behaviour can be interpreted in terms of a ring profile for the signal intensity in the plane as reflected by the contour plots of fig. 2, where shadows of the double-pulse soliton in the plane are represented for different values of .
The figure clearly indicates an increase of the separation between pulses in the
implying that the radius and curvature of the ring signal in the
plane are fixed by the magnitude of .
The double-pulse structure obtained as well as its ring profile emerging in are reminiscent of dipole and ring-vortex solitons, predicted recently in some nonlocal nonlinear media cohen (); bambi (); fei (). However in these previous studies the physical origins of nonlocalities were generally not well specified, whereas in our context the quadratic form of the optical inhomogeneity has an experimental foundation dist1 (); koechner (); eich (). Also, a previous attempt cari () to model the same problem led to eq.(6). In this previous work the equation was treated perturbatively and so results could not reflect the remarkable aspects underlying the physics of thermal lensing in Kerr media. As end remark, the universality of the IST one-soliton solution (9) can be checked by applying any other exact spectral method to the generating equation, such as the Darboux method nakka () combined with a Lax-pair formalism with non-conserved spectral parameters porse (); khawa1 ().
Part of this work was done at the Abdus Salam International Centre for Theoretical Physics (ICTP) Trieste, Italy. I thank Markus Muller, Matteo Marsili, M. Kiselev and V. Kravtsov for their kind hospitalities.
- (1) W. Koechner, Solid State Laser Engineering (Springer-Berlin, edition, 2006).
- (2) J. D. Foster and L. M. Osterink, J. Appl. Phys. 41, 3656 (1970).
- (3) C. A. Klein, Opt. Eng. 29, 343 (1990).
- (4) H. J. Eichler, A. Hasse, R. Menzel and A. Siemoneit, J. Phys. D 26, 1884 (1993).
- (5) M. Ostermeyer, G. Klemz and R. Menzel, Proc. SPIE 4629, 67 (2002).
- (6) I. Moshe and S. Jackel, J. Opt. Soc. Am. B 22, 1228 (2005).
- (7) G. Machavariani, Y. Lumer, I. Moshe, A. Meir, S. Jackel and N. Davidson, Appl. Opt. 46, 3304 (2007).
- (8) Y. Wang, K. Inoue, H. Kan and S. Wada, J. Phys. D 42, 235108 (2009).
- (9) Y. Wang, H. Kan, T. Ogawa and S. Wada, J. Opt. 13, 015703 (2011).
- (10) E. A. Khazanov, O. V. Kulagin, S. Yoshida, D. B. Tanner and D. H. Reitze, IEEE Quantum Electron. 35, 1116 (1999).
- (11) F. Genereux, S. W. Leonard and H. M. Driel, Phys. Rev. B 63, 161101(R) (2001).
- (12) I. Shoji, Y. Sato, S. Kurimura, V. Lupei, T. Taira, A. Ikesue and K. Yoshida, Opt. Lett. 27, 234 (2002).
- (13) W. Koechner, Appl. Opt. 9, 2548(1970).
- (14) A. S. Koujelev and A. E. Dudelzak, Opt. Eng. 47, 085003(2008).
- (15) H. J. Eichler, A. Haase, R. Menzel and A. Siemoneit, J. Phys. D 26, 1884 (1993).
- (16) M. Frede, R. Wilhelm, M. Brendel, C. Fallnich, F. Seifert, B. Willke and K. Danzmann, Optics Express 12, 3581 (2004).
- (17) S. Vasudevan, G. C. K. Chen and B. S. Ahluwalia, Opt. Lett. 33, 2779 (2008).
- (18) M. Andika, G. C. K. Chen and S. Vasudevan, J. Opt. Soc. Am. B 27, 796 (2010).
- (19) S. Wielandy and A. L. Gaeta, Phys. Rev. Lett. 81, 3359(1998).
- (20) A. I. Yakimenko, Y. A. Zaliznyak and Y. Kivshar, Phys. Rev. E 71, 065603 (2005).
- (21) C. Rotschild, O. Cohen, O. Manela, M. Segev and T. Carmon, Phys. Rev. Lett. 95, 213904 (2005).
- (22) S. Skupin, M. Saffman and W. Królikowski, Phys. Rev. Lett. 98, 263902 (2007)
- (23) F. Ye, Y. V. Kartashov, Bambi Hu and L. Torner, Opt. Lett. 34, 2658 (2009).
- (24) A. Gharaati, P. Elahi and S. Cari, Act. Phys. Pol. A 112, 891 (2007)
- (25) R. Balakrishnan, Phys. Rev. B 32 1144 (1985).
- (26) A. M. Dikandé, J. Math. Phys. 49, 073520 (2008).
- (27) R. Sowade, I. Breunig, C. Tulea and K. Buse, Appl. Phys. B 99, 63 (2010).
- (28) F. Ye, B. A. Malomed, Y. He and Bambi Hu, Phys. Rev. A 81, 043816 (2010).
- (29) Z. Xu, L. Li, Z. Li, G. Zhou and K. Nakkeeran, Phys. Rev. E 68, 046605( 2003).
- (30) R. Rhada, V. R. Kumar and K. Porsezian, J. Phys. A 41, 315209 (2008).
- (31) U. Al Khawaya, J. Phys. A 39, 9679 (2006) & A 42, 265206 (2009).