## Passive mode-locking under higher order effects

Theodoros P. Horikis 1 , ∗ and Mark J. Ablowitz 2

1 Department of Mathematics, University of Ioannina, Ioannina 45110, Greece

2 Department of Applied Mathematics, University of Colorado, 526 UCB, Boulder, CO 80309-0526, USA

compiled: August 4, 2014

The response of a passive mode-locking mechanism, where gain and spectral filtering are saturated with the energy and loss saturated with the power, is examined under the presence of higher order effects. These include third order dispersion, self-steepening and Raman gain. The locking mechanism is maintained even with these terms; mode-locking occurs for both the anomalous and normal regimes. In the anomalous regime, these perturbations are found to affect the speed but not the structure of the (locked) pulses. In fact, these pulses behave like solitons of a classical nonlinear Schr¨dinger equation and as such a soliton perturbation theory is o used to verify the numerical observations. In the normal regime, the effect of the perturbations is small, in line with recent experimental observations. The results in the normal regime are verified mathematically using a WKB type asymptotic theory. Finally, bi-solitons are found to behave as dark solitons on top of a stable background and are significantly affected by these perturbations.

OCIS codes: 140.4050 (Mode-locked lasers), 190.5530 (Pulse propagation and temporal solitons),

320.5540 (Pulse shaping)

http://dx.doi.org/10.1364/XX.99.099999

## 1. Introduction

Ultrashort pulses are used frequently in modern applications ranging from femtochemistry and medical imaging to micro-machining and optical communications. Consequently the study and modeling of such pulses is an important area of research. Kerr-lens mode locking (KLM) is a common technique used to generate these pulses. In fact, following the discovery of KLM mode-locked (ML) lasers have revolutionized the field of ultrafast science and their performance has led to their widespread use.

Mode-locked lasers operate under the requirement that sufficient amounts of gain and loss are present, otherwise, they are subject to radiation-mode instability. In addition, passive mode-locking generally utilizes fast saturable absorbers which are effectively simulated by KLM. Saturable absorbers are commonly modeled by a transfer function and as such loss is placed periodically during the evolution and not as part of a distributed model, i.e., an evolution equation. However, the qualitative features are the same in both cases [1] thus indicating that distributive models are very good descriptions of modes in mode-locked lasers.

Modeling these lasers must take into account their physical properties. For this, the distributive powerenergy saturation (PES) model was recently introduced [2]. This system is a variant of the well known nonlinear Schr¨dinger (NLS) equation with additional terms o to account for gain and spectral filtering, saturated with

∗ Corresponding author: horikis@uoi.gr

energy, and loss, saturated with power. It is a generalization of the commonly used master-equation [3, 4] which is obtained as a limiting case. Mode-locking is much more robust in the PES model than in the masterequation. Indeed, while in the master-equation modelocking is achieved only for a narrow range of parameters [5, 6], the PES equation has only one key requirement for mode-locking: sufficient gain. Another important feature of the PES equation is that it leads to physically relevant results for both the anomalous and normal dispersive regimes.

Solitons in ML lasers are fundamentally different in the two dispersive regimes. Pulse formation, in the anomalous regime, is typically dominated by the interplay between dispersion and nonlinearity. Suitable gain media and an effective saturable absorber are required for initiation of pulsed operation. In contrast, pulses found in the normal regime are positively chirped throughout the cavity [7]. They are comprised of an approximately parabolic temporal amplitude profile near the peak of the pulse with a transition to steep decay [7]. Pulses in the normal regime have a much larger temporal width than those in the anomalous regime.

Here we study the response of a ML laser under higher order effects using the PES equation. These effects include third-order dispersion (TOD), self-steepening and Raman gain. Previous studies, in the anomalous dispersion regime, include variants of the master-equation [8, 9] and generalized Ginzburg-Landau type equations [10-12]. In the normal regime recent experiments [13] indicate that similariton (self-similar) evolution in these fiber lasers results in localized modes, i.e. solitons, which

are broad in the time domain; they are robust against perturbations. The effect of TOD on self-similar evolution and parabolic pulses was studied in Refs. [14-16].

