Building on the development of a Hermite-Legendre analysis of one-dimensional gravitating collisionless systems, we present a technique for determining the steady states of such systems. This provides an important component for understanding the physics involved in the relaxation of these kinds of systems. As the dark matter structures in the universe should have traits in common with these systems, insight into this relaxation can provide clues to larger astrophysical questions. For large perturbation strengths, we determine physically motivated parameter ranges for the simplest families of steady states as well as their stability. We also demonstrate that any set of initial conditions in the linear regime can be resolved into unique time-independent and time-dependent modes. Combinations of time-independent modes then describe the steady state of any system linearly perturbed from equilibrium. These results highlight the importance of initial conditions over relaxation mechanisms in the evolution of these systems.

In an effort to better understand collisionless relaxation processes in gravitational systems, we investigate one-dimensional models. Taking advantage of a Hermite-Legendre expansion of relevant distribution functions, we present analytical and numerical behaviors of Maxwell-Boltzmann entropy. In particular, we modestly perturb systems about a separable-solution equilibrium and observe their collisionless evolution to a steady state. We verify the time-independence of fine-grained entropy in these systems before turning our attention to the behavior of coarse-grained entropy. We also verify that there is no analogue to the collisional H-theorem for these systems. Competing terms in the second-order coarse-grained entropy make it impossible to guarantee continuously increasing entropy. However, over dynamical time-scales the coarse-grained entropy generally increases, with small oscillations occurring. The lack of substantive differences between the entropies in test-particle and self-gravitating cases suggests that phase mixing, rather than violent relaxation associated with potential changes, more significantly drives the coarse-grained entropy evolution. The effects of violent relaxation can be better quantified through analysis of energy distributions rather than phase-space distributions.

The current paradigm for understanding galaxy formation in the universe depends on the existence of self-gravitating collisionless dark matter. Modeling such dark matter systems has been a major focus of astrophysicists, with much of that effort directed at computational techniques. Not surprisingly, a comprehensive understanding of the evolution of these self-gravitating systems still eludes us, since it involves the collective nonlinear dynamics of many-particle systems interacting via long-range forces described by the Vlasov equation. As a step towards developing a clearer picture of collisionless self-gravitating relaxation, we analyze the linearized dynamics of isolated one-dimensional systems near thermal equilibrium by expanding their phase space distribution functions f(x,v) in terms of Hermite functions in the velocity variable, and Legendre functions involving the position variable. This approach produces a picture of phase-space evolution in terms of expansion coefficients, rather than spatial and velocity variables. We obtain equations of motion for the expansion coefficients for both test-particle distributions and self-gravitating linear perturbations of thermal equilibrium. N-body simulations of perturbed equilibria are performed and found to be in excellent agreement with the expansion coefficient approach over a time duration that depends on the size of the expansion series used.