Using a numerically exact technique we study spin transport and the evolution of spin-density excitation profiles in a disordered spin-chain with long-range interactions, decaying as a power-law, r −α with distance and α < 2. Our study confirms the prediction of recent theories that the system is delocalized in this parameters regime. Moreover we find that for α > 3/2 the underlying transport is diffusive with a transient super-diffusive tail, similarly to the situation in clean long-range systems. We generalize the Griffiths picture to long-range systems and show that it captures the essential properties of the exact dynamics.Introduction.-Many-body localization (MBL) extends the notion of Anderson localization to interacting systems [1]. For local interactions, its existence is well established theoretically [2, 3] and experimentally in one-dimensional systems [4-6] (see [7] for a recent review), but there is evidence of localization also in twodimensional systems [8][9][10][11][12][13]. For long-range interactions the fate of MBL is less clear. Some studies suggest that the many-body localization is stable for α > 2d [14-17], and some claim that the system is delocalized for all α in the thermodynamic limit [17][18][19]. Finite size systems of size L are claimed to exhibit an effective many-body-like localization transition at a critical disorder which scales with the system size (as a power-law for α < 2d), and diverges in the thermodynamic limit [14,16,17,19,20]. Understanding the dynamics of disordered systems with long-range interactions is of great importance to a number of physical systems, such as nuclear spins [21], dipoledipole interactions of vibrational modes [22][23][24], Frenkel excitons [25], nitrogen vacancy centers in diamond [26][27][28][29][30] and polarons [31]. Long range interactions are also common in atomic and molecular systems, where interactions can be dipolar [32][33][34][35][36][37], van der Waals like [32,38], or even of variable range [39][40][41][42]. Some aspects of the dynamics in such systems were studied numerically in Ref. [43], analytically in Ref. [44] and experimentally in Ref. 45, however spin transport in such systems was not considered.The delocalized phase of one-dimensional systems with local interactions, shows subdiffusive transport [46][47][48][49][50], accompanied by sublinear growth of the entanglement entropy [51][52][53] and intermediate statistics of eigenvalue spacing [54]. Anomalous transport is commonly explained by rare insulating regions, which effectively suppress transport in one-dimensional systems. This mechanism is known as the Griffith's picture [48,55,56] (see Ref.[57] for a recent review and also Ref.[58] were rare regions were introduced externally). In dimensions higher than one the Griffiths picture predicts diffusion, since rare regions can be circumvented [56], however approximate numerical studies [8] as also recent experiments [10,11] suggest that at least for short to inter-arXiv:1911.07857v1 [cond-mat.dis-nn]