arxiv:physics/ v1 28 Aug 2002

Transcription

1 Second Harmonic X-wave Generation Claudio Conti Λ and Stefano Trillo National Institute for the Physics of Matter, INFM- Roma Tre, Via della Vasca Navale 84, Rome - Italy (Dated: January 28, 2003) arxiv:physics/ v1 28 Aug 2002 Abstract In the process of optical second-harmonic generation an X-wave, travelling at the group velocity of the fundamental frequency, spontaneously forms. Its parameters, as well as the angular spectrum, are directly correlated to the group velocity mismatch and to the wavevector mismatch between the two harmonics. A novel mechanism for spatio-temporal localized waves generation is thus proposed. Λ Electronic address: c.conti@ele.uniroma3.it 1

2 According to some authors (see (Recami, 1998)) the first prediction of localized solutions to the linear Maxwell equations dates back to the work of Stratton in (Stratton, 1941) Since then a huge amount of research work has been produced in this field, starting from the article by Durning et al. in (Durnin et al., 1987) Experimental results have confirmed the existence of these waves, (Mugnai et al., 2000; Saari and Reivelt, 1997; Sonajalg et al., 1997) and important theoretical advances have been obtained: from generalization of the propagation invariant solutions (Gutierrez et al., 2000; Piestun and Shamir, 1998), to the description of practical realizations with finite aperture (Shaarawi, 1997; Stepanishen, 1998; Stepanishen and Sun, 1997; Ziolkowski et al., 1993), to the analysis of dispersive media (Porras, 2001; Sonajalg and Saari, 1996; Zamboni-Rached et al., 2002) and waveguides (Zamboni-Rached et al., 2001a). A lot of attention has been dedicated to the superluminal propagation(recami, 1998) and significant possible applications, as femtosecond microscopy and communication systems(lu and He, 1999) have been evisaged. (Saari, 2001) The most simple localized wave, solution to the wave equation in a dispersive media with refractive index n(ω), can be expressed as a superposition of Bessel beams: E(r;z;t)= Z 1 0 S(Ω)J 0 (»r)exp[i(k z z Ωt)]dΩ, (1) where» = fi(ω)sin[ (Ω)] and k z = fi(ω)cos[ (Ω)] are, respectively, the transversal and the longitudinal wavenumbers, with fi(ω) = (Ω=c)n(Ω) and c the vacuum lightvelocity.(sonajalg and Saari, 1996) Eq.(1) describes a non-diffracting-non-dispersive beam if the X-wave dispersion relation k z (Ω) = ff + fi X Ω, (2) is satisfied. Indeed, in such away, the group velocity dispersion k 00 z is zero (together with higher order derivatives of k z ) and the intensity only depends on the combination t fi X z, i.e. 2

