We present simulations of bar-unstable stellar discs in which the bars thicken into box/peanut shapes. Detailed analysis of the evolution of each model revealed three different mechanisms for thickening the bars. The first mechanism is the well-known buckling instability, the second is the vertical excitation of bar orbits by passage through the 2:1 vertical resonance, and the third is a gradually increasing fraction of bar orbits trapped into this resonance. Since bars in many galaxies may have formed and thickened long ago, we have examined the models for fossil evidence in the velocity distribution of stars in the bar, finding a diagnostic to discriminate between a bar that had buckled from the other two mechanisms.

An earlier paper presented the potentially significant discovery that disturbances in simplified simulations of a stellar disc model that was predicted to be stable in linear theory grew to large amplitude over a long period of time. The ultimate appearance of true instabilities was attributed to non-linear scattering by a succession of collective waves excited by shot noise from the finite number of particles. The paper concluded that no finite number of particles, however large, could mimic a smooth disc. As this surprising finding has been challenged as an artifact of the numerical scheme employed, we here present a new calculation of the same model using a different grid geometry that confirms the original behaviour.

We argue that self-excited instabilities are the cause of spiral patterns in simulations of unperturbed stellar discs. In previous papers, we have found that spiral patterns were caused by a few concurrent waves, which we claimed were modes. The superposition of a few steadily rotating waves inevitably causes the appearance of the disc to change continuously, and creates the kind of shearing spiral patterns that have been widely reported. Although we have found that individual modes last for relatively few rotations, spiral activity persists because fresh instabilities appear, which we suspected were excited by the changes to the disc caused by previous disturbances. Here we confirm our suspicion by demonstrating that scattering at either of the Lindblad resonances seeds a new groove-type instability. With this logical gap closed, our understanding of the behaviour in the simulations is almost complete. We believe that our robust mechanism is a major cause of spiral patterns in the old stellar discs of galaxies, including the Milky Way where we have previously reported evidence for resonance scattering in the recently released Gaia data.

The inner disc of the local group galaxy M33 appears to be in settled rotational balance, and near IR images reveal a mild, large-scale, two-arm spiral pattern with no strong bar. We have constructed N-body models that match all the extensive observational data on the kinematics and surface density of stars and gas in the inner part of M33. We find that currently favoured models are unstable to the formation of a strong bar of semi-major axis 2 < a_B < 3 kpc within 1 Gyr, which changes the dynamical properties of the models to become inconsistent with the current, apparently well-settled, state. The formation of a bar is unaffected by how the gas component is modelled, by increasing the mass of the nuclear star cluster, or by making the dark matter halo counter-rotate, but it can be prevented by either reducing the mass-to-light ratio of the stars to Upsilon_V ~ 0.6 or Upsilon_K ~ 0.23 in solar units or by increasing the random motions of the stars. Also a shorter and weaker bar results when the halo is rigid and unresponsive. However, all three near-stable models support multi-arm spirals, and not the observed large-scale bi-symmetric spiral. A two-arm spiral pattern could perhaps be tidally induced, but such a model would require an implausibly low mass disc to avoid a bar and there is no visible culprit. Thus the survival of the current state of this exceptionally well-studied galaxy is not yet understood. We also suspect that many other unbarred galaxies present a similar puzzle.

We compare the distribution in position and velocity of nearby stars from the Gaia DR2 radial velocity sample with predictions of current theories for spirals in disc galaxies. Although the rich substructure in velocity space contains the same information, we find it more revealing to reproject the data into action-angle variables, and we describe why resonant scattering would be more readily identifiable in these variables. We compute the predicted changes to the phase space density, in multiple different projections, that would be caused by a simplified isolated spiral pattern, finding widely differing predictions from each theory. We conclude that the phase space structure present in the Gaia data shares many of the qualitative features expected in the transient spiral mode model. We argue that the popular picture of apparently swing-amplified spirals results from the superposition of a few underlying spiral modes.

The distributions of stars, gas, and dark matter in disk galaxies provide important constraints on galaxy formation models, particularly on small spatial scales (<1kpc). We have designed the RSS Imaging spectroscopy Nearby Galaxy Survey (RINGS) to target a sample of 19 nearby spiral galaxies. For each of these galaxies, we are obtaining and modeling Halpha and H1~21~cm spectroscopic data as well as multi-band photometric data. We intend to use these models to explore the underlying structure and evolution of these galaxies in a cosmological context, as well as whether the predictions of LCDM are consistent with the mass distributions of these galaxies. In this paper, we present spectroscopic imaging data for 14 of the RINGS galaxies observed with the medium spectral resolution Fabby-Perot etalon on the Southern African Large Telescope. From these observations, we derive high spatial resolution line of sight velocity fields of the Halpha line of excited hydrogen, as well as maps and azimuthally averaged profiles of the integrated Halpha and [NII] emission and oxygen abundances. We then model these kinematic maps with axisymmetric models, from which we extract rotation curves and projection geometries for these galaxies. We show that our derived rotation curves agree well with other determinations and the similarity of the projection angles with those derived from our photometric images argues against these galaxies having intrinsically oval disks.

We present high-quality fluid dynamical simulations of isothermal gas flows in a rotating barred potential. We show that a large quantity of gas is driven right into the nucleus of a model galaxy when the potential lacks a central mass concentration, but the inflow stalls at a nuclear ring in comparison simulations that include a central massive object. The radius of the nuclear gas ring increases linearly with the mass of the central object. We argue that bars drive gas right into the nucleus in the early stages of disk galaxy formation, where a nuclear star cluster and perhaps a massive black hole could be created. The process is self-limiting, however, because inflow stalls at a nuclear ring once the mass of gas and stars in the nucleus exceeds ~1% of the disk mass, which shuts off rapid growth of the black hole. We briefly discuss the relevance of these results to the seeding of massive black holes in galaxies, the merger model for quasar evolution, and the existence of massive black holes in disk galaxies that lack a significant classical bulge.

MHONGOOSE is a deep survey of the neutral hydrogen distribution in a representative sample of 30 nearby disk and dwarf galaxies with HI masses from 10^6 to ~10^{11} M_sun, and luminosities from M_R ~ -12 to M_R ~ -22. The sample is selected to uniformly cover the available range in log(M_HI). Our extremely deep observations, down to HI column density limits of well below 10^{18} cm^{-2} - or a few hundred times fainter than the typical HI disks in galaxies - will directly detect the effects of cold accretion from the intergalactic medium and the links with the cosmic web. These observations will be the first ever to probe the very low-column density neutral gas in galaxies at these high resolutions. Combination with data at other wavelengths, most of it already available, will enable accurate modelling of the properties and evolution of the mass components in these galaxies and link these with the effects of environment, dark matter distribution, and other fundamental properties such as halo mass and angular momentum. MHONGOOSE can already start addressing some of the SKA-1 science goals and will provide a comprehensive inventory of the processes driving the transformation and evolution of galaxies in the nearby universe at high resolution and over 5 orders of magnitude in column density. It will be a Nearby Galaxies Legacy Survey that will be unsurpassed until the advent of the SKA, and can serve as a highly visible, lasting statement of MeerKAT's capabilities.

We report idealized simulations that mimic the growth of galaxy disks embedded in responsive halos and bulges. The disks manifested an almost overwhelming tendency to form strong bars that we found very difficult to prevent. We found that fresh bars formed in growing disks after we had destroyed the original, indicating that bar formation also afflicts continued galaxy evolution, and not just the early stages of disk formation. This behavior raises still more insistently the previously unsolved question of how some galaxies avoid bars. Since our simulations included only collisionless star and halo particles, our findings may apply to gas-poor galaxies only; however the conundrum persists for the substantial unbarred fraction of those galaxies. Our original objective was to study how internal dynamics rearranged the distribution of mass in the disk as a generalization of our earlier study with rigid spherical components. With difficulty, we were able to construct some models that were not strongly influenced by bars, and found that halo compression, and angular momentum exchange with the disk did not alter our earlier conclusion that spiral activity is largely responsible for creating smooth density profiles and rotation curves.

We show that the exponential growth rate of a bar in a stellar disk is substantially greater when the disk is embedded in a live halo than in a rigid one having the same mass distribution. We also find that the vigor of the instability in disk-halo systems varies with the shape of the halo velocity ellipsoid. Disks in rigid halos that are massive enough to be stable by the usual criteria, quickly form bars in isotropic halos and much greater halo mass is needed to avoid a strong bar; thus stability criteria derived for disks in rigid halos do not apply when the halo is responsive. The study presented here is of an idealized family of models with near uniform central rotation and that lack an extended halo; we present more realistic models with extended halos in a companion paper. The puzzle presented by the absence of strong bars in some galaxies having gently rising inner rotation curves is compounded by the results presented here.

We have obtained high-spatial-resolution spectrophotometric data on several nearby spiral galaxies with the Southern African Large Telescope (SALT) Fabry-P\'erot interferometer on the Robert Stobie Spectrograph (RSS) as a part of the RSS Imaging spectroscopy Nearby Galaxy Survey (RINGS). We have successfully reduced two tracks of Fabry-P\'erot data for the galaxy NGC 2280 to produce a velocity field of the H-alpha line of excited hydrogen. We have modeled these data with the DiskFit modeling software and found these models to be in excellent agreement both with previous measurements in the literature and with our lower-resolution HI velocity field of the same galaxy. Despite this good agreement, small regions exist where the difference between the H-alpha and HI velocities is larger than would be expected from typical dispersions. We investigate these regions of high velocity difference and offer possible explanations for their existence.

