Currently there is great interest in the outcomes and astrophysical implications of mergers of double degenerate binaries. In a commonly adopted approximation, the components of such binaries are represented by polytropes with an index n=3/2. We present detailed comparisons of stellar mass-transfer and merger simulations of polytropic binaries that have been carried out using two very different numerical algorithms --- a finite-volume "grid" code and a smoothed-particle hydrodynamics (SPH) code. We find that there is agreement in both the ultimate outcomes of the evolutions and the intermediate stages if the initial conditions for each code are chosen to match as closely as possible. We find that even with closely matching initial setups, the time it takes to reach a concordant evolution differs between the two codes because the initial depth of contact cannot be matched exactly. There is a general tendency for SPH to yield higher mass transfer rates and faster evolution to the final outcome. We also present comparisons of simulations calculated from two different energy equations: in one series we assume a polytropic equation of state and in the other series an ideal gas equation of state. In the latter series of simulations an atmosphere forms around the accretor which can exchange angular momentum and cause a more rapid loss of orbital angular momentum. In the simulations presented here, the effect of the ideal equation of state is to de-stabilize the binary in both SPH and grid simulations, but the effect is more pronounced in the grid code.

We present a numerical method for the study of double white dwarf (DWD) binary systems at the onset of super-Eddington mass transfer. We incorporate the physics of ideal inviscid hydrodynamical flow, Newtonian self-gravity, and radiation transport on a three-dimensional uniformly rotating cylindrical Eulerian grid. Care has been taken to conserve the key physical quantities such as angular momentum and energy. Our new method conserves total energy to a higher degree of accuracy than other codes that are presently being used to model mass-transfer in DWD systems. We present the results of verification tests and we simulate the first 20+ orbits of a binary system of mass ratio q = 0.7 at the onset of dynamically unstable direct impact mass transfer. The mass transfer rate quickly exceeds the critical Eddington limit by many orders of magnitude, and thus we are unable to model a trans-Eddington phase. It appears that radiation pressure does not significantly effect the accretion flow in the highly super-Eddington regime. An optically thick common envelope forms around the binary within a few orbits. Although this envelope quickly exceeds the spatial domain of the computational grid, the fraction of the common envelope that exceeds zero gravitational binding energy is extremely small, suggesting that radiation-driven mass loss is insignificant in this regime. It remains to be seen whether simulations that capture the trans-Eddington phase of such flows will lead to the same conclusion or show that substantial material gets expelled.

We test a new "hybrid" scheme for simulating dynamical fluid flows in which cylindrical components of the momentum are advected across a rotating Cartesian coordinate mesh. This hybrid scheme allows us to conserve angular momentum to machine precision while capitalizing on the advantages offered by a Cartesian mesh, such as a straightforward implementation of mesh refinement. Our test focuses on measuring the real and imaginary parts of the eigenfrequency of unstable axisymmetric modes that naturally arise in massless polytropic tori having a range of different aspect ratios, and quantifying the uncertainty in these measurements. Our measured eigenfrequencies show good agreement with the results obtained from the linear stability analysis of Kojima (1986) and from nonlinear hydrodynamic simulations performed on a cylindrical coordinate mesh by Woodward et al. (1994). When compared against results conducted with a traditional Cartesian advection scheme, the hybrid scheme achieves qualitative convergence at the same or, in some cases, much lower grid resolutions and conserves angular momentum to a much higher degree of precision. As a result, this hybrid scheme is much better suited for simulating astrophysical fluid flows, such as accretion disks and mass-transferring binary systems.

