Uncovering CDM halo substructure with tidal streams

Models for the formation and growth of structure in a cold dark matter dominated universe predict that galaxy halos should contain significant substructure. Studies of the Milky Way, however, have yet to identify the expected few hundred sub-halos with masses greater than about 10^6 Msun. Here we propose a test for the presence of sub-halos in the halos of galaxies. We show that the structure of the tidal tails of ancient globular clusters is very sensitive to heating by repeated close encounters with the massive dark sub-halos. We discuss the detection of such an effect in the context of the next generation of astrometric missions, and conclude that it should be easily detectable with the GAIA dataset. The finding of a single extended cold stellar stream from a globular cluster would support alternative theories, such as self-interacting dark matter, that give rise to smoother halos.Comment: 7 pages, 7 figures, submitted to MNRA

Hierarchical formation of bulgeless galaxies II: Redistribution of angular momentum via galactic fountains

Within a fully cosmological hydrodynamical simulation, we form a galaxy which rotates at 140 km/s, and is characterised by two loose spiral arms and a bar, indicative of a Hubble Type SBc/d galaxy. We show that our simulated galaxy has no classical bulge, with a pure disc profile at z=1, well after the major merging activity has ended. A long-lived bar subsequently forms, resulting in the formation of a secularly-formed "pseudo" bulge, with the final bulge-to-total light ratio B/T=0.21. We show that the majority of gas which loses angular momentum and falls to the central region of the galaxy during the merging epoch is blown back into the hot halo, with much of it returning later to form stars in the disc. We propose that this mechanism of redistribution of angular momentum via a galactic fountain, when coupled with the results from our previous study which showed why gas outflows are biased to have low angular momentum, can solve the angular momentum/bulgeless disc problem of the cold dark matter paradigm.Comment: 9 Pages, 10 Figures, accepted MNRAS version. Comments welcom

An Advanced, Three-Dimensional Plotting Library for Astronomy