This posting announces public availability of version 1.2 of the DiskFit software package developed by the authors, which may be used to fit simple non-axisymmetric models either to images or to velocity fields of disk galaxies. Here we give an outline of the capability of the code and provide the link to downloading executables, the source code, and a comprehensive on-line manual. We argue that in important respects the code is superior to rotcur for fitting kinematic maps and to galfit for fitting multi-component models to photometric images.

I present empirical measurements of the rate of relaxation in N-body simulations of stable spherical systems and distinguish two separate types of relaxation: energy diffusion that is largely independent of particle mass, and energy exchange between particles of differing masses. While diffusion is generally regarded as a Fokker-Planck process, it can equivalently be viewed as the consequence of collective oscillations that are driven by shot noise. Empirical diffusion rates scale as N^{-1} in inhomogeneous models, in agreement with Fokker-Planck predictions, but collective effects cause relaxation to scale more nearly as N^{-1/2} in the special case of a uniform sphere. I use four different methods to compute the gravitational field, and a 100-fold range in the numbers of particles in each case. I find the rate at which energy is exchanged between particles of differing masses does not depend at all on the force determination method, but I do find the energy diffusion rate is marginally lower when a field method is used. The relaxation rate in 3D is virtually independent of the method used because it is dominated by distant encounters; any method to estimate the gravitational field that correctly captures the contributions from distant particles must also capture their statistical fluctuations and the collective modes they drive.

This is a brief rebuttal to arXiv:1502.03821, which claims to provide the first observational proof of dark matter interior to the solar circle. We point out that this result is not new, and can be traced back at least a quarter century.

We use the Made-to-Measure method to construct N-body realizations of self-similar, triaxial ellipsoidal halos having cosmologically realistic density profiles. Our implementation parallels previous work with a few numerical refinements, but we show that orbital averaging is an intrinsic feature of the force of change equation and argue that additional averaging or smoothing schemes are redundant. We present models having the Einasto radial mass profile that range from prolate to strongly triaxial. We use a least-squares polynomial fit to the expansion coefficients to obtain an analytical representation of the particle density from which we derive density contours and eccentricity profiles more efficiently than by the usual particle smoothing techniques. We show that our N-body realizations both retain their shape in unconstrained evolution and recover it after large amplitude perturbations.

We show that spiral activity can erase pronounced features in disk galaxy rotation curves. We present simulations of growing disks, in which the added material has a physically motivated distribution, as well as other examples of physically less realistic accretion. In all cases, attempts to create unrealistic rotation curves were unsuccessful because spiral activity rapidly smoothed away features in the disk mass profile. The added material was redistributed radially by the spiral activity, which was itself provoked by the density feature. In the case of a ridge-like feature in the surface density profile, we show that two unstable spiral modes develop, and the associated angular momentum changes in horseshoe orbits remove particles from the ridge and spread them both inwards and outwards. This process rapidly erases the density feature from the disk. We also find that the lack of a feature when transitioning from disk to halo dominance in the rotation curves of disk galaxies, the so called "disk-halo conspiracy", could also be accounted for by this mechanism. We do not create perfectly exponential mass profiles in the disk, but suggest that this mechanism contributes to their creation.

Precise measurements of gas kinematics in the disk of a spiral galaxy can be used to estimate its mass distribution. The Southern African Large Telescope (SALT) has a large collecting area and field of view, and is equipped with a Fabry-Perot interferometer that can measure gas kinematics in a galaxy from the Halpha line. To take advantage of this capability, we have constructed a sample of 19 nearby spiral galaxies, the RSS Imaging and Spectroscopy Nearby Galaxy Survey (RINGS), as targets for detailed study of their mass distributions and have collected much of the needed data. In this paper, we present velocity maps produced from Halpha Fabry-Perot interferometry and HI aperture synthesis for one of these galaxies, NGC 2280, and show that the two velocity measurements are generally in excellent agreement. Minor differences can mostly be attributed to the different spatial distributions of the excited and neutral gas in this galaxy, but we do detect some anomalous velocities in our Halpha velocity map of the kind that have previously been detected in other galaxies. Models produced from our two velocity maps agree well with each other and our estimates of the systemic velocity and projection angles confirm previous measurements of these quantities for NGC 2280.

Athanassoula (2014) has claimed that measurements of the ratio of corotation radius to bar length in galaxies do not place a constraint on the disk to halo mass ratio. Her conclusion was based on a series of simulations published by Athanassoula et al. (2013). Here we show that these results are, in fact, consistent with previous work on the slow down of bars due to dynamical friction because gas inflow rearranges the disk mass before the bar settles. It therefore remains true that a bar rotating sufficiently fast that corotation is not far beyond the bar end requires a near maximum disk.

This posting announces public availability of the GALAXY software package developed by the author over the past 40 years. It is a highly efficient code for the evolution of (almost) isolated, collisionless stellar systems, both disk-like and ellipsoidal. In addition to the N-body code galaxy, which offers eleven different methods to compute the gravitational accelerations, the package also includes sophisticated set-up and analysis software. This paper gives an outline of the contents of the package and provides links to the source code and a comprehensive on-line manual. While not as versatile as tree codes, the particle-mesh methods in this package are shown, for certain restricted applications, to be between 50 and 200 times faster than a widely-used tree code.

We present evidence that recurrent spiral activity, long manifested in simulations of disk galaxies, results from the super-position of a few transient spiral modes. Each mode lasts between 5 and 10 rotations at its corotation radius where its amplitude is greatest. The scattering of stars as each wave decays takes place over narrow ranges of angular momentum, causing abrupt changes to the impedance of the disk to subsequent traveling waves. Partial reflections of waves at these newly created features, allows new standing-wave instabilities to appear that saturate and decay in their turn, scattering particles at new locations, creating a recurring cycle. The spiral activity causes the general level of random motion to rise, gradually decreasing the ability of the disk to support further activity unless the disk contains a dissipative gas component from which stars form on near-circular orbits. We also show that this interpretation is consistent with the behavior reported in other recent simulations with low mass-disks.

We utilize Young's algorithm to model the adiabatic compression of the dark matter haloes of galaxies in the THINGS survey to determine the relationship between the halo fit to the rotation curve and the corresponding primordial halo prior to compression. Young's algorithm conserves radial action and angular momentum, resulting in less halo compression than more widely utilized approximations. We find that estimates of the parameters of the NFW haloes fit to the current dark matter distribution systematically overestimate the concentration and underestimate the virial velocity of the corresponding primordial halo. It is the latter that is predicted by dark matter simulations; so accounting for compression is a necessary step for evaluating whether massive galaxies are consistent with dark matter-only simulations. The inferred primordial haloes broadly follow the c-V200 relation expected in a LCDM cosmogony, but often scatter to lower concentrations. We are unable to obtain fits at all for those galaxies whose current dark matter haloes are poorly described by the NFW form. We thus find a mixed bag: some galaxies are reasonably well described by adiabatic compression within a primordial NFW halo, while others require additional mechanisms that reduce the density of dark matter below the primordial initial condition.

Disk galaxies evolve over time through processes that may rearrange both the radial mass profile and the metallicity distribution within the disk. This review of such slow changes is largely, though not entirely, restricted to internally-driven processes that can be distinguished from evolution driven by galaxy interactions. It both describes our current understanding of disk evolution, and identifies areas where more work is needed. Stellar disks are heated through spiral scattering, which increases random motion components in the plane, while molecular clouds redirect some fraction of the random energy into vertical motion. The recently discovered process of radial migration at the corotation resonance of a transient spiral mode does not alter the underlying structure of the disk, since it neither heats the disk nor causes it to spread, but it does have a profound effect on the expected distribution of metallicities among the disk stars. Bars in disks are believed to be major drivers of secular evolution through interactions with the outer disk and with the halo. Once the material that makes up galaxy disks is converted into stars, their overall angular momentum distribution cannot change by much, but that of the gas is generally far more liable to rearrangement, allowing rings and pseudo-bulges to form. While simulations are powerful tools from which we have learned a great deal, those of disks may suffer from collisional relaxation that requires some results to be interpreted with caution.

I use N-body simulations with two mass species of particles to demonstrate that disk galaxy simulations are subject to collisional relaxation at a higher rate than is widely assumed. Relaxation affects the vertical thickness of the disk most strongly, and drives the velocity ellipsoid to a moderately flattened shape similar to that observed for disk stars in the solar neighborhood. The velocity ellipsoid in simulations with small numbers of particles quickly approaches this shape, but shot noise also dominates the in-plane behavior. Simulations with higher, but reachable, numbers of particles relax slowly enough to be considered collisionless, allowing the in-plane dispersions to rise due to spiral activity without heating the vertical motions. Relaxation may have affected many previously published simulations of the formation and evolution of galaxy disks.

