Negative-frequency dispersive wave generation in quadratic media

Negative-frequency dispersive wave generation in quadratic media

Matteo Conforti, Niclas Westerberg, Fabio Baronio, Stefano Trillo, and Daniele Faccio CNISM, Dipartimento di Ingegneria dell’Informazione, Università di Brescia, Via Branze 38, 25123 Brescia, Italy
School of Engineering and Physical Sciences, SUPA, Heriot-Watt University, Edinburgh EH14 4AS, United Kingdom
Dipartimento di Ingegneria, Università di Ferrara, Via Saragat 1, 44122 Ferrara, Italy
July 16, 2019

We show that the extremely blue-shifted dispersive wave emitted in Kerr media owing to the coupling with the negative-frequency branch [Phys. Rev. Lett. 108, 253901 (2012)] can be observed in quadratic media via second-harmonic generation. Not only such phenomenon is thus independent on the specific nonlinear mechanism, but it is shown to occur regardless of the fact that the process is pumped by a pulse which exhibits soliton-like features or, viceversa, undergoes wave-breaking. A simple unified formula gives the frequencies of the emitted dispersive waves in both cases.

pacs:, 42.65.Re, 52.35.Tc

Introduction. Solitons emit resonant radiation (RR) owing to a universal mechanism of phase-matching with linear waves ruled by perturbing higher-order dispersive terms. Well-known examples range from fiber Wai (); AK95 () to Langmuir plasma Karpman97 () or water wave Karpman98 () solitons. During the last decade, optical fibers offered the unique opportunity to deepen the physics of RR Skryabin03 (); BSY04 (); SG07 (), with important applicative fall-out in supercontinuum generation DudleyRMP06 (), where RR is responsible for broadening the spectrum over the blue-shifted (normally dispersive) region Cristiani04 (). More recently, the field was significantly advanced by important results recognizing the role of RR in turbulence transport Rumpf09 (), the observation of RR in different settings encompassing tapered Stark11 () and noble-gas-filled photonic crystal fibers Joly11 (); Saleh11 (), slow-light waveguides Colman12 (), spatial diffraction in arrays Tran12 (), and second-harmonic generation (SHG) Bache10 (); Zhou12 (). Importantly, it was also shown that, in Kerr media, new frequencies can be generated owing to the coupling with the negative-frequency part of the spectrum, a process termed negative-frequency resonant radiation (NRR) rubino12 (); rubino12b (). With reference to this new phenomenon, the aim of this letter is twofold: (i) to assess the universal nature of NRR by showing that it can be predicted to occur also via pure quadratic nonlinearities, under experimentally viable conditions of SHG; (ii) to generalize the concept of radiation by showing that in fact one does not necessarily need a soliton-like excitation, since RR and NRR can be produced also in the opposite regime where the nonlinearity, instead of compensating the action of group-velocity dispersion (GVD), is such to strongly enforce it leading to wave-breaking (shock formation) Rothenberg89 (). This regime investigated recently in experiments performed in the spatial domain Hoefer06 (); Wan07 () features dispersive shock wave (DSW), as predicted in a seminal work by Gurevich and Pitaevskii GP74 ().

Resonant radiation. Let us first explain the general origin of the NRR from a different perspective as compared with the analysis of Ref. rubino12 (); rubino12b (). We consider an intense pump at carrier frequency , characterized by a complex envelope , which travels with characteristic group-velocity . When such a pulse travels in a nonlinear medium without experiencing significant dispersive effects, its total electric field, which is by definition a real quantity, can be written as , where is the retarded time and , where is the linear wave-number in the moving frame, is the nonlinear correction (phase-shift) due to the nonlinearity, and for the field to be real. Also linear waves (radiation) at frequency can be expressed in terms of positive- and negative-frequency content through the real field . Here is determined by the full dispersive relationship of the material in terms of the real index , neglecting losses for simplicity. Upon substitution , we cast the radiation in the form , where stands for the wave-number in a frame co-moving with the pump envelope . Focusing on the positive frequency content of the radiation, its amplitude starts to grow due to resonant transfer of energy from the pump at the phase-matching frequency such that , which gives the well-known condition for RR Wai (); AK95 (); SG07 (); Colman12 (); Bache10 (). Conversely, what it is usually not recognized is the fact that can be phase-matched also with the negative-frequency content of the pump at a different frequency such that the condition is fulfilled. By summarizing, both the phase-matched frequencies can be obtained by solving a unified set of two equations


We point out that, owing to symmetry, Eqs. (1) have also solutions , i.e. the image frequencies required by the Hermitian symmetry of the radiation field. Moreover, we arrive at Eq. (1) also when starting from the negative frequency content of the radiation .

