Abstract. We study the behavior of perturbations of small nonlinear Dirac standing waves. We assume that the linear Dirac operator of reference H = Dm + V has only two double eigenvalues and that degeneracies are due to a symmetry of H (theorem of Kramers). In this case, we can build a small four-dimensional manifold of stationary solutions tangent to the first eigenspace of H. Then we assume that a resonance condition holds, and we build a center manifold of real codimension 8 around each stationary solution. Inside this center manifold any H s perturbation of stationary solutions, with s > 2, stabilizes towards a standing wave. We also build center-stable and center-unstable manifolds, each one of real codimension 4. Inside each of these manifolds, we obtain stabilization towards the center manifold in one direction of time, while in the other, we have instability. Eventually, outside all of these manifolds, we have instability in the two directions of time. For localized perturbations inside the center manifold, we obtain a nonlinear scattering result. Introduction. We study the asymptotic stability of stationary solutions of a time-dependent nonlinear Dirac equation.A localized stationary solution of a given time-dependent equation represents a bound state of a particle. Like Rañada [39], we call it a particle-like solution (PLS). Many works have been devoted to the proof of the existence of such solutions for a wide variety of equations. Although their stability is a crucial problem (in particular, in a numerical computation or experiment), less attention has been devoted to this issue.In this paper, we deal with the problem of stability of small PLSs of the following nonlinear Dirac equation:where ∇F is the gradient of F : C 4 → R for the standard scalar product of R 8 . Here, D m is the usual Dirac operator (see Thaller [48]where m ∈ R * + , α = (α 1 , α 2 , α 3 ), and β are C 4 Hermitian matrices satisfying whereIn (0.1), V is the external potential field, and F : C 4 → R is a nonlinearity with the following gauge invariance:Some additional assumptions on F and V will be made in what follows. Nonlinearity with no potential arises in some Dirac models introduced by physicists to model either extended particles with self-interaction or particles in space-time with geometrical structure. In the latter case, physicists have shown that a relativistic theory sometimes imposes a fourth order nonlinear potential (i.e., a cubic nonlinearity) such as the square of a quadratic form on C 4 ; see Rañada [39] and the references therein. We added a potential with special features to ensure the existence of small stationary solutions, whose existence and stability are easier to study.Stationary solutions (PLSs) of (0.1) take the form ψ(t, x) = e −iEt φ(x), where φ satisfiesWe show that there exists a manifold of small solutions to (0.3) tangent to the first eigenspace of D m + V (see Proposition 1.1 below).Concerning the asymptotic stability in the Schrödinger equation, the question has been solved in several cases. Fo...