Modulational instability in optical fibers with randomly-kicked normal dispersion

G. Dujardin,1 A. Armaroli,2 S. Rota Nodari,3 A. Mussot,2 A. Kudlinski,2 S. Trillo,4 M. Conforti,2, ∗ and S. De Bi`evre1, † 1Univ. Lille, CNRS UMR 8524 - Laboratoire Paul Painlev´e,Inria, F-59000 Lille, France 2Univ. Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers Atomes et Mol´ecules,F-59000 Lille, France 3Institut de Math´ematiquesde Bourgogne (IMB), CNRS, UMR 5584, Universit´eBourgogne Franche Comt´e,F-21000 Dijon, France 4Department of Engineering, University of Ferrara, I-44122 Ferrara, Italy We study modulational instability (MI) in optical fibers with random group velocity dispersion (GVD) generated by sharply localized perturbations of a normal GVD fiber that are either randomly or periodically placed along the fiber and that have random strength. This perturbation leads to the appearance of low frequency MI side lobes that grow with the strength of the perturbations, whereas they are faded by randomness in their position. If the random perturbations exhibit a finite average value, they can be compared with periodically perturbed fibers, where Arnold tongues appear. In that case, increased randomness in the strengths of the variations tends to affect the Arnold tongues less than increased randomness in their positions.

I. INTRODUCTION fibers [30]. An attempt to consider a GVD perturbed by a Gaussian noise with a finite correlation length was re- The combined effect of nonlinearity and group veloc- ported in [31], but the analysis was not conclusive, and ity dispersion (GVD) may lead to the destabilization of the problem was solved with numerical simulations only the stationary states (plane or continuous waves) of a [32]. given physical system. This phenomenon, known under The question then arises for which type of random the name of modulational instability (MI), consists in GVD processes MI occurs and how the characteristics the exponential growth of small harmonic perturbations of such instabilities depend on the statistical properties of a continuous wave [1]. MI has been pioneered in the of the process. In such generality, the question seems 60s in the context of fluid mechanics [2, 3], electromag- however out of reach and it is consequently of inter- netic waves [4] as well as in plasmas [5], and it has been est to study the problem in a class of random fibers observed in nonlinear fiber optics in the 80s [6]. In uni- that is both experimentally accessible and theoretically form fibers, MI arises for anomalous (negative) GVD, tractable. Our focus here will therefore be on homoge- but it may also appear for normal GVD if polarization neous fibers with a normal GVD perturbed by a set of [7], higher order modes [8] or higher order dispersion are random “kicks”. More specifically, the fibers we consider, considered [9]. A different kind of MI related to a para- described in more detail in Section II, have a GVD given metric resonance mechanism emerges when the disper- by sion or the nonlinearity of the fiber are periodically mod-   X z − Zn ulated [10–13]. Many studies were published to address β2(z) = β2,ref + ∆β2 λnδ . (1) Zref generalizations to high-order dispersion [14, 15], birefrin- n∈Z gence [16], fiber cavities [17–23], and the nonlinear stage of MI [24]. In our previous work [33], we investigated such fibers with The effect of a random variation of GVD on MI has a periodic modulation of the fiber dispersion induced by been studied extensively [25–29] for the particular case a Dirac comb, in which the Zn = nZref are periodically where the GVD is perturbed by a Gaussian white noise, placed along the fiber and the strengths λn of the delta which is explicitly solvable. Under these conditions a de- functions are all equal. Such GVD was shown to be well formation of the conventional MI gain profile due to the approximated in experiments by a periodic series of short random perturbation was found when the unperturbed gaussian-like pulses. In this work, the Zn are chosen to be random points along the fiber, and the λn are in- arXiv:2103.02261v1 [physics.optics] 3 Mar 2021 fiber has an anomalous dispersion. In the case of nor- mal dispersion, the generation of MI sidebands as the dependent and identically distributed centered random result of the random perturbation was reported as well. variables and Zref is a reference length. We will limit our White noise however, which implies arbitrarily high vari- analysis to perturbations of fibers with a normal GVD ations of GVD on arbitrarily small scales, constitutes an because in that case the unperturbed fibers show no MI, idealization that does not always provide a relevant mod- and consequently any MI observed in the randomly per- eling of the randomness that may occur in physical GVD turbed fibers is entirely due to the randomness. The perturbations of the GVD considered in Eq. (1) result in a non-autonomous nonlinear Schr¨odingerequa- tion (NLSE) (see Eq. (2)) determining the evolution of ∗ [email protected] the wave profile as a function of the longitudinal coor- † [email protected] dinate z along the fiber in which the perturbation can 2 be interpreted as a succession of “kicks” taking place at finds points z = Zn. These systems therefore bear an analogy ! β2(z) 2 to the paradigmatic problem of the kicked rotor in clas- 0 − 2 ω ∂zxˆ(z, ω) = xˆ(z, ω). (3) sical and quantum mechanics [34–36], which is why refer β2(z) 2 2 ω + 2γP 0 to them as randomly GVD kicked fibers. We will show that MI occurs in such fibers through Note that this is, for each ω, a non-autonomous lin- a mechanism that is familiar from Anderson localization ear Hamiltonian dynamical system in a two-dimensional theory [37, 38]. The MI gain can indeed be computed phase plane with canonical coordinates (ˆq, pˆ). We wish in terms of a random product of transfer matrices in to study the (in)stability of its fixed point at the ori- the same way as the localization length of the station- ginq ˆ = 0 =p ˆ, as a function of the frequency ω and ary eigenfunctions of the random Schr¨odingeroperator of the properties of β2(z). Since for general β2(z), its in the Anderson model. We are then able to analyze explicit solution cannot be computed analytically, this how the properties of this gain depend on the features is not straightforward. We will concentrate on random of the random process β2(z) and on the frequency of the fiber profiles, as detailed below. harmonic perturbation of the continuous wave. Note that if γ = 0, Eq. (3) is reminiscent of a har- The rest of the paper is organized as follows. In Sec- β2(z) 2 monic oscillator with random frequency k = 2 ω tion II we precisely describe the randomly kicked fibers (also called multiplicative noise) for which a vast liter- under study, then we derive a general expression for the ature exists [30, 39–43]. The focus of the present work is mean MI gain in Section III. In Section IV, we develop to study a “generalised” random oscillator including the a perturbative estimation of the mean MI gain in the term 2γP ) which accounts for the nonlinear effects, where case where the random perturbations of the GVD vanish the random frequency is modeled by the non-stationary on average (λn = 0) and compare it to numerical re- stochastic processes described by Eq. (1). The study of sults. We establish in this manner the existence of MI at stationary colored noise, where classic perturbative tech- low frequencies of which we characterize the properties. niques apply [30], will be the subject of a future work We finally compare the MI gain with the gain computed [44]. from the solution of the nonlinear Schr¨odingerequation The modulational instability of the fiber is expressed and observe a good correspondence. In Section VI we in terms of the sample MI gain G(ω), defined as follows: consider MI in a randomly kicked homogeneous fiber of normal GVD in which the random perturbation does not 1 G(ω) = lim ln kxˆ(z, ω)k, (4) vanish on average (λn 6= 0). The situation is very dif- z→+∞ z ferent since, depending on the nature of the point pro- where k · k designates the euclidean norm. Here sam- cess Zn determining the positions of the kicks along the fiber, there may or there may not be a remnant of Arnold ple stands for a single realization of the random per- tongues, a signature of MI in periodic fibers. In Section V turbation. When G(ω) > 0, this indicates kxˆ(z, ω)k ' we show how GVD kicked fibers can approximate fibers exp(G(ω)z), meaning that the stationary solution is un- with a white noise GVD. Conclusions are drawn in Sec- stable for perturbations with frequency ω. One is inter- tion VII. ested in establishing for which ω, if any, this occurs, and how large G(ω) is in that case. Dispersion-kicked fibers are characterized by the ex- pression of β2(z) given in Eq. (1). Here δ is a sharply R II. MODULATIONAL INSTABILITY IN peaked positive function satisfying δ(z)dz = 1; in R RANDOMLY KICKED FIBERS our theoretical analysis below, we will take δ to be a Dirac delta function; λn are independent, identically dis- We consider the NLSE tributed real random variables with mean λ ≥ 0, and Zref > 0 is a characteristic length associated to the ran- 1 dom sequence of points Zn. We will write i∂ u − β (z)∂2u + γ|u|2u = 0, (2) z 2 2 t λn = λ + εδλn, (5) where γ > 0 is the fiber nonlinear coefficient, and with δλn = 0 and ε > 0 a dimensionless parameter. β2(z) its GVD. We are interested in the modula- We think of this as a fiber with irregularities in its di- √tional (in)stability of the stationary solution u0(z) = ameter of random area, giving rise to effective disper- P exp(iP γz) of Eq. (2). We consider a perturbation sion kicks |λn|∆β2Zref , placed at the points Zn along the of u0(z) in the form u(z, t) = [v(z, t)+1]u0(z), where the fiber. These fibers can be physically fabricated by means perturbation v(z, t) satisfies |v|  1. Writing v = q + ip, of the state-of-the-art fiber-drawing techniques. Indeed, with q and p real functions, inserting this expression some examples of uniformly dispersion-kicked fibers has into Eq. (2), and retaining only the linear terms, we been reported in [33] with a period Zref = 10 m and a obtain a linear system for q and p. Writing x(z, t) = kick width w = 0.14 m and relative large kick strength (q(z, t), p(z, t)) andx ˆ(z, ω) = √1 R x(z, t)e−iωt dt one 2π max |β2(z)|/β2,ref ≈ 35. 3

FIG. 1. Realization of a simple random walk fiber (left panel) and a Poisson fiber (central panel) with zero-mean kick ¯ ¯ strength λ√= 0, and of a simple random walk fiber with λ = 1 (right panel). The kicks are modeled with Gaussian pulses 2 2 2 δ(z) = 1/ 2πw exp(−z / 2w ) of width w = Zref /20. The position of the kicks is indicated by the vertical dotted lines; λn are uniform random variables in [−1, 1].

To compute the MI gain in such fibers, we proceed Hence it follows that for large ω, both Kn and Ln are as follows. We first note that, due to the presence of rotation matrices and consequently their product Φn is ± the delta functions, the left and right limitsx ˆn (ω) = also a rotation matrix. Consequently, for large ω, the MI ± ± (ˆq(Zn , ω), pˆ(Zn , ω)) of the solutionx ˆ(z, ω) at Zn are dif- gain tends to zero. + − ferent and are related byx ˆn (ω) = Knxˆn (ω), where the We finally describe the models we will consider for the random matrix Kn is defined as random positions Zn of the delta-functions, and for their strengths λn. For the Zn, we set ω2 ω2 ! cos(∆β2Zref λn 2 ) − sin(∆β2Zref λn 2 ) K = 2 2 . (6) n ω ω Z0 = 0, ∀n ∈ N,Zn+1 = Zn + Zref jn, (11) sin(∆β2Zref λn 2 ) cos(∆β2Zref λn 2 ) where Zref > 0 and jn is a sequence of independent and On the other hand, for Zn < z < Zn+1, (3) is au- identically distributed positive random variables with tonomous and straightforwardly solved; the solution values in + and with probability density ρ(j ). We − + R n xˆ(z, ω) is smooth in this range. One findsx ˆn+1 = Lnxˆn , assume that where now Z +∞   cos(k∆Zn) −µ sin(k∆Zn) hjni = jρ(j)dj = 1, (12) Ln = −1 , (7) 0 µ sin(k∆Zn) cos(k∆Zn) so that h∆Zni = h(Zn+1 − Zn)i = Zref > 0. Hence with Zn+1 = Zn + ∆Zn is a random walk with drift. We will β2,ref principally consider two cases. First, introducing a new β β  ω2 k2 = 2,ref ω2 2,ref ω2 + 2γP , µ = 2 . (8) parameter ε > 0, we consider 2 2 k Z jn = (1 + εZ δjn) , (13) Note that k is real for all ω when β2,ref > 0 (normal dispersion or defocusing NLSE), and that it is imaginary where now δjn is a sequence of independent, identically for small ω when β2,ref < 0 (anomaluous dispersion or and uniformly distributed random variables in (−1, 1) focusing NLSE). To sum up, we can now describe the and εZ ∈ [0, 1]. In this random walk model the mean − − evolution of this system between z = Zn and z = Zn+1 position of Zn is nZref and their variance as follows. For all ω ∈ and n ∈ R N 2 2 2 2 h(Zn − nZref ) i = nεZ Zref h(δj) i, − − − xˆ = LnKnxˆ = Φnxˆ . (9) n+1 n n grows with n. The increments are independent and iden- − tically distributed and their variance is given by Considering an initial conditionx ˆ0 = − − 2 − (ˆq(Z0 , ω), pˆ(Z0 , ω)) ∈ R with kxˆ0 k = 1, one then finds 2 2 2 2 h(∆Zn − Zref ) i = Zref εZ h(δj) i. 1 1 − G(ω) = lim ln kΦnΦn−1 ... Φ1xˆ0 k. (10) We will refer to this as the simple random walk model. Zref n→+∞ n Second, we will consider the Poisson model where Note that the ∆Zn = Zref jn are independent and identically dis- tributed with exponential density 2γP −4 0 ≤ µ ≤ 1, µ ' 1 − 2 + O(ω ), lim µ = 1. β2,refω ω→+∞ ρ(j) = exp(−j). (14) 4

− 2 3 − 2 The Zn are now the arrival times of a Poisson process Since kxˆn k ≤ kXnk1 ≤ 2 kxˆn k it readily follows that −1 with parameter Zref . We will refer to this as the Poisson fiber. 1 1 G(ω) = lim ln kXnk1, When all λn = 0 the above fibers are homogeneous. 2Zref n→+∞ n We will concentrate here on the defocusing regime, in which no modulational instability occurs when λn = 0. where now k(x, y, z)k1 = |x| + |y| + |z|. From Jensen’s Some examples of realizations of the random walk and inequality, it follows that of the Poisson fibers are reported in Fig. 1, where the 1 1 kicks are modeled by sharp Gaussian functions. Here and G(ω) = lim ln kXnk1 2Zref n→+∞ n in all numerical examples we take ∆β2 = β2,ref = γ = (17) P = Z = 1. 1 1 ref ≤ lim ln kXnk1 = G2(ω). In Section III and IV we consider random kick 2Zref n→+∞ n strengths λn that vanish on average. Our focus is there- fore on the question: what kind of random inhomo- We will refer to G2(ω) as the mean MI gain. It is worth geneities of the GVD can produce MI in an otherwise to point out that this corresponds to one half the growth modulationally stable homogeneous fiber? rate of the average power of the perturbations. In Section VI we then briefly discuss the case λ > As we now show, G2(ω) is larger than the sample 0: in that situation, MI occurs in the form of Arnold MI gain G(ω) itself but can be more readily computed. For that purpose, first note that, since |qˆ(Z−)ˆp(Z−)| ≤ tongues in the limiting case when εZ = 0, resulting in a n n 1 ((ˆq(Z−))2 + (ˆp(Z−))2), and since qˆ(Z−)ˆp(Z−) ≤ periodic GVD. We will investigate the stability of these 2 n n n n Arnold tongues under the random perturbations of the − − |qˆ(Zn )ˆp(Zn )|, we have Zn and/or λn which occur when εZ 6= 0 and/or ε 6= 0. 3 kXnk1 ≤ kXnk1 ≤ kXnk1. III. THE MEAN MI GAIN OF A RANDOMLY 2 KICKED FIBER Since the Mn are identically distributed and mutually in- n dependent, one has moreover Xn = Mn Xn−1 = M X0, Since the random walk process Zn in Eq. (12) has where independent and identically distributed increments the  2 2  Furstenberg theorem [45] (see [38] for a textbook treat- an 2anbn bn ment) asserts that the limit in Eq. (10) exists for almost M = ancn andn + bncn bndn (18) 2 2 every realization of the fiber, that is, for amost every cn 2cndn dn choice of the λn and the jn; it is in addition strictly pos- itive and independent of the realization. As a result, in is independent of n. Hence such random fibers, there is always MI at all values of ω. 1 1 However, it is notoriously difficult to obtain analytical G2(ω) = lim ln kXnk1 expressions for G(ω) as function of the model parame- 2Zref n→+∞ n ters, and hence to assess the strength of the sample MI 1 1 n 1 = lim ln kM X0k1 = ln |x|, gain G(ω). We will therefore follow [28] and introduce 2Zref n→+∞ n 2Zref a suitable mean MI gain based on moments of Eq. (9) (19) which we will refer to as G2(ω). We recall that no useful information can be obtained where x is the eigenvalue of M with largest modulus. from the evolution of first-order moments of Eq. (9). The The matrix M depends on the parameters of the second-order moments have instead to be computed. Let model, in particular on ω, the laws of λn, jn, ε and εZ . It follows that G (ω) can be computed from the spec-   2 an bn trum of M, which is easily determined numerically. The Φn ≡ LnKn = . cn dn explicit analytical computation of closed formulas for its spectrum remains however complicated. As we show be- It is then straightforward to check that for all n ∈ N, low, a perturbative treatment yields an analytic expres-

 − 2   − 2  sion provided that the perturbation is not too large. qˆ(Zn+1) qˆ(Zn ) − − − − Dropping the index n on λn,Zn, we introduce Xn+1 := qˆ(Zn+1)ˆp(Zn+1) = Mn qˆ(Zn )ˆp(Zn ) , − 2 − 2 2 pˆ(Z ) pˆ(Zn ) ω n+1 θ = ∆β Z λ , (20) (15) 2 ref 2 where which corresponds to the phase acquired by the pertur-  2 2  an 2anbn bn bationv ˆ(ω) at each kick, and write Mn = ancn andn + bncn bndn . (16) 2 2 2 2 cn 2cndn dn M = M+ cos θ + M0 sin(2θ) + M− sin θ, (21) 5

FIG. 2. Comparison of numerically computed MI gains in randomly kicked homogeneous fibers with λ = 0; εZ and ε are indicated in each panel. The random kicks are located at Zn as in Eq. (13) (random walk). The δλn and δjn are taken to be uniform in [−1, 1]. The numerically simulated sample MI gain G(ω) (dashed blue line), the mean MI gain G2(ω) computed directly from the largest eigenvalue of M (dotted red line), the perturbative approximation of the mean MI gain (dash-dot black line) and the perturbations’ growth rates (black stars) calculated from numerical solutions of the NLSE are shown. where IV. ZERO AVERAGE KICK AMPLITUDE cos2(k∆Z) −µsin(2k∆Z) µ2sin2(k∆Z) sin(2k∆Z) µ sin(2k∆Z) We consider in this section randomly kicked fibers as M+ =  2µ cos(2k∆Z) − 2  ,  2  in Eq. (1) with λ as in Eq. (5) and λ = 0 and with Z a sin (k∆Z) sin(2k∆Z) cos2(k∆Z) n n µ2 µ random process as in Eq. (11), so that the random fiber can be seen as a perturbation of a homogeneous defocus- µ2sin2(k∆Z) µsin(2k∆Z) cos2(k∆Z) ing fiber (see Fig. 1). We recall the latter is known to µ sin(2k∆Z) sin(2k∆Z) be modulationally stable. The MI gains of such random M− =  − 2 − cos(2k∆Z) 2µ  ,  2  fibers with Zn as in Eq. (13) and Eq. (14) are illustrated 2 sin(2k∆Z) sin (k∆Z) cos (k∆Z) − µ µ2 in Fig. 2 (for the random walk model) and in Fig. 3 (for the Poisson model) for various parameter values ε and and εZ as indicated. It is worth reminding that the sample MI gain is deterministic if calculated at infinite z. In the M0 = numerics we calculated an approximation of the sample  µ sin(2k∆Z) 2 2 2 µ sin(2k∆Z)  − 2 µ sin (k∆Z) − cos (k∆Z) 2 MI gain from Eq. (10) for n = 250 kicks and averaged cos(2k∆Z) 1 sin(2k∆Z) cos(2k∆Z)  2 −(µ + µ ) 2 − 2  . over 500 realizations (dashed blue curves in the figures).  2  sin(2k∆Z) cos2(k∆Z) − sin (k∆Z) − sin(2k∆Z) In order to check that the sample MI gain correctly pre- 2µ µ2 2µ dicts the growth rate of the perturbations we computed From the independence of the random variables it follows the latter from a numerical solution of the NLSE using that the same procedure (black stars in the figures). One ob- serves a good agreement. 2 2 M = M +cos θ + M 0sin 2θ + M −sin θ. (22) One notices that the random perturbation produces The expressions of the coefficients of the matrix M are modulational instabilities that we now further analyze, reported in Appendix B; the (in)stability depends on its for a given distribution of the ∆Zn and the λn, and at spectrum. Note that all information on the randomness fixed ω. This can be done perturbatively in ε. Indeed, in the strengths of the kicks is contained in θ, whereas the since λ = 0, we have λ = εδλ, which is small. We can randomness in the spacings between the kicks is encoded therefore compute the largest eigenvalue of M pertur- batively in ε; its logarithm will yield the mean MI gain in the matrices M0,M±. For a homogeneneous randomly kicked fiber with av- G2(ω). The details of the computation are given in Ap- erage normal GVD and vanishing mean kicking strength pendix A. (λ = 0) we resort to a perturbative analysis, presented To set up the perturbation problem, we proceed as below in Sec. IV. follows. For ε = 0, M = M + so that the first ingredient 6

FIG. 3. Same as Fig. 2, with Zn as in Eq. (14) (Poisson model). The parameter ε is indicated in the panels.

we need are the eigenvalues and eigenvectors of M +. It Recalling that λ = εδλ (See Eq. (5), with λ = 0), is easy to check that M + has the eigenvalue x0 = 1. The we consider the case where the probability distribution corresponding right and left eigenvectors are ν(δλ) of δλ satisfies ν(−δλ) = ν(δλ) so that not only λ = ε R sν(s)ds = 0 but also sin(2θ) = 0, since sin θ is an (0) 2
T (0) 2
T ϕ = µ 0 1 and ψ = 1 0 µ , (23) odd function of δλ. Hence, we can write

T with ψ(0) ϕ(0) = 2µ2. M = M + + ∆Mη, (28) We always suppose ω 6= 0 so that µ 6= 0. The two √ 2
remaining eigenvalues of M + are easily determined to be where η = cos(2θ) − 1 = 2π νˆ(ε∆β2Zref ω ) − νˆ(0) 1 and ∆M = (M + − M −). Z +∞ 2 Finally, for small εω2, x± = exp(±i2k∆Z) = exp(±i2kZref j)ρ(j)dj 0 √ 1 = 2πρˆ(∓2kZ ). η ' − (ε∆β Z ω2)2δλ2. (29) ref 2 2 ref (24)

Since the exp(i2kZref j) lie on the unit circle, clearly The eigenvalue xη of M in Eq. (28) that emanates from |x±| ≤ 1. To apply nondegenerate perturbation theory, x0 = 1 can be expanded as we need that x± 6= 1. To investigate this condition in (1) 2 (2) the two models that we investigate here for the Zn, let xη ≈ 1 + ηx + η x , (30) us first consider the random walk model Eq. (13), and assume the distribution of δj is given by a density σ, so where that 2 2 (1) (1 − µ ) Z x = − 2 , ±i2kZref ±2ikZref εZ s 4µ x± = exp(±i2k∆Z) = e e σ(s)ds. 1 2S(2S − 1) + S2 (25) x(2) = (1 − µ4)2 2 , (31) 16µ4 4S2 + S2 We can then compute x± perturbatively in εZ as follows. 2 Noting that, when εZ = 0, x± = exp(±i2kZref ) and with remembering that hδji = 0, one sees 2 1    2 2 3  S = sin (k∆Z) = 1 − cos(2k∆Z) ,S2 = sin(2k∆Z). x± = exp(±i2kZref ) 1 − 2(kZref εZ ) hδj i + O(εZ ) , 2 (26) (32) so that, with increasing εZ , the eigenvalues x± move ra- The corresponding mean MI gain is then, using Eq. (19) dially inward, towards the origin. Hence |x±| < 1 and, 1 1 a fortiori, x± 6= x0 = 1 for εZ 6= 0. For a Poisson fiber, (1) 2 (2) G2(ω) = ln xη ≈ ln(1+ηx +η x ). (33) where the Zn form a Poisson process, the eigenvalues x± 2Zref 2Zref can be computed explicitly: To lowest order in ε, we therefore find that, approxi- 1 mately, for small ω, x± = , (27) 1 ∓ i2kZref 2 2 2 2 1 (1 − µ ) 1 ∆β2 Zref γP 2 2 G2(ω) ≈ |η| 2 ≈ (εω) δλ . (34) so that |x±| < 1 for all ω 6= 0. 2 4µ 4 β2,ref We can therefore use, in the above cases, a non- degenerate perturbation expansion to compute x0 as a The mean MI gain G2(ω), computed by evaluating the function of ε. We will establish that x0 is real and a spectral radius of M numerically, as well as its approxi- growing function of ε, giving rise to a strictly positive mation from Eq. (33), are shown in Fig. 2 for the simple mean MI gain. random walk model and in Fig. 3 for the Poisson model. 7

There is, in both cases, a MI side lobe. The perturba- partially destroyed when εZ is nonzero, but small. It tive approximation works well for the Poisson model, for completely disappears when εZ approaches its maximal all ω in the range considered (see dashed black curves in possible value, which is 1, as can be seen in the third Fig. 3 ). It does however not capture the vanishing of column of the figure. the MI that occurs at a specific value of ω for εZ = 0 or Since for the homogeneous fiber (β2(z) = β2,ref, εZ = small in the random walk model (see Fig. 2). There is in- 0 = ε) there is no MI at all, it is clear that all MI is cre- deed a marked difference between the shape of these side ated by the randomness. Note however that, whereas in- lobes, depending on which of the two random processes creased fluctuations in the kicking strengths λn increases are chosen for the Zn, that we now explain. the MI, increased fluctuations in the Zn tends to decrease For that purpose, first consider the leftmost column it. of Fig. 2. There εZ = 0, which means the Zn = nZref In Fig. 3 the MI for the Poisson model is illustrated. are distributed periodically along the fiber. The GVD One sees that, as for the simple random walk model, is nevertheless not periodic, since the kick strengths λn there is an MI side lobe starting at low frequencies, but are random. One notices on the figure that, in that case, the frequencies ω` do now no longer play a special role. both the sample MI gain and the mean MI gain show The gain of this side lobe is comparable, in width and a characteristic zero at a precise value of ω. This phe- height, for the same value ε = 0.5 of the strength of the nomenon can be explained as follows. If ω = ω` is such kicks, as in the random walk fiber with εZ = 0.9. Note that kZref = π` for some ` ∈ Z, then the matrices Ln that, in the latter, this means successive Zn can be close, in Eq. (7) are all equal to the identity matrix, and one as in the Poisson fiber. The perturbative treatment of immediately sees from Eq. (10) that G(ω`) = 0, for all ε. the previous section reproduces the mean MI gain quite Indeed the propagation through the constant dispersion accurately. segments of the fiber does not change the perturbation and the kicks act as random rotations, which only change the phase of the perturbationv ˆ, giving as a results a van- V. APPROXIMATING WHITE NOISE GVD ishing MI gain. Using Eq. (8), one sees this corresponds WITH RANDOMLY KICKED GVD to the specific values of ω given by ω = ω`, where Previous work on MI in random fibers has concentrated s  2  π` 2 on GVD perturbed by white noise [26–28]. As pointed ω2 = (γP )2 + − γP . (35) ` β  Z  out above, white noise is not necessarily physically per- 2,ref ref tinent since it requires arbitrarily large variations of the GVD over arbitrarily short distances. In this section we One furthermore readily checks that the eigenvalues of will show that the kicked fibers considered here can, in M in this case are given by 1, cos θ, sin θ, which are less an appropriate parameter regime determined below, and than one in absolute value. Therefore the mean MI gain for sufficiently low frequencies, produce a similar MI gain G (ω ) also vanishes for all ε at ω = ω . Again, this is ap- 2 ` ` as a white noise GVD. Fig. 4 illustrates our findings. parent from Fig. 2 at ω = ω . Interestingly enough, the 1 For a kicked GVD as in (1), with λ = 0, the two-point frequencies determined from Eq. (35), fulfill the paramet- function is ric resonance condition kZref = π` [12, 24, 33]. These fre- 0 quencies correspond to the location of the tips of Arnold (β2(z) − β2,ref)(β2(z ) − β2,ref) = tongues for any periodic fiber with period Zref and whose ! 2 2 X 0 average GVD over one period equals β2,ref. (∆β2Zrefε) δλ ρn(z) δ(z − z ). Quite surprisingly, for kicks with random amplitudes n and zero-mean, the situation is reversed and the system Here the average is taken with respect to λ and Z and becomes stable under perturbations precisely at these n n ρ (z) is the probability distribution function of Z . Note same frequencies. n n that the fiber is therefore delta-correlated in z, but it is It is finally clear from Fig. 2 that the perturbative not stationary, since P ρ (z) is not constant. Never- treatment of the mean MI gain does not function well n n theless, one finds when ω approaches ω1 or ω2. This is as expected since, when ω = ω`, the three eigenvalues of M + coincide: X −1 lim ρn(z) = Zref , x± = 1 = x0. The nondegenerate perturbation theory z→+∞ n used above does then not apply. When ω 6= ω`, this degeneracy is lifted and the perturbation theory yields so that it does becomes stationary for large z. To show increasingly good results as εZ increases, even at ω = ω`. this, one can proceed as follows. Introducing the counting As can be seen in the second column of Fig. 2, both the function N(Z) = ]{n | Zn ≤ Z}, one easily sees that sample and mean MI gains are still diminished in the N(Z) ' Z . On the other hand, N(Z) = R Z n(z)dz, Zref 0 neighbourhood of ω1 and ω2 when εZ is small, and this P where n(z) = δ(z − Zm). So, since is well captured by the perturbative analysis above. This m phenomenon can be seen as a remnant of the underlying dN dN n(z) = (z), n(z) = (z), periodic structure of the random points Zn that is only dz dz 8

FIG. 4. Comparison between the mean MI gain of the white noise model (black dots) and the mean MI gain of the random walk (left and center panels)√ and Poisson models (right panel). MI gain for white noise model is calculated from Eq. (29) form 2 2 [25] with A = 1 (A ≡ P in our notation) and σ = 0.1. Parameters ε and εZ as indicated and Zref as in (37) with δλ = 1/3.

P it follows from n(z) = m ρm(z) that white noise model. Third, it would be a challenge to make fibers with a considerably larger value of ε since X dN 1 they correspond to a small value of Zref . This means lim ρn(z) = lim n(z) = (z) = . z→+∞ z→+∞ dz Zref that one would need to be able to put the sharp peaks n and dips in the fiber diameter very closely together. In conclusion, for z large, one has

0 2 0 (β2(z) − β2,ref)(β2(z ) − β2,ref) ≈ 2σ δ(z − z ), VI. NONZERO AVERAGE KICK AMPLITUDE where In this section we briefly discuss the effect on MI of 2σ2 = (∆β ε)2Z δλ2. (36) randomly placed kicks with random strengths λn that 2 ref are not vanishing on average, so that λ 6= 0. We will One recognizes here the two-point function of a Gaussian separately consider the two cases where the Zn form a white noise. In a kicked fiber, the fluctuations take place simple random walk or a Poisson process, and that lead to different phenomena. on a length scale comparable to Zref and have a strength We first consider the simple random walk, as in proportional to ε. This suggests that, if Zref is small and ε large, with a scaling given by Eq. (13), an example of which is reported in Fig. 1, rightmost panel. In the first column of Fig. 5 the numer- 2 2 2 2 Zref = 2σ /(∆β2 ε δλ ) (37) ically computed sample and mean MI gains are displayed for such a fiber for which εZ = 0, so that Zn = nZref , then the kicked fiber will be statistically close to a white i.e., the delta kicks are placed periodically. The λn are noise fiber and the resulting MI will therefore be similar chosen as in Eq. (5) with average λ = 1. Note that, if in in both fibers. This is indeed illustrated in Fig 4. The addition ε = 0, then this fiber is actually periodic: this is mean MI gain of a white noise fiber is plotted there (black the situation in the top left panel of the figure. It is well circles) using an explicit formula for this gain obtained known that periodic fibers display MI through so-called in [25]. It is compared to the mean MI gain of randomly Arnold tongues [12, 16, 46–48] and this has been experi- kicked fibers, with parameters ε and εZ as indicated and mentally shown in periodically kicked fibers [13, 33]. The with Zref as in (37), computed from the largest eigen- first such Arnold tongue can be observed in the top left value of M. As suggested by the above argument, for panel of Fig. 5. In the other panels of the first column, sufficiently large ε, the MI gain of the kicked fibers con- ε 6= 0 so that the kick strengths are now random, while verges to the one of the white noise fiber. The agreement the Zn remain periodic. The corresponding fibers can is best for small ω, a reflection of the fact that the low therefore be viewed as random perturbations of a periodic frequency perturbations are less sensitive to the rapid fiber. It appears from these data that the fluctuating kick variations of the white noise. strengths only weakly affect the position and strength of We make two further comments. First, the value of the Arnold tongues. In the other columns of Fig. 5, on σ2 = 0.1 is chosen in the numerics because it is the right the other hand, MI gains are shown for fibers for which order of magnitude for the fibers used in the experiments there is randomness in both kick positions Zn and kick described in [33] in which both ε and Zref are of order 1. strengths λn. It appears from these data that randomiz- Note that this means that the mean spacing between the ing the positions of the kicks has a much stronger effect kicks in these fibers is of the same order of magnitude as than randomizing their strengths. As the randomness in- −1 the nonlinear length ZNL = (γP ) . Second, for higher creases, one observes a marked decrease in the MI gain, but intermediate values of ε, the MI lobe of the randomly with a widening of the Arnold tongue, while its position kicked fibers can be higher and wider than the one of the is less affected. Note that the MI is considerably larger in 9

FIG. 5. Sample and mean MI gains for a simple random walk process and λ = 1 and δλn uniform in [−1, 1]. Color code is the same as Fig. 2 these fibers with λ = 1 than with those having λ = 0 as that, if the above condition is satisfied, then θ = nπ in can be seen by comparing the vertical scale of Fig. 5 with Eq. (20). Hence Eq. (21) implies M = M+ and conse- that of Fig. 2 and Fig. 3. In the first case, we are deal- quently M = M +. But we saw that the eigenvalues of ing with random perturbations of a periodic fiber, which M + are 1, x±, with |x±| ≤ 1. Hence the mean MI gain displays MI due to parametric resonance. In the second, also vanishes. the MI is on the contrary generated by the randomness. Let us stress that this phenomenon is not limited to In Fig. 6 the numerically computed sample and mean the Poisson fiber, for which it can be observed in Fig. 6. MI gains are displayed for a Poisson fiber with values In fact, for the simple random walk fiber it is visible in of ε as indicated. For ε = 0 one observes the identical the top panels of Fig. 5. The phenomenon also appears vanishing of these gains at ω ≈ 2.5. This is actually a in periodically kicked fibers where more general phenomenon, not related to the distribution   of the kicks, that is easily understood as follows. Let us X z − nZref β2(z) = β2,ref + ∆β2λ δ , consider a defocusing fiber with, for all n, λn = λ > 0 and Z n ref with arbitrary values for Zn. They could in particular be periodic, quasiperiodic, or random. If we now choose ωn so that the spatial average of the GVD is as Z 2 1 Z ref ωn β = β (z)dz = β + ∆β λ. nπ = λ∆β2Zref , (38) 2,av 2 2,ref 2 2 Zref 0 then we immediately see from (6) that Kn = ±I2. In It is then well known that the Arnold tongues for such a 0 other words, for these values of ω, the kicks have no ef- fiber occur at values ωn defined by fect on the linearized solution of the equation of motion. 0 2 One therefore has Φn = ±Ln and so the sample MI gain ω β Z n = p(γP Z )2 + (nπ)2 − γP Z . is equal to the sample MI gain of the unperturbed fiber. 2,av ref 2 ref ref Since the latter is defocusing, the√ sample MI gain van- ishes. When n = 1, one finds ω = 2π ≈ 2.5. That the Comparing this to Eq. (38), one observes that, for all n, 0 0 mean MI gain must also vanish follows from observing ωn < ωn. Note that both ωn and ωn depend on λ. 10

FIG. 6. Sample and mean MI gains for a Poisson model with λ = 1. Color code is the same as Fig. 2

Condition (38) is complementary to condition (35). produced through parametric resonance in fibers with a The former arises when the spacing between the kicks periodic GVD is considerable larger than the MI gain is constant, whereas the second appears when the am- obtained from random perturbations of the fibers. plitude of the kicks is. The two conditions show that whenever either the position or the amplitude of the kicks are deterministic, particular values of the frequency exist ACKNOWLEDGMENTS where the gain vanishes. The work was supported in part the by the French government through the Programme Investissement VII. DISCUSSION AND CONCLUSIONS d’Avenir with the Labex CEMPI (Grant ANR-11-LABX- 0007-01) and the I-SITE ULNE ( Grant ANR-16-IDEX- We reported on the modulation instability phe- 0004 ULNE, projects VERIFICO, EXAT, FUNHK) nomenon in optical fibers where the GVD is a random managed by the Agence Nationale de la Recherche. process. Most studies of such random fibers have con- The work was also supported by the Nord-Pas-de-Calais centrated on the case when a homogeneous GVD is per- Regional Council and the European Regional Develop- turbed by a stationary white noise. In that case vari- ment Fund through the Contrat de Projets Etat-R´egion´ ous methods exist to compute the MI gain. White noise (CPER) and IRCICA. is however very particular and not always adequate to model physically relevant scenarios. In this paper we in- vestigated the behavior of the MI gain in random fibers Appendix A: Perturbative analysis for which the GVD is of a very different nature. More specifically, we considered a class of experimentally re- The eigenvalue xη of M in Eq. (28) that emanates alizable [33] random fibers in which a constant normal from 1 can be expanded as GVD is perturbed by a sequence of delta kicks with ran- (1) 2 (2) 3 dom positions and amplitudes. xη = 1 + ηx + η x + O(η ), (A1) The main result of our analysis is fourfold. First, we T show that in this situation low frequency MI lobes al- and the corresponding eigenvectors of M and M , ways arise as a result of such random perturbations, in- (0) (1) 2 (2) 3 dependently of the statistical distribution of the position ϕη = ϕ + ηϕ + η ϕ + O(η ), and amplitudes of the kicks. We trace the occurrence of ψ = ψ(0) + ηψ(1) + η2ψ(2) + O(η3), MI in random fibers to a mechanism familiar from the η study of Anderson localization and of the occurrence of where ϕ(0), ψ(0) are defined in Eq. (23). Expanding the positive Lyapounov exponents in chaotic dynamical sys- eigenvalue equation in powers of η leads in the usual man- tems. Second, we show that the specific shape of these ner to the first correction to the eigenvalue: side lobes does depend on these statistical properties and we provide expressions to determine them. In particular, T 1 ψ(0) (M − M )ϕ(0) we find that if either the positions or the amplitudes of (1) + − x = T the kicks are deterministic, then at special frequency val- 2 ψ(0) ϕ(0) ues the MI gain is identically zero. At these same points, 1 (1 − µ2)2 = 2µ2 − (1 + µ4) = − . the MI remains suppressed when the random fluctuations 4µ2 4µ2 of the control parameter remain small enough. Third, we show that the randomly kicked fibers considered be- Note that this correction does not at all depend on the have, in a suitable parameter regime that we identify, distribution of the Zn nor of the λn. as a fiber with a white noise GVD. Finally, we have ob- To obtain a satisfactory expression for the mean MI served that for comparable parameter regimes, the MI gain G2(ω), we need to obtain the second order correction 11 x(2). For that purpose, we need the first correction ϕ(1) following averages of trigonometric functions. For the 3 to the eigenvector corresponding to xη, which is given by terms involving the kicks we have

(1) 1 h T (1)  T (1) i 1 Z 1  ω2  ϕ = ψ ϕ ϕ+ + ψ ϕ ϕ− 2 2 T + − cos θ = cos ∆β2Zref (λ + εx) dx ψ+ϕ+ 2 −1 2 2
2
where 1 1 cos λ∆β2Zref ω sin ε∆β2Zref ω = + 2 , T (0) 2 2 ε∆β2Zref ω T (1) ψ±(∆M)ϕ (0)T (1) ψ±ϕ = , ψ ϕ = 0. (A2) 1 − x±

T 2
2
Here ϕ± are the eigenvectors of M + and of M + corre- 2 1 1 cos λ∆β2Zref ω sin ε∆β2Zref ω sin θ = − 2 , sponding to the eigenvalues x±, given by 2 2 ε∆β2Zref ω

2
T ϕ± = −µ ±iµ 1 , and T ψ = 1 ±2iµ −µ2
, (A3) ± 2
2
sin λ∆β2Zref ω sin ε∆β2Zref ω T 2 sin 2θ = 2 . with ψ±ϕ± = −4µ . Finally, the second order correction ε∆β2Zref ω to the eigenvalue is For the four terms involving ∆Z, we have for the ran- ψ(0)T ∆MR∆Mϕ(0) x(2) = , (A4) dom walk model ψ(0)T ϕ(0) 2 1 1 sin(2kZref εZ ) with cos (k∆Z) = + cos(2kZref ) , 2 2 2kZref εZ  T T  T 1 ϕ+ψ+ ϕ−ψ− 2 ϕ+ψ+ R = T + = T Re . ψ+ϕ+ 1 − x+ 1 − x− ψ+ϕ+ 1 − x+ 2 1 1 sin(2kZref εZ ) sin (k∆Z) = − cos(2kZref ) , Hence 2 2 2kZref εZ 2 ψ(0)T ∆Mϕ ψT ∆Mϕ(0) x(2) = Re + + (0)T (0) T ψ ϕ ψ+ϕ+ 1 − x+ sin(2kεZ Zref ) (0)T T (0) cos (2k∆Z) = cos(2kZref ) , 1 ψ M ϕ ψ M ϕ 2kεZ Zref = Re − + + − . 2ψ(0)T ϕ(0)ψT ϕ 1 − x + + + and (A5)

(0)T sin(2kεZ Zref ) Here we used in the last line that ψ M +ϕ+ = 0 since sin (2k∆Z) = sin(2kZref ) . (0)T T (0) 2kεZ Zref ψ ϕ+ = 0 and similarly ψ+M +ϕ = 0. One finds

ψT M ϕ = (1 − µ4), ψT M ϕ(0) = (1 − µ4)x . For the four terms involving ∆Z, we have for the Pois- 0 − + + − + son model Hence 1 1 1 cos2 (k∆Z) = + , 1  1  2 2 (2) 4 2 2 2 1 + 4k Zref x = 4 (1 − µ ) 1 + Re 16µ x+ − 1 1 2S(2S − 1) + S2 = (1 − µ4)2 2 , (A6) 4 2 2 2 1 1 1 16µ 4S + S2 sin (k∆Z) = − 2 2 , 2 2 1 + 4k Zref where

2 1   S = sin (k∆Z) = 1 − cos(2k∆Z) ,S2 = sin(2k∆Z). 1 2 cos (2k∆Z) = 2 2 , (A7) 1 + 4k Zref

and Appendix B: Computation of M 2kZref sin (2k∆Z) = 2 2 . We describe here how to calculate the coefficients of 1 + 4k Zref the 3 × 3-matrix M in Eq. (22). They depend on the 12

[1] V. Zakharov and L. Ostrovsky, Modulation instability: S. Trillo, Parametric instabilities in modulated fiber ring The beginning, Physica D: Nonlinear Phenomena 238, cavities, Optics Letters 41, 5027 (2016). 540 (2009). [19] F. Copie, M. Conforti, A. Kudlinski, A. Mussot, and [2] T. B. Benjamin and J. E. Feir, The disintegration of wave S. Trillo, Competing Turing and Faraday Instabilities in trains on deep water Part 1. Theory, Journal of Fluid Longitudinally Modulated Passive Resonators, Physical Mechanics 27, 417 (1967). Review Letters 116, 143901 (2016). [3] V. E. Zakharov, Stability of periodic waves of finite am- [20] F. Copie, M. Conforti, A. Kudlinski, S. Trillo, and plitude on the surface of a deep fluid, Journal of Applied A. Mussot, Modulation instability in the weak disper- Mechanics and Technical Physics 9, 86 (1968). sion regime of a dispersion modulated passive fiber-ring [4] V. Bespalov and V. Talanov, Filamentary Structure of cavity, Optics Express 25, 11283 (2017). Light Beams in Nonlinear Liquids, Soviet Journal of Ex- [21] F. Copie, M. Conforti, A. Kudlinski, S. Trillo, and perimental and Theoretical Physics Letters 3, 307 (1966). A. Mussot, Dynamics of Turing and Faraday instabilities [5] Y. H. Ichikawa, T. Suzuki, and T. Taniuti, Modulation in a longitudinally modulated fiber-ring cavity, Optics Instability of Electron Plasma Wave, Journal of the Phys- Letters 42, 435 (2017). ical Society of Japan 34, 1089 (1973). [22] F. Copie, M. Conforti, A. Kudlinski, A. Mussot, F. Bian- [6] K. Tai, A. Hasegawa, and A. Tomita, Observation of calana, and S. Trillo, Instabilities in passive dispersion os- modulational instability in optical fibers, Physical Re- cillating fiber ring cavities, The European Physical Jour- view Letters 56, 135 (1986). nal D 71, 133 (2017). [7] A. L. Berkhoer and V. E. Zakharov, Self Excitation of [23] F. Bessin, F. Copie, M. Conforti, A. Kudlinski, A. Mus- Waves with Different Polarizations in Nonlinear Media, sot, and S. Trillo, Real-Time Characterization of Period- Soviet Journal of Experimental and Theoretical Physics Doubling Dynamics in Uniform and Dispersion Oscillat- 31, 486 (1970). ing Fiber Ring Cavities, Physical Review X 9, 041030 [8] R. H. Stolen, J. E. Bjorkholm, and A. Ashkin, Phase- (2019). matched three-wave mixing in silica fiber optical waveg- [24] M. Conforti, A. Mussot, A. Kudlinski, S. Rota Nodari, uides, Applied Physics Letters 24, 308 (1974). G. Dujardin, S. De Bi´evre,A. Armaroli, and S. Trillo, [9] S. B. Cavalcanti, J. C. Cressoni, H. R. da Cruz, and Heteroclinic Structure of Parametric Resonance in the A. S. Gouveia-Neto, Modulation instability in the region Nonlinear Schr¨odingerEquation, Physical Review Let- of minimum group-velocity dispersion of single-mode op- ters 117, 013901 (2016), arXiv:1606.01658. tical fibers via an extended nonlinear Schr¨odinger equa- [25] F. K. Abdullaev, S. S. S. A. Darmanyan, A. Kobyakov, tion, Physical Review A 43, 6162 (1991). and F. Lederer, Modulational instability in optical fibers [10] N. J. Smith and N. Doran, Modulational instabilities in with variable dispersion, Physics Letters A 220, 213 fibers with periodic dispersion management., Optics let- (1996). ters 21, 570 (1996). [26] F. K. Abdullaev, S. A. Darmanyan, S. Bischoff, and [11] M. Droques, A. Kudlinski, G. Bouwmans, G. Martinelli, M. P. Sørensen, Modulational instability of electromag- and A. Mussot, Experimental demonstration of modula- netic waves in media with varying nonlinearity, Journal tion instability in an optical fiber with a periodic disper- of the Optical Society of America B 14, 27 (1997). sion landscape., Optics letters 37, 4832 (2012). [27] F. K. Abdullaev and J. Garnier, Modulational instability [12] A. Armaroli and F. Biancalana, Tunable modulational of electromagnetic waves in birefringent fibers with peri- instability sidebands via parametric resonance in period- odic and random dispersion, Physical Review E 60, 1042 ically tapered optical fibers, Optics Express 20, 25096 (1999). (2012). [28] J. Garnier and F. K. Abdullaev, Modulational instability [13] A. Mussot, M. Conforti, S. Trillo, F. Copie, and induced by randomly varying coefficients for the nonlin- A. Kudlinski, Modulation instability in dispersion oscil- ear Schr¨odingerequation, Physica D: Nonlinear Phenom- lating fibers, Advances in Optics and Photonics 10, 1 ena 145, 65 (2000). (2018). [29] M. Chertkov, I. Gabitov, and J. Moeser, Pulse confine- [14] M. Droques, A. Kudlinski, G. Bouwmans, G. Martinelli, ment in optical fibers with random dispersion, Proceed- A. Mussot, A. Armaroli, and F. Biancalana, Fourth-order ings of the National Academy of Sciences 98, 14208 dispersion mediated modulation instability in dispersion (2001). oscillating fibers, Optics Letters 38, 3464 (2013). [30] V. N. Kampen, Stochastic Processes in Physics and [15] A. Armaroli and F. Biancalana, Suppression and split- Chemistry (Elsevier, 2007). ting of modulational instability sidebands in periodically [31] M. Karlsson, Four-wave mixing in fibers with randomly tapered optical fibers because of fourth-order dispersion, varying zero-dispersion wavelength, Journal of the Opti- Optics Letters 39, 4804 (2014). cal Society of America B 15, 2269 (1998). [16] A. Armaroli and F. Biancalana, Vector modulational in- [32] M. Farahmand and M. de Sterke, Parametric amplifica- stability induced by parametric resonance in periodically tion in presence of dispersion fluctuations, Optics Express tapered highly birefringent optical fibers, Physical Re- 12, 136 (2004). view A 87, 063848 (2013). [33] S. Rota Nodari, M. Conforti, G. Dujardin, A. Kudlinski, [17] M. Conforti, A. Mussot, A. Kudlinski, and S. Trillo, Mod- A. Mussot, S. Trillo, and S. De Bi`evre,Modulational in- ulational instability in dispersion oscillating fiber ring stability in dispersion-kicked optical fibers, Phys. Rev. A cavities, Optics letters 39, 4200 (2014), arXiv:1405.6041. 92 (2015). [18] M. Conforti, F. Copie, A. Mussot, A. Kudlinski, and [34] B. V. Chirikov, A universal instability of many- 13