The isolation, simple apparent structure, and low luminosity of the nearby spiral galaxy NGC 6503 make it an ideal candidate for an in-depth kinematic and photometric study. We introduce a new publicly available code, DiskFit, that implements procedures for fitting non-axisymmetries in either kinematic or photometric data. We use DiskFit to analyze new Halpha and CO velocity field data as well as HI kinematics from Greisen et al. to search for non-circular motions in the disc of NGC 6503. We find NGC 6503 to have remarkably regular gas kinematics that are well-described by rotation only. We also use DiskFit and a new Ks-band image of NGC 6503 to constrain photometric models of the disc, bar and bulge. We find the galaxy to be photometrically dominated by the disc. We find NGC 6503 to contain a faint bar and an exponential bulge which together contribute only ~5% of the total galaxy light. The combination of our kinematic and photometric DiskFit models suggest that NGC 6503 contains a weak, end-on bar that may have produced its Type II surface brightness profile but is unlikely to be responsible for its strong sigma-drop.

We describe the current plans for a spectroscopic survey of millions of stars in the Milky Way galaxy using the Guo Shou Jing Telescope (GSJT, formerly the Large Area Multi-Object Spectroscopic Telescope - LAMOST). The survey will obtain spectra for 2.5 million stars brighter than $r<19$ during dark/grey time, and 5 million stars brighter than $r<17$ or $J<16$ on nights that are moonlit or have low transparency. The survey will begin in fall of 2012, and will run for at least four years. The telescope design constrains the optimal declination range for observations to $10^\circ<\delta<50^\circ$, and site conditions lead to an emphasis on stars in the direction of the Galactic anticenter. The survey is divided into three parts with different target selection strategies: disk, anticenter, and spheroid. The resulting dataset will be used to study the merger history of the Milky Way, the substructure and evolution of the disks, the nature of the first generation of stars through identification of the lowest metallicity stars, and star formation through study of open clusters and the OB associations. Detailed design of the LEGUE survey will be completed after a review of the results of the pilot survey in summer 2012.

We describe a general target selection algorithm that is applicable to any survey in which the number of available candidates is much larger than the number of objects to be observed. This routine aims to achieve a balance between a smoothly-varying, well-understood selection function and the desire to preferentially select certain types of targets. Some target-selection examples are shown that illustrate different possibilities of emphasis functions. Although it is generally applicable, the algorithm was developed specifically for the LAMOST Experiment for Galactic Understanding and Exploration (LEGUE) survey that will be carried out using the Chinese Guo Shou Jing Telescope. In particular, this algorithm was designed for the portion of LEGUE targeting the Galactic halo, in which we attempt to balance a variety of science goals that require stars at fainter magnitudes than can be completely sampled by LAMOST. This algorithm has been implemented for the halo portion of the LAMOST pilot survey, which began in October 2011.

We use samples of local main-sequence stars to show that the radial gradient of [Fe/H] in the thin disk of the Milky Way decreases with mean effective stellar temperature. Many of these stars are visiting the solar neighborhood from the inner and outer Galaxy. We use the angular momentum of each star about the Galactic center to determine the guiding center radius and to eliminate the effects of epicyclic motion, which would otherwise blur the estimated gradients. We interpret the effective temperature as a proxy for mean age, and conclude that the decreasing gradient is consistent with the predictions of radial mixing due to transient spiral patterns. We find some evidence that the trend of decreasing gradient with increasing mean age breaks to a constant gradient for samples of stars whose main-sequence life-times exceed the likely age of the thin disk.

The origin of spiral patterns in galaxies is still not fully understood. Similar features also develop readily in N-body simulations of isolated cool, collisionless disks, yet even here the mechanism has yet to be explained. In this series of papers, I present a detailed study of the origin of spiral activity in simulations in the hope that the mechanism that causes the patterns is also responsible for some of these features galaxies. In this first paper, I use a suite of highly idealized simulations of a linearly stable disk that employ increasing numbers of particles. While the amplitudes of initial non-axisymmetric features scale as the inverse square-root of the number of particles employed, the final amplitude of the patterns is independent of the particle number. I find that the amplitudes of non-axisymmetric disturbances grow in two distinct phases: slow growth occurs when the relative overdensity is below ~2%, but above this level the amplitude rises more rapidly. I show that all features, even of very low amplitude, scatter particles at the inner Lindblad resonance, changing the distribution of particles in the disk in such a way as to foster continued growth. Stronger scattering by larger amplitude waves provokes a vigorous instability that is a true linear mode of the modified disk.

We present a study of the extent to which the Sellwood & Binney radial migration of stars is affected by their vertical motion about the midplane. We use both controlled simulations in which only a single spiral mode is excited, as well as slightly more realistic cases with multiple spiral patterns and a bar. We find that rms angular momentum changes are reduced by vertical motion, but rather gradually, and the maximum changes are almost as large for thick disc stars as for those in a thin disc. We find that particles in simulations in which a bar forms suffer slightly larger angular momentum changes than in comparable cases with no bar, but the cumulative effect of multiple spiral events still dominates. We have determined that vertical action, and not vertical energy, is conserved on average during radial migration.

This chapter reviews theoretical work on the stellar dynamics of galaxy disks. All the known collective global instabilities are identified, and their mechanisms described in terms of local wave mechanics. A detailed discussion of warps and other bending waves is also given. The structure of bars in galaxies, and their effect on galaxy evolution, is now reasonably well understood, but there is still no convincing explanation for their origin and frequency. Spiral patterns have long presented a special challenge, and ideas and recent developments are reviewed. Other topics include scattering of disk stars and the survival of thin disks.

We extract a sample of disc stars within 200 pc of the Sun from the RAVE and SDSS surveys. Distances are estimated photometrically and proper motions are from ground-based data. We show that the velocity-space substructure first revealed in the Geneva-Copenhagen sample is also present in this completely independent sample. We also evaluate action-angle variables for these stars and show that the Hyades stream stars in these data are again characteristic of having been scattered at a Lindblad resonance. Unfortunately, analysis of such local samples can determine neither whether it is an inner or an outer Lindblad resonance, nor the multiplicity of the pattern.

N-body simulations of disc galaxies that display recurrent transient spiral patterns are comparatively easy to construct, but are harder to understand. In this paper, I summarise the evidence from such experiments that the spiral patterns result from a recurrent spiral instability cycle. Each wave starts as rapidly growing, small-amplitude instability caused by a deficiency of particles at a particular angular momentum. The resulting large-amplitude wave creates, through resonant scattering, the conditions needed to precipitate a new instability.

The rate of internally-driven evolution of galaxy discs is strongly affected by the lifetimes of the spiral patterns they support. Evolution is much faster if the spiral patterns are recurrent short-lived transients rather than long-lived, quasi-steady features. As rival theories are still advocated based on these two distinct hypotheses, I review the evidence that bears on the question of the lifetimes of spiral patterns in galaxies. Observational evidence from external galaxies is frustratingly inconclusive, but the velocity distribution in the solar neighbourhood is more consistent with the transient picture. I present simulations of galaxy models that have been proposed to support quasi-steady, two-arm spiral modes that in fact evolve quickly due to multi-arm instabilities. I also show that all simulations to date manifest short-lived patterns, despite claims to the contrary. Thus the transient hypothesis is favoured by both numerical results and the velocity distribution in the solar neighbourhood.

Stars in the solar neighbourhood do not have a smooth distribution of velocities. Instead, the distribution of velocity components in the Galactic plane manifests a great deal of kinematic substructure. Here I present an analysis in action-angle variables of the Geneva-Copenhagen survey of ~14,000 nearby F & G dwarfs with distances and full space motions. I show that stars in the so-called "Hyades stream" have both angle and action variables characteristic of their having been scattered at an inner Lindblad resonance of a rotating disturbance potential. This discovery seems to favour spiral patterns as recurrent, short-lived instabilities.

After a short review of the principal theories of spiral structure in galaxies, I describe two new developments. First, it now seems clear that linear theory cannot yield a full description for the development of spiral patterns because N-body simulations suggest that non-linear effects manifest themselves at a relative overdensity of ~2%, which is well below the believed spiral amplitudes in galaxies. Second, I summarize the evidence that some stars in the solar neighborhood have been scattered at an inner Lindblad resonance. This evidence strongly supports a picture of spirals as a recurring cycle of transient instabilities, each caused by resonant scattering by a previous wave.

High-quality velocity maps of galaxies frequently exhibit signatures of non-circular streaming motions. We here apply the software tool, "velfit" recently proposed by Spekkens & Sellwood, to five representative galaxies from the THINGS sample. We describe the strengths and weaknesses of the tool, and show that it is both more powerful and yields results that are more easily interpreted than the commonly used procedure. We demonstrate that it can estimate the magnitudes of forced non-circular motions over a broad range of bar strengths from a strongly barred galaxy, through cases of mild bar-like distortions to placing bounds on the shapes of halos in galaxies having extended rotation curves. We identify mild oval distortions in the inner parts of two dwarf galaxies, NGC 2976 and NGC 7793, and show that the true strength of the non-axisymmetric gas flow in the strongly barred galaxy NGC 2903 is revealed more clearly in our fit to an optical Halpha map than to the neutral hydrogen data. The method can also yield a direct estimate of the ellipticity of a slowly-rotating potential distortion in the flat part of a rotation curve, and we use our results to place tight bounds on the possible ellipticity of the outer halos of NGC 3198 and NGC 2403.

