Faint counts as a function of morphological type in a hierarchical merger model

The unprecedented resolution of the refurbished Wide Field and Planetary Camera 2 (WFPC2) on the Hubble Space Telescope (HST) has led to major advances in our understanding of galaxy formation. The high image quality in the Medium Deep Survey and Hubble Deep Field has made it possible, for the first time, to classify faint distant galaxies according to morphological type. These observations have revealed a large population of galaxies classed as irregulars or which show signs of recent merger activity. Their abundance rises steeply with apparent magnitude, providing a likely explanation for the large number of blue galaxies seen at faint magnitudes. We demonstrate that such a population arises naturally in a model in which structure forms hierarchically and which is dynamically dominated by cold dark matter. The number counts of irregular, spiral and elliptical galaxies as a function of magnitude seen in the HST data are well reproduced in this model.We present detailed predictions for the outcome of spectroscopic follow-up observations of the HST surveys. By measuring the redshift distributions of faint galaxies of different morphological types, these programmes will provide a test of the hierarchical galaxy formation paradigm and might distinguish between models with different cosmological parameters.Comment: 5 pages, 3 postscript figures included. To be published as a Letter in Monthly Notices of the RAS. Postscript version available at http://star-www.dur.ac.uk/~cmb/counts.htm

GALAXY DYNAMICS IN CLUSTERS

We use high resolution simulations to study the formation and distribution of galaxies within a cluster which forms hierarchically. We follow both dark matter and baryonic gas which is subject to thermal pressure, shocks and radiative cooling. Galaxy formation is identified with the dissipative collapse of the gas into cold, compact knots. We examine two extreme representations of galaxies during subsequent cluster evolution --- one purely gaseous and the other purely stellar. The results are quite sensitive to this choice. Gas-galaxies merge efficiently with a dominant central object while star-galaxies merge less frequently. Thus, simulations in which galaxies remain gaseous appear to suffer an ``overmerging'' problem, but this problem is much less severe if the gas is allowed to turn into stars. We compare the kinematics of the galaxy population in these two representations to that of dark halos and of the underlying dark matter distribution. Galaxies in the stellar representation are positively biased (\ie over-represented in the cluster) both by number and by mass fraction. Both representations predict the galaxies to be more centrally concentrated than the dark matter, whereas the dark halo population is more extended. A modest velocity bias also exists in both representations, with the largest effect, $\sigma_{gal}/\sigma_{DM} \simeq 0.7$, found for the more massive star-galaxies. Phase diagrams show that the galaxy population has a substantial net inflow in the gas representation, while in the stellar case it is roughly in hydrostatic equilibrium. Virial mass estimators can underestimate the true cluster mass by up to a factor of 5. The discrepancy is largest if only the most massive galaxies are used, reflecting significant mass segregation.Comment: 30 pages, self-unpacking (via uufiles) postscript file without figures. Eighteen figures (and slick color version of figure 3) and entire paper available at ftp://oahu.physics.lsa.umich.edu/groups/astro/fews Total size of paper with figures is ~9.0 Mb uncompressed. Submitted to Ap.J

Dynamical Evolution of Globular Clusters in Hierarchical Cosmology

We probe the evolution of globular clusters that could form in giant molecular clouds within high-redshift galaxies. Numerical simulations demonstrate that the large and dense enough gas clouds assemble naturally in current hierarchical models of galaxy formation. These clouds are enriched with heavy elements from earlier stars and could produce star clusters in a similar way to nearby molecular clouds. The masses and sizes of the model clusters are in excellent agreement with the observations of young massive clusters. Do these model clusters evolve into globular clusters that we see in our and external galaxies? In order to study their dynamical evolution, we calculate the orbits of model clusters using the outputs of the cosmological simulation of a Milky Way-sized galaxy. We find that at present the orbits are isotropic in the inner 50 kpc of the Galaxy and preferentially radial at larger distances. All clusters located outside 10 kpc from the center formed in the now-disrupted satellite galaxies. The spatial distribution of model clusters is spheroidal, with a power-law density profile consistent with observations. The combination of two-body scattering, tidal shocks, and stellar evolution results in the evolution of the cluster mass function from an initial power law to the observed log-normal distribution.Comment: 5 pages, proceedings of IAU 246 "Dynamical Evolution of Dense Stellar Systems", eds. Vesperini, Giersz, Sill

