Search

When starting an N-body simulation of an isolated galaxy, it is desirable to select particles from a distribution function to ensure that the model is in equilibrium. Random sampling from a DF is widely used, but results in a set of particles that differs by shot noise from that intended. This paper presents a method to reduce sampling noise that has been developed by the author in a many collaborations over a number of years. The technique has been partly described in past papers, though the ideas have not previously been gathered together, nor have its advantages been clearly demonstrated in past work. Of course, sampling errors can also be reduced by a brute force increase in the number of particles, but methods to achieve the same effect with fewer particles have obvious advantages. Here we not only describe the method, but also present three sets of simulations to illustrate the practical advantages of reducing sampling error. The improvements are not dramatic, but are clearly worth having., Comment: Accepted to appear in MNRAS. 10 pages, 5 figures, two minor revisions

The problem of how some disk galaxies avoid forming bars remains unsolved. Many galaxy models having reasonable properties continue to manifest vigorous instabilities that rapidly form strong bars and no widely-accepted idea has yet been advanced to account for how some disk galaxies manage to avoid this instability. It is encouraging that not all galaxies formed in recent cosmological simulations possess bars, but the dynamical explanation for this result is unclear. The unstable mode that creates a bar is understood as a standing wave in a cavity that reflects off the disk center and the corotation radius, with amplification at corotation. Here we use simulations to address one further idea that may perhaps inhibit the feedback loop and therefore contribute to stability, which is to make the disk center dynamically hot and/or to taper away mass from the inner disk, which could be masked by a bulge. Unfortunately, we find that neither strategy makes much difference to the global stability of the disk in the models we have tried. While deep density cutouts do indeed prevent feedback through the center, they still reflect incoming waves and thereby provoke a slightly different instability that again leads to a strong bar., Comment: 22 pages, 15 figures, to appear in Apj. Comments welcome

This paper continues a series reporting different aspects of the behaviour of disc galaxy simulations that support spiral instabilities. The focus in this paper is to demonstrate how linear spiral instabilities saturate and decay, and how the properties of the disc affect the limiting amplitude of the spirals. Once again, we employ idealized models that each possess a single instability that we follow until it has run its course. Remarkably, we find a tight correlation between the growth rate of the mode and its limiting amplitude, albeit from only six simulations. We show that non-linear orbit deflections near corotation cause the mode to saturate, and that the more time available in a slowly-growing mode creates the critical deflections at lower amplitude. We also find that scattering at the inner Lindblad resonance is insignificant until after the mode has saturated. Our objective in this series of papers, which we believe we have now achieved, has been to develop a convincing and well-documented account of the physical behaviour of the spiral patterns that have been observed in simulations by others, and by ourselves, for many decades. Understanding the simulations is an important step towards the greater objective, which is to find observational evidence from galaxies that could confront the identified mechanism., Comment: 13 pages, 13 figures, to appear in MNRAS, corrected Fig 13 only

Recent 21 cm line observations of the ultra-diffuse galaxy AGC~114905 indicate a rotating disc largely supported against gravity by orbital motion, as usual. Remarkably, this study has revealed that the form and amplitude of the HI rotation curve is completely accounted for by the observed distribution of baryonic matter, stars and neutral gas, implying that no dark halo is required. It is surprising to find a DM-free galaxy for a number of reasons, one being that a bare Newtonian disk having low velocity dispersion would be expected to be unstable to both axi- and non-axisymmetric perturbations that would change the structure of the disc on a dynamical timescale, as has been known for decades. We present $N$-body simulations of the DM-free model, and one having a low-density DM halo, that confirm this expectation: the disc is chronically unstable to just such instabilities. Since it is unlikely that a galaxy that is observed to have a near-regular velocity pattern would be unstable, our finding calls into question the suggestion that the galaxy may lack, or have little, dark matter. We also show that if the inclination of this near face-on system has been substantially overestimated, the consequent increased amplitude of the rotation curve would accommodate a halo massive enough for the galaxy to be stable., Comment: 10 pages, to appear in MNRAS, minor revisions requested by the referee

Any perturbation to a disc galaxy that creates a misalignment between the planes of the inner and outer disc, will excite a slowly evolving bending wave in the outer disc. The torque from the stiff inner disc drives a retrograde, leading-spiral bending wave that grows in amplitude as it propagates outward over a period of several Gyr. This behaviour creates warps that obey the rules established from observations, and operates no matter what the original cause of the misalignment between the inner and outer disc. The part of the disc left behind by the outwardly propagating wave is brought into alignment with the inner disc. Here we confirm that mild warps in simulations of disc galaxies can be excited by shot noise in the halo, as was recently reported. We show that the quadrupole component of the noise creates disc distortions most effectively. Bending waves caused by shot noise in carefully constructed equilibrium simulations of isolated galaxies are far too mild to be observable, but perturbations from halo substructure and galaxy assembly must excite larger amplitude bending waves in real galaxies., Comment: 8 pages, 10 figs, revised version submitted to MNRAS

