We consider one-dimensional infinite chains of harmonic oscillators with stochastic perturbations and long-range interactions which have polynomial decay rate x −θ , x → ∞, θ > 1, where x ∈ Z is the interaction range. We prove that if 2 < θ ≤ 3, then the time evolution of the macroscopic thermal energy distribution is superdiffusive and governed by a fractional diffusion equation with exponent 3 7−θ , while if θ > 3, then the exponent is 3 4 . The threshold is θ = 3 because the derivative of the dispersion relation diverges as k → 0 when θ ≤ 3.