Existence and stability of multihumped femtosecond solitons

The propagation of bound soliton pairs in nonlinear photonic crystal fibers has recently been experimentally observed. The system may be modeled by a generalized nonlinear Schrodinger equation (GNLSE) which includes higher intrapulse Raman Scattering, self-steepening and higher order dispersion. Here, we find multihumped pulses as result of an accelerating similarity reduction of a GNLSE containing the intrapulse Raman scattering. Numerical simulations of the suitable GNLSE using these solutions as input showed that they are not stable, however, they may be related with the experimentally observed bound pairs since they propagate steadily for distances compared to the ones observed.