A leading formation scenario for R Coronae Borealis (RCB) stars invokes the merger of degenerate He and CO white dwarfs (WD) in a binary. The observed ratio of 16O/18O for RCB stars is in the range of 0.3-20 much smaller than the solar value of ~500. In this paper, we investigate whether such a low ratio can be obtained in simulations of the merger of a CO and a He white dwarf. We present the results of five 3-dimensional hydrodynamic simulations of the merger of a double white dwarf system where the total mass is 0.9 Mdot and the initial mass ratio (q) varies between 0.5 and 0.99. We identify in simulations with $q\lesssim0.7$ a feature around the merged stars where the temperatures and densities are suitable for forming 18O. However, more 16O is being dredged-up from the C- and O-rich accretor during the merger than the amount of 18O that is produced. Therefore, on a dynamical time scale over which our hydrodynamics simulation runs, a 16O/18O ratio of ~2000 in the "best" case is found. If the conditions found in the hydrodynamic simulations persist for 10^6 seconds the oxygen ratio drops to 16 in one case studied, while in a hundred years it drops to ~4 in another case studied, consistent with the observed values in RCB stars. Therefore, the merger of two white dwarfs remains a strong candidate for the formation of these enigmatic stars.

We have developed a self-consistent-field technique similar to the one described by Hachisu, Eriguchi, & Nomoto (1986b) that can be used to construct detailed force-balanced models of synchronously rotating, double white dwarf (DWD) binaries that have a wide range of total masses, mass ratios, and separations. In addition to providing a computational tool that can be used to provide quiet initial starts for dynamical studies of the onset of mass transfer in DWD systems, we show that this SCF technique can be used to construct model sequences that mimic the last portion of the detached inspiral phase of DWD binary evolutions, and semi-detached model sequences that mimic a phase of conservative mass transfer.

We present preliminary results from recent high-resolution double-degenerate merger simulations with the Smooth Particle Hydrodynamics (SPH) technique. We put particular emphasis on verification and validation in our effort and show the importance of details in the initial condition setup for the final outcome of the simulation. We also stress the dynamical importance of including shocks in the simulations. These results represent a first step toward a suite of simulations that will shed light on the question whether double-degenerate mergers are a viable path toward type 1a supernovae. In future simulations, we will make use of the capabilities of the NuGrid collaboration in post-processing SPH particle trajectories with a complete nuclear network to follow the detailed nuclear reactions during the dynamic merger phase.

We investigate the influence of magnetic fields upon the dynamics of and resulting gravitational waves from a binary neutron star merger in full general relativity coupled to ideal magnetohydrodynamics (MHD). We consider two merger scenarios, one where the stars begin with initially aligned poloidal magnetic fields and one with no magnetic field. Both mergers result in a strongly differentially rotating object. In comparison to the non-magnetized scenario, the aligned magnetic fields delay the final merger of the two stars. During and after merger we observe phenomena driven by the magnetic field, including Kelvin-Helmholtz instabilities in shear layers, winding of the field lines, and transition from poloidal to toroidal fields. These effects not only produce electromagnetic radiation, but also can have a strong influence on the gravitational waves. Thus, there are promising prospects for studying such systems with both types of waves.

We model two mergers of orbiting binary neutron stars, the first forming a black hole and the second a differentially rotating neutron star. We extract gravitational waveforms in the wave zone. Comparisons to a post-Newtonian analysis allow us to compute the orbital kinematics, including trajectories and orbital eccentricities. We verify our code by evolving single stars and extracting radial perturbative modes, which compare very well to results from perturbation theory. The Einstein equations are solved in a first order reduction of the generalized harmonic formulation, and the fluid equations are solved using a modified convex essentially non-oscillatory method. All calculations are done in three spatial dimensions without symmetry assumptions. We use the \had computational infrastructure for distributed adaptive mesh refinement.

This paper has been removed by arXiv admin because it was an erroneous duplicate of astro-ph/9411031.

