Exact numerical solutions for dark waves on the discrete nonlinear Schrödinger equation

In this paper we study numerically existence and stability of exact dark waves on the (nonintegrable) discrete nonlinear Schrödinger equation for a finite one-dimensional lattice. These are solutions that bifurcate from stationary dark modes with constant background intensity and zero intensity at a site, and whose initial state translates exactly one site each period of the internal oscillations. We show that exact dark waves are characterized by an oscillatory background whose wavelength is closely related with the velocity. Faster dark waves require smaller wavelengths. For slow enough velocity dark waves are linearly stable, but when trying to continue numerically a solution towards higher velocities bifurcations appear, due to rearrangements in the oscillatory tail in order to make possible a decreasing of the wavelength. However, in principle, one might control the stability of an exact dark wave adjusting a phase factor which plays the role of a discreteness parameter. In addition, we also study the regimes of existence and stability for stationary discrete gray modes, which are exact solutions with phase-twisted constant-amplitude background and nonzero minimum intensity. Also such solutions develop envelope oscillations on top of the homogeneous background when continued into moving phase-twisted solutions.