We demonstrate that the chaotic nature of N-body systems can lead to macroscopic variations in the evolution of collisionless simulations containing rotationally supported discs. The unavoidable stochasticity that afflicts all simulations generally causes mild differences between the evolution of similar models but, in order to illustrate that this is not always true, we present a case that shows extreme bimodal divergence. The divergent behaviour occurs in two different types of code and is independent of all numerical parameters. We identify and give explicit illustrations of several sources of stochasticity, and also show that macroscopic variations in the evolution can originate from differences at the round-off error level. We obtain somewhat more consistent results from simulations in which the halo is set up with great care compared with those started from more approximate equilibria, but we have been unable to eliminate diverging behaviour entirely because the main sources of stochasticity are intrinsic to the disc. We show that the divergence is only temporary and that halo friction is merely delayed, for a substantial time in some cases. We argue that the delays are unlikely to arise in real galaxies, and that our results do not affect dynamical friction constraints on halo density. Stochastic variations in the evolution are inevitable in all simulations of disc-halo systems, irrespective of how they were created, although their effect is generally far less extreme than we find here. The possibility of divergent behaviour complicates comparison of results from different workers.

The first phase of stellar evolution in the history of the universe may be Dark Stars, powered by dark matter heating rather than by fusion. Weakly interacting massive particles, which are their own antiparticles, can annihilate and provide an important heat source for the first stars in the universe. This and the previous contribution present the story of Dark Stars. In this second part, we describe the structure of Dark Stars and predict that they are very massive ($\sim 800 M_\odot$), cool (6000 K), bright ($\sim 10^6 L_\odot$), long-lived ($\sim 10^6$ years), and probable precursors to (otherwise unexplained) supermassive black holes. Later, once the initial dark matter fuel runs out and fusion sets in, dark matter annihilation can predominate again if the scattering cross section is strong enough, so that a Dark Star is born again.

The first phase of stellar evolution in the history of the universe may be Dark Stars, powered by dark matter heating rather than by fusion. Weakly interacting massive particles, which are their own antiparticles, can annihilate and provide an important heat source for the first stars in the universe. This and the following contribution present the story of Dark Stars. In this first part, we describe the conditions under which dark stars form in the early universe: 1) high dark matter densities, 2) the annihilation products get stuck inside the star, and 3) dark matter heating wins over all other cooling or heating mechanisms.

We test the proposal by El-Zant et al that the dark matter density of halos could be reduced through dynamical friction acting on heavy baryonic clumps in the early stages of galaxy formation. Using N-body simulations, we confirm that the inner halo density cusp is flattened to 0.2 of the halo break radius by the settling of a single clump of mass \ga 0.5% of the halo mass. We also find that an ensemble of 50 clumps each having masses \ga 0.2% can flatten the cusp to almost the halo break radius on a time scale of \sim9 Gyr, for an NFW halo of concentration 15. We summarize some of the difficulties that need to be overcome if this mechanism is to resolve the apparent conflict between the observed inner densities of galaxy halos and the predictions of LCDM.

The first phase of stellar evolution in the history of the universe may be Dark Stars, powered by dark matter heating rather than by fusion. Weakly interacting massive particles, which are their own antiparticles, can annihilate and provide an important heat source for the first stars in the the universe. This talk presents the story of these Dark Stars. We make predictions that the first stars are very massive ($\sim 800 M_\odot$), cool (6000 K), bright ($\sim 10^6 L_\odot$), long-lived ($\sim 10^6$ years), and probable precursors to (otherwise unexplained) supermassive black holes. Later, once the initial DM fuel runs out and fusion sets in, DM annihilation can predominate again if the scattering cross section is strong enough, so that a Dark Star is born again.

I show that the predicted densities of the inner dark matter halos in LCDM models of structure formation appear to be higher than estimates from real galaxies and constraints from dynamical friction on bars. This inconsistency would not be a problem for the LCDM model if physical processes that are omitted in the collisionless collapse simulations were able to reduce the dark matter density in the inner halos. I review the mechanisms proposed to achieve the needed density reduction.

The first stars in the universe form inside $\sim 10^6 M_\odot$ dark matter (DM) haloes whose initial density profiles are laid down by gravitational collapse in hierarchical structure formation scenarios. During the formation of the first stars in the universe, the baryonic infall compresses the dark matter further. The resultant dark matter density is presented here, using an algorithm originally developed by Young to calculate changes to the profile as the result of adiabatic infall in a spherical halo model; the Young prescription takes into account the non-circular motions of halo particles. The density profiles obtained in this way are found to be within a factor of two of those obtained using the simple adiabatic contraction prescription of Blumenthal et al. Our results hold regardless of the nature of the dark matter or its interactions and rely merely on gravity. If the dark matter consists of weakly interacting massive particles, which are their own antiparticles, their densities are high enough that their annihilation in the first protostars can indeed provide an important heat source and prevent the collapse all the way to fusion. In short, a ``Dark Star'' phase of stellar evolution, powered by DM annihilation, may indeed describe the first stars in the universe.

Spiral patterns are important agents of galaxy evolution. In this review, I describe how the redistribution of angular momentum by recurrent transient spiral patterns causes the random speeds of stars to rise over time, metallicity gradients to be reduced, and drives large-scale turbulence in the disk, which could be important for galactic dynamos. I also outline a possible mechanism for the recurrence of spiral instabilities and supporting evidence from solar neighborhood kinematics. Finally, I confirm that cloud scattering alone would predict the local velocity ellipsoid to be flattened, contrary to long-held expectations.

The predicted central densities of dark matter halos in LCDM models exceed those observed in some galaxies. Weinberg & Katz argue that angular momentum transfer from a rotating bar in the baryonic disk can lower the halo density, but they also contend that N-body simulations of this process will not reveal the true continuum result unless many more than the usual numbers of particles are employed. Adopting their simplified model of a rotating rigid bar in a live halo, I have been unable to find any evidence to support their contention. I find that both the angular momentum transferred and the halo density change are independent of the number of particles over the range usually employed up to that advocated by these authors. I show that my results do not depend on any numerical parameters, and that field methods perform equally with grid methods. I also identify the reasons that the required particle number suggested by Weinberg & Katz is excessive. I further show that when countervailing compression by baryonic settling is ignored, moderate bars can reduce the mean density of the inner halo by 20% - 30%. Long, massive, skinny bars can reduce the mean inner density by a factor ~10. The largest density reductions are achieved at the expense of removing most of the angular momentum likely to reside in the baryonic component. Compression of the halo by baryonic settling must reduce, and may even overwhelm, the density reduction achievable by bar friction.

The question whether the dark matter halo density in the centers of galaxies could be changed through interactions with a rotating bar in the baryonic disk is of considerable current interest. While N-body simulations have been used to address this question, it has also been claimed that results from such simulations cannot be trusted. Based on a perturbative treatment of resonant exchanges between orbits and a rotating perturbation, Weinberg & Katz contend that N-body simulations of this process will not reveal the continuum result unless many more than the usual numbers of particles are employed. Here I report a study designed to examine their contention, finding results that show no dependence on the number of particles over the range usually employed up to that advocated by these authors. I show that my results are independent of all numerical parameters, and that field methods perform equally with grid methods in this respect. I also identify the reasons that the required particle number suggested by Weinberg & Katz is excessive.

We present new observations of the strongly-barred galaxy NGC 1365, including new photometric images and Fabry-Perot spectroscopy, as well as a detailed re-analysis of the neutral hydrogen observations from the VLA archive. We find the galaxy to be at once remarkably bi-symmetric in its I-band light distribution and strongly asymmetric in the distribution of dust and in the kinematics of the gas in the bar region. The velocity field mapped in the H-alpha line reveals bright HII regions with velocities that differ by 60 to 80 km/s from that of the surrounding gas, which may be due to remnants of infalling material. We have attempted hydrodynamic simulations of the bar flow to estimate the separate disk and halo masses, using two different dark matter halo models and covering a wide range of mass-to-light ratios (Upsilon) and bar pattern speeds (Omega_p). None of our models provides a compelling fit to the data, but they seem most nearly consistent with a fast bar, corotation at sim 1.2r_B, and Upsilon_I simeq 2.0 +- 1.0, implying a massive, but not fully maximal, disk. The fitted dark halos are unusually concentrated, a requirement driven by the declining outer rotation curve.

We give a brief description of a new model for non-circular motions in disk galaxy velocity fields, that does not invoke epicycles. We assume non-circular motions to stem from a bar-like or oval distortion to the potential, as could arise from a triaxial halo or a bar in the mass distribution of the baryons. We apply our model to the high-quality CO and Halpha kinematics of NGC 2976 presented by Simon et al. 2003; it fits the data as well as their model with unrealistic radial flows, but yields a steeper rotation curve. Our analysis and other evidence suggests that NGC 2976 hosts a bar, implying a large baryonic contribution to the potential and thus limiting the allowed dark matter halo density.