We present numerical simulations of dynamically unstable mass transfer in a double white dwarf binary with initial mass ratio, q = 0.4. The binary components are approximated as polytropes of index n = 3/2 and the initially synchronously rotating, semi-detached equilibrium binary is evolved hydrodynamically with the gravitational potential being computed through the solution of Poisson's equation. Upon initiating deep contact in our baseline simulation, the mass transfer rate grows by more than an order of magnitude over approximately ten orbits, as would be expected for dynamically unstable mass transfer. However, the mass transfer rate then reaches a peak value, the binary expands and the mass transfer event subsides. The binary must therefore have crossed the critical mass ratio for stability against dynamical mass transfer. Despite the initial loss of orbital angular momentum into the spin of the accreting star, we find that the accretor's spin saturates and angular momentum is returned to the orbit more efficiently than has been previously suspected for binaries in the direct impact accretion mode. To explore this surprising result, we directly measure the critical mass ratio for stability by imposing artificial angular momentum loss at various rates to drive the binary to an equilibrium mass transfer rate. For one of these driven evolutions, we attain equilibrium mass transfer and deduce that effectively q_crit has evolved to approximately 2/3. Despite the absence of a fully developed disk, tidal interactions appear effective in returning excess spin angular momentum to the orbit.

Using a three-dimensional nonlinear hydrodynamic code, we examine the dynamical stability of more than twenty self-gravitating, compressible, ellipsoidal fluid configurations that initially have the same velocity structure as Riemann S-type ellipsoids. Our focus is on ``adjoint'' configurations, in which internal fluid motions dominate over the collective spin of the ellipsoidal figure; Dedekind-like configurations are among this group. We find that, although some models are stable and some are moderately unstable, the majority are violently unstable toward the development of $m=1$, $m=3$, and higher-order azimuthal distortions that destroy the coherent, $m=2$ bar-like structure of the initial ellipsoidal configuration on a dynamical time scale. The parameter regime over which our models are found to be unstable generally corresponds with the regime over which incompressible Riemann S-type ellipsoids have been found to be susceptible to an elliptical strain instability \citep{LL96}. We therefore suspect that an elliptical instability is responsible for the destruction of our compressible analogs of Riemann ellipsoids. The existence of the elliptical instability raises concerns regarding the final fate of neutron stars that encounter the secular bar-mode instability and regarding the spectrum of gravitational waves that will be radiated from such systems.

Detached, inspiraling and semi-detached, mass-transferring double white dwarf (DWD) binary systems are both expected to be important sources for the proposed space-based gravitational-wave detector, LISA. The mass-radius relationship of individual white dwarf stars in combination with the constraints imposed by Roche geometries permit us to identify population boundaries for DWD systems in LISA's ``absolute'' amplitude-frequency diagram. With five key population boundaries in place, we are able to identify four principal population sub-domains, including one sub-domain that identifies where progenitors of Type Ia supernovae will reside. Given one full year of uninterrupted operation, LISA should be able to measure the rate at which the gravitational-wave frequency $f$ and, hence, the orbital period is changing in the highest frequency subpopulation of our Galaxy's DWD systems. We provide a formula by which the distance to each DWD system in this subpopulation can be determined; in addition, we show how the masses of the individual white dwarf stars in mass-transferring systems may be calculated.

A one-armed spiral instability has been found to develop in differentially rotating stellar models that have a relatively stiff, $n=1$ polytropic equation of state and a wide range of rotational energies. This suggests that such instabilities can arise in neutron stars that are differentially, although not necessarily rapidly, rotating. The instability seems to be directly triggered by the presence of a corotation resonance inside the star. Our analysis also suggests that a resonant cavity resulting from a local minimum in the radial vortensity profile of the star plays an important role in amplifying the unstable mode. Hence, it appears as through this instability is closely related to the so-called ``Rossby wave instability'' \citep{LLCN99} that has been found to arise in accretion disks. In addition to the one-armed ($m=1$) spiral mode, we have found that higher-order ($m = 2$ and $m=3$) nonaxisymmetric modes also can become unstable if corotation points that resonate with the eigenfrequencies of these higher-order modes also appear inside the star. The growth rate of each mode seems to depend on the location of its corotation radius with respect to the vortensity profile (or on the depth of its corotation radius inside the vortensity well). The existence of such instabilities makes the stability criterion for differentially rotating neutron stars non-unique. Also, the gravitational-waves emitted from such unstable systems generally will not have a monochromatic frequency spectrum.