Here we focus on the equation

$$i \frac { \partial u } { \partial z } + \frac { d _ { 0 } } { 2 } \frac { \partial ^ { 2 } u } { \partial t ^ { 2 } } + | u | ^ { 2 } u & = \frac { i g u } { 1 + \epsilon E } + \frac { i \tau } { 1 + \epsilon E } \frac { \partial ^ { 2 } u } { \partial t ^ { 2 } } - \frac { i l u } { 1 + \delta P } \\ & + i \beta \frac { \partial ^ { 3 } u } { \partial t ^ { 3 } } - i \gamma \frac { \partial ( | u | ^ { 2 } u ) } { \partial t } + R u \frac { \partial ( | u | ^ { 2 } ) } { \partial t }$$

where g represents gain saturated with the energy E = ∫ ∞ -∞ | u | 2 dt , τ is spectral filtering, l is the gain saturated with power P = | u | 2 , β represents TOD, γ selfsteepening and R the Raman gain. The parameters /epsilon1 &gt; 0 and δ &gt; 0 correspond to a measure of the saturation energy and power of the system, respectively. Unless otherwise stated all parameters, except d 0 , are positive. When β = γ = R = 0 this equation is the PES model. Finally, other linear and higher order nonlinear dissipative terms, such as terms proportional to iu and i u | | 2 u are included in Eq. (1) as limiting cases in the expansion for the saturable terms and will not be considered here.

## 2. The anomalous dispersion regime

Here we consider Eq. (1) with d 0 &gt; 1 and we start our analysis by analyzing the effect of TOD. To do so, we fix the parameters such that τ = l = 0 1, . d 0 = /epsilon1 = δ = 1 and β = 0 01 while . γ = R = 0 and vary the gain parameter g . A unit gaussian is evolved under Eq. (1) with the resulting evolution shown in Fig. 1.

As seen here, the initial gaussian undergoes a locking procedure and subsequently locks to a pulse of constant height/shape. After that happens the effect of TOD shifts the soliton with constant velocity (seen as a straight line for the soliton center in the top of Fig. 1). Interestingly the resulting mode is a solution of the unperturbed NLS equation. Furthermore, when the gain is increased the angle of the displacement becomes more acute. This is reminiscent of the response of a single soliton of the classical NLS equation to TOD. Indeed, recall that under TOD the NLS soliton is slowed down and as a result the soliton peak is shifted by an amount that increases linearly with distance [17, 18].

Similarly, we repeat the analysis for β = R = 0 and γ = 0 01 so as to study the response to self-steepening. . The resulting evolution is depicted in Fig. 2.

The similarities between Figs. 1 and 2 are remarkable. Indeed, after mode-locking occurs the pulse is shifted linearly in a way that is similar to the way a classical NLS soliton propagates under this effect.

Finally, we repeat the analysis for β = γ = 0 and R = 0 01 so as to . study the response to the Raman gain. The resulting evolution is depicted in Fig. 3.

The phenomenon is repeated only now the displacement is sharper. Again, recall that Raman gain on solitons of the classical NLS equation gives rise to the so-

Fig. 1. (Color online) Mode-locking under TOD with (a) g = 0 2 and (b) . g = 0 3. . The bottom figures show the pulse at a distance z = 150 (a sufficient distance of propagation after mode-locking has occurred). The dashed lines correspond to the relative classical NLS solitons.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0041v1-figure-86633020-e5bf-4ab5-b7dc-3dccae1124fe.jpg)

Fig. 2. (Color online) Mode-locking under self-steepening with (a) g = 0 2 and (b) . g = 0 3. . Profiles of mode-locked pulses resemble those of Fig. 1 (bottom).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0041v1-figure-88bc5d0b-11e2-4a3b-bcdf-d57ec2381c05.jpg)

called self-frequency shift [17, 18]. Note the size of the pulse amplitude is largely unaffected.

It is also important to stress that this locking behavior strongly depends on the saturation terms. If these terms/effects are absent more complicated dynamics can occur [12].