We present a new procedure to fit non-axisymmetric flow patterns to 2-D velocity maps of spiral galaxies. We concentrate on flows caused by bar-like or oval distortions to the total potential that may arise either from a non-axially symmetric halo or a bar in the luminous disk. We apply our method to high-quality CO and Halpha data for the nearby, low-mass spiral NGC 2976 previously obtained by Simon et al., and find that a bar-like model fits the data at least as well as their model with large radial flows. We find supporting evidence for the existence of a bar in the baryonic disk. Our model suggests that the azimuthally averaged central attraction in the inner part of this galaxy is larger than estimated by these authors. It is likely that the disk is also more massive, which will limit the increase to the allowed dark halo density. Allowance for bar-like distortions in other galaxies may either increase or decrease the estimated central attraction.

We present a non-parametric method for decomposition of the light of disk galaxies into disk, bulge and bar components. We have developed and tested the method on a sample of 68 disk galaxies for which we have acquired I-band photometry. The separation of disk and bar light relies on the single assumption that the bar is a straight feature with a different ellipticity and position angle from that of the projected disk. We here present the basic method, but recognise that it can be significantly refined. We identify bars in only 47% of the more nearly face-on galaxies in our sample. The fraction of light in the bar has a broad range from 1.3% to 40% of the total galaxy light. If low-luminosity galaxies have more dominant halos, and if halos contribute to bar stability, the luminosity functions of barred and unbarred galaxies should differ markedly; while our sample is small, we find only a slight difference of low significance.

Using two-dimensional velocity maps and I-band photometry, we have created mass models of 40 spiral galaxies using the Milgrom relation (the basis of modified Newtonian dynamics, or MOND) to complement previous work. A Bayesian technique is employed to compare several different dark matter halo models to Milgrom and Newtonian models. Pseudo-isothermal dark matter halos provide the best statistical fits to the data in a majority of cases, while the Milgrom relation generally provides good fits as well. We also find that Milgrom models give mass-to-light ratios that roughly correlate with galaxy color, as predicted by stellar population models. A subsample of galaxies in the Hydra cluster follow a tight relation between mass-to-light and color, but one that is significantly different from relations found in previous studies. Ruling out the Milgrom relation with rotational kinematics is difficult due to systematic uncertainties in the observations as well as underlying model assumptions. We discuss in detail two galaxies for which the Milgrom relation appears to fail and find that relaxing the assumption of constant stellar mass-to-light ratio can maintain Milgrom models' viability.

Some 30% of disc galaxies have a pronounced central bar feature in the disc plane and many more have weaker features of a similar kind. Kinematic data indicate that the bar constitutes a major non-axisymmetric component of the mass distribution and that the bar pattern tumbles rapidly about the axis normal to the disc plane. The observed motions are consistent with material within the bar streaming along highly elongated orbits aligned with the rotating major axis. A barred galaxy may also contain a spheroidal bulge at its centre, spirals in the outer disc and, less commonly, other features such as a ring or lens. Mild asymmetries in both the light and kinematics are quite common. We review the main problems presented by these complicated dynamical systems and summarize the effort so far made towards their solution, emphasizing results which appear secure. (Truncated)

Recent ideas for the origin and persistence of the warps commonly observed in disc galaxies have focused on cosmic infall. We present N-body simulations of an idealized form of cosmic infall onto a disc galaxy and obtain a warp that closely resembles those observed. The inner disc tilts remarkably rigidly, indicating strong cohesion due to self-gravity. The line of nodes of the warp inside $R_{26.5} ~ 4.5 \Rd$ is straight, while that beyond $R_{26.5}$ generally forms a loosely-wound, leading spiral in agreement with Briggs's rules. We focus on the mechanism of the warp and show that the leading spiral arises from the torques from the misaligned inner disc and its associated inner oblate halo. The fact that the line of nodes of most warps forms a leading spiral might imply that the disc mass is significant in the centre. If the line of nodes can be traced to very large radii in future observations, it may reveal information on the mass distribution of the outer halo. The warp is not strongly damped by the halo because the precession rate of the inner disc is slow and the inner halo generally remains aligned with the inner disc. Thus even after the imposed quadrupolar perturbation is removed, the warp persists for a few Gyrs, by which time another infall event can be expected.

It is well-established that strong bars rotating in dense halos generally slow down as they lose angular momentum to the halo through dynamical friction. Angular momentum exchanges between the bar and halo particles take place at resonances. While some particles gain and others lose, friction arises when there is an excess of gainers over losers. This imbalance results from the generally decreasing numbers of particles with increasing angular momentum, and friction can therefore be avoided if there is no gradient in the density of particles across the major resonances. Here we show that anomalously weak friction can occur for this reason if the pattern speed of the bar fluctuates upwards. After such an event, the density of resonant halo particles has a local inflexion created by the earlier exchanges, and bar slowdown can be delayed for a long period; we describe this as a metastable state. We show that this behavior in purely collisionless N-body simulations is far more likely to occur in methods with adaptive resolution. We also show that the phenomenon could arise in nature, since bar-driven gas inflow could easily raise the bar pattern speed enough to reach the metastable state. Finally, we demonstrate that mild external, or internal, perturbations quickly restore the usual frictional drag, and it is unlikely therefore that a strong bar in a galaxy having a dense halo could rotate for a long period without friction.

It has been known for some time that rotating bars in galaxies slow due to dynamical friction against the halo. However, recent attempts to use this process to place constraints on the dark matter density in galaxies and possibly also to drive dark matter out of the center have been challenged. This paper uses simplified numerical experiments to clarify several aspects of the friction mechanism. I explicitly demonstrate the Chandrasekhar scaling of the friction force with bar mass, halo density, and halo velocity dispersion. I present direct evidence that exchanges between the bar and halo orbits at major resonances are responsible for friction and study both individual orbits and the net changes at these resonances. I also show that friction alters the phase space density of particles in the vicinity of a major resonance, which is the reason the magnitude of the friction force depends on the prior evolution. I demonstrate that bar slow down can be captured correctly in simulations having modest spatial resolution and practicable numbers of particles. Subsequent papers in this series delineate the dark matter density that can be tolerated in halos of different density profiles.

A bar rotating in a pressure-supported halo generally loses angular momentum and slows down due to dynamical friction. Valenzuela & Klypin report a counter-example of a bar that rotates in a dense halo with little friction for several Gyr, and argue that their result invalidates the claim by Debattista & Sellwood that fast bars in real galaxies require a low halo density. We show that it is possible for friction to cease for a while should the pattern speed of the bar fluctuate upward. The reduced friction is due to an anomalous gradient in the phase-space density of particles at the principal resonance created by the earlier evolution. The result obtained by Valenzuela & Klypin is probably an artifact of their adaptive mesh refinement method, but anyway could not persist in a real galaxy. The conclusion by Debattista & Sellwood still stands.

The initial radial density profiles of dark matter halos are laid down by gravitational collapse in hierarchical structure formation scenarios and are subject to further compression as baryons cool and settle to the halo centers. We here describe an explicit implementation of the algorithm, originally developed by Young, to calculate changes to the density profile as the result of adiabatic infall in a spherical halo model. Halos with random motion are more resistant to compression than are those in which random motions are neglected, which is a key weakness of the simple method widely employed. Young's algorithm results in density profiles in excellent agreement with those from N-body simulations. We show how the algorithm may be applied to determine the original uncompressed halos of real galaxies, a step which must be computed with care in order to enable a confrontation with theoretical predictions from theories such as LCDM.

We model the mass distributions of 40 high surface brightness spiral galaxies inside their optical radii, deriving parameters of mass models by matching the predicted velocities to observed velocity maps. We use constant mass-to-light disk and bulge models, and we have tried fits with no halo and with three different halo density profiles. The data require a halo in most, but not all, cases, while in others the best fit occurs with negligible mass in the luminous component, which we regard as unphysical. All three adopted halo profiles lead to fits of about the same quality, and our data therefore do not constrain the functional form of the halo profile. The halo parameters display large degeneracies for two of the three adopted halo functions, but the separate luminous and dark masses are better constrained. However, the fitted disk and halo masses vary substantially between the adopted halo models, indicating that even high quality 2-D optical velocity maps do not provide significant constraints on the dark matter content of a galaxy. We demonstrate that data from longslit observations are likely to provide still weaker constraints. We conclude that additional information is needed in order to constrain the separate disk and halo masses in a galaxy.

Newtonian mechanics indicates that galaxies and galaxy clusters are much more massive than we would have guessed from their luminosities, with the discrepancy being generally attributed to dark matter halos. An alternative hypothesis is that accelerations in very weak gravitational fields are larger than predicted by Newton's laws, and there is no need for dark matter. Even though we do not currently have a satisfactory theory associated with this rival hypothesis, we can ask whether any observational tests could rule it out or prefer it over the dark matter hypothesis. Current evidence suggests that neither hypothesis enjoys a decisive advantage over the other. If dark matter turns out to be the correct interpretation however, then theories of galaxy formation face some quite severe fine-tuning problems.