Hydrodynamical simulations of semi-detached, polytropic binary stars are presented in an effort to study the onset and stability of dynamical mass transfer events. Initial, synchronously rotating equilibrium models are constructed using a self-consistent-field technique and then evolved with an Eulerian hydrodynamics code in a fully self-consistent manner. We describe code improvements introduced over the past few years that permit us to follow dynamical mass-transfer events through more than 30 orbits. Mass-transfer evolutions are presented for two different initial configurations: A dynamically unstable binary with initial mass ratio (donor/accretor) $q_0 = 1.3$ that leads to a complete merger in $\sim 10$ orbits; and a double-degenerate binary with initial mass ratio $q_0 = 0.5$ that, after some initial unstable growth of mass transfer, tends to separate as the mass-transfer rate levels off.

A three-dimensional, Newtonian hydrodynamic technique is used to follow the postbounce phase of a stellar core collapse event. For realistic initial data we have employed post core-bounce snapshots of the iron core of a 20 solar mass star. The models exhibit strong differential rotation but have centrally condensed density stratifications. We demonstrate for the first time that such postbounce cores are subject to a so-called low-T/|W| nonaxisymmetric instability and, in particular, can become dynamically unstable to an m=1 - dominated spiral mode at T/|W| ~ 0.08. We calculate the gravitational wave emission by the instability and find that the emitted waves may be detectable by current and future GW observatories from anywhere in the Milky Way.

We have modelled the nonlinear development of the secular bar-mode instability that is driven by gravitational radiation-reaction (GRR) forces in rotating neutron stars. In the absence of any competing viscous effects, an initially uniformly rotating, axisymmetric $n=1/2$ polytropic star with a ratio of rotational to gravitational potential energy $T/|W| = 0.181$ is driven by GRR forces to a bar-like structure, as predicted by linear theory. The pattern frequency of the bar slows to nearly zero, that is, the bar becomes almost stationary as viewed from an inertial frame of reference as GRR removes energy and angular momentum from the star. In this ``Dedekind-like'' state, rotational energy is stored as motion of the fluid in highly noncircular orbits inside the bar. However, in less than 10 dynamical times after its formation, the bar loses its initially coherent structure as the ordered flow inside the bar is disrupted by what appears to be a purely hydrodynamical, short-wavelength, ``shearing'' type instability. The gravitational waveforms generated by such an event are determined, and an estimate of the detectability of these waves is presented.

Nonlinear evolution of the gravitational radiation (GR) driven instability in the r-modes of neutron stars is studied by full numerical 3D hydrodynamical simulations. The growth of the r-mode instability is found to be limited by the formation of shocks and breaking waves when the dimensionless amplitude of the mode grows to about three in value. This maximum mode amplitude is shown by numerical tests to be rather insensitive to the strength of the GR driving force. Upper limits on the strengths of possible nonlinear mode--mode coupling are inferred. Previously unpublished details of the numerical techniques used are presented, and the results of numerous calibration runs are discussed.

We describe computational tools that have been developed to simulate dynamical mass transfer in semi-detached, polytropic binaries that are initially executing synchronous rotation upon circular orbits. Initial equilibrium models are generated with a self-consistent field algorithm; models are then evolved in time with a parallel, explicit, Eulerian hydrodynamics code with no assumptions made about the symmetry of the system. Poisson's equation is solved along with the equations of ideal fluid mechanics to allow us to treat the nonlinear tidal distortion of the components in a fully self-consistent manner. We present results from several standard numerical experiments that have been conducted to assess the general viability and validity of our tools, and from benchmark simulations that follow the evolution of two detached systems through five full orbits (up to approximately 90 stellar dynamical times). These benchmark runs allow us to gauge the level of quantitative accuracy with which simulations of semi-detached systems can be performed using presently available computing resources. We find that we should be able to resolve mass transfer at levels $\dot{M} / M > few x 10^-5$ per orbit through approximately 20 orbits with each orbit taking about 30 hours of computing time on parallel computing platforms.

