Gravitational instability and affine dynamics of gaseous astrophysical discs

We develop several aspects of the theory of gaseous astrophysical discs in which the gravity of the disc makes a significant contribution to its structure and dynamics. We show how the internal gravitational potential can be expanded in powers of the aspect ratio of the disc (or of a structure within it) and separated into near and far contributions. We analyse the hydrostatic vertical structure of a wide family of disc models, both analytically and numerically, and show that the near contribution to the internal gravitational potential energy can be written in an almost universal form in terms of the surface density and scaleheight. We thereby develop an affine model of the dynamics of (generally non-hydrostatic) self-gravitating discs in which this contribution to the energy acts as a gravitational pressure in the plane of the disc. This combines with and significantly reinforces the gas pressure, allowing us to define an enhanced effective sound speed and Toomre stability parameter Q for self-gravitating discs. We confirm that this theory fairly accurately reproduces the onset of axisymmetric gravitational instability in discs with resolved vertical structure. Among other things, this analysis shows that the critical wavelength is on the order of twenty times the scaleheight, helping to justify the validity of the affine model. The weakly nonlinear theory also typically exhibits subcritical behaviour, with equilibrium solutions of finite amplitude being found in the linearly stable regime Q > 1 for adiabatic exponents less than 1.50.