A Monte Carlo Method for Calculating Lynden-Bell Equilibrium in Self-Gravitating Systems

We present a Monte Carlo approach that allows us to easily implement Lynden-Bell (LB) entropy maximization for an arbitrary initial particle distribution. The direct maximization of LB entropy for an arbitrary initial distribution requires an infinite number of Lagrange multipliers to account for the Casimir invariants. This has restricted studies of Lynden-Bell's violent relaxation theory to only a very small class of initial conditions of a very simple waterbag form, for which the entropy maximization can be performed numerically. In the present approach, an arbitrary initial distribution is discretized into density levels which are then evolved using an efficient Monte Carlo algorithm towards the final equilibrium state. A comparison is also made between the LB equilibrium and explicit Molecular Dynamics simulations. We find that for most initial distributions, relaxation is incomplete and the system is not able to reach the state of maximum LB entropy. In particular, we see that the tail of the stationary particle distribution is very different from the one predicted by the theory of violent relaxation, with a hard cutoff instead of an algebraic decay predicted by LB's theory.