The evolution of a neutron-star r-mode driven unstable by gravitational radiation (GR) is studied here using numerical solutions of the full non-linear fluid equations. The amplitude of the mode grows to order unity before strong shocks develop which quickly damp the mode. In this simulation the star loses about 40% of its initial angular momentum and 50% of its rotational kinetic energy before the mode is damped. The non-linear evolution causes the fluid to develop strong differential rotation which is concentrated near the surface and especially near the poles of the star.

In an effort to better understand the formation and evolution of barred galaxies, we have examined the properties of equatorial orbits in the effective potential of one model of a rapidly rotating, steady-state gas-dynamical bar that has been constructed via a self-consistent hydrodynamical simulation. Using a ``Restriction Hypothesis'' to determine initial conditions, we find that a significant fraction of orbits in this potential are quasi-ergodic and that regular orbits have a ``bowtie'' shape in contrast to the more typical x1 orbits. This bowtie orbit should give a boxy-peanut shape to such systems.

Compact astrophysical objects that rotate rapidly may encounter the dynamical ``bar instability.'' The bar-like deformation induced by this rotational instability causes the object to become a potentially strong source of gravitational radiation. We have carried out a set of long-duration simulations of the bar instability with two Eulerian hydrodynamics codes. Our results indicate that the remnant of this instability is a persistent bar-like structure that emits a long-lived gravitational radiation signal.

Warped H I gas layers in the outer regions of spiral galaxies usually display a noticeably twisted structure. This structure almost certainly arises primarily as a result of differential precession in the H I disk as it settles toward a preferred orientation in an underlying dark halo potential well that is not spherically symmetric. In an attempt to better understand the structure and evolution of these twisted, warped disk structures, we have adopted the "twist-equation" formalism originally developed by Petterson (1977) to study accretion onto compact objects. Utilizing more recent treatments of this formalism, we have generalized the twist-equation to allow for the treatment of non-Keplerian disks and from it have derived a steady-state structure of twisted disks that develops from free precession in a nonspherical, logarithmic halo potential. We have used this steady-state solution to produce H I maps of five galaxies (M83, NGC 300, NGC 2841, NGC 5033, NGC5055), which match the general features of the observed maps of these galaxies quite well. In addition, the model provides an avenue through which the kinematical viscosity of the H I disk and the quadrupole distortion of the dark halo in each galaxy can be quantified. This generalized equation can also be used to examine the time-evolutionary behavior of warped galaxy disks.

The orbital separation of compact binary stars will shrink with time due to the emission of gravitational radiation. This inspiralling phase of a binary system's evolution generally will be very long compared to the system's orbital period, but the final coalescence may be dynamical and driven to a large degree by hydrodynamic effects, particularly if there is a critical separation at which the system becomes dynamically unstable toward merger. Indeed, if weakly relativistic systems (such as white dwarf-white dwarf binaries) encounter a point of dynamical instability at some critically close separation, coalescence may be entirely a classical, hydrodynamic process. Therefore, a proper investigation of this stage of binary evolution must include three-dimensional hydrodynamic simulations. We have constructed equilibrium sequences of synchronously rotating, equal-mass binaries in circular orbit with a single parameter - the binary separation - varying along each sequence. Sequences have been constructed with various polytropic as well as realistic white dwarf and neutron star equations of state. Using a Newtonian, finite-difference hydrodynamics code, we have examined the dynamical stability of individual models along these equilibrium sequences. Our simulations indicate that no points of instability exist on the sequences we analyzed that had relatively soft equations of state (polytropic sequences with polytropic index $n=1.0$ and 1.5 and two white dwarf sequences). However, we did identify dynamically unstable binary models on sequences with stiffer equations of state ($n=0.5$ polytropic sequence and two neutron star sequences). We thus infer that binary systems with soft equations of state are not driven to merger by a dynamical instability.