Spirals in galaxies have long been thought to be caused by gravitational instability in the stellar component of the disk, but discerning the precise mechanism had proved elusive. Tidal interactions, and perhaps bars, may provoke some spiral responses, but spirals in many galaxies must be self-excited. We survey the relevant observational data and aspects of disk dynamical theory. The origin of the recurring spiral patterns in simulations of isolated disk galaxies has recently become clear and it is likely that the mechanism is the same in real galaxies, although evidence to confirm this supposition is hard to obtain. As transient spiral activity increases random motion, the patterns must fade over time unless the disk also contains a dissipative gas component. Continuing spiral activity alters the structure of the disks in other ways: reducing metallicity gradients and flattening rotation curves are two of the most significant. The overwhelming majority of spirals in galaxies have two- or three-fold rotational symmetry, indicating that the cool, thin disk component is massive. Spirals in simulations of halo-dominated disks instead manifest many arms, and consequently do not capture the expected full spiral-driven evolution. We conclude by identifying areas where further work is needed., Comment: 50 pages, 9 figures, to appear in vol 60 of ARAA (2022) - this version includes revisions to the figures made by the copy-editor - Supplementary video available at http://www.physics.rutgers.edu/~sellwood/supp_material.html

In order to address the question of whether spiral disturbances in galaxy discs are gravitationally coupled to the halo, we conduct simulations of idealized models of disc galaxies. We compare growth rates of spiral instabilities in identical mass models in which the halo is held rigid or is represented by particles drawn from an equilibrium distribution function. We examine cases of radial and azimuthal bias in the halo velocity ellipsoid in one of our models, and an isotropic velocity distribution in both. We find at most marginal evidence for an enhanced growth rate of spiral modes caused by a halo supporting response. We also find evidence for very mild dynamical friction between the spiral disturbance and the halo. We offer an explanation to account for the different behaviour between spiral modes and bar modes, since earlier work had found that bar instabilities became significantly more vigorous when a responsive halo was substituted for an equivalent rigid mass distribution. The barely significant differences found here justify the usual simplifying approximation of a rigid halo made in studies of spiral instabilities in galaxies., Comment: 7 pages, 4 figures, accepted to appear in MNRAS

Assessing the likelihood that the rotation curve of a galaxy matches predictions from galaxy formation simulations requires that the uncertainties in the circular speed as a function of radius derived from the observational data be statistically robust. Few uncertainties presented in the literature meet this requirement. In this paper we present a new standalone tool, makemap, that estimates the fitted velocity at each pixel from Gauss-Hermite fits to a 3D spectral data cube, together with its uncertainty obtained from a modified bootstrap procedure. We apply this new tool to neutral hydrogen spectra for 18 galaxies from the THINGS sample, and present new velocity maps with uncertainties. We propagate the estimated uncertainties in the velocity map into our previously-described model fitting tool DiskFit to derive new rotation curves. The uncertainties we obtain from these fits take into account not only the observational errors, but also uncertainties in the fitted systemic velocity, position of the rotation centre, inclination of the galaxy to the line of sight, and forced non-circular motion. They are therefore much better-defined than values that have previously been available. Our estimated uncertainties on the circular speeds differ from previous estimates by factors ranging up to of five, being smaller in some cases and larger in others. We conclude that kinematic models of well-resolved HI datasets vary widely in their precision and reliability, and therefore potentially in their value for comparisons with predictions from cosmological galaxy formation simulations., Comment: 31 pages, 22 figures, accepted to appear in MNRAS. A complete version with full resolution figures is posted at https://www.physics.rutgers.edu/~sellwood/makemap.pdf (v2 fixed typo in Vsys for N2403)