## 2.A. Perturbation theory

Now we shall explain the above features using soliton perturbation theory valid for small β , γ , R . This ap-

Fig. 3. (Color online) Mode-locking under Raman gain with (a) g = 0 2 and (b) . g = 0 3. . Profiles of mode-locked pulses resemble those of Fig. 1 (bottom).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0041v1-figure-7e8889ea-ed5d-413b-94c8-a92b7fa68fdc.jpg)

proach leads to very good qualitative results since when β = γ = R = 0 general initial conditions, such as unit gaussians, will evolve to solitons of the unperturbed classical NLS equation [19]. To do this, we consider a solution ansatz for Eq. (1) of the form

$$u ( z, t ) = \eta ( z ) \text{sech} [ \eta ( z ) ( t - t _ { 0 } ( z ) ) ] e ^ { i \phi ( z ) - i \sigma ( z ) t }$$

Under the effect of the perturbed PES the soliton parameters evolve according to [17, 18, 20]

$$\frac { d \eta } { d z } & = \frac { 2 g \eta } { 1 + 2 \epsilon \eta } - \frac { 2 \tau ( \eta ^ { 2 } + 3 \sigma ^ { 2 } ) } { 3 ( 1 + 2 \epsilon \eta ) } \eta - \frac { 2 l \tanh ^ { - 1 } \left ( \frac { \sqrt { \delta } \eta } { \sqrt { 1 + \delta \eta ^ { 2 } } } \right ) } { \sqrt { \delta + \delta ^ { 2 } \eta ^ { 2 } } } \\ d \sigma \, & = \, 4 \tau \eta ^ { 2 } \sigma \, \quad \, 8 \, _ { \ D \dots } 4 \, \quad \, \text{fig. 4} }$$

$$\frac { d \sigma } { d z } = - \frac { 4 \tau \eta ^ { 2 } \sigma } { 3 ( 1 + 2 \epsilon \eta ) } - \frac { 8 } { 1 5 } R \eta ^ { 4 } & ( 2 b ) \quad \begin{array} { c } \text{Fig.} \\ \text{$value} \\ l = t \end{array}$$

$$\frac { d t _ { 0 } } { d z } = - \sigma + \beta ( \eta ^ { 2 } + 3 \sigma ^ { 2 } ) + \gamma \eta ^ { 2 } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 0 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdot \cdots } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \sum _ { \substack { 1 \\ 1 } } ^ { \cdot \cdots } \\ \end{pmatrix}$$

$$\stackrel { \ddots } { \frac { d \phi } { d z } } = \beta \sigma ( 3 \eta ^ { 2 } + \sigma ^ { 2 } ) + \gamma \eta ^ { 2 } + \frac { 1 } { 2 } ( \eta ^ { 2 } - \sigma ^ { 2 } ) + t _ { 0 } \frac { d \sigma } { d z } \ \ ( 2 d ) \quad \ \ \sigma \ a _ { \sigma \ a _ { \dots } }$$

When the pulse has mode-locked, the zeroes of the right hand side of Eq. (2a)-(2b) give the stationary solutions for η and σ . We denote the mode locked value of η and σ as η ML and σ ML , respectively. As such when the equilibrium has been reached

$$\sigma \equiv \sigma _ { M L } = - \frac { 2 R ( 1 + 2 \epsilon \eta ) \eta ^ { 2 } } { 5 \tau }$$

this results in the following equation for η

$$\max \text{rows in one moving equation } \eta \\ \frac { d \eta } { d z } = \frac { 2 g \eta } { 1 + 2 \epsilon \eta } - \frac { 2 \left [ 2 5 \tau ^ { 2 } \eta ^ { 2 } + 1 2 ( 2 \epsilon R \eta + R ) ^ { 2 } \eta ^ { 4 } \right ] } { 7 5 ( \tau + 2 \epsilon \tau \eta ) } \\ - \frac { 2 l \tanh ^ { - 1 } \left ( \frac { \sqrt { \delta } \eta } { \sqrt { 1 + \delta \eta ^ { 2 } } } \right ) } { \sqrt { \delta + \delta ^ { 2 } \eta ^ { 2 } } }$$