Quadratic media. In the following, we address the open question as to whether NRR can be observed in a quadratic medium. Our aim is to show that, in such media, Eqs. (1) accurately predict both the RR and NRR frequencies even if the pump pulse is not strictly invariant (strictly non-dispersive), provided one is able to accurately estimate its velocity . In fact, a deviation from the ideal non-dispersive behavior of the pump is even beneficial since the growth of RR and NRR become significant when the pump undergoes a strong spectral broadening thereby seeding the phase-matched radiation frequencies. When nonlinearity and GVD act so as to mutually balance each other, this requires to operate with pulses which exhibit compression, i.e. higher-order solitons. However we go further by showing that also in the opposite regime, where a pulse experience strong temporal broadening, RR and NRR generation of comparable magnitude can be emitted. In order to demonstrate this, since the radiation detunings can be extremely large, we employ a description based on numerical integration of the Unidirectional Pulse Propagation Equation (UPPE2) implemented for anisotropic media Conforti11 (). The latter is suitable to describe ultra-broadband propagation, not relying on the separation of spectral envelopes around the carriers Bache10 (); reviewSHG ()), whose validity breaks down in the regime considered here. Furthermore this approach allows us to account for the full (all orders) dispersion through the Sellmeier equations which characterize any specific material nikogosyan ().

Figure 1: (Color online) Time domain electric field nota1 () evolution in a BBO crystal (ordinary polarization). (a) Linear scale. The black line is the trajectory of a pulse propagating with group velocity . Inset: snapshot at the point maximal compression compared with the input. (b) Same in log scale (dB). Here the dashed black and dashed white lines follow the peak of the RR ( nm) and the NRR ( nm), respectively.
Figure 2: (Color online) (a) Graphical solution of Eq.(1): the intersections between the blue curve and the red and green horizontal lines, standing for and , respectively, give the radiation frequencies ( nm) and ( nm). Here ( nm), (as derived from Fig. 1), and (negligible nonlinear shift). (b) Color level plot of the evolution of the electric field spectrum (ordinary polarization, log scale).

As a first example we consider a medium with anomalous GVD (), with SHG occurring in the regime of high negative mismatch , which results into an effective focusing Kerr nonlinearity, supporting solitary wave propagation reviewSHG (). These conditions can be realized e.g. in a BaBO (BBO) crystal at carrier wavelength nm. By exploiting type I () SHG in a crystal with orientation angles and (quadratic nonlinear coefficients pm/V, and pm/V nikogosyan ()), we obtain from Sellmeier equations a phase mismatch and a GVD psm. We consider the propagation of an ordinarily polarized hyperbolic secant pulse with fs duration. We show typical results obtained for a soliton number , where and (input peak intensity TW/cm in vacuum). Here is the effective Kerr nonlinear index due to cascading, being the vacuum impedance. Figure 1 shows the time domain evolution of the ordinarily polarized electric field nota1 (). Radiation starts to be emitted at the activation length mm, where the maximal pulse compression and spectral broadening are achieved [see Fig. 2(b)]. After this stage soliton fission occurs with the two constituent solitons separating asymptotically. The temporal evolution in log-scale reported in Fig. 1(b) clearly shows that the emitted radiation, which is slower, possesses two distinct branches traveling at different velocities, which turn out to correspond to the RR and the NRR dispersive waves. Indeed the central frequency of these two branches found from the spectral evolution in Fig. 2(b) to be ( nm) and ( nm), are accurately described by Eq. (1). The graphical solution of this equation displayed in Fig. 2(a), shows indeed that such values of and are obtained as the intersection of the dispersion curve with the wavenumber of the positive- [] and negative-frequency [] components of the pump pulse, respectively. We emphasize that, while the nonlinear phase shift turns out to negligible in this case, it is of paramount importance to accurately estimate the pump velocity around the activation length where the radiation is emitted, since the curve is dramatically affected by even small errors in the value of . Here we extract the correct value of from the time domain evolution, finding , which correctly describes the pulse velocity at its maximal compression as shown by the black line in Fig. 1(a).

Figure 3: (Color online) Color level plot of time domain evolution of the electric field nota1 () (ordinary polarization): (a) linear scale, early stage ( mm). Inset: snapshot at the point of shock formation compared with the input. The black line stands for the shock velocity ; (b) log (dB) scale, long range evolution ( mm). The dashed black and dashed white correspond to the RR and the NRR, respectively.
Figure 4: (Color online) As in Fig. 2 for the case of wave-breaking in GaSe described in Fig. 3. In (a) the graphical solutions yielding ( nm) and ( nm) are obtained with the full expression of (solid horizontal lines), while for comparison the dashed lines stand for the corresponding quantities calculated with .