We present a study of the spiral responses in a stable disc galaxy model to co-orbiting perturbing masses that are evenly spaced around rings. The amplitudes of the responses, or wakes, are proportional to the masses of the perturbations, and we find that the response to a low-mass ring disperses when it is removed -- behaviour that is predicted by linear theory. Higher mass rings cause nonlinear changes through scattering at the major resonances, provoking instabilities that were absent before the scattering took place. The separate wake patterns from two rings orbiting at differing frequencies, produce a net response that is an apparently shearing spiral. When the rings have low mass, the evolution of the simulation is both qualitatively and quantitatively reproduced by linear superposition of the two separate responses. We argue that apparently shearing transient spirals in simulations result from the superposition of two or more steadily rotating patterns, each of which is best accounted for as a normal mode of the non-smooth disc., Comment: 14 pages, 19 figures and one animation. The animation may be viewed at http://www.physics.rutgers.edu/~sellwood/animated-fig.html

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., Comment: 18 pages, 17 figures, accepted to appear in MNRAS, comments welcome

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., Comment: Accepted to appear in MNRAS. 4 pages, 2 figures

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., Comment: 18 pages, 15 figures, minor revisions in response to a thoughtful referee report

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., Comment: 15 pages, 8 figures, to appear in MNRAS. A nymber of revisions from v1

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., Comment: Revised version accepted to appear in MNRAS. Some significant improvements. A full resolution version of Fig 4 is available from http://www.physics.rutgers.edu/~sellwood/mult_res.pdf

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., Comment: Accepted to appear in AJ (28 pages, 22 Figures)

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., Comment: 11 pages, 6 figures, accepted to appear in ApJ

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., Comment: To be published in "MeerKAT Science: On the Pathway to the SKA". Proceedings of Science. Workshop held 25-27 May, 2016 Stellenbosch, South Africa

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., Comment: 13 pages, 7 figures, accepted to appear in ApJ

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., Comment: 8 pages, 4 figures. Revised version submitted to ApJ

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., Comment: 8 pages, 8 figures, Conference proceeding for SALT Science Conference 2015

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., Comment: 3 pages, the package is available at http://www.physics.rutgers.edu/~spekkens/diskfit/

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., Comment: 7 pages, 4 figures, accepted to appear in MNRAS. Very minor changes in proofs

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., Comment: Revised version submitted to MNRAS. 14 pages, 7 figures - removed extra figure file

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., Comment: Revised version submitted to AJ. 14 pages, 11 figures

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., Comment: Minor corrections in proofs and updates to references, 16 pages, 16 figures

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., Comment: 2 pages

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., Comment: This is a brief summary, the full manual is at http://www.physics.rutgers.edu/~sellwood/manual.pdf and the source code is at http://www.physics.rutgers.edu/galaxy

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., Comment: Includes proof corrections. To appear in ApJ, 13 pages, 11 figures. A version with full resolution figures is available at http://www.physics.rutgers.edu/~sellwood/modes.pdf

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., Comment: 13 pages, 23 figures, accepted for publication in MNRAS

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., Comment: Minor revisions to become as it will appear in the journal

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., Comment: Accepted to ApJ Letters April 23, uses emulateapj.sty

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., Comment: Accepted for publication in MNRAS; 15 pages, 15 figures, 3 tables; The DiskFit code can be found at http://www.physics.rutgers.edu/~spekkens/diskfit/

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., Comment: 17 pages, 7 figures, accepted for publication in RAA

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., Comment: 19 pages, 1 figure, accepted for publication in RAA

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., Comment: Accepted to appear in ApJ, 8 pages, 8 figures

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., Comment: 11 pages, 10 figures, revised after referee's report, accepted to appear ApJ, minor textual revisions

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., Comment: 23 pages, 25 figures, accepted to appear in MNRAS

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., Comment: Included copy editor's corrections. Accepted to appear in MNRAS. 9 pages, 11 figures, uses mn2e LaTeX

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., Comment: Conference article from 1989 posted so as to allow access to the full text through ADS

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., Comment: Revised version submitted to MNRAS, 11 pages, 9 figures, uses MN macros

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., Comment: Chapter accepted to appear in Planets, Stars and Stellar Systems, vol 5, ed G. Gilmore. 32 pages, 17 figures. Includes minor corrections made in proofs. Uses emulateapj.sty

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., Comment: Invited review presented at "Evolution of Planetary and Stellar Systems", Prato, Italy. 6 pages, 7 figures, uses emulateapj

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., Comment: 12 pages, 10 figures to appear in MNRAS. Minor revisions from original version

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., Comment: 12 pages, 8 color figures, to appear in MNRAS. Version 2 of the software can be downloaded from http://www.physics.rutgers.edu/~spekkens/velfit/

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., Comment: To appear in MNRAS, 21 pages 30 figures. Minor changes from v1: we stress that divergent evolution is only temporary, and provide more evidence for chaos. A tar file (346MB) containing the starting coordinates of all the models shown in Figure 4, together with a README explanation is available at http://www.physics.rutgers.edu/~sellwood/SDmodels.tar

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., Comment: 7 pages, 6 figures, uses emulateapj

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., Comment: 5 pages, Conference Proceeding for IAU Symposium 255: Low-Metallicity Star Formaion: From the First Stars to Dwarf Galaxies

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., Comment: 12 pages, 5 figures, to appear in "The Galaxy Disk in Cosmological Context" eds Andersen, Bland-Hawthorn & Nordstrom