Dissipative nonlinear structures in optical fiber systems III:dissipative solitons via nonlinear mapping

In the third and final talk on dissipative structures in fiber applications, we discuss mathematical techniques that can be used to characterize modern laser systems that consist of several discrete elements. In particular, we use a nonlinear mapping technique to evaluate high power laser systems where significant changes in the pulse evolution per cavity round trip is observed. We demonstrate that dissipative soliton solutions might be effectively described using this Poincaré mapping approach.

A perturbative analysis of dispersion-managed solitons

We study soliton solutions of the path-averaged propagation equation governing the transmission of dispersion-managed (DM) optical pulses in the (practical) limit when residual dispersion and nonlinearity only slightly affect the pulse dynamics over one compensation period. In the case of small dispersion map strengths, the averaged pulse dynamics is governed by a perturbed form of the nonlinear Schrödinger equation; applying a perturbation theory – elsewhere developed – based on inverse scattering theory, we derive an analytic expression for the envelope of the DM soliton. This expression correctly predicts the power enhancement arising from the dispersion management. Theoretical results are verified by direct numerical simulations.

One-dimensional optical wave turbulence:experiment and theory

We present a review of the latest developments in one-dimensional (1D) optical wave turbulence (OWT). Based on an original experimental setup that allows for the implementation of 1D OWT, we are able to show that an inverse cascade occurs through the spontaneous evolution of the nonlinear field up to the point when modulational instability leads to soliton formation. After solitons are formed, further interaction of the solitons among themselves and with incoherent waves leads to a final condensate state dominated by a single strong soliton. Motivated by the observations, we develop a theoretical description, showing that the inverse cascade develops through six-wave interaction, and that this is the basic mechanism of nonlinear wave coupling for 1D OWT. We describe theory, numerics and experimental observations while trying to incorporate all the different aspects into a consistent context. The experimental system is described by two coupled nonlinear equations, which we explore within two wave limits allowing for the expression of the evolution of the complex amplitude in a single dynamical equation. The long-wave limit corresponds to waves with wave numbers smaller than the electrical coherence length of the liquid crystal, and the opposite limit, when wave numbers are larger. We show that both of these systems are of a dual cascade type, analogous to two-dimensional (2D) turbulence, which can be described by wave turbulence (WT) theory, and conclude that the cascades are induced by a six-wave resonant interaction process. WT theory predicts several stationary solutions (non-equilibrium and thermodynamic) to both the long- and short-wave systems, and we investigate the necessary conditions required for their realization. Interestingly, the long-wave system is close to the integrable 1D nonlinear Schrödinger equation (NLSE) (which contains exact nonlinear soliton solutions), and as a result during the inverse cascade, nonlinearity of the system at low wave numbers becomes strong. Subsequently, due to the focusing nature of the nonlinearity, this leads to modulational instability (MI) of the condensate and the formation of solitons. Finally, with the aid of the probability density function (PDF) description of WT theory, we explain the coexistence and mutual interactions between solitons and the weakly nonlinear random wave background in the form of a wave turbulence life cycle (WTLC).

Ion-acoustic envelope modes in a degenerate relativistic electron-ion plasma

A self-consistent relativistic two-fluid model is proposed for one-dimensional electron-ion plasma dynamics. A multiple scales perturbation technique is employed, leading to an evolution equation for the wave envelope, in the form of a nonlinear Schrödinger type equation (NLSE). The inclusion of relativistic effects is shown to introduce density-dependent factors, not present in the non-relativistic case - in the conditions for modulational instability. The role of relativistic effects on the linear dispersion laws and on envelope soliton solutions of the NLSE is discussed.

Discrete breathers in a two-dimensional hexagonal Fermi-Pasta-Ulam lattice