We use a free-energy minimization approach to describe the secular and dynamical instabilities as well as the bifurcations along equilibrium sequences of rotating, self-gravitating fluid systems. Our approach is fully nonlinear and stems from the Landau-Ginzburg theory of phase transitions. Here we examine higher than 2nd-harmonic disturbances applied to Maclaurin spheroids, the corresponding bifurcating sequences, and their relation to nonlinear fission processes. The triangle and ammonite sequences bifurcate from the two 3rd-harmonic neutral points on the Maclaurin sequence while the square and one-ring sequences bifurcate from two of the three known 4th harmonic neutral points. In the other three cases, secular instability does not set in at the corresponding bifurcation points because the sequences stand and terminate at higher energies relative to the Maclaurin sequence. There is no known bifurcating sequence at the point of 3rd-harmonic dynamical instability. Our nonlinear approach easily identifies resonances between the Maclaurin sequence and various multi-fluid-body sequences that cannot be detected by linear stability analyses. Resonances appear as first-order phase transitions at points where the energies of the two sequences are nearly equal but the lower energy state belongs to one of the multi-fluid-body sequences.

In Papers I and II, we have used a free-energy minimization approach that stems from the Landau-Ginzburg theory of phase transitions to describe in simple and clear physical terms the secular and dynamical instabilities as well as the bifurcations along equilibrium sequences of rotating, self-gravitating fluid systems. Here we investigate the secular and dynamical 3rd-harmonic instabilities that appear first on the Jacobi sequence of incompressible zero- vorticity ellipsoids. Poincare found a bifurcation point on the Jacobi sequence where a 3rd-harmonic mode becomes neutral. A sequence of pear-shaped equilibria branches off at this point but stands at higher energies. Therefore, the Jacobi ellipsoids remain secularly and dynamically stable. Cartan found that dynamical 3rd-harmonic instability also sets in at the Jacobi-pear bifurcation. We find that Cartan's instability leads to differentially rotating objects and not to uniformly rotating pear-shaped equilibria. We demonstrate that the pear-shaped sequence exists at higher energies and at higher rotation relative to the Jacobi sequence. The Jacobi ellipsoid undergoes fission on a secular time scale and a short-period binary is produced. The classical fission hypothesis of binary-star formation of Poincare and Darwin is thus feasible.

We use a free-energy minimization approach to describe the secular and dynamical instabilities as well as the bifurcations along equilibrium sequences of rotating, self-gravitating fluid systems. Our approach is fully nonlinear and stems from the Ginzburg-Landau theory of phase transitions. In this paper, we examine fourth-harmonic axisymmetric disturbances in Maclaurin spheroids and fourth-harmonic nonaxisymmetric disturbances in Jacobi ellipsoids. These two cases are very similar in the framework of phase transitions. Irrespective of whether a nonlinear first-order phase transition occurs between the critical point and the higher turning point or an apparent second-order phase transition occurs beyond the higher turning point, the result is fission (i.e. ``spontaneous breaking'' of the topology) of the original object on a secular time scale: the Maclaurin spheroid becomes a uniformly rotating axisymmetric torus and the Jacobi ellipsoid becomes a binary. The presence of viscosity is crucial since angular momentum needs to be redistributed for uniform rotation to be maintained. The phase transitions of the dynamical systems are briefly discussed in relation to previous numerical simulations of the formation and evolution of protostellar systems.

