# Optical Kerr spatiotemporal dark extreme waves

###### Abstract

We study the existence and propagation of multidimensional dark non-diffractive and non-dispersive spatiotemporal optical wave-packets in nonlinear Kerr media. We report analytically and confirm numerically the properties of spatiotemporal dark lines, X solitary waves and lump solutions of the nonlinear Schrödinger equation (NLSE). Dark lines, X waves and lumps represent holes of light on a continuous wave background. These solitary waves are derived by exploiting the connection between the NLSE and a well-known equation of hydrodynamics, namely the Kadomtsev-Petviashvili (KP) equation. This finding opens a novel path for the excitation and control of spatiotemporal optical solitary and rogue waves, of hydrodynamic nature.

Optical Kerr spatiotemporal dark extreme waves

Keywords: Nonlinear optics, self-action effects, kerr effect, solitons

## 1 Introduction

Techniques to shape and control the propagation of electromagnetic radiation are of paramount importance in many fields of basic science and applied research, such as atomic physics, spectroscopy, communications, material processing, and medicine [1, 2]. Among these, of particular interest are methods that produce localized and distortion-free wave-packets, i.e. free from spatial and temporal spreading due to diffraction or material group velocity dispersion (GVD), respectively [3]. There are two main strategies to achieve propagation invariant electromagnetic wave packets [4]. The first methodology is based on the spatiotemporal synthesis of a special input wave, so that diffractive and dispersive effects compensate for each other upon linear propagation in the material. Building blocks of these linear light bullets are Bessel beams and their linear combinations, along with Airy pulses [5]. These waveforms enable the generation spatiotemporal invariant packets such as the Airy-Bessel beams, and the so-called X-waves, obtained by a linear superposition of Bessel beams with different temporal frequencies. The second approach involves the generation of solitary waves, that exploit the nonlinear (quadratic or cubic) response of the material for compensating diffractive and dispersive wave spreading [6, 7]. Although successfully exploited in propagation models, that describe for example temporal solitons in optical fibers and spatial solitons in slab waveguides, in more than one dimension spatiotemporal solitons have so far largely eluded experimental observation, owing to their lack of stability associated with the presence of modulation instability (MI), collapse and filamentation.

Here, we overview our recent contributions to the field of non-diffractive and non-dispersive wave-packets in Kerr media [8, 9, 10], by deriving analytically and confirming numerically the existence and propagation of novel multidimensional dark non-diffractive and non-dispersive spatiotemporal solitons propagating in i) self-focusing and normal dispersion Kerr media, and in ii) self-defocusing and anomalous dispersion Kerr media. The analytical dark solitary solutions are derived by exploiting the connection between the NLSE and the Kadomtsev-Petviashvili (KP) equation [11], a well-known equation of hydrodynamics. Our results extend and confirm the connection between nonlinear wave propagation in optics and hydrodynamics, that was first established in the ’s [12, 13, 14].

## 2 Optical NLSE–hydrodynamic KP mapping

In the presence of group-velocity dispersion and one-dimensional diffraction, the dimensionless time-dependent paraxial wave equation in cubic Kerr media reads as [15]:

(1) |

namely the (2+1)D, or more precisely (1+1+1)D, NLSE, where stands for the complex wave envelope, and represent the retarded time, in the frame traveling at the natural group-velocity, and the spatial transverse coordinate, respectively, and is the longitudinal propagation coordinate. Each subscripted variable in Eq. (1) stands for partial differentiation. are normalized real constants that describe the effect of dispersion, diffraction and Kerr nonlinearity, respectively.

We refer (1) to as elliptic NLSE if , and hyperbolic NLSE if . In the case of weak nonlinearity, weak diffraction and slow modulation, the dynamics of optical NLSE dark envelopes may be related to the hydrodynamic KP variable as follows [8, 9]:

(2) |

where stands for a background continuous wave amplitude, represents a small amplitude variation, say with and the order one background ; ; satisfies the KP equation,

(3) |

where , with , (see [8] for further details).

Of interest in the optical context, the elliptic anomalous dispersion and self-defocusing regime [8] (, ) leads to the KP-I regime, while the hyperbolic normal dispersion and self-focusing regime [9] (, ), leads to the KP-II regime.

