2.3. Gravitational stability of flat systems
Lin asked [Woltjer in 1961]: What are the circumstances that
would be needed for either one or both of the stellar and inter
stellar parts of a supposedly smooth galactic disk to remain
gravitationally stable against all large scale disturbances?
Toomre 1964, p.1217
The importance of collective effects in our Galaxy was first
clearly pointed out by
Toomre (1964).
He showed that in the
disk the stellar motions are sufficiently coherent to make it
almost vulnerable to collapse. He also pointed out that the
scale on which this would occur is quite large.
Kalnajs 1971, p.275
As we have seen, Safronov already raised the question of gravitational instability in flat rotating systems, aiming at the breakup of a protoplanetary cloud into detached rings. Toomre, interested in basically smoother objects like galaxies, turned in 1961 to a rather close, although opposite in accent, topic, and by the summer of 1963 he prepared an article "On the gravitational stability of a disk of stars" (Toomre 1964, hereinafter T64).
The paper started with the general presentation of the problem as it was then seen.
"The wellknown instabilities of those Maclaurin spheroids whose rotational flattening exceeds a certain fairly moderate value suggest that the other sufficiently flattened, rotating, and selfgravitating systems might in some sense likewise be unstable. At any rate, these instabilities have been often cited as a likely reason why one does not observe elliptical galaxies exceeding a certain degree of oblateness. It is only when we turn to consider what are now thought to be the distributions of all but the youngest stars in the disks of the ordinary (as opposed to the barred) spiral galaxies that this classical result suggests a serious dilemma: How is it conceivable, in spite of these or analogous instabilities, that so much of the fainter stellar matter within such galaxies  and certainly the S0 galaxies  should today appear distributed relatively evenly over disks with something like a tentoone flattening?" (T64, p.1217)
The detailed study of the problem was preceded by a primary, qualitative stability estimate.
A rotating thin cold disk, in an approximate equilibrium between gravity and centrifugal forces acting on each mass element, is prevented from general contraction, still not from fragmentation. Smallsize clumpings arise everywhere in such a disk, and then collapse, their gravity taking excess over rotation. But if largersized, they do not go as these two factors counteract each other. The demarcation length scale L_{T} proves plain comeasurable with the disk radius R. Thus the cold model, for all specifications it may have, is clearly unstable. ^{47} The part played by random motions is best visualized with an immovable sheet model. There instability is avoided if stars (other mass elements), having an rms velocity ñ, cross a clumping zone in a time not exceeding that needed for an e fold amplitude growth as registered in the cold case. Hence the largest yet ungrowing disturbance is found on an L_{J} c^{2} / Gµ_{0} scale, which is essentially the Jeans stability criterion. Now, letting the sheet rotate, one sees the two characteristic scales, L_{T} and L_{J}, be present (Fig.5). L_{J} gets closer to L_{T} for higher velocity dispersions, until they coincide at c's as high  in the order of magnitude  as the rotational velocity, thus meaning full stabilization against this sort of disturbances.
The strict analysis of axisymmetric disturbances to a razorthin disk, performed in T64, supported these rough estimates. In the cold case, it led to a local dispersion relation
(5a) 
or
(5b) 
linking the wave frequency in units of , = / , with a critical wavenumber
(6) 
the one to determine the shortest wavelength _{T} 2 / k_{T} of ungrowing (^{2} 0) disturbances (Fig.6). ^{48} The hotdisk analysis detected the minimum radial velocity dispersion at which the system is still resistant against all axisymmetric disturbances (Fig.7): ^{49}^{,}^{50}
(7) 
The realtominimum velocitydispersion ratio
(8) 
thus got a local diskstability parameter. ^{51}^{,}^{52} In a marginally stable state Q = 1, disturbances of _{0} 0.55 _{T} proved most unpliant and barely suppressible. Our solar neighborhood would have such a _{0} 5  8 kpc, but if some Q 1  1.5 were not preferred empirically, implying a certain stability reserve. Of course, "it was as yet impossible to rule out instabilities altogether", but should any actually be present, they would not do with scales responding to the challenging 2kpc spacings, as these "must almost certainly be judged as stable". This "is important as an argument against any suggestion that the existing spiral structure in this Galaxy might be the result of collective stellar instabilities" of the sort considered (T64, p.1236).
Figure 6. (left) The dispersion relation
curve for radial oscillations and tightly wrapped spiral waves in a
cold disk.

