At low temperature, a quasi-one-dimensional ensemble of atoms with attractive interaction forms a bright soliton. When exposed to a weak and smooth external potential, the shape of the soliton is hardly modified, but its center-of-mass motion is affected. We show that in a spatially correlated disordered potential, the quantum motion of a bright soliton displays Anderson localization. The localization length can be much larger than the soliton size and could be observed experimentally.PACS numbers: 03.75. Lm,72.15.Rn,05.30.Jp At zero temperature, cold atoms interacting attractively in a one dimensional (1D) system tend to cluster together, forming a bright soliton. Explicit solutions of the many-body problem can be found in some cases, for example for contact interactions [1]. Using external potentials, it has been experimentally shown how to put solitons in motion [2]. What happens to a soliton exposed to a disordered potential? If the potential is strong, it will destroy the soliton. If it is sufficiently weak and smooth not to perturb the soliton shape, one expects the soliton to undergo multiple scattering, diffusive motion and possibly Anderson localization [3]. Indeed, propagation of waves in a disordered potential is profoundly affected by Anderson localization. Multiple scattering on random defects yields exponentially localized density profiles and a suppression of the usual diffusive transport associated with incoherent wave scattering [4]. In 1D Anderson localization is a ubiquitous phenomenon [5], which has been recently observed for cold atomic matter waves [6]. It is important to understand how it is modified when interactions between particles are taken into account.We consider a Bose-Einstein condensate in a quasi-1D geometry. Within mean-field theory, it is described by the Gross-Pitaevskii energy functionalin units of E 0 = 4mω 2 ⊥ a 2 , l 0 = /2|a|mω ⊥ , and t 0 = /4a 2 mω 2 ⊥ for energy, length and time, respectively. Here, ω ⊥ denotes the transverse harmonic confinement frequency, a the atomic s-wave scattering length, and µ the chemical potential. The cases of repulsive and attractive atomic interaction are covered by g = ±1.The dynamics is extremely different in both cases, so that we first discuss the case of attractive interaction, g = −1. The ground state of (1) is the bright soliton [7] φ 0 (z − q) = N 2ξnormalized to the total number of particles N . The chemical potential is µ = −N 2 /8 and the soliton width is ξ = 2/N. This ground-state solution has an arbitrary center-of-mass (CM) position q and an arbitrary global phase θ that spontaneously break the translational and the U (1) gauge symmetry of the energy functional (1), respectively. These degrees of freedom appear as zeroenergy modes of Bogoliubov theory, and their quantum dynamics requires special attention [8,9]. The energy functional (1) is no longer translation invariant when a potential term dzV (z)|φ| 2 is added. If V (z) is sufficiently weak and smooth, the soliton shape remains unchanged to lowest order in V , and on...