Clearly, η has at least one equilibrium at η = 0. To understand its stability consider the first derivative of η z with respect to η evaluated at η = 0,

$$\frac { d } { d \eta } \left ( \frac { d \eta } { d z } \right ) \Big | _ { \eta = 0 } = 2 ( g - l ) \\ \intertext { I, \eta = 0 \text{ is an unstable equilibrium} } \cdot \sum \dots \sum \dots$$

When g &gt; l, η = 0 is an unstable equilibrium while for g &lt; l, η = 0 is stable. The phase portrait for various values of g &lt; l indicates that η z is a monotonic decreasing function in η and so η → 0 as z →∞ for any initial condition η &gt; 0. Physically this corresponds to loss overtaking gain and the pulse decays. For g &gt; l , η z is initially positive and as η increases the filtering term becomes dominant and η z takes on negative values. Thus, there must exist another equilibrium. Plotting the phase portrait for various values of g for g &gt; l shows there is a single stable equilibrium for η &gt; 0, as shown in Fig. 4.

η

Fig. 4. Phase portrait of the amplitude equation for several values of g illustrating the single stable equilibrium. Here l = 0 1. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0041v1-figure-f5385219-c88a-4e40-942f-83b26c67b291.jpg)

/negationslash

As η approaches its stable equilibrium, η ML = 0, and σ approaches its stable equilibrium, σ ML , the system of equations (2) tend to

$$\frac { d \eta } { d z } \Big | _ { \eta = \eta _ { M L } } = \frac { d \sigma } { d z } \Big | _ { \sigma = \sigma _ { M L } } = 0 \\ \text{ans the solution tends to a particular } \Big | _ { \quad }.$$

which means the solition tends to a particular NLS soliton of speed zero and height/width determined by η ML , thus suggesting the mode-locking capabilities of Eq. (1). In fact, η ML is an attractor that any arbitrary initial condition will eventually converge to, even when it is far way from its soliton solution [2, 19].

After η = η ML and mode-locking has been completed, we consider Eq. (2c)-(2d), which, interestingly, are independent of g, τ, l . This now explains the above dynamics. Indeed, when η = η ML (and the mode-locking mechanism has been completed) we get by direct integration

of Eq. (2c)