Laboratory searches for the detection of gravitational waves have focused on the detection of burst signals emitted during a supernova explosion, but have not resulted in any confirmed detections. An alternative approach has been to search for continuous wave (CW) gravitational radiation from the Crab pulsar. In this paper, we examine the possibility of detecting CW gravitational radiation from pulsars and show that nearby millisecond pulsars are generally much better candidates. We show that the minimum strain h_c ~ 10E-26 that can be detected by tuning an antenna to the frequency of the milli- second pulsar PSR 1957+20, with presently available detector technology, is orders of magnitude better than what has been accomplished so far by observing the Crab pulsar, and within an order of magnitude of the maximum strain that may be produced by it. In addition, we point out that there is likely to be a population of rapidly rotating neutron stars (not necessarily radio pulsars) in the solar neighborhood whose spindown evolution is driven by gravitational radiation. We argue that the projected sensitivity of modern resonant detectors is sufficient to detect the subset of this population that lies within 0.1 kpc of the sun.

We have previously introduced the parameter `alpha' as an indicator of stability to m=2 nonaxisymmetric modes in rotating, self-gravitating, axisymmetric, gaseous and stellar systems. This parameter can be written as a function of the total rotational kinetic energy, the total gravitational potential energy, and as a function of the topology/connectedness and the geometric shape of a system. Here we extend the stability criterion to nonaxisymmetric equilibrium systems, such as ellipsoids and elliptical disks and cylinders. We test the validity of this extension by considering predictions for previously published, gaseous and stellar, nonaxisymmetric models. The above formulation and critical values account accurately for the stability properties of m=2 modes in gaseous Riemann S-type ellipsoids (including the Jacobi and Dedekind ellipsoids) and elliptical Riemann disks as well as in stellar elliptical Freeman disks and cylinders: all these systems are dynamically stable except for the stellar elliptical Freeman disks that exhibit a relatively small region of m=2 dynamical instability.

We analyze previous results on the stability of uniformly and differentially rotating, self-gravitating, gaseous and stellar, axisymmetric systems to derive a new stability criterion for the appearance of toroidal, m=2 Intermediate (I) and bar modes. In the process, we demonstrate that the bar modes in stellar systems and the m=2 I-modes in gaseous systems have many common physical characteristics and only one substantial difference: because of the anisotropy of the stress tensor, dynamical instability sets in at lower rotation in stellar systems. This difference is reflected also in the new stability criterion. The new stability parameter "alpha" is formulated first for uniformly rotating systems and is based on the angular momentum content rather than on the energy content of a system. For stability of stellar systems "alpha" = 0.254-0.258, while it is in the range of 0.341-0.354 for gaseous systems. For uniform rotation, one can write "alpha" as a function of the total (rotational) kinetic and gravitational energies, and of a parameter which is characteristic of the topology/connectedness and the geometric shape of a system. "Alpha" can be extended to and calculated for a variety of differentially rotating, gaseous and stellar, axisymmetric disk and spheroidal models whose equilibrium structures and stability characteristics are known. We also estimate "alpha" for gaseous toroidal models and for stellar disk systems embedded in an inert or responsive "halo". We find that the new stability criterion holds equally well for all these previously published axisymmetric models.

A free-energy minimization approach is used to address the secular & dynamical instabilities & the bifurcations along sequences of rotating, self-gravitating fluid and stellar systems. Our approach stems from the Landau-Ginzburg theory of phase transitions. We focus on the Maclaurin sequence of oblate spheroids & on the effects of second-harmonic disturbances. Second-order phase transitions appear on the Maclaurin sequence also at the points of dynamical instability & of bifurcation of the Dedekind sequence. Distinguishing characteristic of each second-order phase transition is the (non)conservation of an integral of motion (e.g. circulation) which determines the appearance of the transition. Circulation is not conserved in stellar systems because the stress-tensor gradient terms that appear in the Jeans equations of motion include viscosity-like off-diagonal terms of the same order of magnitude as the conventional pressure gradient terms. This explains why the Jacobi bifurcation is a point of dynamical instability in stellar systems but only a point of secular instability in viscous fluids. The second-order phase transitions are discussed in relation to the dynamical instability of stellar systems, the lambda-transition of liquid He-4, the second-order phase transition in superconductivity & the mechanism of spontaneous symmetry breaking.