Planetary Nebulae Kinematics in M31

We present kinematics of 135 planetary nebulae in M31 from a survey covering 3.9 square degrees and extending out to 15 kpc from the southwest major axis and more than 20 kpc along the minor axis. The majority of our sample, even well outside the disk, shows significant rotational support (mean line-of-sight velocity 116 km/s). We argue that these PN belong to the outer part of M31's large de Vaucouleurs bulge. Only five PN have velocities clearly inconsistent with this fast rotating bulge. All five may belong to tidal streams in M31's outer halo. One is projected on the Northern Spur, and is counter-rotating with respect to the disk there. Two are projected along the major axis at X=-10 kpc and have M32-like velocities; they could be debris from that galaxy. The remaining two halo PN are located near the center of the galaxy and their velocities follow the gradient found by Ibata et al. (2004), implying that these PN could belong to the Southern Stream. If M31 has a non-rotating, pressure-supported halo, we have yet to find it, and it must be a very minor component of the galaxy.Comment: accepted to ApJ; main body of paper is 36 pages, including 14 figure

Physical Bias of Galaxies From Large-Scale Hydrodynamic Simulations

We analyze a new large-scale ($100h^{-1}$Mpc) numerical hydrodynamic simulation of the popular $\Lambda$CDM cosmological model, including in our treatment dark matter, gas and star-formation, on the basis of standard physical processes. The method, applied with a numerical resolution of $<200h^{-1}$kpc (which is still quite coarse for following individual galaxies, especially in dense regions), attempts to estimate where and when galaxies form. We then compare the smoothed galaxy distribution with the smoothed mass distribution to determine the "bias" defined as $b\equiv (\delta M/M)_{gal}/(\delta M/M)_{total}$ on scales large compared with the code numerical resolution (on the basis of resolution tests given in the appendix of this paper). We find that (holding all variables constant except the quoted one) bias increases with decreasing scale, with increasing galactic age or metallicity and with increasing redshift of observations. At the $8h^{-1}$Mpc fiducial comoving scale bias (for bright regions) is 1.35 at $z=0$ reaching to 3.6 at $z=3$, both numbers being consistent with extant observations. We also find that $(10-20)h^{-1}$Mpc voids in the distribution of luminous objects are as observed (i.e., observed voids are not an argument against CDM-like models) and finally that the younger systems should show a colder Hubble flow than do the early type galaxies (a testable proposition). Surprisingly, little evolution is found in the amplitude of the smoothed galaxy-galaxy correlation function (as a function of {\it comoving} separation). Testing this prediction vs observations will allow a comparison between this work and that of Kauffmann et al which is based on a different physical modelingmethod.Comment: in press, ApJ, 26 latex pages plus 7 fig

A G1-like globular cluster in NGC 1023

The structure of a very bright (MV = -10.9) globular cluster in NGC 1023 is analyzed on two sets of images taken with the Hubble Space Telescope. From careful modeling of King profile fits to the cluster image, a core radius of 0.55+/-0.1 pc, effective radius 3.7+/-0.3 pc and a central V-band surface brightness of 12.9+/-0.5 mag / square arcsec are derived. This makes the cluster much more compact than Omega Cen, but very similar to the brightest globular cluster in M31, G1 = Mayall II. The cluster in NGC 1023 appears to be very highly flattened with an ellipticity of about 0.37, even higher than for Omega Cen and G1, and similar to the most flattened clusters in the Large Magellanic Cloud.Comment: 14 pages, 3 figures, 1 table. Accepted for AJ, Oct 200

