Large self-gravitating stellar systems share with correlated liquids in condensed matter physics a pattern of hierarchical density variations. While it takes the microscopic time resolution to discern the correlated dynamics of the critical opalescence, characteristic astronomical times hide fluctuational dynamics of stellar liquids, where, governed by interstellar correlations, denser clusters of stars assemble and disperse. For a semi-isolated galaxy, these dynamical fluctuations are dense globular clusters. For a galaxy cluster, these dynamical fluctuations are the member galaxies, elliptical ones in the interior. Bright over-density fluctuations, galaxies, are exhibits of only a small fraction of stars found in a cosmic stellar liquid, the dark matter. Here I report a fluctuational gravitational collapse as a property of a self-gravitating system in the virial equilibrium.

The most probable state of an infinite self-gravitating gas in the dynamical equilibrium is defined by `gravitational haziness', a parameter representing many-body effects and formally like the temperature in the case of thermal equilibrium. A kinetic equation is constructed using a concept of statistical equipartition of the virial among subsystems of the self-gravitating gas. A closed equation for the gravitational potential is conjectured as a special property of the kinetic equation. An equilibrium particle distribution function in the phase space, an analog of the Maxwell-Boltzmann weight, and a galaxy equation of state are found for all `gravitational haziness'. The first law of a `hazydynamics' (thermodynamics) states that the total mass of an astronomical stellar collection is the sum of the Archimedes displaced mass and an excess `gobbled' mass determined by the `gravitational haziness' and history.