dimensional oscillator systems, Physics Reports 52, 263 plicative noise, Physica A: Statistical Mechanics and its (1979). Applications 352, 309 (2005). [35] C. Tian, A. Altland, and M. Garst, Theory of the an- [42] K. Mallick and P.-E. Peyneau, Phase diagram of derson transition in the quasiperiodic kicked rotor, Phys. the random frequency oscillator: The case of Orn- Rev. Lett. 107, 074101 (2011). stein–Uhlenbeck noise, Physica D 221, 72 (2006). [36] G. Casati, B. V. Chirikov, F. M. Izraelev, and J. Ford, [43] F. J. Poulin and G. R. Flierl, The stochastic Mathieu’s Stochastic behavior of a quantum pendulum under a peri- equation, Proceedings of the Royal Society A: Mathemat- odic perturbation, in Stochastic behavior in classical and ical, Physical and Engineering Sciences 464, 1885 (2008). quantum Hamiltonian systems (Volta Memorial Conf., [44] A. Armaroli, M. Conforti, G. Dujardin, A. Kudlinski, Como, 1977), Lecture Notes in Phys., Vol. 93 (Springer, A. Mussot, S. Trillo, and S. De Bi`evre,From parametric Berlin-New York, 1979) pp. 334–352. to stochastic resonances in nonlinear Schr¨odingerequa- [37] P. W. Anderson, Absence of Diffusion in Certain Random tion with random dispersion, In preparation. Lattices, Physical Review 109, 1492 (1958). [45] H. Furstenberg, Noncommuting random products, Trans. [38] P. Bougerol and J. Lacroix, Products of random matri- Amer. Math. Soc. 108, 377 (1963). ces with applications to Schr¨odingeroperators, Progress [46] L. Landau and E. Lifshitz, Mechanics, third edition ed. in Probability and Statistics, Vol. 8 (Birkh¨auserBoston, (Butterworth-Heinemann, Oxford, 1976). Inc., Boston, MA, 1985) pp. xii+283. [47] V. I. Arnol’d, Mathematical Methods of Classical Me- [39] N. G. Van Kampen, Stochastic differential equations, chanics, 2nd ed. (Springer, New York Heidelberg Dor- Physics Reports 24, 171 (1976). drecht London, 2013). [40] R. V. Bobryk and A. Chrzeszczyk, Colored-noise-induced [48] A. H. Nayfeh and D. T. Mook, Nonlinear oscillations parametric resonance, Physica A: Statistical Mechanics (Wiley, 1979). and its Applications 316, 225 (2002). [41] M. Gitterman, Classical harmonic oscillator with multi-