Lieb-Robinson-type bounds are reported for a large class of classical Hamiltonian lattice models. By a suitable rescaling of energy or time, such bounds can be constructed for interactions of arbitrarily long range. The bound quantifies the dependence of the system's dynamics on a perturbation of the initial state. The effect of the perturbation is found to be effectively restricted to the interior of a causal region of logarithmic shape, with only small, algebraically decaying effects in the exterior. A refined bound, sharper than conventional Lieb-Robinson bounds, is required to correctly capture the shape of the causal region, as confirmed by numerical results for classical long-range XY chains. We discuss the relevance of our findings for the relaxation to equilibrium of long-range interacting lattice models.In many nonrelativistic lattice systems, and despite the absence of Lorentz covariance, physical effects are mostly restricted to a causal region, often in the shape of an effective "light cone," with only tiny effects leaking out to the exterior. The technical tool, known as the LiebRobinson bound [1,2], to quantify this statement in a quantum mechanical context is an upper bound on the norm of the commutator ½O A ðtÞ; O B ð0Þ, where O A ð0Þ and O B ð0Þ are operators supported on the subspaces of the Hilbert space corresponding to nonoverlapping regions A and B of the lattice. The importance of such a bound lies in the fact that a multitude of physically relevant results can be derived from it. Examples are bounds on the creation of equal-time correlations [3], on the transmission of information [4], and on the growth of entanglement [5], the exponential spatial decay of correlations in the ground state of a gapped system [6], or a Lieb-Schultz-Mattis theorem in higher dimensions [7]. Experimental observations related to Lieb-Robinson bounds have also been reported [8].The original proof by Lieb and Robinson [1] requires interactions of finite range. An extension to power-lawdecaying long-range interactions has been reported in Refs. [3,6]. In this case the effective causal region is no longer cone shaped, and the spatial propagation of physical effects is not limited by a finite group velocity [9]. For "strong long-range interactions," i.e., when the interaction potential decays proportionally to 1=r α with an exponent α smaller than the lattice dimension d, the theorems in [3,6] do not apply and no Lieb-Robinson-type results are known. This fact nicely fits into the larger picture that, for α ≤ d, the behavior of long-range interacting systems often differs substantially from that of short-range interacting systems. Examples of such differences include nonequivalent equilibrium statistical ensembles and negative response functions [10], or the occurrence of quasistationary states whose lifetimes diverge with the system size [11,12]. The latter is a dynamical phenomenon, and it has been conjectured in [13] that some of its properties are universal and in some way connected to Lieb-Robinson bounds.In mos...