Without loss of generality, we may set the following constraints to the coefficients of Eq. (1), ; moreover, we fix (thus ). Note that, with the previous relations among its coefficients, in the case (, ), Eq. (3) reduces to the standard KP-I form: . Whereas in the case (, ), Eq. (3) reduces to the standard KP-II form . Moreover, the imposed constraints to the coefficients of Eq. (1) also fix the scaling between the dimensionless variables in Eq. (1) and the corresponding real-world quantities . The longitudinal scaling factor turns out to be , where is the usual nonlinear length associated with the intensity of the background and , being the Kerr nonlinear index and the vacuum wavenumber. The “transverse” scales read as and , where and are the group-velocity dispersion and the linear refractive index, respectively.

## 3 Normal dispersion and self-focusing regime: nonlinear lines and X-waves

At first, we consider the case of normal dispersion and self-focusing nonlinearity [9, 10]. We proceed to consider the existence and propagation of NLSE dark line solitary waves, which are predicted by the existence of KP-II bright line solitons [16, 17]. When considering the small amplitude regime, a formula for an exact line bright soliton of Eq. (3) can be expressed as follows [16, 17]: where rules the amplitude and width of the soliton, is the angle measured from the axis in the counterclockwise, is the velocity in -direction. Notice that is of order . Moreover we obtain The analytical spatiotemporal envelope intensity profile of a NLSE dark line solitary wave is given by the mapping (2) exploiting the KP bright soliton expression. The intensity dip of the dark line solitary wave is , the velocity in the -direction. We numerically verified the accuracy of the analytically predicted dark line solitary waves of the NLSE. To this end, we made use of a standard split-step Fourier technique, which is commonly adopted in the numerical solution of the NLSE (1). Figure 1 shows the numerical spatiotemporal envelope intensity profile of a NLSE dark line solitary wave, which corresponds to the predicted analytical dynamics.

As can be seen from the images, the numerical solutions of the NLSE show an excellent agreement with the analytical approximate NLSE solitary solutions.

In the long wave context, the KP-II equation admits complex soliton solutions, mostly discovered and demonstrated in the last decade, which may describe non-trivial web patterns generated under resonances of line-solitons [16, 17]. Here, we consider the resonances of four line solitons, which give birth to the so-called O-type bright X-shaped two-soliton solution of the KP-II (the name O-type is due to the fact that this solution was originally found by using the Hirota bilinear method). When considering the small amplitude regime, the formula of the O-type solution of Eq. (3) can be expressed as follows, where the function is given by with are small real positive parameters which are related to the amplitude, width and the angle of the O-type X-soliton solutions. The corresponding (2+1)D NLSE dark X solitary wave is directly given through the mapping Eq. (2), by exploiting the soliton expression for .

We numerically verified the accuracy of the analytically predicted O-type dark X solitary wave of the NLSE. Fig. 2 shows the profile of the numerical solution of the hyperbolic NLSE at , and at . In this particular example we have chosen , . Specifically, Fig. 2 illustrates a solitary solution which describes the X-interaction of four dark line solitons. The maximum value of the dip in the interaction region is . Asymptotically, the solution reduces to two line dark waves for and two for , with intensity dips and characteristic angles , measured from the axis. Numerical simulations and analytical predictions are in excellent agreement. We estimate the error between the asymptotic formula and the X solitary wave in the numerics to be lower than

## 4 Anomalous dispersion and self-defocusing regime: dark lumps

Next, we consider the case of anomalous dispersion and self-defocusing nonlinearity [8]. We proceed to verify the existence of (2+1)D NLSE dark-lump solitary waves, as predicted by the solutions of KP-I through Eq.(2) (see [18] for details). When considering the small amplitude regime (), a form of KP lump-soliton solution of Eq. (3) can be expressed as . The parameter rules the amplitude/width and velocity properties of the KP lump soliton. The lump peak amplitude in the plane is ; the velocity in the -direction is . Moreover, The analytical spatiotemporal envelope intensity profile of a NLSE dark solitary wave is given by the mapping (2), which exploits the KP bright lump expression. Then, we numerically verified the accuracy of the analytically predicted dark lumps solitary waves of the NLSE. Figure 3 shows the numerical spatio-temporal envelope intensity profile of a NLSE dark lump solitary wave in the - plane (), at the input and after the propagation distance , for . In the numerics, the initial dark NLSE profile, of KP-I lump origin, propagates stably in the -direction, with virtually negligible emission of dispersive waves, with the predicted velocity , and intensity dip of . Thus, the predicted theoretical dark lump solitary waves of Eq. (2) are well confirmed by numerical simulations.

