We study the spreading of density-density correlations in Bose-Hubbard models after a quench of the interaction strength, using time-dependent variational Monte Carlo simulations. It gives access to unprecedented long propagation times and to dimensions higher than one. In both one and two dimensions, we find ballistic light-cone spreading of correlations and extract accurate values of the light-cone velocity in the superfluid regime. We show that the spreading of correlations is generally supersonic, with a light-cone propagating faster than sound modes but slower than the maximum group velocity of density excitations, except at the Mott transition, where all the characteristic velocities are equal. Further, we show that in two dimensions the correlation spreading is highly anisotropic and presents nontrivial interference effects.
