It can be shown that any sufficiently fast decaying smooth solution will eventually split into a finite superposition of solitons travelling to the right plus a decaying dispersive part travelling to the left. This was first observed by Zabusky & Kruskal (1965) and can be rigorously proven using the nonlinear steepest descent analysis for oscillatory Riemann–Hilbert problems.