Hamiltonian systems with long-range interactions give rise to long lived out of equilibrium macroscopic states, so-called quasi-stationary states. We show here that, in a suitably generalized form, this result remains valid for many such systems in the presence of dissipation. Using an appropriate mean-field kinetic description, we show that models with dissipation due to a viscous damping or due to inelastic collisions admit "scaling quasi-stationary states", i.e., states which are quasi-stationary in rescaled variables. A numerical study of one dimensional self-gravitating systems confirms both the relevance of these solutions, and gives indications of their regime of validity in line with theoretical predictions. We underline that the velocity distributions never show any tendency to evolve towards a Maxwell-Boltzmann form.

One dimensional versions of cosmological N-body simulations have been shown to share many qualitative behaviours of the three dimensional problem. They can resolve a large range of time and length scales, and admit exact numerical integration. We use such models to study how non-linear clustering depends on initial conditions and cosmology. More specifically, we consider a family of models which, like the 3D EdS model, lead for power-law initial conditions to self-similar clustering characterized in the strongly non-linear regime by power-law behaviour of the two point correlation function. We study how the corresponding exponent \gamma depends on the initial conditions, characterized by the exponent n of the power spectrum of initial fluctuations, and on a single parameter \kappa controlling the rate of expansion. The space of initial conditions/cosmology divides very clearly into two parts: (1) a region in which \gamma depends strongly on both n and \kappa and where it agrees very well with a simple generalisation of the so-called stable clustering hypothesis in three dimensions, and (2) a region in which \gamma is more or less independent of both the spectrum and the expansion of the universe. We explain the observed location of the boundary in (n, \kappa) space dividing the "stable clustering" region from the "universal" region. We compare and contrast our findings to results in three dimensions, and discuss in particular the light they may throw on the question of "universality" of non-linear clustering in this context.

Studies of a class of infinite one dimensional self-gravitating systems have highlighted that, on the one hand, the spatial clustering which develops may have scale invariant (fractal) properties, and, on the other, that they display "self-similar" properties in their temporal evolution. The relevance of these results to three dimensional cosmological simulations has remained unclear. We show here that the measured exponents characterizing the scale-invariant non-linear clustering are in excellent agreement with those derived from an appropriately generalized "stable-clustering" hypothesis. Further an analysis in terms of "halos" selected with a friend-of-friend algorithm reveals that such structures are, statistically, virialized across the range of scales corresponding to scale-invariance. Thus the strongly non-linear clustering in these models is accurately described as a virialized fractal structure, very much in line with the "clustering hierarchy" which Peebles originally envisaged qualitatively as associated with stable clustering. If transposed to three dimensions these results would imply, notably, that cold dark matter halos (or even subhalos) are 1) not well modeled as smooth objects, and 2) that the supposed "universality" of their profiles is, like apparent smoothness, an artefact of poor numerical resolution.

We formalize a classification of pair interactions based on the convergence properties of the {\it forces} acting on particles as a function of system size. We do so by considering the behavior of the probability distribution function (PDF) P(F) of the force field F in a particle distribution in the limit that the size of the system is taken to infinity at constant particle density, i.e., in the "usual" thermodynamic limit. For a pair interaction potential V(r) with V(r) \rightarrow \infty) \sim 1/r^a defining a {\it bounded} pair force, we show that P(F) converges continuously to a well-defined and rapidly decreasing PDF if and only if the {\it pair force} is absolutely integrable, i.e., for a > d-1, where d is the spatial dimension. We refer to this case as {\it dynamically short-range}, because the dominant contribution to the force on a typical particle in this limit arises from particles in a finite neighborhood around it. For the {\it dynamically long-range} case, i.e., a \leq d-1, on the other hand, the dominant contribution to the force comes from the mean field due to the bulk, which becomes undefined in this limit. We discuss also how, for a \leq d-1 (and notably, for the case of gravity, a=d-2) P(F) may, in some cases, be defined in a weaker sense. This involves a regularization of the force summation which is generalization of the procedure employed to define gravitational forces in an infinite static homogeneous universe. We explain that the relevant classification in this context is, however, that which divides pair forces with a > d-2 (or a < d-2), for which the PDF of the {\it difference in forces} is defined (or not defined) in the infinite system limit, without any regularization. In the former case dynamics can, as for the (marginal) case of gravity, be defined consistently in an infinite uniform system.

The dynamics of infinite, asymptotically uniform, distributions of self-gravitating particles in one spatial dimension provides a simple toy model for the analogous three dimensional problem. We focus here on a limitation of such models as treated so far in the literature: the force, as it has been specified, is well defined in infinite point distributions only if there is a centre of symmetry (i.e. the definition requires explicitly the breaking of statistical translational invariance). The problem arises because naive background subtraction (due to expansion, or by "Jeans' swindle" for the static case), applied as in three dimensions, leaves an unregulated contribution to the force due to surface mass fluctuations. Following a discussion by Kiessling, we show that the problem may be resolved by defining the force in infinite point distributions as the limit of an exponentially screened pair interaction. We show that this prescription gives a well defined (finite) force acting on particles in a class of perturbed infinite lattices, which are the point processes relevant to cosmological N-body simulations. For identical particles the dynamics of the simplest toy model is equivalent to that of an infinite set of points with inverted harmonic oscillator potentials which bounce elastically when they collide. We discuss previous results in the literature, and present new results for the specific case of this simplest (static) model starting from "shuffled lattice" initial conditions. These show qualitative properties (notably its "self-similarity") of the evolution very similar to those in the analogous simulations in three dimensions, which in turn resemble those in the expanding universe.