Still, the linear theory developed could not lay claim to very much. So it did not elucidate the cause of stellar disk heating, it even could not show any definitely what was to become with primary condensations appearing in a tentatively cold disk in one or two revolutions already. "It must not be presumed that such initial clumpings would necessarily have led to the formation of any permanent irregularities", Toomre noticed. "On the contrary, it seems much more likely that the bulk of the stars involved in any given (generally nonaxisymmetric) instability [...] would eventually have dispersed themselves upon emerging from the opposite sides of the aggregation and upon experiencing the shearing effect of differential rotation".
"It follows that an initially unstable disk of stars should probably have undergone not just one but several successive generations of instabilities, after each of which the system would have been left somewhat less unstable than it was previously. In particular, it seems likely that before very many rotation periods had elapsed, the disk would have approached a new equilibrium state that was again fairly regular and quite possibly axisymmetric, but in which the random velocities at the various radii had become  and would henceforth remain  about equal to the minimum values needed for complete stability: (T64, p.1237)." ^{53}
Besides, since the total gravitational energy of the disk would have had to be the same during its evolution (the virial theorem), "the said redistribution of stars could not simply have consisted of an overall contraction, but would have had to entail a contraction perhaps in the inner parts of the disk jointly with a net expansion of the outer portions" (T64, p.1237)  as it was already seen by LyndenBell (1960b) from the gasdynamical viewpoint.
As regards nonaxisymmetric disturbances, it was pointed out in T64 that because of the specific action of the Coriolis force those are restrained even more effectively than radial disturbances, thus requiring no addition for c_{r, min}. However, Toomre remarked, a question that his discussion left "completely unanswered" was "to what extent a similar amount of random motion [Q = 1] might affect the character of the most extensive nonaxisymmetric disturbances, in particular those which ought to determine whether or not a given disk might prefer to develop into a barlike structure" (T64, p.1235). ^{54}
^{47} Toomre got this estimate by the fall of 1961 and was struck with the fact that nothing had ever been said on the thing just shocking with its as simply derivable inference that cold disks be prone to violent instabilities. (Toomre) Back.
^{48} Analyzing axisymmetric disturbances to a flattened rotating cloud, Safronov (1960a, b) did not solve the Poisson equation. He was guided by the notice that short radial waves find adequate the cylindric approximation for a torus (ring). But the cylinder is the sum of `rods', or elementary cylinders whose individual gravity is given by a simple formula, so that the business is just to integrate in infinite limits the elementary contributions over longitudal and transversal variables xand z. There Safronov was not perfect, however. His gently stratified cloud turned a stiff 2h thick plate as he took his introduced density function _{0}(z) out of integration over z. His subsequent integration over xwas in an interval of ± /4; that, he argued, ensured a predominant contribution to the perturbed force (which is qualitatively true). Had he integrated in infinite limits, and first  most trivially  over x, the gravity term in his Eqn (3) would have become 2 Gk _{0}(z / h)e^{kz} dz, and with the exponential factor serving as a thickness correction he would have accurately managed with any density profile  and, most obviously, would have found that in the zerothickness limit that factor simplifies to unity, the integral just gives the surface density µ_{0}, so that the gravity term converts into  2 G µ_{0} k, the form in which it was presented soon by Toomre (1964) in frames of `regular' methods of the potential theory. Back.
^{49} To solve the Vlasov kinetic equation, Toomre used the characteristics method that for some threedimensional purposes had already served LyndenBell (1962), who in his turn cited the original source (Bernstein 1958) where that method had genuinely helped with the general disperion relation for the mathematically similar problem with a Maxwellian plasma in a magnetic field. Back.
^{50} Because of a technical error in Toomre's analysis, this minimum value was initially overestimated by 20%. Not so little if one considers that the difference in c_{r, min} for star and gas disk models (the latter case admits a much simpler analysis) reaches 7% only. It is this "substantial error" which was detected in 1963 by Kalnajs (cf. Sect. 2.4), as reported frankly in T64 (p.1233). Back.
^{51} Formally, the `Q parameter' (8) was introduced by Julian and Toomre (1966). Back.
^{52} This quantitative analysis refines the above view of disk stabilization as it shows via Eqs (6) and (7) that locally the result is attained already once L_{J} / L_{T} = (3.36/2 )^{2} 0.286 (0.25 in a gas disk). Back.
^{53} Asked to reminisce on how he had originally understood those dispersion velocities "about equal" to the needed minimum in the new equilibrium state  on whether or not this was a factual suggestion of marginal stability of our stellar disk, or some extra amount was yet permitted for its stability  Toomre has responded: "It is hard for me to reconstruct from this vantage point what exactly I meant or hoped by that statement. Probably I was mostly just trying to rationalize the surprising fact which I had then unearthed that the minimum theoretically needed c_{r, min} and the observed amounts seemed to agree so well within their considerable uncertainties, meaning within a factor of 1.5 or thereabouts, rather than some 2 or 3 or 4 [...] From about 1966 onwards, I was surely of the opinion that any Q less than about 1.5 here was highly suspect, if not downright ludicrous, because of fierce heating of cooler disks by their embedded gas complexes. But that came a little later. In 1964 my views were no doubt more permissive toward Q = 1.0". (Toomre) Back.
^{54} Real progress in the study of this problem first came half a decade later. Back.