Chemo-dynamical evolution of Globular Cluster Systems

We studied the relation between the ratio of rotational velocity to velocity dispersion and the metallicity (/\sigma_{v}-metallicity relation) of globular cluster systems (GCS) of disk galaxies by comparing the relation predicted from simple chemo-dynamical models for the formation and evolution of disk galaxies with the observed kinematical and chemical properties of their GCSs. We conclude that proto disk galaxies underwent a slow initial collapse that was followed by a rapid contraction and derive that the ratio of the initial collapse time scale to the active star formation time scale is \sim 6 for our Galaxy and \sim 15 for M31. The fundamental formation process of disk galaxies was simulated based on simple chemo-dynamical models assuming the conservation of their angular momentum. We suggest that there is a typical universal pattern in the /\sigma_{v}-metallicity relation of the GCS of disk galaxies. This picture is supported by the observed properties of GCSs in the Galaxy and in M31. This relation would deviate from the universal pattern, however, if large-scale merging events took major role in chemo-dynamical evolution of galaxies and will reflect the epoch of such merging events. We discuss the properties of the GCS of M81 and suggest the presence of past major merging event.Comment: 25 pages, 8 figures, Accepted for publication in the Astrophysical Journa

A Comparison of Semi-Analytic and Smoothed Particle Hydrodynamics Galaxy Formation

We compare the statistical properties of galaxies found in two different models of hierarchical galaxy formation: the semi-analytic model of Cole et al. and the smoothed particle hydrodynamics (SPH) simulations of Pearce et al. Using a `stripped-down' version of the semi-analytic model which mimics the resolution of the SPH simulations and excludes physical processes not included in them, we find that the two models produce an ensemble of galaxies with remarkably similar properties, although there are some differences in the gas cooling rates and in the number of galaxies that populate halos of different mass. The full semi-analytic model, which has effectively no resolution limit and includes a treatment of star formation and supernovae feedback, produces somewhat different (but readily understandable) results. Agreement is particularly good for the present-day global fractions of hot gas, cold dense (i.e. galactic) gas and uncollapsed gas, for which the SPH and stripped-down semi-analytic calculations differ by at most 25%. In the most massive halos, the stripped-down semi-analytic model predicts, on the whole, up to 50% less gas in galaxies than is seen in the SPH simulations. The two techniques apportion this cold gas somewhat differently amongst galaxies in a given halo. This difference can be tracked down to the greater cooling rate in massive halos in the SPH simulation compared to the semi-analytic model. (abridged)Comment: 19 pages, 13 figures, to appear in MNRAS. Significantly extended to explore galaxy progenitor distributions and behaviour of models at high redshift

The construction of non-spherical models of quasi-relaxed stellar systems

Spherical models of collisionless but quasi-relaxed stellar systems have long been studied as a natural framework for the description of globular clusters. Here we consider the construction of self-consistent models under the same physical conditions, but including explicitly the ingredients that lead to departures from spherical symmetry. In particular, we focus on the effects of the tidal field associated with the hosting galaxy. We then take a stellar system on a circular orbit inside a galaxy represented as a "frozen" external field. The equilibrium distribution function is obtained from the one describing the spherical case by replacing the energy integral with the relevant Jacobi integral in the presence of the external tidal field. Then the construction of the model requires the investigation of a singular perturbation problem for an elliptic partial differential equation with a free boundary, for which we provide a method of solution to any desired order, with explicit solutions to two orders. We outline the relevant parameter space, thus opening the way to a systematic study of the properties of a two-parameter family of physically justified non-spherical models of quasi-relaxed stellar systems. The general method developed here can also be used to construct models for which the non-spherical shape is due to internal rotation. Eventually, the models will be a useful tool to investigate whether the shapes of globular clusters are primarily determined by internal rotation, by external tides, or by pressure anisotropy.Comment: AASTeX v5.2, 37 pages with 2 figures, accepted for publication in The Astrophysical Journa