We consider a two-dimensional Fermi-Pasta-Ulam (FPU) lattice with hexagonal symmetry. Using asymptotic methods based on small amplitude ansatz, at third order we obtain a eduction to a cubic nonlinear Schr{\"o}dinger equation (NLS) for the breather envelope. However, this does not support stable soliton solutions, so we pursue a higher-order analysis yielding a generalised NLS, which includes known stabilising terms. We present numerical results which suggest that long-lived stationary and moving breathers are supported by the lattice. We find breather solutions which move in an arbitrary direction, an ellipticity criterion for the wavenumbers of the carrier wave, symptotic estimates for the breather energy, and a minimum threshold energy below which breathers cannot be found. This energy threshold is maximised for stationary breathers, and becomes vanishingly small near the boundary of the elliptic domain where breathers attain a maximum speed. Several of the results obtained are similar to those obtained for the square FPU lattice (Butt \& Wattis, {\em J Phys A}, {\bf 39}, 4955, (2006)), though we find that the square and hexagonal lattices exhibit different properties in regard to the generation of harmonics, and the isotropy of the generalised NLS equation.

Discrete breathers in a two-dimensional Fermi-Pasta-Ulam lattice

Using asymptotic methods, we investigate whether discrete breathers are supported by a two-dimensional Fermi-Pasta-Ulam lattice. A scalar (one-component) two-dimensional Fermi-Pasta-Ulam lattice is shown to model the charge stored within an electrical transmission lattice. A third-order multiple-scale analysis in the semi-discrete limit fails, since at this order, the lattice equations reduce to the (2+1)-dimensional cubic nonlinear Schrödinger (NLS) equation which does not support stable soliton solutions for the breather envelope. We therefore extend the analysis to higher order and find a generalised $(2+1)$-dimensional NLS equation which incorporates higher order dispersive and nonlinear terms as perturbations. We find an ellipticity criterion for the wave numbers of the carrier wave. Numerical simulations suggest that both stationary and moving breathers are supported by the system. Calculations of the energy show the expected threshold behaviour whereby the energy of breathers does {\em not} go to zero with the amplitude; we find that the energy threshold is maximised by stationary breathers, and becomes arbitrarily small as the boundary of the domain of ellipticity is approached.

Spectral approach for the soliton and periodic solutions of the nonlinear wave equation

A spectral method that obtains the soliton and periodic solutions to the nonlinear wave equation is presented. The results show that the nonlinear group velocity is a function of the frequency shift as well as of the soliton power. When the frequency shift is a function of time, a solution in terms of the Jacobian elliptic function is obtained. This solution is periodic in nature, and, to generate such an optical pulse train, one must simultaneously amplitude- and frequency-modulate the optical carrier. Finally, we extend the method to include the effect of self-steepening.

Asymptotic soliton train solutions of Kaup-Boussinesq equations

Asymptotic soliton trains arising from a 'large and smooth' enough initial pulse are investigated by the use of the quasiclassical quantization method for the case of Kaup-Boussinesq shallow water equations. The parameter varying along the soliton train is determined by the Bohr-Sommerfeld quantization rule which generalizes the usual rule to the case of 'two potentials' h(0)(x) and u(0)(x) representing initial distributions of height and velocity, respectively. The influence of the initial velocity u(0)(x) on the asymptotic stage of the evolution is determined. Excellent agreement of numerical solutions of the Kaup-Boussinesq equations with predictions of the asymptotic theory is found. (C) 2003 Elsevier B.V. All rights reserved.

Asymptotic soliton train solutions of the defocusing nonlinear Schrodinger equation

Asymptotic behavior of initially large and smooth pulses is investigated at two typical stages of their evolution governed by the defocusing nonlinear Schrodinger equation. At first, wave breaking phenomenon is studied in the limit of small dispersion. A solution of the Whitham modulational equations is found for the case of dissipationless shock wave arising after the wave breaking point. Then, asymptotic soliton trains arising eventually from a large and smooth initial pulse are studied by means of a semiclassical method. The parameter varying along the soliton train is calculated from the generalized Bohr-Sommerfeld quantization rule, so that the distribution of eigenvalues depends on two functions-intensity rho(0)(x) of the initial pulse and its initial chirp v(0)(x). The influence of the initial chirp on the asymptotic state is investigated. Excellent agreement of the numerical solution of the defocusing NLS equation with predictions of the asymptotic theory is found.

A family of interesting exact solutions of the sine-Gordon equation

By using AKNS [Phys. Rev. Lett. 31 (1973) 125] system and introducing the wave function, a family of interesting exact solutions of the sine-Gordon equation are constructed. These solutions seem to be some soliton, kink, and anti-kink ones respectively for the different choice of the spectrum, whereas due to the interaction between two traveling-waves they have some properties different from usual soliton, kink, and anti-kink solutions.

Modeling relativistic soliton interactions in overdense plasmas: A perturbed nonlinear Schrodinger equation framework

We investigate the dynamics of localized solutions of the relativistic cold-fluid plasma model in the small but finite amplitude limit, for slightly overcritical plasma density. Adopting a multiple scale analysis, we derive a perturbed nonlinear Schrodinger equation that describes the evolution of the envelope of circularly polarized electromagnetic field. Retaining terms up to fifth order in the small perturbation parameter, we derive a self-consistent framework for the description of the plasma response in the presence of localized electromagnetic field. The formalism is applied to standing electromagnetic soliton interactions and the results are validated by simulations of the full cold-fluid model. To lowest order, a cubic nonlinear Schrodinger equation with a focusing nonlinearity is recovered. Classical quasiparticle theory is used to obtain analytical estimates for the collision time and minimum distance of approach between solitons. For larger soliton amplitudes the inclusion of the fifth-order terms is essential for a qualitatively correct description of soliton interactions. The defocusing quintic nonlinearity leads to inelastic soliton collisions, while bound states of solitons do not persist under perturbations in the initial phase or amplitude

On asymptotic solutions of integrable wave equations

Asymptotic 'soliton train' solutions of integrable wave equations described by inverse scattering transform method with second-order scalar eigenvalue problem are considered. It is shown that if asymptotic solution can be presented as a modulated one-phase nonlinear periodic wavetrain, then the corresponding Baker-Akhiezer function transforms into quasiclassical eigenfunction of the linear spectral problem in weak dispersion limit for initially smooth pulses. In this quasiclassical limit the corresponding eigenvalues can be calculated with the use of the Bohr Sommerfeld quantization rule. The asymptotic distributions of solitons parameters obtained in this way specify the solution of the Whitham equations. (C) 2001 Elsevier B.V. B.V. All rights reserved.

Black soliton in a quasi-one-dimensional trapped fermion-fermion mixture

Employing a time dependent mean-field-hydrodynamic model we study the generation of black solitons in a degenerate fermion-fermion mixture in a cigar-shaped geometry using variational and numerical solutions. The black soliton is found to be the first stationary vibrational excitation of the system and is considered to be a nonlinear continuation of the vibrational excitation of the harmonic oscillator state. We illustrate the stationary nature of the black soliton, by studying different perturbations on it after its formation.

Soliton-cuspon interaction for the Camassa-Holm equation

The interaction of different kinds of solitary waves of the Camassa-Holm equation is investigated. We consider soliton-soliton, soliton-cuspon and cuspon-cuspon interactions. The description of these solutions had previously been shown to be reducible to the solution of an algebraic equation. Here we give explicit examples, numerically solving these algebraic equations and plotting the corresponding solutions. Further, we show that the interaction is elastic and leads to a shift in the position of the solitons or cuspons. We give the analytical expressions for this shift and represent graphically the coupled soliton-cuspon, soliton-soliton and cuspon-cuspon interactions.

Numerical and variational solutions of the dipolar Gross-Pitaevskii equation in reduced dimensions

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)