We present numerical results within the one-dimensional disordered Hubbard model for several characteristic indicators of the many-body localization (MBL). Considering traditionally studied charge disorder (i.e., the same disorder strength for both spin orientations) we find that even at strong disorder all signatures consistently show that while charge degree of freedom is non-ergodic, the spin is delocalized and ergodic. This indicates the absence of the full MBL in the model that has been simulated in recent cold-atom experiments. Full localization can be restored if spin-dependent disorder is used instead. 71.27.+a, 71.30.+h, 71.10.Fd Introduction.-The many-body localization (MBL) is a phenomenon whereby an interacting many-body system localizes due to disorder, proposed [1, 2] in analogy to the Anderson localization of noninteracting particles [3, 4]. The MBL physics has attracted a broad attention of theoreticians. Yet, it has so far been predominantly studied within the prototype model, i.e., the one-dimensional (1D) model of interacting spinless fermions with random potentials, equivalent to the anisotropic spin-1/2 Heisenberg chain with random local fields. Emerging from these studies are main hallmarks of the MBL state of the system: a) the Poisson many-body level statistics [5][6][7][8][9], in contrast to the Wigner-Dyson one for normal ergodic systems, b) vanishing of d.c. transport at finite temperatures T > 0, including the T → ∞ limit [10][11][12][13][14][15][16][17], c) logarithmic growth of the entanglement entropy [18][19][20], as opposed to linear growth in generic systems, d) an existence of a set of local integrals of motion [21][22][23][24], and e) a non-ergodic time evolution of (all) correlation functions and of quenched initial states [25][26][27][28][29]. Because of these unique properties, the MBL can be used, e.g. to protect quantum information [30,31]. For more detailed review see Refs. [32,33].The experimental evidence for the MBL comes from recent experiments on cold atoms in optical lattices [34][35][36][37] and ion traps [38]. In particular, for strong disorders, experiments reveal non-ergodic decay of the initial density profile in uncoupled [34] and coupled [36] 1D fermionic chains, as well as the vanishing of d.c. mobility in a 3D disordered lattice [35]. In contrast to most numerical studies, being based on the spinless fermion models, the cold-atom experiments simulate a disordered Hubbard model. The latter has been much less investigated theoretically [34,39,40], whereby results show that density imbalance might be non-ergodic at strong disorder [34,39], in accordance with experiments [34,36].The essential difference with respect to the interacting spinless model is that Hubbard model has two local degrees of freedom: charge (density) and spin. The aim of this Letter is to present numerical evidence that in the case of a (charge) potential disorder and finite repulsion U > 0 (as e.g. realized in the cold-atom experiments), both degrees behave qualitatively different. In...