$$t _ { 0 } ( z ) & = t _ { 0 } ( 0 ) + \left \{ \left ( \beta + \gamma + \frac { 2 R } { 5 \tau } \right ) \eta _ { M L } ^ { 2 } \quad \begin{array} { c } \text{$\text{$\eta_{ML}$} \\ \text{$\text{has}$} \\ \text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$
$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$                                                                                                                                                                                                         \text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{	$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text	$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\test$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text
$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text`$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$	\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\tt$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$}\end{array} \right \}$$

The equations for the equilibria give now more insight as to the way each term affects the propagation. For instance, in the absence of spectral filtering there still exists a stable (nonzero) equilibrium for η but not for σ . While the amplitude is affected only in its magnitude the other quantities change in a more prominent way. Indeed, by direct integration of Eqs. (2b) and (2c) we obtain at η = η ML

$$\sigma ( z ) & = - \frac { 8 } { 1 5 } R \eta _ { M L } ^ { 4 } z + \sigma ( 0 ) \\ t _ { 0 } ( z ) & = \frac { 6 4 } { 2 2 5 } R ^ { 2 } \beta \eta _ { M L } ^ { 8 } z ^ { 3 } - \frac { 4 } { 1 5 } R \eta _ { M L } ^ { 4 } [ 6 \beta \sigma ( 0 ) - 1 ] z ^ { 2 } \\ & + [ ( \beta + \gamma ) \eta _ { M L } ^ { 2 } + 3 \beta \sigma ^ { 2 } ( 0 ) - 3 \sigma ( 0 ) ] z + t _ { 0 } ( 0 ) \\ \text{which suggests that the trajectory the soliton follows}$$

which suggests that the trajectory the soliton follows is no longer linear. Similarly when the Raman gain is not considered the only equilibrium is found for σ ≡ 0. Clearly when τ = R = 0 the soliton moves with constant velocity equal to ( β + γ η ) 2 ML ( σ ≡ 0). Also, note that as supported by the numerical computations the final amplitude of the pulse is not affected by the parameters β, γ .

The increase of η with g also explains the pronounced difference in the speed of the pulse when g is increased. In Fig. 5 we depict the dependance of the amplitude to the gain parameter, which attests to that. In fact, Eq.

η

Fig. 5. (Color online) The amplitude dependance on gain based on Eq. (2a).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0041v1-figure-6b4d51b3-9946-4595-89ad-8489bbdd6ab0.jpg)

(3) clearly suggests that the trajectories are highly dependent on powers of η ML . Hence, any increase on this amplitude is enhanced (when η ML &gt; 1 or g &gt; 0 22) mak-. ing the response to the perturbations more pronounced. This analysis also suggests that nonlinear phase shifts cannot compensate for TOD (or self-steepening) in this regime, unlike in the normal regime [21].

Thus, the mode-locking mechanism can now be summarized as follows: (i) solve Eqs. (2a)-(2b), to determine η ML and σ ML ; (ii) assume the mode-locking mechanism has been completed; (iii) treat the pulses as classical NLS solitons; and (iv) use Eq. (3) and (2d) to determine the adiabatic change of the rest of the soliton parameters.

## 3. The normal dispersion regime

We turn our attention to the normal regime ( d 0 = -1). The main difference these pulses exhibit from their anomalous regime counterparts is that they are highly chirped, very wide (almost parabolic) pulses in the time domain; the modes depend critically on the values of g, τ, l .

Recent experiments indicate that localized modessolitons in the normal regime are robust. They are not affected by perturbations [13] such as TOD and self steepening; this is also in agreement with the observation that nonlinear phase shifts (in terms of Raman gain here) can compensate for TOD or self-steepening [21]. To demonstrate this phenomena, we evolve a unit gaussian under Eq. (1), with all parameters as before ( τ = l = 0 1, . β = γ = R = 0 01) only now with . g = 1 5. . The reason for this increase in gain is the nature/size of these pulses (cf. Fig. 6). Since these pulses are much wider one expects that more energy (gain) is required for pulses to lock to noticeable amplitudes. Furthermore, in Fig. 6 (bottom) we depict the resulting pulse and its corresponding phase. The dashed line indicates the soliton of the unperturbed (i.e. β = γ = R = 0) PES equation. Indeed the experimental observations [7] on the shape and phase of the pulses are met. We also confirm subsequent observations that unlike the solitons of the anomalous regime these wide localized pulses are resistant to perturbations even when allowed to propagate for greater distances.

## 3.A. Asymptotic theory

Following the theory developed in Ref. [22] we write the solution of Eq. (1) in the form ψ z, t ( ) = A exp( iµz + iθ ), where A = A t ( ) and θ = θ t ( ) are the pulse amplitude and phase respectively. Since these pulses are slowly varying in t , with large phase, we can introduce a slow time scale in the equation and using perturbation theory, the soliton system can be reduced to simpler ordinary differential equations for the amplitude and the phase of the pulse. Substituting into the equation and equating

Fig. 6. (Color online) Mode-locking in the normal regime (a) under the unperturbed PES equation with β = γ = R = 0 and (b) under the PES with the higher order effects included. The bottom figure depicts the final pulse. The dashed line is the soliton of the unperturbed PES equation.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0041v1-figure-4bb22dfb-591b-4d27-ad75-3329879cf708.jpg)

real and imaginary parts we get

$$& - \mu A - \frac { d _ { 0 } } { 2 } ( A _ { t t } - A \theta _ { t } ^ { 2 } ) + A ^ { 3 } = - \frac { \tau } { 1 + \epsilon E } ( 2 A _ { t } \theta _ { t } + A \theta _ { t t } ) \\ & + \beta ( A \theta _ { t t } - 3 A _ { t t } \theta _ { t } - 3 A _ { t } \theta _ { t t } - A \theta _ { t t t } ) + \gamma A ^ { 3 } \theta _ { t } + 2 R A ^ { 2 } A _ { t }$$

$$+ & \beta ( A \theta _ { t t } - 3 A _ { t t } \theta _ { t } - 3 A _ { t } \theta _ { t t t } ) + \gamma A ^ { 3 } \theta _ { t } + 2 R A ^ { 2 } A _ { t } \\ & - \frac { d _ { 0 } } { 2 } ( 2 A _ { t } \theta _ { t } + A \theta _ { t t } ) = \frac { g } { 1 + \epsilon E } A + \frac { \tau } { 1 + \epsilon E } ( A _ { t t } + A \theta _ { t } ^ { 2 } ) \\ & - \frac { l } { 1 + \delta A ^ { 2 } } A + \beta ( 3 A _ { t } \theta _ { t } ^ { 2 } - 3 A \theta _ { t } \theta _ { t t } + A _ { t t t } ) - 3 \gamma A ^ { 2 } A _ { t }$$

where we have replaced d 0 →-d 0 &lt; 0 since we are in the normal dispersion regime. We then take the characteristic time length of the pulse to be such that we can define a scaling in the independent variable of the form ε = T/ε or T = εt so that R = R εt ( ), θ t = O (1) (i.e. θ is large) and θ tt = O ε ( ) where ε /lessmuch 1. Furthermore in the above simulations β , γ and R are all O ε ( ). Then Eq. (4a) becomes

$$- \mu A + \frac { d _ { 0 } } { 2 } A \theta _ { t } ^ { 2 } + A ^ { 3 } = \frac { d _ { 0 } } { 2 } \varepsilon ^ { 2 } A _ { T T } - \frac { \tau } { 1 + \epsilon E } ( 2 \varepsilon A _ { T } \theta _ { t } + A \theta _ { t t } ) \quad \text{lutio} \quad \text{--aged}$$

while the leading order equation is

$$\theta _ { t } ^ { 2 } = \frac { 2 } { d _ { 0 } } ( \mu - A ^ { 2 } )$$

Using the same argument and this newly derived Eq. (5) we then find that Eq. (4b) to leading order reads

$$A _ { t } & = - \text{sgn} ( t ) \frac { \sqrt { 2 ( \mu - A ^ { 2 } ) / d _ { 0 } } } { 3 A ^ { 2 } - 2 \mu } \\ & \quad \times \left ( \frac { g } { 1 + \epsilon E } - \frac { 4 \tau ( \mu - A ^ { 2 } ) / d _ { 0 } } { 1 + \epsilon E } - \frac { l } { 1 + \delta A ^ { 2 } } \right ) \quad ( 6 ) \\ \text{This is now a nonlocal first order differential equation}$$

This is now a nonlocal first order differential equation for A = A t ( ), since E = ∫ ∞ -∞ A 2 dt . We also note that imposing θ t ( t = 0) = 0 ⇒ µ ≈ A 2 (0).

To remove the nonlocality another condition is needed and it is based on the singular points of Eq. (6). Recall, A t ( ) is a decaying function in t ∈ [0 , + ∞ ) and µ ≈ A 2 (0) ≥ A t ( ) 2 . Thus there exists a point in t such that the denominator in the equation becomes zero. To remove the singularity we require that the numerator of the equation is also zero at the same point which leads to

$$1 + \epsilon E = \frac { 1 } { l } \left ( g - \frac { 2 \tau } { 3 d _ { 0 } } \mu \right ) \left ( 1 + \frac { 2 \delta } { 3 } \mu \right ) \\ \psi = \left ( \lambda \right ) \cdot \left ( \lambda \right ) \cdot \right ) \cdot \right ) \cdot \right )$$

Thus Eq. (6) is now a first order equation which can be solved by standard numerical methods and analyzed by phase plane methods. The resulting solutions from the PES equation and the reduced equations are compared in Fig. 7.

Fig. 7. (Color online) Solutions of the complete PES equation and the reduced system.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0041v1-figure-14092df1-db5e-4ad7-a06d-7f605f54ec70.jpg)

## 3.B. Higher order states

Remarkably, in the normal regime the PES also exhibits higher-order solutions in term of bi-solitons [23]. These are pairs of regular solitons whose peak amplitudes have a π difference in phase and an appropriate separation. They, also, differ significantly from the higher-order solutions of the classical NLS and the dispersion managed solitons since they do not exhibit any oscillatory

or breathing behavior as they propagate in the cavity. A typical such mode is given in Fig. 8. The dashed line corresponds to the relative 'ground state', i.e. the localized pulse of the unperturbed PES equation.

Fig. 8. (Color online) A typical bisoliton of the PES. Dashed lines correspond to the soliton and phase of the localized pulse of the unperturbed PES equation.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0041v1-figure-05368f53-fb97-4d53-b1c0-fbaaa9a3f057.jpg)

Higher order states and interaction of these structures have been discussed in Ref. [24]. It is, however, interesting to see their response to TOD and Raman gain. In order to generate these modes we use an antisymmetric initial condition such as u (0 , t ) = t exp( -t 2 ). The resulting evolution, under Eq. (1), is shown in Fig. 9.

Fig. 9. (Color online) Bisoliton evolution (a) under the PES equation; (b) under Eq. (1).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0041v1-figure-fa4f1387-950a-4081-9c21-4e7089b7f68c.jpg)

Under these higher order effects the dip characterizing the bi-soliton, which separates the two solitons, behaves as a dark soliton on a nearly constant background and starts moving until it disappears; see Fig. 9. This is reminiscent of the dark solitons generated in Bose-Einstein condensates [25]. We will further explore this elsewhere.

## 4. Conclusions

To conclude, we have studied the response of a PES model with additional higher order terms: third order dispersion, self-steepening and Raman gain. In the anomalous dispersion regime, it is found that the PES

equation provides the mode-locking mechanism and after that the resulting pulses propagate in a manner similar to solitons of a perturbed classical NLS equation. Perturbative soliton theory confirms numerical observations. On the other hand, in the normal dispersion regime, ground state pulses are more robust and are largely unaffected by these perturbing influences as also observed in experiments. However, this is not the case with higher order modes which are unstable when perturbed.

## Acknowledgments

This research was partially supported by the United States Air Force Office of Scientific Research (USAFOSR) under grant FA9550-12-1-0207.

## References

- [1] M. J. Ablowitz and T. P. Horikis, 'Solitons and spectral renormalization methods in nonlinear optics,' Eur. Phys. J. Special Topics 173 , 147-166 (2009).
- [2] M. J. Ablowitz, T. P. Horikis, and B. Ilan, 'Solitons in dispersion-managed mode-locked lasers,' Phys. Rev. A 77 , 033814 (2008).
- [3] H. A. Haus, 'Theory of mode locking with a fast saturable absorber,' J. Appl. Plys. 46 , 3049-3058 (1975).
- [4] H. A. Haus, J. G. Fujimoto, and E. P. Ippen, 'Analytic theory of additive pulse and kerr lens mode locking,' IEEE J. Quant. Elec. 28 , 2086-2096 (1992).
- [5] J. N. Kutz, 'Mode-locked soliton lasers,' SIAM Rev. 48 , 629-678 (2006).
- [6] T. Kapitula, J. N. Kutz, and B. Sandstede, 'Stability of pulses in the master mode-locking equation,' J. Opt. Soc. Am. B 19 , 740746 (2002).
- [7] F. O. Ilday, J. R. Buckley, W. G. Clark, and F. W. Wise, 'Self-similar evolution of parabolic pulses in a laser,' Phys. Rev. Lett. 92 , 213901 (2004).
- [8] H. A. Haus, J. D. Moores, and L. E. Nelson, 'Effect of third-order dispersion on passive mode locking,' Opt. Lett. 18 , 51-53 (1993).
- [9] Q. Lin and I. Sorokina, 'High-order dispersion effects in solitary mode-locked lasers: side-band generation,' Opt. Commun. 153 , 285-288 (1998).
- [10] V. L. Kalashnikov, A. Fern´ndez, a and A. Apolonski, 'High-order dispersion in chirped-pulse oscillators,' Opt. Express 16 , 4206-4216 (2008).
- [11] L. Song, L. Li, Z. Li, and G. Zhou, 'Effect of third-order dispersion on pulsating, erupting and creeping solitons,' Opt. Commun. 249 , 301-309 (2005).
- [12] E. Sorokin, N. Tolstik, V. L. Kalashnikov, and I. T. Sorokina, 'Chaotic chirped-pulse oscillators,' Opt. Express 21 , 29567-29577 (2013).
- [13] B. Oktem, C. ¨ lg¨d¨r, U u u and F. O. Ilday, 'Solitonsimilariton fibre laser,' Nature Photonics 4 , 307-311 (2010).
- [14] J. M. Dudley, C. Finot, D. J. Richardson, and G. Millot, 'Self-similarity in ultrafast nonlinear optics,' Nature Physics 3 , 597-603 (2007).
- [15] B. G. Bale and S. Boscolo, 'Impact of third-order fibre dispersion on the evolution of parabolic optical pulses,' J. Opt. 12 , 015202 (2010).
- [16] S. Zhang, G. Zhao, A. Luo, and Z. Zhang, 'Thirdorder dispersion role on parabolic pulse propagation in

- dispersion-decreasing fiber with normal group-velocity dispersion,' Appl. Phys. B 94 , 227-232 (2009).
- [17] A. Hasegawa and Y. Kodama, Solitons in Optical Communications (Claredo Press, 1995).
- [18] Y. S. Kivshar and G. P. Agrawal, Optical Solitons (Academic Press, 2003).
- [19] M. J. Ablowitz and T. P. Horikis, 'Pulse dynamics and solitons in mode-locked lasers,' Phys. Rev. A 78 , 011802(R) (2008).
- [20] M. J. Ablowitz, T. P. Horikis, S. D. Nixon, and Y. Zhu, 'Asymptotic analysis of pulse dynamics in mode-locked lasers,' Stud. Appl. Math 122 , 411-425 (2009).
- [21] S. Zhou, L. Kuznetsova, A. Chong, and F. W. Wise, 'Compensation of nonlinear phase shifts with third-
- order dispersion in short-pulse fiber amplifiers,' Opt. Express 13 , 4869-4877 (2005).
- [22] M. J. Ablowitz and T. P. Horikis, 'Solitons in normally dispersive mode-locked lasers,' Phys. Rev. A 79 , 063845 (2009).
- [23] A. Chong, W. H. Renninger, and F. W. Wise, 'Observation of antisymmetric dispersion-managed solitons in a mode-locked laser,' Opt. Lett. 33 , 1717-1719 (2008).
- [24] M. J. Ablowitz, T. P. Horikis, and S. D. Nixon, 'Soliton strings and interactions in mode-locked lasers,' Opt. Comm. 282 , 4127-4135 (2009).
- [25] D. J. Frantzeskakis, 'Dark solitons in atomic boseeinstein condensates: from theory to experiments,' J. Phys. A: Math. Theor. 43 , 213001 (2010).
Theodoros P. Horikis Mark J. Ablowitz
Passive mode-locking under higher order effects
The response of a passive mode-locking mechanism, where gain and spectral
filtering are saturated with the energy and loss saturated with the power, is
examined under the presence of higher order effects. These include third order
dispersion, self-steepening and Raman gain. The locking mechanism is maintained
even with these terms; mode-locking occurs for both the anomalous and normal
regimes. In the anomalous regime, these perturbations are found to affect the
speed but not the structure of the (locked) pulses. In fact, these pulses
behave like solitons of a classical nonlinear Schrodinger equation and as such
a soliton perturbation theory is used to verify the numerical observations. In
the normal regime, the effect of the perturbations is small, in line with
recent experimental observations. The results in the normal regime are verified
mathematically using a WKB type asymptotic theory. Finally, bi-solitons are
found to behave as dark solitons on top of a stable background and are
significantly affected by these perturbations.