In the second case we consider the opposite sign of dispersion, namely normal GVD, yet with the same sign of mismatch. Under these conditions, the cascading nonlinearity does not compensate for GVD-induced temporal broadening but rather enforces it. We have carefully chosen the operating conditions to work in a regime where the nonlinearity initially dominates over the dispersion (i.e. weakly dispersing regime), which is characteristic of the formation of DSW Rothenberg89 (); Hoefer06 (); Wan07 (); Conforti12 (). In this regime, we find viable conditions for the observation of NRR e.g. in Gallium Selenide (GaSe), which is characterized by a large nonlinear coefficient pm/V nikogosyan (). By operating at central wavelength nm with type I () SHG in a crystal oriented at angles and , we find from the Sellmeier formulas, a phase mismatch , and a GVD  psm. We show in Figs. 3-4 the outcome of our simulations obtained from a ordinarily polarized input pulse with fs duration (FWHMI), gaussian shape, and input peak intensity GWcm (in vacuum). In this case is several orders of magnitude larger than , and the dynamics is essentially dominated by the nonlinearity. The latter is responsible for the pulse temporal broadening and steepening shown in Fig. 3(a). In particular steepening is found to occur on the trailing edge until a gradient catastrophe leads to the formation of a shock wave (maximal steep front) at mm [see inset in Fig. 3(a)]. Here the dynamics is essentially different from Kerr media where two symmetric shocks are formed over the leading and trailing edges Rothenberg89 (), as also confirmed recently with reference to spatial dynamics Hoefer06 (); Wan07 (). The reason is that, in SHG the repeated up- () and down-conversion () gives rise, not only to the well known effective Kerr nonlinearity, but also to a leading-order steepening term Ilday04 (), owing to the group-velocity mismatch, which induces the shock to be asymmetric Conforti12 (). Whenever the group velocity at fundamental frequency is sufficiently larger than that at second-harmonic, this term dominates and leads to shock formation on the trailing edge Ilday04 (); Conforti12 (). Once formed, the shock front travels with a characteristic velocity which we estimate numerically to be , while it develops fast oscillations due to the GVD [these occur, in this case, on a small scale due to the absence of a pulse background Conforti12 () and hence are not visible in Fig. 3(a)]. The shock formation is accompanied by an abrupt spectral broadening [see Fig. 4(b)] and the consequent emission of radiation. The latter is emitted along two branches, as shown by the temporal evolution in log-scale displayed in Fig. 3(b). The different speeds of these branches arise from their different frequencies, which are found from the spectrum in Fig. 4(b) to be ( nm) and ( nm). Also in this case we denoted such frequencies as RR and NRR since they agree perfectly well with the values obtained by the graphical solution of Eq. (1), illustrated in Fig. 4(a). Again, it is crucial to have an accurate estimate of the pump velocity in the proximity of the radiation emission distance (). However, in this case the larger nonlinear coefficient results in a non-negligible . Indeed we predict the correct values of and by evaluating as arising from from local self-phase modulation, i.e. , being the location of the shock. We point out that, unlike the previous example, here the two RR and NRR frequencies are red- and blue-shifted, respectively, with respect to the second-harmonic frequency, whereas Fig. 4(b) shows also a weak component at that corresponds to the slowest wave in Fig. 3b, i.e. non-phase-matched third harmonic generation from the cascaded process .

Finally, we point out that nonlinearities can compete with quadratic ones due to the high intensities (especially in BBO case) involved. However, both BBO and GaSe exhibit a focusing Kerr nonlinear index, which simply results in a lowered intensity threshold for radiation emission without any significant change to the dynamics illustrated above. This is confirmed by additional simulations (not reported) where we account also for the intrinsic Kerr nonlinearity.

Conclusions. In summary, we have demonstrated the generality of the NRR phenomenon, by predicting that it can be observed in quadratic media, under different scenarios that involve pumping either with soliton-like pulses or, in the opposite regime, when pulses undergo wave-breaking. This represents a substantial step forward towards the understanding and management of ultrafast cascading nonlinearities for producing broadband emission (supercontinuum) in standard crystals.

