We study numerically the nonlinear dynamics of a shear banding interface in two dimensional planar shear flow, within the non-local Johnson Segalman model. Consistent with a recent linear stability analysis, we find that an initially flat interface is unstable with respect to small undulations for sufficiently small ratio of the interfacial width ℓ to cell length Lx. The instability saturates in finite amplitude interfacial fluctuations. For decreasing ℓ/Lx these undergo a non equilibrium transition from simple travelling interfacial waves with constant average wall stress, to periodically rippling waves with a periodic stress response. When multiple shear bands are present we find erratic interfacial dynamics and a stress response suggesting low dimensional chaos.PACS numbers: 47.50.+d, 47.20.-k, 36.20.-r. Complex fluids such as polymers, liquid crystals and surfactant solutions have mesoscopic structure that is readily perturbed by flow [1]. For example, wormlike surfactant micelles with lengths of the order of microns can be induced to stretch, disentangle and entangle, and break or increase in length. Their mechanical response is therefore highly non-Newtonian, with shear flows inducing normal stresses (e.g. σ xx − σ yy ), and with the shear stress σ xy being a nonlinear function of applied shear ratė γ. Recent work on such fluids has led to a fairly consistent picture of "shear banding": coexistence in shear flow of viscously thicker (nascent) and thinner (flow-induced) bands of material flowing at different local shear rates, for an overall average imposed shear rate. This phenomenon can be described by constitutive models for which the shear stress is a non-monotonic function of shear rate for homogeneous flow. This leads to a separation into bands of differing shear rate that coexist at a common shear stress [2,3]. Although most studies have assumed a flat interface between the bands and (with a few exceptions [4]) predicted a timeindependent banded state, an accumulating body of data has demonstrated that the average shear stress, and the banding interface, can fluctuate [5,6,7,8,9,10,11,19]. An important question is then whether these fluctuations resemble small amplitude capillary waves stabilised by surface tension, or whether they arise from an underlying instability, stabilised at large amplitude by nonlinearities. In this Letter we give strong evidence supporting the latter scenario, via the first theoretical study of the nonlinear dynamics of a shear banding interface.The model -The generalised Navier Stokes equation for a viscoelastic material in a Newtonian solvent of viscosity η and density ρ is:where v(r) is the velocity field. The pressure P is determined by incompressibility, ∇ · v = 0. The viscoelastic stress Σ(r) evolves according to the non-local ("diffusive") Johnson Segalman (DJS) model [12,13] (with plateau modulus G and relaxation time τ . D and Ω are the symmetric and antisymmetric parts of the velocity gradient tensor, (∇v) αβ ≡ ∂ α v β . For a = 1 and ℓ = 0 this model reduces t...