More than two thirds of disk galaxies are barred to some degree. Many today harbor massive concentrations of gas in their centers, and some are known to possess supermassive black holes (SMBHs) and their associated stellar cusps. Previous theoretical work has suggested that a bar in a galaxy could be dissolved by the formation of a mass concentration in the center, although the precise mass and degree of central concentration required is not well-established. We report an extensive study of the effects of central masses on bars in high-quality N-body simulations of galaxies. We have varied the growth rate of the central mass, its final mass and degree of concentration to examine how these factors affect the evolution of the bar. Our main conclusions are: (1) Bars are more robust than previously thought. The central mass has to be as large as several percent of the disk mass to completely destroy the bar on a short timescale. (2) For a given mass, dense objects cause the greatest reduction in bar amplitude, while significantly more diffuse objects have a lesser effect. (3) The bar amplitude always decreases as the central mass is grown, and continues to decay thereafter on a cosmological time-scale. (4) The first phase of bar-weakening is due to the destruction by the CMC of lower-energy, bar-supporting orbits, while the second phase is a consequence of secular changes to the global potential which further diminish the number of bar-supporting orbits. We provide detailed phase-space and orbit analysis to support this suggestion. Thus current masses of SMBHs are probably too small, even when dressed with a stellar cusp, to affect the bar in their host galaxies. The molecular gas concentrations found in some barred galaxies are also too diffuse to affect the amplitude of the bar significantly.

Any connection between central activity and the large-scale dynamics of disk galaxies requires an efficient mechanism to remove angular momentum from the orbiting material. The only viable means of achieving inflow from kiloparsec scales is through gravitational stresses created by bars and/or mergers. The inflow of gas in bars today appears to stall at a radius of few hundred parsec, however, forming a nuclear ring. Here we suggest that bars in the early Universe may have avoided this problem, and propose that the progenitors of central supermassive black holes (SMBHs) are created by gas that is driven deep into the centers of galaxies by bars in the early stages of disk formation. The coincidence of the QSO epoch with galaxy formation, the short lifetimes of QSOs, and the existence of SMBHs in the centers of most bright galaxies are all naturally accounted for by disk dynamics in this model. The progenitor SMBHs are the seeds for brighter QSO flares during galaxy mergers. We present a new study of bar weakening by central mass concentrations, which shows that bars are less easily destroyed than previously thought. An extremely massive and compact central mass can, however, dissolve the bar, which creates a pseudo-bulge component in the center of the disk.

Many barred galaxies today harbor massive concentrations of gas in their centers, and some are known to possess supermassive black holes (SBHs) and their associated stellar cusps. Previous theoretical work has suggested that a bar in a galaxy could be dissolved by the formation of a mass concentration in the center, although the precise mass and degree of central concentration required is not well-established. Here we report an extensive study of the effects of central masses on bars in high-quality N-body simulations of galaxies. We study models containing both strong and weak bars, and quantify the change in bar amplitude as a mass concentration is grown in their centers. We find that bars are more robust than previously thought. We have varied the growth rate of the central mass, its final mass and the degree of concentration of the mass. Our main conclusions are: (1) the central mass has to be as large as several percent of the disk mass to completely destroy the bar; (2) for a given mass, objects whose scale-length is a few pc or less cause the greatest reduction in bar amplitude, while significantly more diffuse objects have a lesser effect; and (3) the bar amplitude always decreases as the central mass is grown, and continues to decay thereafter on a cosmological timescale. Thus the masses of SBHs are probably too small, even when dressed with a stellar cusp, to affect the bar amplitude significantly.

Self-consistent bars that form in galaxies embedded within cuspy halos are unable to flatten the cusp. Short bars form in models with quasi-flat rotation curves. They lose angular momentum to the halo through dynamical friction, but the continuous concentration of mass within the disk as the bar grows actually compresses the halo further, overwhelming any density reduction due to the modest angular momentum transfer to the halo. Thus the Weinberg-Katz proposed solution to the non-existence of the predicted cuspy halos from CDM simulations would seem to be unworkable. I also find that the concerns over the performance of N-body codes raised by these authors do not apply to the methods used here.

We investigate the impact of nonaxisymmetric structure on estimates of galaxy inclinations and position angles. A new minimization technique is used to obtain estimates of inclination and position angle from a global fit to either photometric or kinematic data. We discuss possible systematic uncertainties which are much larger than statistical uncertainties. Our investigation reveals that systematic uncertainties associated with fitting photometric data dominate the formal statistical uncertainties. For our sample of inclined galaxies, we estimate that nonaxisymmetric features introduce inclination and position angle uncertainties of approximately 5 degrees, on average. The magnitudes of these uncertainties weaken the arguments for intrinsically elliptical galaxy disks.

Transient spiral waves of moderate amplitude cause substantial changes to the angular momenta of many stars in a galaxy disk. Stars near to corotation are affected most strongly: for a wave of ~20% overdensity, the rms change for particles near to corotation is ~18% of their initial angular momenta. Yet these changes produce no increase in random motion near corotation, while the same wave causes mild heating near its Lindblad resonances. The direction of radial migration depends both on initial location relative to corotation and on azimuth. Streamlines of cold gas become very significantly distorted before they cross, suggesting that gas is also mixed radially by this mechanism. Such continuous churning of the disk has profound implications for abundance gradients in galaxies. Furthermore, it is quite possible that the Sun has migrated significantly from its radius of birth, consistent with the claims of Wielen et al.

We show that spiral waves in galaxy discs churn the stars and gas in a manner that largely preserves the overall angular momentum distribution and leads to little increase in random motion. Changes in the angular momenta of individual stars are typically as large as ~50% over the lifetime of the disk. The changes are concentrated around the corotation radius for an individual spiral wave, but since transient waves with a wide range of pattern speeds develop in rapid succession, the entire disk is affected. This behaviour has profound consequences for the metallicity gradients with radius in both stars and gas. The ISM is also stirred by the same mechanism. We find observational support for stirring, propose a simple model for the distribution of stars over metallicity and age, and discuss other possible consequences.

Dynamical mass discrepancies in galaxies have two possible explanations: the existence of large amounts of dark matter or the breakdown of Newtonian gravity. True dark matter halos of galaxies could differ in several respects from apparent halos inferred from misinterpreted modified gravity, and we discuss two observational tests to probe these differences. Modified gravity gives generic, geometrical predictions for the shapes of apparent halos, and also requires that the inferred excess mass strongly correlates with the visible baryonic mass. Thus modified gravity is more easily falsified, but we find both hypotheses to be consistent with the meager data on galaxy and cluster halos; dark matter is not required by current data. We strongly encourage the acquisition of more data, which could ultimately eliminate one or other of these possibilities.

I give a very brief review of aspects of internal dynamics that affect the global shape of a galaxy, focusing on triaxiality, bars and warps. There is general agreement that large central masses can destroy triaxial shapes, but recent simulations of this process seem to suffer from numerical difficulties. Central black holes alone are probably not massive enough to destroy global triaxiality, but when augmented by star and gas concentrations in barred galaxies, the global shape may be affected. Even though we do not understand the origin of bars in galaxies, they are very useful as probes of the dark matter density of the inner halo. Finally, I note that dynamical friction acts to reduce a misalignment between the spin axes of the disk and halo, producing a nice warp in the outer disk which has many of the properties of observed galactic warps.

We demonstrate the existence of long-lived, large-amplitude eccentric modes in massive, annular particle disks in orbit about a central mass. The lopsided modes we have found precess slowly in the prograde direction at a rate which increases with disk mass and decreases with the amplitude of the distortion. The lopsidedness generally survives for as long as we run the calculations and may last indefinitely; we found no significant decay in over 700 disk-particle orbits in one case. These strongly non-linear modes are readily created using a number of different perturbing rules, but we find no evidence for linear instabilities in unperturbed disks. Our results suggest that Tremaine's eccentric disk model for the nucleus of M31 may be viable, but none of our models matches all aspects of the observational data.

The success of the Lambda-CDM model on large scales does not extend down to galaxy scales. We list a dozen problems of the dark matter hypothesis, some of which arise in specific models for the formation of structure in the universe, while others are generic and require fine tuning in any dark matter theory. Modifications to the theory, such as adding properties to the DM particles beyond gravitational interactions, or simply a better understanding of the physics of galaxy formation, may resolve some problems, but a number of conspiracies and correlations are unlikely to yield to this approach. The alternative is that mass discrepancies result from of a non-Newtonian law of gravity, a hypothesis which avoids many of the more intractable problems of dark matter. A modified law of gravity is not without formidable difficulties of its own, but it is no longer obvious that they are any more daunting than those facing DM.

We report a dynamical determination of the separate contributions of disk and dark halo masses to the rotation curve of a spiral galaxy. We use fluid-dynamical models of gas flow in the barred galaxy NGC 4123 to constrain the dynamical properties of the galaxy: disk M/L, bar pattern speed, and the central density and scale radius of the dark halo. We derive a realistic barred potential directly from the light distribution. For each model we assume a value of the stellar M/L and a bar pattern speed Omega_p and add a dark halo to fit the rotation curve. We then compute the gas flow velocities with a 2-D gas dynamical code, and compare the model flow patterns to a 2-D velocity field derived from Fabry-Perot observations. The strong shocks and non-circular motions in the observed gas flow require a high stellar M/L and a fast-rotating bar. Models with I-band disk M/L of 2.0 -- 2.5 h_75, or 80 -- 100% of the maximum disk value, are highly favored. The corotation radius of the bar must be <= 1.5 times the bar semi-major axis. These results contradict some recent claimed ``universal'' galaxy disk/halo relations, since NGC 4123 is of modest size (rotation curve maximum 145 km/sec, and V_flat = 130 km/sec) yet is quite disk-dominated. The dark halo of NGC 4123 is less concentrated than favored by current models of dark halos based on cosmological simulations. Since some 30% of bright disk galaxies are strongly barred and have dust lanes indicating shock morphology similar to that of NGC 4123, it is likely that they also have high stellar M/L and low density halos. We suggest that luminous matter dominates inside the optical radius R_25 of high surface brightness disk galaxies.