Funding from MIUR (grant PRIN 2009P3K72Z) and EPSRC (grant EP/.J00443X/1) is gratefully acknowledged.


  • (1) P. K. A. Wai, C. R. Menyuk, Y. C. Lee, and H. H. Chen, Opt. Lett. 11, 464 (1986); ibidem, 12, 628 (1987); P. K. A. Wai, H. H. Chen, and Y. C. Lee, Phys. Rev. A 41, 426 (1990); H. H. Kuehl and C. Y. Zhang, Phys. Fluids B 2, 889 (1990).
  • (2) N. Akhmediev and M. Karlsson, Phys. Rev. A 51, 2602 (1995).
  • (3) V. I. Karpman and H. Schamel, Phys. Plasmas 4, 120 (1997).
  • (4) V. I. Karpman, Phys. Rev. E 58, 5070 (1998).
  • (5) D. V. Skryabin, F. Luan, J. C. Knight, and P. St. J. Russell, Science 301, 1705 (2003).
  • (6) F. Biancalana, D. V. Skryabin, and A. V. Yulin, Phys. Rev. E 70, 016615 (2004).
  • (7) D. V. Skryabin and A. V. Gorbach, Rev. Mod. Phys. 82, 1287 (2010).
  • (8) J. M. Dudley, G. Genty, and S. Coen, Rev. Mod. Phys.  78, 1135 (2006).
  • (9) I. Cristiani, R. Tediosi, L. Tartara, and V. Degiorgio, Opt. Express 12, 124 (2004).
  • (10) B. Rumpf, A. C. Newell, and V. E. Zakharov, Phys. Rev. Lett. 103, 074502 (2009).
  • (11) S. P. Stark, A. Podlipensky, and P. St. J. Russell, Phys. Rev. Lett. 106, 083903 (2011).
  • (12) N.Y. Joly, J. Nold, W. Chang, P. Hölzer, A. Nazarkin, G. K. L. Wong, F. Biancalana, and P. St. J. Russell, Phys. Rev. Lett. 106, 203901 (2011).
  • (13) M. F. Saleh, W. Chang, Hölzer, A. Nazarkin, J.C. Travers, N. Y. Joly, P. St. J. Russell, and F. Biancalana, Phys. Rev. Lett. 107, 203902 (2011).
  • (14) P. Colman, S. Combrié, G. Lehoucq, A. de Rossi, and S. Trillo, Phys. Rev. Lett. 109, 093901 (2012).
  • (15) T. X. Tran and F. Biancalana, arXiV 1210.5204.
  • (16) M. Bache, O. Bang, B. B. Zhou, J. Moses, J. and F. W. Wise, Phys. Rev. A 82, 063806 (2010). C. R. Phillips, Carsten Langrock, J. S. Pelc, M. M. Fejer, I. Hartl, and Martin E. Fermann, Opt. Express 19, 18754 (2011).
  • (17) B. B. Zhou, A. Chon, F. W. Wise, and M. Bache, Phys. Rev. Lett. 109, 043902 (2012).
  • (18) E. Rubino, J. McLenaghan, S. C. Kehr, F. Belgiorno, D. Townsend, S. Rohr, C. E. Kuklevicz, U. Leonhardt, F. König, and D. Faccio, Phys. Rev. Lett. 108, 253901 (2012).
  • (19) E. Rubino, A. Lotti, F. Belgiorno, S.L. Cacciatori, A. Couairon, U. Leonhardt, D. Faccio, Scient. Rep., in press (2012), arXiv:1211.0256v1
  • (20) J. E. Rothenberg and D. Grischkowsky, Phys. Rev. Lett. 62, 531 (1989).
  • (21) M. A. Hoefer, M. J. Ablowitz, I. Coddington, E. A. Cornell, P. Engels, and V. Schweikhard, Phys. Rev. A74, 023623 (2006).
  • (22) W. Wan, S. Jia, and J. W. Fleischer, Nature Phys. 3, 46 (2007); N. Ghofraniha, C. Conti, G. Ruocco, S. Trillo, Phys. Rev. Lett. 99, 043903 (2007).
  • (23) A.V. Gurevich and L.P. Pitaevskii, Sov. Phys. JETP 38, 291 (1974);
  • (24) M. Conforti, F. Baronio, and C. De Angelis, Phys. Rev. A 81, 053841 (2010); J. Opt. Soc. Am. B 28, 1231 (2011).
  • (25) A. V. Buryak, P. Di Trapani, D. V. Skryabin, and S. Trillo, Phys. Rep. 370, 63 (2002).
  • (26) D. N. Nikogosyan, Nonlinear Optical Crystals: A Complete Survey, Springer (Berlin, 2005)
  • (27) F. Ö. Ilday, K. Beckwitt, Y.-F. Chen, H. Lim, and F. W. Wise, J. Opt. Soc. Am. B 21, 376 (2004).
  • (28) M. Conforti, F. Baronio, and S. Trillo, Opt. Lett. 37, 1082 (2012).
  • (29) We display the inverse Fourier transform of the positive frequency components of the spectrum. This allows us to directly visualize the envelope, eliminating the fast oscillations of the carrier while retaining the full broadband features of the spectrum.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description