We remark that the KP-I equation admits other types of lump solutions which have several peaks with the same amplitude in the asymptotic stages . We call such lump solution multi-pole lump. Here we show that D NLSE can also support such lump solutions. We consider multi-pole lump solution with two peaks, which is expressed as: , where , and , . define the dislocation; are arbitrary complex parameters. The analytical spatiotemporal envelope intensity profile of a NLSE dark solitary wave is again given by the mapping (2).

Figure 4 shows the initial spatio-temporal envelope intensity profile of a two-peaked NLSE dark lump in the plane, along with the numerically computed profiles after propagation distances , for (). In particular, Fig. 4 depicts the scattering interaction of the two-peaked waves: two dark lumps approach each other along the -axis, interact, and subsequently recede along the -axis. These solutions exhibit anomalous (nonzero deflection angles) scattering due to multi-pole structure in the wave function of the inverse scattering problem. We remark that the numerical result of NLSE dynamics is in excellent agreement with analytical dark solitary solution Eq. (2) with KP-I multi-pole lump solution.

## 5 Instabilities and Experimental Feasibility

Let us discuss the important issue of the stability of the predicted dark line, X solitary waves and lumps. Two instability factors may affect the propagation of these waves. The first one is the modulation instability (MI) of the continuous wave background. In the case of normal dispersion and self-focusing, , , MI is of the conical type [19]. Generally speaking, MI can be advantageous to form X waves from arbitrary initial conditions both in the absence or in the presence of the background. However, for sufficiently long propagation distances the MI of the CW background may compete and ultimately destroy the propagation of dark solitary waves and their interactions. In the case of anomalous dispersion and in the self-defocusing regime, , , MI is absent, thus lumps are not affected by MI.

The second mechanism is related to the transverse instability of the line solitons that compose the asymptotic state of the X wave. We point out that such instability is known to occur for the NLSE, despite the fact that line solitons are transversally stable in the framework of the KP-II (unlike those of the KP-I) [11]. However, in our simulations of the NLSE, these transverse instabilities never appear, since they are extremely long-range, especially for shallow solitons. In fact, we found that the primary mechanism that affects the stability of dark line and X solitary waves is the MI of the CW background.

Let us briefly discuss a possible experimental setting in nonlinear optics for the observation of cubic spatiotemporal solitary wave dynamics of hydrodynamic origin. As to (2+1)D spatiotemporal dynamics, one may consider optical propagation in a planar glass waveguide (e.g., see the experimental set-up of Ref. [20]), or a quadratic lithium niobate crystal, in the regime of large phase-mismatch, which mimics an effective Kerr nonlinear regime (e.g., see the experimental set-up of Ref. [21]). As far as the (2+1)D spatial dynamics is concerned, one may consider using a CW Ti:sapphire laser pulse propagating in a nonlinear medium composed of atomic-rubidium vapor (e.g., see the experimental set-up of Ref. [14]), or a bulk quadratic lithium niobate crystal, again in the regime of large phase-mismatch (e.g., see the experimental set-up of Ref. [22, 23]).

## 6 Conclusions

We have analytically predicted a new class of dark solitary wave solutions, that describe non-diffractive and non-dispersive spatiotemporal localized wave packets propagating in optical Kerr media. We numerically confirmed the existence of nonlinear lines, X-waves, lumps and peculiar scattering interactions of the solitary waves of the (2+1)D NLSE. The key novel property of these solutions is that their existence and interactions are inherited from the hydrodynamic soliton solutions of the well known KP equation. Our findings open a new avenue for research in spatiotemporal extreme nonlinear optics. Given that deterministic rogue and shock wave solutions, so far, have been essentially restricted to (1+1)D models, future research on multidimensional spatiotemporal nonlinear waves will lead to a substantial qualitative enrichment of the landscape of extreme wave phenomena.

## Acknowledgments

We acknowledge the financial support of the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 740355).

## References

