Topological Control of Extreme Waves

TL;DR

This paper demonstrates the first experimental topological control of extreme nonlinear waves, enabling transitions between different wave states by shaping wave genus through nonlinear Schrödinger equation solutions.

Contribution

It introduces a novel topological control method for nonlinear waves, experimentally observing supervised transitions between wave states with different genera.

Findings

First experimental observation of topological control of extreme waves.

Controlled transition from dispersive shock to rogue waves.

Shaping wave genus via time-dependent nonlinearity in a photorefractive medium.

Abstract

From optics to hydrodynamics, shock and rogue waves are widespread. Although they appear as distinct phenomena, new theories state that transitions between extreme waves are allowed. However, these have never been experimentally observed because of the lack of control strategies. We introduce a new concept of nonlinear wave topological control, based on the one-to-one correspondence between the number of wave packet oscillating phases and the genus of toroidal surfaces associated with the nonlinear Schr\"odinger equation solutions by the Riemann theta function. We prove it experimentally by reporting the first observation of supervised transitions between extreme waves with different genera, like the continuous transition from dispersive shock to rogue waves. Specifically, we use a parametric time-dependent nonlinearity to shape the asymptotic wave genus. We consider the box problem in a focusing Kerr-like photorefractive medium and tailor time-dependent propagation coefficients, as nonlinearity and dispersion, to explore each region in the state-diagram and include all the dynamic phases in the nonlinear wave propagation. Our result is the first example of the topological control of integrable nonlinear waves. This new technique casts light on dispersive shock waves and rogue wave generation, and can be extended to other nonlinear phenomena, from classical to quantum ones. The outcome is not only important for fundamental studies and control of extreme nonlinear waves, but can be also applied to spatial beam shaping for microscopy, medicine and spectroscopy, and to the broadband coherent light generation.

Equations

$$ıϵ{\partial }_{\zeta }\psi +\frac{{ϵ}^2}{2}{\partial }_{\xi }^2\psi +|\psi {|}^2\psi =0,$$