The non-circular streaming motions in barred galaxies are sensitive to the mass of the bar and can be used to lift the degeneracy between disk and dark matter halo encountered when fitting axisymmetric rotation curves of disk galaxies. In this paper, we present photometric and kinematic observations of NGC 4123, a barred galaxy of modest size (V_rot = 130 km/sec, L = 0.7 L_*), which reveal strong non-circular motions. The bar has straight dust lanes and an inner Lindblad resonance. The disk of NGC 4123 has no sign of truncation out to 10 scale lengths, and star-forming regions are found well outside R_25. A Fabry-Perot H-alpha velocity field shows velocity jumps of >100 km/sec at the location of the dust lanes within the bar, indicating shocks in the gas flow. VLA observations yield the velocity field of the H I disk. Axisymmetric mass models yield good fits to the rotation curve outside the bar regionfor disk I-band M/L of 2.25 or less, and dark halos with either isothermal or power-law profiles can fit the data well. In a companion paper, we model the full 2-D velocity field, including non-circular motions, to determine the stellar M/L and the mass of the dark halo.

The paper claimed that significant radial migration of stars in a stellar disk like that of the Milky Way could not occur. We now think that while the treatment of the effects of molecular clouds was correct, the paper seriously underestimated the ability of spiral arms to shift the radii of stars that corotate with them. Consequently, it is likely that significant radial migration_is_ possible.

The central densities of dark matter (DM) halos are much lower than predicted in cold DM models of structure formation. Confirmation that they have cores with a finite central density would allow us to rule out many popular types of collisionless particle as candidates for DM. Any model that leads to cusped halos (such as cold DM) is already facing serious difficulties on small scales and hot DM models have been excluded. Here I show that fermionic warm DM is inconsistent with the wide range of phase space densities in the DM halos of well-observed nearby galaxies.

We show that bars in galaxy models having halos of moderate density and a variety of velocity distributions all experience a strong drag from dynamical friction unless the halo has large angular momentum in the same sense as the disk. The frictional drag decreases the bar pattern speed, driving the corotation point out to distances well in excess of those estimated in barred galaxies. The halo angular momentum required to avoid strong braking is unrealistically large, even when rotation is confined to the inner halo only. We conclude, therefore, that bars are able to maintain their observed high pattern speeds only if the halo has a central density low enough for the disk to provide most of the central attraction in the inner galaxy. We present evidence that this conclusion holds for all bright galaxies.

We confirm that a high rate of shear at the center is able to stabilize an entire stellar disk against bar-forming (m=2) modes, irrespective of the dark halo density. Our simulations of unstable power-law disks also yield the growth rate and pattern speed of the dominant lop-sided (m=1) mode in close agreement with that predicted by linear theory in our cleanest case. The one-armed modes, which dominate models with extensive disks, can be controlled by tapering the outer disk. Two-armed instabilities of hard-centered disks are more difficult to identify because they are easily overwhelmed by particle noise. Nevertheless, we have detected the predicted m=2 modes in simulations with very large numbers of particles. Such bisymmetric instabilities in these disks are provoked only by sharp edges and are therefore easily eliminated. We have constructed a cool disk model with an almost flat rotation curve and a quasi-exponential density profile that is unambiguously stable. The halo in this stable model has a large core radius, with the disk and bulge providing almost all the rotational support in the inner parts.

Having once thought that we understood how some galaxies were barred and had difficulty accounting for the absence of bars in others, it now seems that we have the opposite problem. Most real galaxies have centres dense enough to inhibit bar formation, even if they have massive discs. This is also true for barred galaxies, which therefore could not have acquired their bars as a result of a self-excited, global instability. There is a hint from the morphology of galaxies in the Hubble Deep Fields that the growth of bars might be a slow or secular process. Here I discuss possible mechanisms that could form bars long after the disc is assembled.

A long-standing controversy in studies of spiral structure has concerned the lifetimes of individual spiral patterns. Much theoretical work has sought quasi-stationary spiral modes while N-body simulations have consistently displayed recurrent, short-lived patterns. The simulations manifest a recurrent cycle of true instabilities related to small-scale features in the angular momentum distribution of particles, with the decay of each instability seeding the growth of the next. Data from the recent Hipparcos mission seem to offer support for the recurrent transient picture.

Longitude-velocity (l-V) diagrams of H I and CO gas in the inner Milky Way have long been known to be inconsistent with circular motion in an axisymmetric potential. Several lines of evidence suggest that the Galaxy is barred, and gas flow in a barred potential could be consistent with the observed ``forbidden'' velocities and other features in the data. We compare the H I observations to l-V diagrams synthesized from 2-D fluid dynamical simulations of gas flows in a family of barred potentials. The gas flow pattern is very sensitive to the parameters of the assumed potential, which allows us to discriminate among models. We present a model that reproduces the outer contour of the H I l-V diagram reasonably well; this model has a strong bar with a semimajor axis of 3.6 kpc, an axis ratio of approximately 3:1, an inner Lindblad resonance (ILR), and a pattern speed of 42 km/s/kpc, and matches the data best when viewed from 34\deg to the bar major axis. The behavior of the models, combined with the constraint that the shocks in the Milky Way bar should resemble those in external barred galaxies, leads us to conclude that wide ranges of parameter space are incompatible with the observations. In particular we suggest that the bar must be fairly strong, must have an ILR, and cannot be too end-on, with the bar major axis at 35\deg +/- 5\deg to the line of sight. The H I data exhibit larger forbidden velocities over a wider longitude range than are seen in molecular gas; this important difference is the reason our favored model differs so significantly from other recently proposed models.

A model that accounts for the brief flaring of QSOs in the early stages of galaxy formation is proposed. I argue that a bar must develop early in the life of nearly every galaxy and that gas to create and fuel the QSO is driven into the center of the galaxy by the bar. The QSO lifetime, and the mass of its central engine, are also controlled by large-scale dynamics, since the fuel supply is shut off after a short period by the development of an inner Lindblad resonance. This resonance causes the gas inflow along the bar to stall at a distance of a few hundred parsecs from the center. The ILR develops as a result of previous inflow, making quasar activity self-limiting. The bars are weakened and can be destroyed by the central mass concentration formed in this way.

In this review, I discuss just three aspects of the stability and evolution of galactic discs. (1) I first review our understanding of the bar instability and how it can be controlled. Disc galaxies in which the orbital speed does not decrease much towards the centre have no difficulty avoiding bars, even when dark matter makes an insignificant contribution to the inner part of the rotation curve. (2) I then briefly discuss interactions between disturbances in the discs of galaxies and the spherical components, which generally exert a damping effect through dynamical friction. The fact that bars in real galaxies appear to rotate quite rapidly, seems to require dark matter halos to have large, low-density cores. (3) In the remainder of the article, I consider the theory of spiral structure. The new development here is that the distribution function for stars in the Solar neighbourhood, as measured by HIPPARCOS, is far less smooth than most theoretical work had previously supposed. The strong variations in the values of the \DF over small ranges in angular momentum have the appearance of having been caused by scattering at Lindblad resonances with spiral patterns. This result, if confirmed when the radial velocity data become available, supports the picture of spiral patterns as dynamical instabilities driven by substructure in the \DF. The details of how decaying patterns might seed conditions for a new instability remain unclear, and deserve fresh attention.

While the total interior mass of a galaxy is reasonably well determined by a good rotation curve, the relative contributions from disk, bulge and halo are only weakly constrained by one-dimensional data. Barred galaxies are intrinsically more complicated, but provide much tighter constraints on the disk masses and support the idea that most of the mass in the inner parts of bright galaxies is in their stars. There appears to be no systematic difference in dark matter content between barred and unbarred galaxies, consistent with the theoretical result that the global stability of galaxies with dense centers does not depend on their halo fraction. The rotation curve shapes of lower luminosity and low-surface-brightness galaxies, on the other hand, indicate significant mass in the DM halo even near their centers. We find that most DM halos appear to have large cores, inconsistent with the predictions from cosmological simulations. We also show that such large-core halos can result from compression by disk infall of physically reasonable initial halos. Maximum disks, while apparently required by the data, do seem to present some puzzles; most notably they re-open the old disk-halo ``conspiracy'' issue and incorrectly predict that surface brightness should be a second parameter in the Tully-Fisher relation.