3 fi X = k 0 z is the inverse velocity of the wave. Given the spectrum S(Ω), ff and fi X are the two parameters of the localized wave. In the case of vacuum, when ff = 0 the spatio-temporal (or angular) spectrum has an X-shape (» 2 = Ω2 c 2 (1 c 2 fi 2 X)) and, due to the invariant propagation, the far-field coincides with the near-field; thus the wave retains its X-shape in the physical space (x;y;z;t) and it is superluminal: cjfi X j < 1 (in the Fourier-transformable stationary regime no general theorem about causality and Maxwell equations is violated (Levi-Civita, 1931), as discussed extensively in the literature (Lunardi, 2001; de Oliveira et al., 2001; Recami et al., 2000; Rodrigues et al., 2001). When ff 6= 0 the X in the spectrum is curved and subluminal propagation is allowed. (Salo and Salomaa, 2001) In general these localized waves have infinite energy, and only approximations can be realized in practice (as for non-diffracting Bessel beams). Finite energy (not propagation invariant) solutions have been obtained.(de Oliveira and Rodrigues, 2001; Zamboni-Rached et al., 2001b) In the case of dispersive media the situation is much more complicated (the general angular spectra is given by» 2 = fi(ω) 2 (fi X Ω+ff) 2, encompassing evanescent waves). In (Reivelt and Saari, 2000) methods for generating X-waves in the optical domain have been suggested: linear dispersive optical elements may help to satisfy eq.(2). Experiments were realized by computer holograms. (Sonajalg et al., 1997) Very recently experimental results (Valiulis et al., 2001) and theoretical investigations (Conti et al., 2002a,b; Trillo et al., 2002) have pointed out the relevant rule of X-waves in the field of nonlinear optics. As shown in (Trillo et al., 2002), in Kerr and in quadratic media a plane wave breaks up, due to spatio-temporal modulational instability (the conical emission), with an angular-spectrum gain strongly resembling that of a localized wave. In (Conti et al., 2002b), X-shaped numerical solutions of nonlinear wave equations have been reported 3

4 for the first time. With reference to second harmonic generation (SH-G), these solutions are given by two X-waves, one at the fundamental frequency (FF) and the other at the SH, travelling locked together (they cannot be considered solitary waves, or solitons, in the usual sense, because they continue to exist in the absence of nonlinearity). Since X-waves exist in the linear regime (i.e. in the absence of any coupling between FF and SH), but also as solutions to nonlinear wave equations, it is natural to study the mechanism of their formation in the intermediate stages. We will show in this Letter that X-waves are spontaneously produced by an FF pulsed pump beam generating SH during propagation. We consider the process of SHG in the undepleted pump approximation. For slowly-varying envelopes and paraxial beams, the SH amplitude a obeys the following r 2 2k a + iffla + 2 χa2 e i fiz =0. (3) 2 The parameters in (3) are: the wavevector mismatch fi = k 2 2k 1 with k j = fi(! j ) (j =1; 2), and the inverse SH group velocity fi S = fi 0 (! 2 ), being! 2 the SH and! 1 the FF frequency. We neglect group velocity dispersion (GVD), as being typically small in the experiments, while the group velocity mismatch (GVM), given by fi F fi S (fi F = fi 0 (! 1 )), is dominant. In (3) the source term is given by the squared amplitude A 2 (measured in W=m 2 as jaj 2 ) of the FF beam, weighted by the relevant nonlinear coefficient χ. Note that a practically unavoidable loss term ffl >0 has been included in the model for technical convenience. The results will concern the case ffl! 0 +. For the sake of simplicity, we will neglect the diffraction at FF, i.e. we will consider sufficiently large input FF beams. Consequently A will be an envelope travelling undistorted with group velocity 1=fi F. The initial condition for the SH is assumed to be a(z = 0) = 0, i.e. no SH energy is seeded at the input face of the nonlinear crystal. To integrate equation (3) we make use of the Fourier- Bessel representation: a(r;t;z)= 1 2ß R 1 R A(!;»; z)j 0(»r)e i!t»d»d!. The transformed 4

5 β S z β F z t β F z β S z t FIG. 1 Orientation (and support) of generated SH X-wave for fi F >fi S (top) and fi F <fi S (bottom). source term in (3) is expressed as χψ(!;»)exp(i!fi F z), assuming circular symmetry for the input FF beam. The solution in the space (!;»; z) is found to be: A = χψ ei(!fi F fi)z e i(fis!»2 =2k 2 +iffl)z. (4) (fi F fi S )! +» 2 =2k 2 iffl We will limit the analysis to the case fi F 6= fi S ; when fi F = fi S higher order dispersion terms must be taken into consideration, as it will be discussed in fore-coming publications. The inverse Fourier transform of (4) can be obtained by applying residue theorems (the details will be given elsewhere). For the field a we found: a = 0 when t < min(fi F z;fi S z), or t > max(fi F z;fi S z), and a = i fi t fi Sz fi F fi S jfi F fi S j Z 1 0 X(»)e i» 2 t fi F z 2 k 2 fi F fi S J 0 (»r)»d» (5) elsewhere (see fig. 1), with X = ψ[(» 2 =2k 2 + fi)=(fi F fi S );»]. Hence in the region covered by the SH (travelling at speed 1=fi S ) emitted at the center of the FF beam, a is given by a non-diffracting beam which propagates at the group velocity of the FF. Indeed the intensity jaj 2 only depends on z via the combination t fi F z. Only one-half of the localized wave is generated due to causality: the peak moves at the FF group velocity, while its costituents (the SH beam) moves at their characteristic velocity 1=fi S. As a result, the extension of the wave is equal to jfi F fi S jz. This is similar to what has been theoretically predicted for 5

6 practical finite-aperture realizations of vacuum localized waves. (Ziolkowski et al., 1993) As in that case, it can been shown that the a field is the sum of two X-waves which cancel each other in the non-causal regions, a form of the extintion theorem. (Born and Wolf, 1980) In the source-sink picture developed by Saari (Saari, 2001), only the source cone is present (as physically reasonable). Eq. (5) strongly resembles the vacuum nondiffracting beams developed in (Salo and Salomaa, 2001). However, since we are dealing with a dispersive media, we need to show that (5) belongs to the class given by (1). By applying the paraxial approximation, we have k z = [fi(ω) 2» 2 ] 1=2 ο = fi(ω)»2 =2fi(Ω); next by expanding fi(ω) around! 2, up to the first order dispersion, and neglecting space-time coupling terms (i.e. terms like!» 2 ), we have k z = k 2 + fi S!» 2 =2k 2, with! =Ω! 2. Since eq.(2), has to be satified we obtain the following identity:! =(k 2 fi X! 2 ff)=(fi X fi S ). The latter, used in (1), after some manipulations yelds: E = e 2 e i(k 2 fi X! 2 ff) t fi S z fi X fi S Z 1 0 X(»)J 0 (»r)e i» 2 t fi X z 2k 2 fi X fi S»d» (6) which is exactly (5), with the needed addition of the optical carrier wave e 2 = e i(k 2z! 2 t) and an appropriate redefinition of the spectrum X(»). By comparison we found fi X = fi F and ff = k 2 fi X! 2 fi =2k 1 fi F! 2. The parameters of the X-wave are directly related to those of the SHG process, as could be expected. Therefore we have demonstrated that for SHG eqs. (2) is automatically satisfied and given by k z = fi F!+2k 1. Correspondingly the angular spectrum is given by the equation (» ο = k 2, is the radial angle)! = fi k 2 fi F fi S 2(fi F fi S ) 2, (7) which is an horizontal parabola in the angle-wavelenght plane, with position and convessity 6

7 determined by the the GVM and the wavector mismatch. Its experimental observation is a clear signature of the presence of angular dispersion and of the formation of SH-Xwaves. R 1 After equation (5) we may determine the energy of the SH: E(z) =2ß 1R 1 0 jaj2 rdrdt = χ 2 F FF ( fi;fi F fi S ) jfi F fi S j z, where F FF = R 1 0 jxj2»d». E linearly during propagation, its finitess is a consequence of the X-wave compact support, and hence of causality. To make a practical example we assume the FF amplitude given by a gaussian pulsed beam with waist w F and duration t F : A = A 0 exp[ (t fi F z) 2 =2t 2 F r 2 =2w 2 F ]. In this case F FF can be evaluated explicitly yelding E = χ2 k 2 E 2 F z 2 p e k 2 w F 2 jfi F fi S j 2 p 2t F e ο [1 + erf(ffο)], (8) 2ßt F with E F the FF energy, ο = k 2 w 2 F ( fi fi X ), fi X = k 2 w 2 F (fi F fi S ) 2 =2t 2 F, and ff = t F = p 2jfi F fi S jk 2 w 2 F. Eq. (8) shows that for large j fij the energy conversion is negligibile, while the maximum efficiency is obtained for a particular value of fi (in general different from zero) which can evaluated numerically. Choosing realistic parameters (Valiulis et al., 2001) (w F =10 5, t F =10 13, jfi F fi S j = , k 2 ο = 8ß 106, SI units) ff ο = 1 and the maximum conversion efficiency is attained for fi ο = fi X ο = jfif fi S j=t F, i.e. the inverse spatial separation between FF and SH during t F. At this point it is necessary to discuss the adopted nomenclature. Rigourously speaking, eq. (1) defines the so-called Focus Waves Mode (FWM), a more general class than X-waves, according to the literature.(reivelt and Saari, 2000) However we are dealing with almost monochromatic optical waves, and hence FWM are the so called X-Bessel (Saari and Reivelt, 1997; Sonajalg and Saari, 1996). On the other hand eq. (5) resembles the Focused X-waves (Saari, 2001). Furthermore the causality implies that only one-half of the double cone is generated, and thus they can be identified with the so-called Y-, or Hankel-, waves. (Salo et al., 2000) For the sake of simplicity we used the term X-wave" from the beginning. 7

8 In conclusion we have theoretically predicted that in the process of second harmonic generation an X-wave, travelling at the FF velocity, and with the ff parameter directly connected to the wavevector mismatch, is spontaneously generated. For a given input FF wave, the X-wave can be modified by acting on the parameters ruling the process of SHG, e.g. by varying the temperature and the orientation of the crystal, while the resulting angular spectra (the far field) is given by horizontal parabolas (Eq.(7)). This mechanism appears as a valid alternative to the generation of optical localized-waves by linear dispersive devices, or computer holograms. Our results can also be interpreted in a different way. The fact that the X-waves spontaneously form unreveals their fundamental role. In any ultrafast parametric process (e.g. fourwave mixing or optical parametric amplification) the phase-matching condition produces the correlation between spatial and temporal degrees of freedom at the basis of X-waves. We believe that, by taking into account such a framework, new generations of parametric sources could be developed. Acknowledgments We thanks P. Di Trapani and E. Recami for fruitful discussions. References Born, M., and E. Wolf, 1980, Principles of Optics (Pergamon, New York), 6 edition. Conti, C., P. D. Trapani, S. Trillo, G. Valiulis, S. Minardi, and O. Jedrkiewicz, 2002a, in Quantum Electronics and Laser Science Conference (Optical Society of America, Washinton,D.C.), p. QTuJ6. 8

9 Conti, C., S. Trillo, P. D. Trapani, O. Jedrkiewicz, G. Valiulis, and J. Trull, 2002b, arxiv:physics/ Phys. Rev. Lett. submitted. Durnin, J., J. J. Miceli, and J. H. Eberly, 1987, Phys. Rev. Lett. 58(15), Gutierrez, J. C., M. D. Iturbe-Castillo, and S. Chavez-Cerda, 2000, Opt. Lett. 25, Levi-Civita, T., 1931, Caratteristiche dei sistemi differenziali e propagazione ondosa (Zanichelli, Bologna, Italy). Lu, J. Y., and S. P. He, 1999, Opt. Commun. 161, 187. Lunardi, J. T., 2001, Phys. Lett. A 291, 66. Mugnai, D., A. Ranfagni, and R. Ruggeri, 2000, Phys. Rev. Lett. 84, de Oliveira, E. C., and W. A. Rodrigues, 2001, Phys. Lett. A 291, 367. de Oliveira, E. C., W. A. Rodrigues, D. S. Thober, and A. L. Xavier, 2001, Phys. Lett. A 284, 296. Piestun, R., and J. Shamir, 1998, J. Opt. Soc. Am. A 15, Porras, M. A., 2001, Opt. Lett. 26, Recami, E., 1998, Phys. A 252, 586. Recami, E., F. Fontana, and R. Garavaglia, 2000, Int. Journ. of Mod. Phys. A 15, Reivelt, K., and P. Saari, 2000, J. Opt. Soc. Am. A 17, Rodrigues, W. A., D. S. Thober, and A. L. Xavier, 2001, Phys. Lett. a 284, 217. Saari, P., 2001, arxiv:physics/ Saari, P., and K. Reivelt, 1997, Phys. Rev. Lett. 79, Salo, J., J. Fagerholm, A. T. Friberg, and M. M. Salomaa, 2000, Phys. Rev. E 62(2), Salo, J., and M. M. Salomaa, 2001, ARLO 2, 31. Shaarawi, A. M., 1997, J. Opt. Soc. Am. A 14. 9

10 Sonajalg, H., M. Ratsep, and P. Saari, 1997, Opt. Lett. 22(5), 310. Sonajalg, H., and P. Saari, 1996, Opt. Lett. 21, Stepanishen, P. R., 1998, J. Acous. Soc. Am. 104, Stepanishen, P. R., and J. Sun, 1997, J. Acous. Soc. Am. 102, Stratton, J. A., 1941, Electromagnetic Theory (McGraw-Hill, New York, USA). Trillo, S., C. Conti, P. D. Trapani, O. Jedrkiewicz, J. Trull, G. Valiulis, and G. Bellanca, 2002, Opt. Lett. 27, Valiulis, G., J. Kilius, O. Jedrkiewicz, A. Bramati, S. Minardi, C. Conti, S. Trillo, A. Piskarskas, and P. D. Trapani, 2001, in Quantum Electronics and Laser Science Conference (Optical Society of America, Washington, D.C.), volume 57 of Trends in Optics and Photonics, postdeadline Papers QPD10-1. Zamboni-Rached, M., K. Z. Nobrega, H. E. Hernandez-Figueroa, and E. Recami, 2002, unpublished. Zamboni-Rached, M., E. Recami, and F. Fontana, 2001a, Phys. Rev. E 64, Zamboni-Rached, M., E. Recami, and H. E. Harnandez-Figueroa, 2001b, arxiv:physics/ Ziolkowski, R. W., I. M. Besieris, and A. M. Shaarawi, 1993, J. Opt. Soc. Am. A 10,