We present a new, three-dimensional (3D) plotting library with advanced features, and support for standard and enhanced display devices. The library - S2PLOT - is written in C and can be used by C, C++ and FORTRAN programs on GNU/Linux and Apple/OSX systems. S2PLOT draws objects in a 3D (x,y,z) Cartesian space and the user interactively controls how this space is rendered at run time. With a PGPLOT inspired interface, S2PLOT provides astronomers with elegant techniques for displaying and exploring 3D data sets directly from their program code, and the potential to use stereoscopic and dome display devices. The S2PLOT architecture supports dynamic geometry and can be used to plot time-evolving data sets, such as might be produced by simulation codes. In this paper, we introduce S2PLOT to the astronomical community, describe its potential applications, and present some example uses of the library.Comment: 12 pages, 10 eps figures (higher resolution versions available from http://astronomy.swin.edu.au/s2plot/paperfigures). The S2PLOT library is available for download from http://astronomy.swin.edu.au/s2plo

Resolving the Structure of Cold Dark Matter Halos

We examine the effects of mass resolution and force softening on the density profiles of cold dark matter halos that form within cosmological N-body simulations. As we increase the mass and force resolution, we resolve progenitor halos that collapse at higher redshifts and have very high densities. At our highest resolution we have nearly 3 million particles within the virial radius, several orders of magnitude more than previously used and we can resolve more than one thousand surviving dark matter halos within this single virialised system. The halo profiles become steeper in the central regions and we may not have achieved convergence to a unique slope within the inner 10% of the virialised region. Results from two very high resolution halo simulations yield steep inner density profiles, $\rho(r)\sim r^{-1.4}$. The abundance and properties of arcs formed within this potential will be different from calculations based on lower resolution simulations. The kinematics of disks within such a steep potential may prove problematic for the CDM model when compared with the observed properties of halos on galactic scales.Comment: Final version, to be published in the ApJLetter

High-Redshift Galaxies: Their Predicted Size and Surface Brightness Distributions and Their Gravitational Lensing Probability

Direct observations of the first generation of luminous objects will likely become feasible over the next decade. The advent of the Next Generation Space Telescope (NGST) will allow imaging of numerous galaxies and mini-quasars at redshifts z>5. We apply semi-analytic models of structure formation to estimate the rate of multiple imaging of these sources by intervening gravitational lenses. Popular CDM models for galaxy formation yield a lensing optical depth of about 1% for sources at redshift 10. The expected slope of the luminosity function of the early sources implies an additional magnification bias of about 5, bringing the fraction of lensed sources at z=10 to about 5%. We estimate the angular size distribution of high-redshift disk galaxies and find that most of them are more extended than the resolution limit of NGST, roughly 0.06 arcseconds. We also show that there is only a modest redshift evolution in the mean surface brightness of galaxies at z>2. The expected increase by 1-2 orders of magnitude in the number of resolved sources on the sky, due to observations with NGST, will dramatically improve upon the statistical significance of existing weak lensing measurements. We show that, despite this increase in the density of sources, confusion noise from z>2 galaxies is expected to be small for NGST observations.Comment: 27 pages, 8 PostScript figures (of which two are new), revised version accepted for Ap

Forming Disk Galaxies in Lambda CDM Simulations

We used fully cosmological, high resolution N-body + SPH simulations to follow the formation of disk galaxies with rotational velocities between 135 and 270 km/sec in a Lambda CDM universe. The simulations include gas cooling, star formation, the effects of a uniform UV background and a physically motivated description of feedback from supernovae. The host dark matter halos have a spin and last major merger redshift typical of galaxy sized halos as measured in recent large scale N--Body simulations. The simulated galaxies form rotationally supported disks with realistic exponential scale lengths and fall on both the I-band and baryonic Tully Fisher relations. An extended stellar disk forms inside the Milky Way sized halo immediately after the last major merger. The combination of UV background and SN feedback drastically reduces the number of visible satellites orbiting inside a Milky Way sized halo, bringing it in fair agreement with observations. Our simulations predict that the average age of a primary galaxy's stellar population decreases with mass, because feedback delays star formation in less massive galaxies. Galaxies have stellar masses and current star formation rates as a function of total mass that are in good agreement with observational data. We discuss how both high mass and force resolution and a realistic description of star formation and feedback are important ingredients to match the observed properties of galaxies.Comment: Revised version after the referee's comments. Conclusions unchanged. 2 new plots. MNRAS in press. 20 plots. 21 page

Two-Stream Instability of Counter-Rotating Galaxies

The present study of the two-stream instability in stellar disks with counter-rotating components of stars and/or gas is stimulated by recently discovered counter-rotating spiral and S0 galaxies. Strong linear two-stream instability of tightly-wrapped spiral waves is found for one and two-armed waves with the pattern angular speed of the unstable waves always intermediate between the angular speed of the co-rotating matter ($+\Omega$) and that of the counter-rotating matter ($-\Omega$). The instability arises from the interaction of positive and negative energy modes in the co- and counter-rotating components. The unstable waves are in general convective - they move in radius and radial wavenumber space - with the result that amplification of the advected wave is more important than the local growth rate. For a galaxy of co-rotating stars and counter-rotating stars of mass-fraction $\xi_* < {1\over 2}$, or of counter-rotating gas of mass-fraction $\xi_g < {1\over 2}$, the largest amplification is usually for the one-armed leading waves (with respect to the co-rotating stars). For the case of both counter-rotating stars and gas, the largest amplifications are for $\xi_*+\xi_g \approx {1\over 2}$, also for one-armed leading waves. The two-armed trailing waves usually have smaller amplifications. The growth rates and amplifications all decrease as the velocity spreads of the stars and/or gas increase. It is suggested that the spiral waves can provide an effective viscosity for the gas causing its accretion.Comment: 14 pages, submitted to ApJ. One table and 17 figures can be obtained by sending address to R. Lovelace at [email protected]

The Expected Mass Function for Low Mass Galaxies in a CDM Cosmology: Is There a Problem?

It is well known that the mass function for_halos_ in CDM cosmology is a relatively steep power law for low masses, possibly too steep to be consistent with observations. But how steep is the_galaxy_ mass function? We have analyzed the stellar and gas mass functions of the first massive luminous objects formed in a \Lambda CDM universe, as calculated in the numerical simulation described in Gnedin (2000ab). We found that while the dark matter mass function is steep, the stellar and gas mass functions are flatter for low mass objects. The stellar mass function is consistently flat at the low mass end. Moreover, while the gas mass function follows the dark matter mass function until reionization at z~7, between z=7 and z=4, the gas mass function also flattens considerably at the low mass end. At z=4, the gas and stellar mass functions are fit by a Schechter function with \alpha ~ -1.2 +/- 0.1, significantly shallower than the dark matter halo mass function and consistent with some recent observations. The baryonic mass functions are shallower because (a) the dark matter halo mass function is consistent with the Press-Schechter formulation at low masses n(M) M^-2 and (b) heating/cooling and ionization processes appear to cause baryons to collect in halos with the relationship M_b M_d^4 at low masses. Combining (a) and (b) gives n(M_b) M_b^-5/4, comparable to the simulation results. Thus, the well known observational fact that low mass galaxies are underabundant as compared to expectations from numerical dark matter simulations or Press-Schechter modeling of CDM universes emerges naturally from these results, implying that perhaps no ``new physics'' beyond the standard model is needed.Comment: Submitted to ApJ, 17 pages including 6 figure

Axiomatic approach to radiation reaction of scalar point particles in curved spacetime

Several different methods have recently been proposed for calculating the motion of a point particle coupled to a linearized gravitational field on a curved background. These proposals are motivated by the hope that the point particle system will accurately model certain astrophysical systems which are promising candidates for observation by the new generation of gravitational wave detectors. Because of its mathematical simplicity, the analogous system consisting of a point particle coupled to a scalar field provides a useful context in which to investigate these proposed methods. In this paper, we generalize the axiomatic approach of Quinn and Wald in order to produce a general expression for the self force on a point particle coupled to a scalar field following an arbitrary trajectory on a curved background. Our equation includes the leading order effects of the particle's own fields, commonly referred to as ``self force'' or ``radiation reaction'' effects. We then explore the equations of motion which follow from this expression in the absence of non-scalar forces.Comment: 17 pages, 1 figur