A galaxy disk embedded in a rotating halo experiences a dynamical friction force which causes it to warp when the angular momentum axes of the disk and halo are misaligned. Our fully self-consistent simulations of this process induce long-lived warps in the disk which mimic Briggs's rules of warp behavior. They also demonstrate that random motion within the disk adds significantly to its stiffness. Moreover, warps generated in this way have no winding problem and are more pronounced in the extended \h1 disk. As emphasized by Binney and his co-workers, angular momentum misalignments, which are expected in hierarchical models of galaxy formation, can account for the high fraction of warped galaxies. Our simulations exemplify the role of misaligned spins in warp formation even when the halo density is not significantly flattened.

We propose that massive central objects form in the centers of the bars which must develop in young high-surface density galactic disks. Large-scale dynamics shuts off the growth of the central mass before it reaches $\sim 2%$ of the disk mass at the time, but this mass is sufficient to weaken the bar substantially. Subsequent evolution of the galaxy can either complete the destruction of the bar or cause it to recover, depending upon the angular momentum distribution of later infalling material. We produce massive, fully self-gravitating disks having roughly flat rotation curves which are quite stable. If at least part of the central masses we require constitute the engines of QSOs, then our picture naturally accounts for their redshift dependence, since the fuel supply is shut off by the development of an inner Lindblad resonance. A prediction is that massive objects should not be found in halo dominated galaxies, such as low-luminosity, or low-surface brightness galaxies.

When present, extended disks of neutral hydrogen around spiral galaxies show a remarkably uniform velocity dispersion of approx 6 km/s. Since stellar winds and supernovae are largely absent in such regions, neither the magnitude nor the constancy of this number can be accounted for in the classical picture in which interstellar turbulence is driven by stellar energy sources. Here we suggest that magnetic fields with strengths of a few microgauss in these extended disks allow energy to be extracted from galactic differential rotation through MHD driven turbulence. The magnitude and constancy of the observed velocity dispersion may be understood if its value is Alfvenic. Moreover, by providing a simple explanation for a lower bound to the gaseous velocity fluctuations, MHD processes may account for the sharp outer edge to star formation in galaxy disks.

We study the vertical heating and thickening of galaxy disks due to accretion of small satellites. Our simulations are restricted to axial symmetry, which largely eliminates numerical evolution of the target galaxy but requires the trajectory of the satellite to be along the symmetry axis of the target. We find that direct heating of disk stars by the satellite is not important because the satellite's gravitational perturbation has little power at frequencies resonant with the vertical stellar orbits. The satellite does little damage to the disk until its decaying orbit resonantly excites large-scale disk bending waves. Bending waves can damp through dynamical friction from the halo or internal wave-particle resonances; we find that wave-particle resonances dominate the damping. The principal vertical heating mechanism is therefore dissipation of bending waves at resonances with stellar orbits in the disk. Energy can thus be deposited some distance from the point of impact of the satellite. The net heating from a tightly bound satellite can be substantial, but satellites that are tidally disrupted before they are able to excite bending waves do not thicken the disk.

We use fully self-consistent N-body simulations of barred galaxies to show that dynamical friction from a dense dark matter halo dramatically slows the rotation rate of bars. Our result supports previous theoretical predictions for a bar rotating within a massive halo. On the other hand, low density halos, such as those required for maximum disks, allow the bar to continue to rotate at a high rate. There is somewhat meager observational evidence indicating that bars in real galaxies do rotate rapidly and we use our result to argue that dark matter halos must have a low central density in all high surface brightness disk galaxies, including the Milky Way. Bars in galaxies that have larger fractions of dark matter should rotate slowly, and we suggest that a promising place to look for such candidate objects is among galaxies of intermediate surface brightness.

The radial distribution of mass in a disk galaxy is strongly constrained by its rotation curve. The separate contributions from the individual stellar populations and dark matter (DM) are not easily disentangled, however, especially since there is generally no feature to indicate where the component dominating the central attraction switches from luminous to dark matter. Here I summarize three recent thesis projects at Rutgers University which all suggest that DM has a low density in the inner parts of bright galaxies, and that most of the mass therefore resides in the disk. In addition, I present some preliminary work on the Milky Way.

I compare various popular and unpopular techniques for simulating large collisionless stellar systems. I give a quantitative comparison of the raw cpu times required for five separate codes, including tree codes and basis function expansions, which demonstrate that grid codes are most efficient for large numbers of particles. Since efficiency is only one consideration when choosing a code, I discuss other strengths and weaknesses of the various methods. While some applications may require the maximum possible number of particles, I argue that quiet start techniques can often permit reliable results to be obtained with moderate particle numbers. I suggest a quiet start procedure for spherical stellar systems and show that it leads to a significant reduction in the relaxation rate. The combination of efficient codes and quiet starts allows many galactic dynamical problems to be tackled without the need for supercomputers.

We have examined the stability of a sequence of oblate elliptical galaxy models having the Stackel form suggested by Kuz'min & Kutuzov. We have employed the 2-integral DFs given by Dejonghe & de Zeeuw for which flattened non-rotating models are characterized by counter-streaming motion and are radially cool; we introduce net rotation in some models by changing the sign of the z-component of angular momentum for a fraction of the particles. We have found that all non-rotating and slowly rotating members of this sequence rounder than E7 are stable, and that even maximally rotating models rounder than E4 are stable. In the absence of strong rotation, the most disruptive instability, and the last to be stabilized by increasing thickness, is a lopsided (m=1) mode. This instability appears to be driven by counter-rotation in radially cool models. Its vigour is lessened as rotation is increased, but it remains strong even in models with net angular momentum 90% of that of a maximally rotating model before finally disappearing in maximally rotating models. Strongly rotating models are more unstable to bar-forming modes which afflict maximally rotating models with c/a >~ 0.5, but this mode is quickly stabilized by moderate fractions of counter-rotating particles. Bending instabilities appear not to be very important; they are detectable in the inner parts of the flatter models, but are less vigorous and more easily stabilized than the lop-sided or bar modes in every case. We briefly discuss the possible relevance of the lop-sided instability to the existence of many lop-sided disk galaxies.

We describe fully self-consistent N-body experiments of barred galaxies with massive halos. A rotating bar is braked through dynamical friction with the halo, which occurs on a short time scale when the central density of the halo is high. On the other hand, friction is weak in a model with a central halo density low enough that the disk dominates the rotation curve in the inner parts; this model supports a fast bar for a Hubble time. We conclude that real barred galaxies in which the bar rotates rapidly, which is believed to be the rule, must have close to maximum disks.

Self-gravitating stellar disks with random motion support both exponentially growing and, in some cases, purely oscillatory axisymmetric bending modes, unlike their cold disk counterparts. A razor-thin disk with even a very small degree of random motion in the plane is both unstable and possesses a discrete spectrum of neutral modes, irrespective of the sharpness of the edge. Random motion normal to the disk plane has a stabilizing effect but at the same time allows bending waves to couple to the internal vibrations of the particles, which causes the formerly neutral modes to decay through Landau damping. Focusing first on instabilities, I here determine the degree of random motion normal to the plane needed to suppress global, axisymmetric, bending instabilities in a family of self-gravitating disks. As found previously, bending instabilities are suppressed only when the thickness exceeds that expected from a na\"\i ve local criterion when the degree of pressure support within the disk plane is comparable to, or exceeds, the support from rotation. A modest disk thickness is adequate for the bending stability of most disk galaxies, except perhaps near their centers. The discretization of the neutral spectrum in a zero-thickness disk is due to the existence of a turning point for bending waves in a warm disk, which is absent when the disk is cold. When the disk is given a finite thickness, the discrete neutral modes generally become strongly damped through wave-particle interactions. It is surprising therefore that I find some simulations of warm, stable disks can support (quasi-)neutral, large-scale, bending modes that decay very slowly, if at all.

A strong bar rotating within a massive halo should lose angular momentum to the halo through dynamical friction, as predicted by Weinberg. We have conducted fully self-consistent, numerical simulations of barred galaxy models with a live halo population and find that bars are indeed braked very rapidly. Specifically, we find that the bar slows sufficiently within a few rotation periods that the distance from the centre to co-rotation is more than twice the semi-major axis of the bar. Observational evidence (meagre) for bar pattern speeds seems to suggest that this ratio typically lies between 1.2 to 1.5 in real galaxies. We consider, a number of possible explanations for this discrepancy between theoretical prediction and observation, and conclude that no conventional alternative seems able to account for it.

I review the various mechanisms for creating bars in rotating stellar disks, and conclude that the swing-amplified feed-back loop, which produces rapidly tumbling bars, remains the most probable. The bar continues to evolve after its formation in a number of ways; here I discuss what appears to be inevitable thickening normal to the plane, continued spiral activity in the outer disk and also underscore the increasingly important problem presented by angular momentum loss to the halo. Finally, I examine possible means, excluding interaction, by which bars in galaxies could be destroyed.

The stability of a galaxy model is most easily assessed through N-body simulation. Particle-mesh codes have been widely used for this purpose, since they enable the largest numbers of particles to be employed. We show that the functional expansion technique, originally proposed by Clutton-Brock for other simulation problems, is in fact superior for stability work. For simulations of linear evolution it is not much slower than grid methods using the same number of particles, and reproduces analytical results with much greater accuracy. This success rests on its ability to represent global modes with a modest number of basis functions; grid methods may be more effective for other applications, however. Our conclusions are based on implementations of functional expansion and grid algorithms for disk galaxies.