$$\psi(\xi,0)=\left\{\begin{array}[]{lr}q&\mbox{for }|\xi|<l\\ 0&\mbox{elsewhere}\end{array}\right.,$$

$$\begin{array}[]{l}\zeta=\frac{z}{\epsilon z_{\mathrm{D}}},\;\;\xi=\frac{2x}{W_{0}},\;\;\psi=\frac{A}{\sqrt{I_{0}}},\end{array}$$

$$ϵ=\frac{\lambda }{\pi W_0}\sqrt{\frac{I_{\mathrm{S}}}{2n_0\delta n_0{I}_{0}f(t)}},$$

$$x=x_0\pm \frac{W_0}{2t_0}(t-t_0)=x_0\pm v(t-t_0),$$

$$v=\frac{W_0}{2t_0}=\frac{32\delta n_0{L}^2}{I_{\mathrm{S}}{n}_0{W}_0^2{U}_0\tau }P,$$

$$ı{\partial }_{z}A+\frac{1}{2k}{\nabla }^{2}A+\frac{k}{n_0}\delta n(I)A=0,$$

$$ı{\partial }_{z}A+\frac{1}{2k}{\partial }_x^{2}A+2\rho (t)|A{|}^{2}A=0,$$

$$A(x,0)=\left\{\begin{array}[]{cc}\sqrt{I_{0}}&\mbox{for }|x|\leq\frac{1}{2}W_{0}\\ 0&\mbox{elsewhere}\end{array}\right..$$

$$\psi (\xi ,\zeta =0)=q\exp \left\{-\frac{1}{2}{\left(\frac{\xi }{l}\right)}^{24}\right\}.$$

Full text

I Introduction

In 1967 Gardner, Greene, Kruskal, and Miura developed a mathematical method - the inverse scattering transform (IST) 1967Gardner - disclosing the inner features of nonlinear waves in hydrodynamics, plasma physics, nonlinear optics and many other physical systems 1972Zakharov ; 1974Gardner ; 1974Ablowitz . According to IST, one also predicts the periodical regeneration of the initial status, as in the Fermi-Pasta-Ulam-Tsingou recurrence 2018Pierangeli ; 2018Mussot .

The nonlinear Schrödinger equation (NLSE) 1987Akhmediev is a cornerstone of IST for detailing dispersive phenomena, as dispersive shock waves (DSWs) 2006Hoefer ; 2007Wan ; 2016ElHoefer , rogue waves (RWs) 2007Solli ; 2011Bonatto ; 2014Dudley ; 2016Akhmediev and shape invariant solitons Solitons ; OpticalSolitons ; 2014Bulso . DSWs regularize catastrophic discontinuities by mean of rapidly oscillations 2007Ghofraniha ; 2017Marcucci ; 2007El ; 2012GhofranihaGentilini ; 2013Kartashov . RWs are giant disturbances appearing and disappearing abruptly in a nearly constant background 2010Erkintalo ; 2011Erkintalo ; 2012Chabchoub ; 2013Akhmediev ; 2015Agafontsev ; 2015Kibler ; 2015Pierangeli ; 2015Toenger ; 2016Suret ; 2017Tikan ; 2018Gelash ; 2019Gelash . Solitons are particle-like dispersion-free wave packets that can form complex interacting assemblies, ranging from crystals to gases 1971Zakharov ; Solitons ; OpticalSolitons ; 2005El ; 2018Gelash ; 2019Redor .

DSWs, RWs, and soliton gases (SGs) are related phenomena, and all appear in paradigmatic nonlinear evolutions, as the box problem for the focusing NLSE 2010Wan ; 2014Picozzi ; 2016ElKhamis ; 2018Audo ; 2018Biondini ; 2019Biondini . However, for the box problem in the small-dispersion NLSE, IST becomes unfeasible. In this extreme regime, the problem can be tackled by the so-called finite-gap theory 2016Bertola ; 2016ElKhamis . It turns out that extreme waves are described in terms of one single mathematical entity, the Riemann theta function, and classified by a topological index, the genus $g$ (see Fig. 1). In nonlinear wave theory, $g$ represents the number of oscillating phases, and evolves during light propagation: “single phase” DSWs have $g=1$, RWs have $g\sim 2$ and SGs have $g>>2$. This creates a fascinating connection between extreme waves and topology. Indeed, the same genus $g$ allows a topological classification of surfaces, to distinguish, for examples, a torus and sphere (Fig. 1). The question lies open if this elegant mathematical classification of extreme waves can inspire new applications. Can it modify the basic paradigm by which the asymptotic evolution of a wave is encoded in its initial shape, opening the way to controlling extreme waves, from lasers to earthquakes?

Here, inspired by the topological classification, we propose and demonstrate the use of topological indices to control the generation of extreme waves with varying genera $g$ 2018Audo . We consider the NLSE box problem where, according to recent theoretical results 2016ElKhamis , light experiences various dynamic phases during propagation, distinguished by diverse genera. In particular, for high values of a nonlinearly-scaled propagation distance $\zeta$, one has $g\sim\zeta$. By continuously varying $\zeta$, we can change $g$ and explore all the possible dynamic phases (see Fig. 1, where $\zeta$ is given in terms of the observation time $t$, detailed below). We experimentally test this approach in photorefractive materials, giving evidence of an unprecedent control of nonlinear waves, which allows the first observation of the transition from focusing DSWs to RWs.

II Results

II.1 Time-dependent Spatial Box Problem

We consider the NLSE

[TABLE]

where $\psi=\psi(\xi,\zeta)$ is the normalized complex field envelope, $\zeta$ is the propagation coordinate, $\xi$ is the transverse coordinate and $\epsilon>0$ is the dispersion parameter. We take a rectangular barrier as initial condition

[TABLE]

that is, a box of finite height $q>0$, length $2l>0$ and genus $g=0$. In our work, we fix $q=l=1$. Eq. (1) with (2) is known as the NLSE box problem, or the dam break problem, which exhibits some of the most interesting dynamic phases in nonlinear wave propagation 2016ElKhamis ; 2017Xu . The initial evolution presents the formation of two wave trains counterpropagating that regularize the box discontinuities. These wave trains are single-phase DSWs ($g=1$). Their two wavefronts superimpose in the central part of the box (see Fig. 1a) - occurring at $\zeta=\zeta_{0}:=\frac{l}{2\sqrt{2}q}$ - and generate a breather lattice of genus $g=2$, a two-phase quasi-periodic wave resembling an ensemble of Akhmediev breathers (ABs) 2014Dudley ; 2015Kibler . Since both $\xi-$ and $\zeta-$period increase with $\zeta$, the oscillations at $\xi\simeq 0$ become locally approximated by Peregrine solitons (PSs) 2010Kibler ; 2014Dudley ; 2018Randoux ; 2019Xu . At long propagation distance $\zeta>>\zeta_{0}$, the wave train becomes multi-phase and generates a SG with $g\sim\zeta$.

In Figure 1a, we report the wave dynamics in physical units, as we make specific reference to our experimental realization of the NLSE box problem for spatial optical propagation in photorefractive media (PR). In these materials, the optical nonlinearity is due to the time-dependent accumulation of free carriers that induces a time-varying low-frequency electric field. Through the electro-optical effect, the charge accumulation results into a time-varying nonlinearity. The corresponding time-profile can be controlled by an external applied voltage and the intensity level 1998DelRe ; 2006DelRe ; 2009DelRe . These features enable to experimentally implement our topological control technique. In PR, Eq. (1) describes an optical beam with complex amplitude $A(z,x,t)$ and intensity $I=|A|^{2}$ through the transformation (see Methods)

[TABLE]

with $W_{0}$ the initial beam waist along $x$-direction, $z_{\mathrm{D}}=\frac{\pi n_{0}W_{0}^{2}}{2\lambda}$ the diffraction length, $n=n_{0}+\frac{2\delta n_{0}I}{I_{\mathrm{S}}}f(t)$ the refractive index, $\delta n_{0}>0$ the nonlinear coefficient, $I_{\mathrm{S}}$ the saturation intensity, $I_{0}$ the initial intensity. For PR

[TABLE]

namely, the dispersion is modulated by the time-dependent crystal response function $f(t)=1-\exp\left(-t/\tau\right)$, with the saturation time $\tau$ fixed by the input power and the applied voltage 2018Pierangeli .

II.2 Genus Control

For a given propagation distance $L$ (the length of the photorefractive crystal), the genus of the final state is determined by the detection time $t$, which determines $\epsilon$, $\zeta=\frac{L}{\epsilon z_{\mathrm{D}}}$, and $g$, correspondingly. The genus time-dependence is sketched in Fig. 1a. The output wave profile depends on its genus content, which varies with $t$.

Following the theoretical approach in 2016ElKhamis , the two separatrix equations divide the evolution diagram in Fig. 1a in three different areas: the flat box plateau with genus $g=0$, the lateral counterpropagating DSWs with genus $g=1$, and the RWs after the DSW-collision point (corresponding to the separatrices intersection) with genus $g=2$. The two separatrices (dashed lines in Fig. 1a) have equations

[TABLE]

with $(t_{0},x_{0})$ the DSW-collision point, with $t_{0}\simeq\frac{\tau I_{\mathrm{S}}n_{0}W_{0}^{2}}{64I_{0}\delta n_{0}L^{2}}$ and $x_{0}$ given by the central position of the box. It turns out that the shock velocity is

[TABLE]

proportional to the input power, as experimentally demonstrated and detailed below.

Eqs. (5) express the genus time-dependence for its first three values $g=0,1,2$. It allows designing the waveshape, before the experiment, by associating a specific combination of the topological indices, and to predict the detection time corresponding to the target topology. In other words, by properly choosing the experimental conditions, we can manage to predict the occurrence of a given extreme wave by using the expected genus $g$. According to Eq. (4), we use time $t$ and initial waist $W_{0}$ to vary $\epsilon$. The accessible states are outlined in the phase diagram in Fig. 1b, in terms of $\epsilon$ and $W_{0}$. Choosing $W_{0}=100\mu$m as in Fig. 1a, by varying $t$ one switches from DSWs to RWs, and then to SGs.

II.3 Supervised Transition from Shock to Rogue Waves

The case $W_{0}=140\mu$m is illustrated in Fig. 2a by numerical simulations.The two focusing DSWs and the SG are visible at the beginning and at the end of temporal evolution, respectively (see phase-diagram in Fig. 1b). As soon as an initial super-Gaussian wave (Fig. 2b, see Methods) starts to propagate, two DSWs appear on the beam borders (Fig. 2c) and propagate towards the beam central part (Fig. 2d). Experimental proof of the genuine nonlinear nature of the beam evolution at this regime, not due to modulation instability arising from noise in the central part of the box, is shown in Supplementary Information (Suppl. Fig. 1). When the DSWs superimpose, ABs are generated (Fig. 2e). From the analytical NLSE solutions for the focusing dam break problem 2016ElKhamis , we see that ABs have $\xi$-period increasing with $\zeta$. Moreover, one finds that $\partial_{t}\zeta>0$, therefore the period in the $x$-direction must increase with time, and central peaks appear upon evolution. These peaks are well approximated by PSs, for large $t$, as confirmed by Figs. 2f,g.

The occurrence of RWs in the large box regime is proved also by statistical analysis, illustrated in Supplementary Information (Suppl. Figs. 2h,i).

Figure 3 shows the experimental observation of the controlled dynamics simulated in Fig. 2. Figure 3a sketches the experimental setup, detailed in Methods. A quasi-one-dimensional box-shaped beam propagates in a photorefractive crystal, and the optical intensity distribution is detected at different times. The observations of shock velocities and beam propagation for $W_{0}=140\mu$m are reported in Figs. 3b,c, respectively. In Fig. 3c, we see an initial DSW phase that undergoes into a train of large amplitude waves. In this regime, we identify breather-like structure (ABs, inset in Fig. 3c) that evolves into a SG at large propagation time. The DSW phase is investigated varying the input power. We find a linear increasing behavior of the shock velocity when increasing the power (Fig. 3b), as predicted by Eq. (6). The shock velocity is proportional to the distance between the two counterpropagating DSWs at a fixed time. We measured the width $\Delta x$ of the plateau at time $\bar{t}\sim 30$s. Referring to Eq. (6), we obtain the normalized velocity $\bar{v}=v/v_{0}$, with $v_{0}=L/\bar{t}$.

II.4 Peregrine Solitons Emergence

Figure 4 illustrates the numerically determined dynamics at smaller values of the beam waist ($W_{0}=10\mu$m), a regime in which the generation of single PSs is evident. The intensity profile is reported in Fig. 4a. As shown in Fig. 1b, one needs to carefully choose $W_{0}$ for observing a RWs generation without the DSWs occurrence. For $W_{0}=10\mu$m, the super-Gaussian wave (Fig. 4b) generates a PS (Figs. 4c-e). The following dynamics shows the higher-order PS emergence (Figs. 4f,g), each order with a higher genus.

Figures 5a-g report the experimental results for the case $W_{0}=30\mu$m. Observations of the Peregrine-like soliton generation are shown, both in intensity (Figs. 5a-d) and in phase (Figs. 5e-g). For a small initial waist, a localized wave, well described by the PS (Figs. 5b,d), forms and recurs without a visible wave breaking. This dynamics is in close agreement with simulations in Figs. 4d-g, where the PS is repeatedly destroyed and generated, each time at a higher order. Phase measurements are illustrated in Figs. 5e-g. Each PS has two phase signatures: a longitudinal smooth phase shift of $2\pi$ and a transversal rectangular phase shift profile, with height $\pi$ and basis as wide as the PS width 2018Randoux ; 2019Xu . Such signatures are here both experimentally demonstrated. From Fig. 5e, which shows interference pattern during the first PS occurrence, we obtain the longitudinal phase shift behavior in Fig. 5g, by a cosinusoidal fitting along the central propagation outline. Fig. 5f reports the experimental transversal phase shift profile along $x$. A comparison with the measured interference fringes is also illustrated in the inset, which directly shows the phase jump (topological defect). Stressing the significance of these results is very important, because they are a proof of the topological control: the genus is determined by the input waist and time of detection. Indeed, the longitudinal phase shift represents the transition from genus [math] to genus $2$, whereas the transversal PS phase shift outline unveils the value $g=2$, equal to the number of phase jumps (first from [math] to $\pi$, then again from $\pi$ to [math]). This is summarized in Fig. 5h, which sketches numerical simulations of phase behavior at $W_{0}=10\mu$m, normalized in $[-\pi,\pi]$. Fig. 5h gives a picture of genera changes, PS occurrence and phase discontinuities. The genus is zero and the phase profile is flat until the first PS occurrences. After that, the phase value changes and the phase transverse profile presents two jumps of $\pi$.

The statistical properties of the PS intensity are illustrated in Supplementary Information (Suppl. Figs. 2f,g), and they confirm the occurrence of RWs in the small box regime.

III Discussion

The topological classification of nonlinear beam propagation by the genera of the Riemann theta functions opens a new route to control the generation of extreme waves. We demonstrated the topological control for the focusing box problem in optical propagation in photorefractive media. By using the time-dependent photorefractive nonlinearity, we could design the final state of the wave evolution in a predetermined way and explore all the possible dynamic phases in the nonlinear propagation.

Such a novel control strategy enabled the first observation of the continuous transition from dispersive shock to rogue waves and soliton gases, demonstrating that different extreme wave phenomena are deeply linked, and also that a proper tuning of their topological content in their nonlinear evolution allows transformations from one state to another. The further numerical and experimental analysis reported in Supplementary Information proves that this new control paradigm in third-order media has a broad range of validity, where it is not affected by linear effects, like modulation instability or loss, but its nature is genuinely nonlinear.

In conclusion, our result is the first example of the topological control of integrable nonlinear waves. This new technique casts light on dispersive shock waves and rogue wave generation. It is general, not limited to the photorefractive media, and can be extended to other nonlinear phenomena, from classical to quantum ones. These outcomes are not only important for fundamental studies and control of extreme nonlinear waves, but further developments in the use of topological concepts in nonlinear physics can allow innovative applications for engineering strongly nonlinear phenomena, as in spatial beam shaping for microscopy, medicine and spectroscopy, and coherent supercontinuum light sources for telecommunication.

IV Methods

IV.1 Photorefractive Media

Starting from Maxwell’s equations in a medium with a third-order-nonlinear polarization, in paraxial and slowly varying envelope approximations, one can derive the propagation equation of the complex optical field envelope $A(x,y,z)$:

[TABLE]

with $z$ the longitudinal coordinate, $x,y$ the transverse coordinates and $n=n_{0}+\delta n(I)$ the refractive index, weakly depending on the intensity $I=|A|^{2}$ $\left(\delta n(I)<<n_{0}\right)$.

Eq. (7) is the nonlinear Schrödinger equation (NLSE) and rules laser beam propagation in centrosymmetric Kerr media. For PR, the refractive index perturbation depends also parametrically on time, i.e., $\delta n=\delta n(I,t)$. In fact, the amplitude of the nonlinear self-interaction increases, on average, with the exposure time up to a saturation value, on a slow timescale, typically seconds for peak intensities of a few kWcm*-2* 2009DelRe .

In our centrosymmetric photorefractive crystal, at first approximation $\delta n=\frac{-\delta n_{0}}{\left(1+\frac{I}{I_{\mathrm{S}}}\right)^{2}}f(t)$, with $f(t)$ the response function. $\delta n_{0}$ includes the electro-optic effect coefficient 1998DelRe ; 2006DelRe ; 2009DelRe . For weak intensities $I<<I_{\mathrm{S}}$, we obtain a Kerr-like regime with $\delta n=2\delta n_{0}\frac{I}{I_{\mathrm{S}}}f(t)$, apart from a constant term. We consider the case $\partial_{y}A\sim 0$ (strong beam anisotropy), thus we look for solutions of the $(1+1)$-dimensional NLSE for the envelope $A\sim A(x,z)$:

[TABLE]

with $\rho(t)=\frac{2\pi\delta n_{0}}{\lambda I_{\mathrm{S}}}f(t)$ and the field envelope initial profile

[TABLE]

One obtains Eq. (8) from Eq. (1) through the transformation (3). We stress that, in this case, the dispersion parameter depends on time, as follows from Eq. (4).

IV.2 Numerical Simulations

We solve numerically Eq. (1) by a one-parameter-depending beam propagation method (BPM) with a symmetrized split-step in the code core Datadriven . We use a high-order super-Gaussian initial condition

[TABLE]

For each temporal value, Eq. (1) solutions have different dispersion parameter $\epsilon$ and final value of $\zeta$, because from Eqs. (3) it reads $\zeta_{\mathrm{fin}}=\frac{4L}{\epsilon(t)kW_{0}^{2}}$, where $L$ is the crystal length. In Figs. 2,4, we show the numerical results. The propagation in time considers $\psi(\xi,\zeta_{\mathrm{fin}})$, which corresponds to detections at end of the crystal.

IV.3 Experimental Set-up

A $y$-polarized optical beam at wavelength $\lambda=532$nm from a continuous $80$mW Nd:YAG laser source is focused by a cylindrical lens down to a quasi-one-dimensional beam with waist $U_{0}=15\mu$m along the $y$-direction. The initial box shape is obtained by a mask of tunable width, placed in proximity of the input face of the photorefractive crystal. A sketch of the optical system is shown in Fig. 3a. The beam is launched into an optical quality specimen of $2.1^{(x)}\times 1.9^{(y)}\times 2.5^{(z)}$mm $K_{0.964}Li_{0.036}Ta_{0.60}Nb_{0.40}O_{3}$ (KLTN) with Cu and V impurities ($n_{0}=2.3$). The crystal exhibits a ferroelectric phase transition at the Curie temperature $T_{\mathrm{C}}=284$K. Nonlinear light dynamics are studied in the paraelectric phase at $T=T_{\mathrm{C}}+8$K, a condition ensuring a large nonlinear response and a negligible effect of small-scale disorder 2016Pierangeli . The time-dependent photorefractive response sets in when an external bias field $E$ is applied along $y$ (voltage $V=500\rm{V}$). To have a so-called Kerr-like (cubic) nonlinearity from the photorefractive effect, the crystal is continuously pumped with an $x$-polarized $15$mW laser at $\lambda=633$nm. The pump does not interact with the principal beam propagating along the $z$ axis and only constitutes the saturation intensity $I_{S}$. The spatial intensity distribution is measured at the crystal output as a function of the exposure time $t$ by means of a high-resolution imaging system composed of an objective lens ($NA=0.5$) and a CCD camera at $15$Hz.

In the present case, evolution is studied at a fixed value of $z$ (the crystal output) by varying the exposure time $t$. In fact, the average index change grows and saturates according to a time dependence well defined by the saturation time $\tau\sim 100$s once the input beam intensity, applied voltage, and temperature have been fixed.

V Data Availability

All data are available in this submission.

Acknowledgements

We acknowledge M. Conforti, S. Gentilini, P. G. Grinevich, A. Mussot, P. M. Santini, S. Trillo, and V. E. Zakharov for fruitful conversations on related topics. We thank MD Deen Islam for technical support in the laboratory. The present research was supported by PRIN2015 NEMO project (2015KEZNYM grant), H2020 QuantERA QUOMPLEX (grant number 731473), H2020 PhoQus (grant number 820392), and Sapienza Ateneo (2016 and 2017 programs).

Author Contribution

G.M. and D.P. equally contributed to this work. G.M., R.K.L, and C.C. conceived the idea and the theoretical framework; D.P., E.D., and C.C. conceived its experimental realization. G.M. and C.C. developed the theoretical background. G.M. performed the numerical simulations. D.P. carried out experiments and data analysis. A.J.A. designed and fabricated the photorefractive crystal. All authors discussed the results and wrote the paper.

Competing Interest

The authors declare no competing interests.

Additional Information

Supplementary Information accompanies this paper.

Correspondence and requests for materials should be addressed to G.M.