- [1] R. Boyd, Nonlinear Optics, 3rd ed. (Academic Press, London, 2008).
- [2] A. Weiner, Ultrafast Optics (Wiley, New York, 2009).
- [3] H.E. Hernandez-Figueroa, E. Recami, and M. Zamboni-Rached, Localized Waves (Wiley, New York, 2008).
- [4] D. Majus, G. Tamosauskas, I. Grazuleviciute, N. Garejev, A. Lotti, A. Couairon, D. Faccio, and A. Dubietis, “Nature of Spatiotemporal Light Bullets in Bulk Kerr Media,” Phys. Rev. Lett. 112, 193901 (2014).
- [5] J. Durnin, J.J. Miceli, and J. H. Eberly, “Diffraction-free beams”, Phys. Rev. Lett. 58, 1499 (1987).
- [6] Y. Silberberg, “Collapse of optical pulses,” Opt. Lett. 15, 1282 (1990).
- [7] X. Liu, L. J. Qian, and F.W. Wise, “Generation of optical spatiotemporal solitons,” Phys. Rev. Lett. 82, 4631 (1999).
- [8] F. Baronio, S. Wabnitz, and Y. Kodama, “Optical Kerr Spatiotemporal Dark-Lump Dynamics of Hydrodynamic Origin,” Phys. Rev. Lett. 116, 173901 (2016).
- [9] F. Baronio, S. Chen, M. Onorato, S. Trillo, S. Wabnitz, and Y. Kodama, “Spatiotemporal optical dark X solitary waves,” Opt. Lett. 41, 5571-574 (2016).
- [10] F. Baronio, M. Onorato, S. Chen, S. Trillo, Y. Kodama, and S. Wabnitz, “Optical-fluid dark line and X solitary waves in Kerr media,” Opt. Data Process. and Storage 3, 1-7 (2017).
- [11] B.B. Kadomtsev and V.I. Petviashvili, “On the stability of solitary waves in weakly dispersing media,” Sov. Phys. - Dokl. 15, 539-541 (1970).
- [12] E.A. Kuznetsov and S.K. Turitsyn, “Instability and collapse of solitons in media with a defocusing nonlinearity,” Sov. Phys. JEPT 67, 1583-1588 (1988).
- [13] D.E. Pelinovsky, Y.A. Stepanyants, and Y.S. Kivshar, “Self-focusing of plane dark solitons in nonlinear defocusing media,” Phys. Rev. E 51, 5016 (1995).
- [14] V. Tikhonenko, J. Christou, B. Luther-Davies, and Y. S. Kivshar, “Observation of vortex solitons created by the instability of dark soliton stripes,” Opt. Lett. 21, 1129-1131 (1996).
- [15] C. Conti, S. Trillo, P. Di Trapani, G. Valiulis, A. Piskarskas, O. Jedrkiewicx, and J. Trull, “Nonlinear Electromagnetic X Waves”, Phys. Rev. Lett. 90, 170406 (2003).
- [16] Y. Kodama, “KP Solitons in shallow water,” J. Phys. A: Math. Theor. 43, 434004 (2010).
- [17] W. Li, H. Yeh, and Y. Kodama, “On the Mach reflection of a solitary wave: revisited,” J. Fluid Mech. 672, 326-357 (2011).
- [18] M. J. Ablowitz and H. Segur, Solitons and the inverse scattering transform, SIAM Stud. in Appl. Math. (SIAM, Philadelphia, 1981).
- [19] H. C. Yuen and B. M. Lake, “Instabilities of Waves on Deep Water,” Ann. Rev. Fluid Mech. 12, 303-334 (1980).
- [20] H. S. Eisenberg, R. Morandotti, Y. Silberberg, S. Bar-Ad, D. Ross, and J. S. Aitchison, “Kerr Spatiotemporal Self-Focusing in a Planar Glass Waveguide,” Phys. Rev. Lett. 87, 043902 (2001).
- [21] F. Baronio, C. De Angelis, M. Marangoni, C. Manzoni, R. Ramponi, and G. Cerullo, “Spectral shift of femtosecond pulses in nonlinear quadratic PPSLT Crystals,” Opt. Express 14, 4774-4779 (2006).
- [22] F. Baronio, C. De Angelis, P. Pioger, V. Couderc, A. Barthelemy, “Reflection of quadratic solitons at the boundary of nonlinear media,” Opt. Lett. 29, 986-988 (2004).
- [23] K. Krupa, A. Labruyere, A. Tonello, B. M. Shalaby, V. Couderc, F. Baronio, and A. B. Aceves, “Polychromatic filament in quadratic media: spatial and spectral shaping of light in crystals,” Optica 2, 1058-1064 (2015).