We test and improve the numerical schemes in our smoothed particle hydrodynamics (SPH) code for cosmological simulations, including the pressure-entropy formulation (PESPH), a time-dependent artificial viscosity, a refined timestep criterion, and metal-line cooling that accounts for photoionisation in the presence of a recently refined Haardt \& Madau (2012) model of the ionising background. The PESPH algorithm effectively removes the artificial surface tension present in the traditional SPH formulation, and in our test simulations it produces better qualitative agreement with mesh-code results for Kelvin-Helmholtz instability and cold cloud disruption. Using a set of cosmological simulations, we examine many of the quantities we have studied in previous work. Results for galaxy stellar and HI mass functions, star formation histories, galaxy scaling relations, and statistics of the Ly$\alpha$ forest are robust to the changes in numerics and microphysics. As in our previous simulations, cold gas accretion dominates the growth of high-redshift galaxies and of low mass galaxies at low redshift, and recycling of winds dominates the growth of massive galaxies at low redshift. However, the PESPH simulation removes spurious cold clumps seen in our earlier simulations, and the accretion rate of hot gas increases by up to an order of magnitude at some redshifts. The new numerical model also influences the distribution of metals among gas phases, leading to considerable differences in the statistics of some metal absorption lines, most notably NeVIII.

Many phenomenologically successful cosmological galaxy formation simulations employ kinetic winds to model galactic outflows, a crucial ingredient in obtaining predictions that agree with various observations. Yet systematic studies of how variations in kinetic wind scalings might alter observable galaxy properties are rare. Here we employ GADGET-3 simulations to study how the baryon cycle, stellar mass function, and other galaxy and CGM predictions vary as a function of the assumed outflow speed $v_w$ and the scaling of the mass loading factor $\eta$ with velocity dispersion $\sigma$. We design our fiducial model to reproduce the measured wind properties at 25% of the virial radius from the Feedback In Realistic Environments (FIRE) simulations. We find that a strong dependence of $\eta \sim \sigma^5$ in low mass haloes with $\sigma < 106\ \mathrm{km\ s^{-1}}$ is required to match the faint end of the stellar mass functions at $z > 1$. The wind speed also has a major impact, with faster winds significantly reducing wind recycling and heating more halo gas. Both effects result in less stellar mass growth in massive haloes and impact high ionization absorption in halo gas. We cannot simultaneously match the stellar content at $z=2$ and $z=0$ within a single model, suggesting that an additional feedback source such as AGN might be required in massive galaxies at lower redshifts, but the amount needed depends strongly on assumptions regarding the outflow properties. We run a 50 $\mathrm{Mpc/h}$, $2\times576^3$ simulation with our fiducial parameters and show that it matches a range of star-forming galaxy properties at $z\sim0-2$. In closing, the results from simulations of galaxy formation are much more sensitive to small changes in the feedback implementation than to the hydrodynamic technique.

The Great Andromeda Galaxy (M31) is the nexus of the near-far galaxy evolution connection and a principal data point for near-field cosmology. Due to its proximity (780 kpc), M31 can be resolved into individual stars like the Milky Way (MW). Unlike the MW, we have the advantage of a global view of M31, enabling M31 to be observed with techniques that also apply to more distant galaxies. Moreover, recent evidence suggests that M31 may have survived a major merger within the last several Gyr, shaping the morphology of its stellar halo and triggering a starburst, while leaving the stellar disk largely intact. The MW and M31 thus provide complementary opportunities for in-depth studies of the disks, halos, and satellites of L* galaxies. Our understanding of the M31 system will be transformed in the 2020s if they include wide field facilities for both photometry (HST-like sensitivity and resolution) and spectroscopy (10-m class telescope, >1 sq. deg. field, highly multiplexed, R~ 3000 to 6000). We focus here on the power of these facilities to constrain the past, present, and future merger history of M31, via chemo-dynamical analyses and star formation histories of phase-mixed stars accreted at early times, as well as stars in surviving tidal debris features, M31's extended disk, and intact satellite galaxies that will eventually be tidally incorporated into the halo. This will yield an unprecedented view of the hierarchical formation of the M31 system and the subhalos that built it into the L* galaxy we observe today.

The 3D velocities of M31 and M33 are important for understanding the evolution and cosmological context of the Local Group. Their most massive stars are detected by Gaia, and we use Data Release 2 (DR2) to determine the galaxy proper motions (PMs). We select galaxy members based on, e.g., parallax, PM, color-magnitude-diagram location, and local stellar density. The PM rotation of both galaxies is confidently detected, consistent with the known line-of-sight rotation curves: $V_{\rm rot} = -206\pm86$ km s$^{-1}$ (counter-clockwise) for M31, and $V_{\rm rot} = 80\pm52$ km s$^{-1}$ (clockwise) for M33. We measure the center-of-mass PM of each galaxy relative to surrounding background quasars in DR2. This yields that $({\mu}_{\alpha*},{\mu}_{\delta})$ equals $(65 \pm 18 , -57 \pm 15)$ $\mu$as yr$^{-1}$ for M31, and $(31 \pm 19 , -29 \pm 16)$ $\mu$as yr$^{-1}$ for M33. In addition to the listed random errors, each component has an additional residual systematic error of 16 $\mu$as yr$^{-1}$. These results are consistent at 0.8$\sigma$ and 1.0$\sigma$ with the (2 and 3 times higher-accuracy) measurements already available from Hubble Space Telescope (HST) optical imaging and VLBA water maser observations, respectively. This lends confidence that all these measurements are robust. The new results imply that the M31 orbit towards the Milky Way is somewhat less radial than previously inferred, $V_{\rm tan, DR2+HST} = 57^{+35}_{-31}$ km s$^{-1}$, and strengthen arguments that M33 may be on its first infall into M31. The results highlight the future potential of Gaia for PM studies beyond the Milky Way satellite system.

The Orphan Stream is one of the most prominent tidal streams in the Galactic halo. Using data on red giants, RR Lyrae, and horizontal branch stars from Gaia and other surveys, we determine the proper motion of the Orphan Stream over a path of more than 90 degrees on the sky. We also provide updated tracks for the sky position, distance, and radial velocity of the stream. Our tracks in these latter dimensions mostly agree with previous results. However, there are significant corrections to the earlier distance and latitude tracks as the stream approaches the Galactic disk. Stream stars selected with three-dimensional kinematics display a very tight red giant sequence. Concordantly, we find that applying a proper motion cut removes the most metal-rich stars from earlier spectroscopic samples of stream stars, though a significant dispersion remains indicating a dwarf galaxy origin. The deceleration of the stream towards its leading end suggests a circular velocity of ~200 km/s at a galactocentric radius ~30 kpc, consistent with other independent evidence. However, the track of the stream departs significantly from an orbit; the spatial track does not point along the same direction as the velocity vector, and it exhibits a lateral wiggle that is unlikely to match any reasonable orbit. The low metallicity and small dispersion of the stream in the various coordinates point to a progenitor with a relatively low dynamical mass ~$10^8$ Msun.

Line-of-sight kinematic studies indicate that many Galactic globular clusters have a significant degree of internal rotation. However, three-dimensional kinematics from a combination of proper motions and line-of-sight velocities are needed to unveil the role of angular momentum in the formation and evolution of these old stellar systems. Here we present the first quantitative study of internal rotation on the plane-of-the-sky for a large sample of globular clusters using proper motions from Gaia DR2. We detect signatures of rotation in the tangential component of proper motions for 11 out of 51 clusters at a $>$3-sigma confidence level, confirming the detection reported in Gaia collaboration et al. (2018) for 8 clusters, and additionally identify 11 GCs with a 2-sigma rotation detection. For the clusters with a detected global rotation, we construct the two-dimensional rotation maps and proper motion rotation curves, and we assess the relevance of rotation with respect to random motions ($V/\sigma\sim0.08-0.51$). We find evidence of a correlation between the degree of internal rotation and relaxation time, highlighting the importance of long-term dynamical evolution in shaping the clusters current properties. This is a strong indication that angular momentum must have played a fundamental role in the earliest phases of cluster formation. Finally, exploiting the spatial information of the rotation maps and a comparison with line-of-sight data, we provide an estimate of the inclination of the rotation axis for a subset of 8 clusters. Our work demonstrates the potential of Gaia data for internal kinematic studies of globular clusters and provides the first step to reconstruct their intrinsic three-dimensional structure.

We present Hubble Space Telescope (HST) absolute proper motion (PM) measurements for 20 globular clusters (GCs) in the Milky Way (MW) halo at Galactocentric distances $R_{\rm GC} \approx 10-100$ kpc, with median per-coordinate PM uncertainty 0.06 mas yr$^{-1}$. Young and old halo GCs do not show systematic differences in their 3D Galactocentric velocities, derived from combination with existing line-of-sight velocities. We confirm the association of Arp 2, Pal 12, Terzan 7, and Terzan 8 with the Sagittarius (Sgr) stream. These clusters and NGC 6101 have tangential velocity $V_{\rm tan} > 290$ km s$^{-1}$, whereas all other clusters have $V_{\rm tan} < 200$ km s$^{-1}$. NGC 2419, the most distant GC in our sample, is also likely associated with the Sgr stream, whereas NGC 4147, NGC 5024, and NGC 5053 definitely are not. We use the distribution of orbital parameters derived using the 3D velocities to separate halo GCs that either formed within the MW or were accreted. We also assess the specific formation history of e.g. Pyxis and Terzan 8. We constrain the MW mass via an estimator that considers the full 6D phase-space information for 16 of the GCs from $R_{\rm GC} = 10$ to 40 kpc. The velocity dispersion anisotropy parameter $\beta = 0.609^{+0.130}_{-0.229}$. The enclosed mass $M (<39.5 \rm{kpc}) = 0.61^{+0.18}_{-0.12} \times 10^{12}$ M$_{\odot}$, and the virial mass $M_\rm{vir} = 2.05^{+0.97}_{-0.79} \times 10^{12}$ M$_{\odot}$, are consistent with, but on the high side among recent mass estimates in the literature.

Recent maps of the halo using RR Lyrae from Pan-STARRS1 have clearly depicted the spatial structure of the Sagittarius stream. These maps show the leading and trailing stream apocenters differ in galactocentric radius by a factor of two, and also resolve substructure in the stream at these apocenters. Here we present dynamical models that reproduce these features of the stream in simple Galactic potentials. We find that debris at the apocenters must be dynamically young, in the sense of being stripped off in the last two pericentric passages, while the Sagittarius dwarf is currently experiencing a third passage. The ratio of apocenters is sensitive to both dynamical friction and the outer slope of the Galactic rotation curve. These dependences can be understood with simple regularities connecting the apocentric radii, circular velocities, and orbital period of the progenitor. The effect of dynamical friction on the stream can be constrained using substructure within the leading apocenter. Our ensembles of models are not intended as statistically proper fits to the stream. Nevertheless, out of the range of models we consider, we consistently find the mass within 100 kpc to be $\sim 7 \times 10^{11} \, M_{\odot}$, with a nearly flat rotation curve between 50 and 100 kpc. This points to a more extended Galactic halo than assumed in some current models. As in previous work, we find prolate or triaxial halos ease agreement with the track of the leading stream. We display the behavior of our models in various observational spaces and characterize the substructure expected within the stream. In particular, the young trailing stream visible near trailing apocenter should exhibit a tight trend of velocity with distance separate from the older debris, and we suggest that this will serve as an especially useful probe of the outer Galactic potential.

Stellar tidal streams provide an opportunity to study the motion and structure of the disrupting galaxy as well as the gravitational potential of its host. Streams around the Milky Way are especially promising as phase space positions of individual stars will be measured by ongoing or upcoming surveys. Nevertheless, it remains a challenge to accurately assess distances to stars farther than 10 kpc from the Sun, where we have the poorest knowledge of the Galaxy's mass distribution. To address this we present observations of 32 candidate RR Lyrae stars in the Orphan tidal stream taken as part of the Spitzer Merger History and Shape of the Galactic Halo (SMHASH) program. The extremely tight correlation between the periods, luminosities, and metallicities of RR Lyrae variable stars in the Spitzer IRAC $\mathrm{3.6 \mu m}$ band allows the determination of precise distances to individual stars; the median statistical distance uncertainty to each RR Lyrae star is $2.5\%$. By fitting orbits in an example potential we obtain an upper limit on the mass of the Milky Way interior to 60 kpc of $\mathrm{5.6_{-1.1}^{+1.2}\times 10^{11}\ M_\odot}$, bringing estimates based on the Orphan Stream in line with those using other tracers. The SMHASH data also resolve the stream in line--of--sight depth, allowing a new perspective on the internal structure of the disrupted dwarf galaxy. Comparing with N--body models we find that the progenitor had an initial dark halo mass of approximately $\mathrm{3.2 \times 10^{9}\ M_\odot}$, placing the Orphan Stream's progenitor amongst the classical dwarf spheroidals.

We report the detection of spatially distinct stellar density features near the apocenters of the Sagittarius (Sgr) stream's main leading and trailing arm. These features are clearly visible in a high-fidelity stellar halo map that is based on RR Lyrae from Pan-STARRS1: there is a plume of stars 10 kpc beyond the apocenter of the leading arm, and there is a "spur" extending to 130 kpc, almost 30 kpc beyond the previously detected apocenter of the trailing arm. Such apocenter substructure is qualitatively expected in any Sgr stream model, as stars stripped from the progenitor at different pericenter passages become spatially separated there. The morphology of these new Sgr stream substructures could provide much-needed new clues and constraints for modeling the Sgr system, including the level of dynamical friction that Sgr has experienced. We also report the discovery of a new, presumably unrelated halo substructure at 80 kpc from the Sun and $10^\circ$ from the Sgr orbital plane, which we dub the Outer Virgo Overdensity.

We present a comprehensive analysis of the structural properties and luminosities of the 23 dwarf spheroidal galaxies that fall within the footprint of the Pan-Andromeda Archaeological Survey (PAndAS). These dwarf galaxies represent the large majority of Andromeda's known satellite dwarf galaxies and cover a wide range in luminosity ($-11.6<M_V<-5.8$ or $10^{4.2}< L <10^{6.5} L_\odot$) and surface brightness ($25.1<\mu_0<29.3$ mag/arcsec$^2$). We confirm most previous measurements, but find And XIX to be significantly larger than before ($r_h=3065^{+1065}_{-935}$ pc, $M_V=-10.1^{+0.8}_{-0.4}$) and cannot derive parameters for And XXVII as it is likely not a bound stellar system. We also significantly revise downward the luminosities of And~XV and And~XVI, which are now $M_V\sim-7.5$ or $L\sim10^5 L_\odot$. Finally, we provide the first detailed analysis of Cas II/And XXX, a fairly faint system ($M_V=-8.0^{+0.4}_{-0.3}$) of typical size ($r_h=270\pm50$ pc), located in close proximity to the two bright elliptical dwarf galaxies NGC 147 & 185. Combined with the set of homogeneous distances published in an earlier contribution, our analysis dutifully tracks all relevant sources of uncertainty in the determination of the properties of the dwarf galaxies from the PAndAS photometric catalogue. We further publish the posterior probability distribution functions of all the parameters we fit for in the form of MCMC chains available online; these inputs should be used in any analysis that aims to remain truthful to the data and properly account for covariance between parameters.

The thermal and chemical properties of the hot diffuse intragroup medium (IGrM) provide important constraints on the feedback processes associated with massive galaxy formation and evolution. Here we explore these constraints via a detailed analysis of the global properties of simulated z<3 galaxy groups from a cosmological simulation that includes a well-constrained prescription for stellar/supernovae-powered galactic outflows but no AGN feedback. Our aims are to (a) establish a baseline against which we will compare future models; (b) identify model successes due to stellar/supernovae-powered outflows; and (c) pinpoint features that signal the need for, and constrain the nature of, AGN feedback. Our simulation successfully reproduces key observed z=0 group IGrM properties, including the various X-ray Lx - Tx - entropy scaling relations, for all but the most massive groups. The z<1 redshift evolution of these also agree with observations. Contrary to expectations, the simulated groups' IGrM does not suffer catastrophic cooling. Yet, the z=0 group stellar mass is ~ 2X too large. This is due to the build-up of cold gas in the massive galaxies before they are incorporated inside groups. This in turn indicates that other feedback mechanisms must activate in real galaxies once their stellar masses grow to a few X 10^{10} M_sun. We show that these must be powerful enough to expel a significant fraction of the gas from the galactic halos. Gentle maintenance-mode (quenching) AGN feedback, as seen in galaxy clusters, will not do. Just as importantly, we find that the stellar/supernovae-powered winds are essential for understanding the IGrM metal abundances. Our simulation is able to reproduce the observed relationship between the global IGrM iron and silicon abundance and the group X-ray temperature, and these results ought to be relatively insensitive to the addition of AGN feedback.

In this paper we discuss a method for the generation of mock tidal streams. Using an ensemble of simulations in an isochrone potential where the actions and frequencies are known, we derive an empirical recipe for the evolving satellite mass and the corresponding mass loss rate, and the ejection conditions of the stream material. The resulting stream can then be quickly generated either with direct orbital integration, or by using the action-angle formalism. The model naturally produces streaky features within the stream. These are formed due to the radial oscillation of the progenitor and the bursts of stars emitted near pericenter, rather than clumping at particular oscillation phases as sometimes suggested. When detectable, these streaky features are a reliable diagnostic for the stream's direction of motion and encode other information on the progenitor and its orbit. We show several tests of the recipe in alternate potentials, including a case with a chaotic progenitor orbit which displays a marked effect on the width of the stream. Although the specific ejection recipe may need adjusting when elements such as the orbit or satellite density profile are changed significantly, our examples suggest that model tidal streams can be quickly and accurately generated by models of this general type for use in Bayesian sampling.

We exploit data from the Pan-Andromeda Archaeological Survey (PAndAS) to study the extended structures of M31's dwarf elliptical companions, NGC147 and NGC185. Our wide-field, homogeneous photometry allows to construct deep colour-magnitude diagrams (CMDs) which reach down to $\sim3$ mag below the red giant branch (RGB) tip. We trace the stellar components of the galaxies to surface brightness of $\mu_g \sim 32$ mag arcsec$^{-2}$ and show they have much larger extents ($\sim5$ kpc radii) than previously recognised. While NGC185 retains a regular shape in its peripheral regions, NGC147 exhibits pronounced isophotal twisting due to the emergence of symmetric tidal tails. We fit single Sersic models to composite surface brightness profiles constructed from diffuse light and star counts and find that NGC147 has an effective radius almost 3 times that of NGC185. In both cases, the effective radii that we calculate are larger by a factor of $\sim2$ compared to most literature values. We also calculate revised total magnitudes of $M_g=-15.36\pm0.04$ for NGC185 and $M_g=-16.36\pm0.04$ for NGC147. Using photometric metallicities computed for RGB stars, we find NGC185 to exhibit a metallicity gradient of [Fe/H]$\sim-0.15$ dex/kpc over the radial range 0.125 to 0.5 deg. On the other hand, NGC147 exhibits almost no metallicity gradient, $\sim-0.02$ dex/kpc from 0.2 to 0.6 deg. The differences in the structure and stellar populations in the outskirts of these systems suggest that tidal influences have played an important role in governing the evolution of NGC147.

We present a detailed kinematic analysis of the outer halo globular cluster (GC) system of M31. Our basis for this is a set of new spectroscopic observations for 78 clusters lying at projected distances between Rproj ~20-140 kpc from the M31 centre. These are largely drawn from the recent PAndAS globular cluster catalogue; 63 of our targets have no previous velocity data. Via a Bayesian maximum likelihood analysis we find that GCs with Rproj > 30 kpc exhibit coherent rotation around the minor optical axis of M31, in the same direction as more centrally- located GCs, but with a smaller amplitude of 86+/-17 km s-1. There is also evidence that the velocity dispersion of the outer halo GC system decreases as a function of projected distance from the M31 centre, and that this relation can be well described by a power law of index ~ -0.5. The velocity dispersion profile of the outer halo GCs is quite similar to that of the halo stars, at least out to the radius up to which there is available information on the stellar kinematics. We detect and discuss various velocity correlations amongst subgroups of GCs that lie on stellar debris streams in the M31 halo. Many of these subgroups are dynamically cold, exhibiting internal velocity dispersions consistent with zero. Simple Monte Carlo experiments imply that such configurations are unlikely to form by chance, adding weight to the notion that a significant fraction of the outer halo GCs in M31 have been accreted alongside their parent dwarf galaxies. We also estimate the M31 mass within 200 kpc via the Tracer Mass Estimator, finding (1.2 - 1.6) +/- 0.2 10^{12}M_sun. This quantity is subject to additional systematic effects due to various limitations of the data, and assumptions built in into the TME. Finally, we discuss our results in the context of formation scenarios for the M31 halo.

We examine the statistics of the low-redshift Lyman-alpha forest from smoothed particle hydrodynamic simulations in light of recent improvements in the estimated evolution of the cosmic ultraviolet background (UVB) and recent observations from the Cosmic Origins Spectrograph (COS). We find that the value of the metagalactic photoionization rate required by our simulations to match the observed properties of the low-redshift Lyman-alpha forest is a factor of 5 larger than the value predicted by state-of-the art models for the evolution of this quantity. This mismatch results in the mean flux decrement of the Lyman-alpha forest being underpredicted by at least a factor of 2 (a 10-sigma discrepancy with observations) and a column density distribution of Lyman-alpha forest absorbers systematically and significantly elevated compared to observations over nearly two decades in column density. We examine potential resolutions to this mismatch and find that either conventional sources of ionizing photons (galaxies and quasars) must be significantly elevated relative to current observational estimates or our theoretical understanding of the low-redshift universe is in need of substantial revision.

We examine the star formation histories (SFHs) of galaxies in smoothed particle hydrodynamics (SPH) simulations, compare them to parametric models that are commonly used in fitting observed galaxy spectral energy distributions, and examine the efficacy of these parametric models as practical tools for recovering the physical parameters of galaxies. The commonly used tau-model, with SFR ~ exp(-t/tau), provides a poor match to the SFH of our SPH galaxies, with a mismatch between early and late star formation that leads to systematic errors in predicting colours and stellar mass-to-light ratios. A one-parameter lin-exp model, with SFR ~ t*exp(-t/tau), is much more successful on average, but it fails to match the late-time behavior of the bluest, most actively star-forming galaxies and the passive, "red and dead" galaxies. We introduce a 4-parameter model, which transitions from lin-exp to a linear ramp after a transition time, which describes our simulated galaxies very well. We test the ability of these parametrised models to recover (at z=0, 0.5, and 1) the stellar mass-to-light ratios, specific star formation rates, and stellar population ages from the galaxy colours, computed from the full SPH star formation histories using the FSPS code of Conroy et al. (2009). Fits with tau-models systematically overestimate M/L by ~ 0.2 dex, overestimate population ages by ~ 1-2 Gyr, and underestimate sSFR by ~ 0.05 dex. Fits with lin-exp are less biased on average, but the 4-parameter model yields the best results for the full range of galaxies. Marginalizing over the free parameters of the 4-parameter model leads to slightly larger statistical errors than 1-parameter fits but essentially removes all systematic biases, so this is our recommended procedure for fitting real galaxies.

We reveal the highly structured nature of the Milky Way stellar halo within the footprint of the PAndAS photometric survey from blue main sequence and main sequence turn-off stars. We map no fewer than five stellar structures within a heliocentric range of ~5 to 30 kpc. Some of these are known (the Monoceros Ring, the Pisces/Triangulum globular cluster stream), but we also uncover three well-defined stellar structures that could be, at least partly, responsible for the so-called Triangulum/Andromeda and Triangulum/Andromeda 2 features. In particular, we trace a new faint stellar stream located at a heliocentric distance of ~17 kpc. With a surface brightness of \Sigma_V ~ 32-32.5 mag/arcsec^2, it follows an orbit that is almost parallel to the Galactic plane north of M31 and has so far eluded surveys of the Milky Way halo as these tend to steer away from regions dominated by the Galactic disk. Investigating our follow-up spectroscopic observations of PAndAS, we serendipitously uncover a radial velocity signature from stars that have colors and magnitudes compatible with the stream. From the velocity of eight likely member stars, we show that this stellar structure is dynamically cold, with an unresolved velocity dispersion that is lower than 7.1 km/s at the 90-percent confidence level. Along with the width of the stream (300-650 pc), its dynamics points to a dwarf-galaxy-accretion origin. The numerous stellar structures we can map in the Milky Way stellar halo between 5 and 30 kpc and their varying morphology is a testament to the complex nature of the stellar halo at these intermediate distances.

We present an analysis of the large-scale structure of the halo of the Andromeda galaxy, based on the Pan-Andromeda Archeological Survey (PAndAS), currently the most complete map of resolved stellar populations in any galactic halo. Despite copious substructure, the global halo populations follow closely power law profiles that become steeper with increasing metallicity. We divide the sample into stream-like populations and a smooth halo component. Fitting a three-dimensional halo model reveals that the most metal-poor populations ([Fe/H]<-1.7) are distributed approximately spherically (slightly prolate with ellipticity c/a=1.09+/-0.03), with only a relatively small fraction (42%) residing in discernible stream-like structures. The sphericity of the ancient smooth component strongly hints that the dark matter halo is also approximately spherical. More metal-rich populations contain higher fractions of stars in streams (86% for [Fe/H]>-0.6). The space density of the smooth metal-poor component has a global power-law slope of -3.08+/-0.07, and a non-parametric fit shows that the slope remains nearly constant from 30kpc to 300kpc. The total stellar mass in the halo at distances beyond 2 degrees is 1.1x10^10 Solar masses, while that of the smooth component is 3x10^9 Solar masses. Extrapolating into the inner galaxy, the total stellar mass of the smooth halo is plausibly 8x10^9 Solar masses. We detect a substantial metallicity gradient, which declines from [Fe/H]=-0.7 at R=30kpc to [Fe/H]=-1.5 at R=150kpc for the full sample, with the smooth halo being 0.2dex more metal poor than the full sample at each radius. While qualitatively in-line with expectations from cosmological simulations, these observations are of great importance as they provide a prototype template that such simulations must now be able to reproduce in quantitative detail.

We undertake the first detailed analysis of the stellar population and spatial properties of a diffuse substructure in the outer halo of M31. The South-West Cloud lies at a projected distance of ~100 kpc from the centre of M31, and extends for at least ~50 kpc in projection. We use Pan-Andromeda Archaeological Survey photometry of red giant branch stars to determine a distance to the South-West Cloud of 793 +/- 45 kpc. The metallicity of the cloud is found to be [Fe/H] = -1.3 +/- 0.1. This is consistent with the coincident globular clusters PAndAS-7 and PAndAS-8, which have metallicities determined using an independent technique of [Fe/H] = -1.35 +/- 0.15. We measure a brightness for the Cloud of M_V = -12.1 mag; this is ~75 per cent of the luminosity implied by the luminosity-metallicity relation. Under the assumption that the South-West Cloud is the visible remnant of an accreted dwarf satellite, this suggests that the progenitor object was amongst M31's brightest dwarf galaxies prior to disruption.

We present a generic algorithm to search for dwarf galaxies in photometric catalogs and apply it to the Pan-Andromeda Archaeological Survey (PAndAS). The algorithm is developed in a Bayesian framework and, contrary to most dwarf-galaxy-search codes, makes use of both the spatial and color-magnitude information of sources in a probabilistic approach. Accounting for the significant contamination from the Milky Way foreground and from the structured stellar halo of the Andromeda galaxy, we recover all known dwarf galaxies in the PAndAS footprint with high significance, even for the least luminous ones. Some Andromeda globular clusters are also recovered and, in one case, discovered. We publish a list of the 143 most significant detections yielded by the algorithm. The combined properties of the 39 most significant isolated detections show hints that at least some of these trace genuine dwarf galaxies, too faint to be individually detected. Follow-up observations by the community are mandatory to establish which are real members of the Andromeda satellite system. The search technique presented here will be used in an upcoming contribution to determine the PAndAS completeness limits for dwarf galaxies. Although here tuned to the search of dwarf galaxies in the PAndAS data, the algorithm can easily be adapted to the search for any localised overdensity whose properties can be modeled reliably in the parameter space of any catalog.

M31 has a giant stream of stars extending far to the south and a great deal of other tidal debris in its halo, much of which is thought to be directly associated with the southern stream. We model this structure by means of Bayesian sampling of parameter space, where each sample uses an N-body simulation of a satellite disrupting in M31's potential. We combine constraints on stellar surface densities from the Isaac Newton Telescope survey of M31 with kinematic data and photometric distances. This combination of data tightly constrains the model, indicating a stellar mass at last pericentric passage of log(M_s / Msun) = 9.5+-0.1, comparable to the LMC. Any existing remnant of the satellite is expected to lie in the NE Shelf region beside M31's disk, at velocities more negative than M31's disk in this region. This rules out the prominent satellites M32 or NGC 205 as the progenitor, but an overdensity recently discovered in M31's NE disk sits at the edge of the progenitor locations found in the model. M31's virial mass is constrained in this model to be log(M200) = 12.3+-0.1, alleviating the previous tension between observational virial mass estimates and expectations from the general galactic population and the timing argument. The techniques used in this paper, which should be more generally applicable, are a powerful method of extracting physical inferences from observational data on tidal debris structures.

We present the first kinematic analysis of the far outer halo globular cluster (GC) population in the Local Group galaxy M31. Our sample contains 53 objects with projected radii of ~20-130 kpc, of which 44 have no previous spectroscopic information. GCs with projected radii >30 kpc are found to exhibit net rotation around the minor axis of M31, in the same sense as the inner GCs, albeit with a smaller amplitude of 79 +/-19 km/s. The rotation-corrected velocity dispersion of the full halo GC sample is 106 +/-12 km/s, which we observe to decrease with increasing projected radius. We find compelling evidence for kinematic-coherence amongst GCs which project on top of halo substructure, including a clear signature of infall for GCs lying along the North-West stream. Using the tracer mass estimator, we estimate the dynamical mass of M31 within 200 kpc to be M_M31 = (1.2-1.5) +/- 0.2 x 10^12 M_sun. This value is highly dependent on the chosen model and assumptions within.

Dwarf satellite galaxies are thought to be the remnants of the population of primordial structures that coalesced to form giant galaxies like the Milky Way. An early analysis noted that dwarf galaxies may not be isotropically distributed around our Galaxy, as several are correlated with streams of HI emission, and possibly form co-planar groups. These suspicions are supported by recent analyses, and it has been claimed that the apparently planar distribution of satellites is not predicted within standard cosmology, and cannot simply represent a memory of past coherent accretion. However, other studies dispute this conclusion. Here we report the existence (99.998% significance) of a planar sub-group of satellites in the Andromeda galaxy, comprising approximately 50% of the population. The structure is vast: at least 400 kpc in diameter, but also extremely thin, with a perpendicular scatter <14.1 kpc (99% confidence). Radial velocity measurements reveal that the satellites in this structure have the same sense of rotation about their host. This finding shows conclusively that substantial numbers of dwarf satellite galaxies share the same dynamical orbital properties and direction of angular momentum, a new insight for our understanding of the origin of these most dark matter dominated of galaxies. Intriguingly, the plane we identify is approximately aligned with the pole of the Milky Way's disk and is co-planar with the Milky Way to Andromeda position vector. The existence of such extensive coherent kinematic structures within the halos of massive galaxies is a fact that must be explained within the framework of galaxy formation and cosmology.

Large scale surveys of the prominent members of the Local Group have provided compelling evidence for the hierarchical formation of massive galaxies, revealing a wealth of substructure that is thought to be the debris from ancient and on-going accretion events. In this paper, we compare two extant surveys of the M31-M33 subgroup of galaxies; the Pan-Andromeda Archaeological Survey (PAndAS) of the stellar structure, and a combination of observations of the HI gaseous content, detected at 21cm. Our key finding is a marked lack of spatial correlation between these two components on all scales, with only a few potential overlaps between stars and gas.The paucity of spatial correlation significantly restricts the analysis of kinematic correlations, although there does appear to the HI kinematically associated with the Giant Stellar Stream where it passes the disk of M31. These results demonstrate that that different processes must significantly influence the dynamical evolution of the stellar and HI components of substructures, such as ram pressure driving gas away from a purely gravitational path. Detailed modelling of the offset between the stellar and gaseous substructure will provide a determination of the properties of the gaseous halo of M31 and M33.

We use data from the Pan-Andromeda Archaeological Survey (PAndAS) to search for evidence of an extended halo component belonging to M33 (the Triangulum Galaxy). We identify a population of red giant branch (RGB) stars at large radii from M33's disk whose connection to the recently discovered extended "disk substructure" is ambiguous, and which may represent a "bona-fide" halo component. After first correcting for contamination from the Milky Way foreground population and misidentified background galaxies, we average the radial density of RGB candidate stars over circular annuli centered on the galaxy and away from the disk substructure. We find evidence of a low-luminosity, centrally concentrated component that is everywhere in our data fainter than mu_V ~ 33 mag arcsec^(-2). The scale length of this feature is not well constrained by our data, but it appears to be of order r_exp ~ 20 kpc; there is weak evidence to suggest it is not azimuthally symmetric. Inspection of the overall CMD for this region that specifically clips out the disk substructure reveals that this residual RGB population is consistent with an old population with a photometric metallicity of around [Fe/H] ~ -2 dex, but some residual contamination from the disk substructure appears to remain. We discuss the likelihood that our findings represent a bona-fide halo in M33, rather than extended emission from the disk substructure. We interpret our findings in terms of an upper limit to M33's halo that is a few percent of its total luminosity, although its actual luminosity is likely much less.

In `A Bayesian Approach to Locating the Red Giant Branch Tip Magnitude (PART I),' a new technique was introduced for obtaining distances using the TRGB standard candle. Here we describe a useful complement to the technique with the potential to further reduce the uncertainty in our distance measurements by incorporating a matched-filter weighting scheme into the model likelihood calculations. In this scheme, stars are weighted according to their probability of being true object members. We then re-test our modified algorithm using random-realization artificial data to verify the validity of the generated posterior probability distributions (PPDs) and proceed to apply the algorithm to the satellite system of M31, culminating in a 3D view of the system. Further to the distributions thus obtained, we apply a satellite-specific prior on the satellite distances to weight the resulting distance posterior distributions, based on the halo density profile. Thus in a single publication, using a single method, a comprehensive coverage of the distances to the companion galaxies of M31 is presented, encompassing the dwarf spheroidals Andromedas I - III, V, IX-XXVII and XXX along with NGC147, NGC 185, M33 and M31 itself. Of these, the distances to Andromeda XXIV - XXVII and Andromeda XXX have never before been derived using the TRGB. Object distances are determined from high-resolution tip magnitude posterior distributions generated using the Markov Chain Monte Carlo (MCMC) technique and associated sampling of these distributions to take into account uncertainties in foreground extinction and the absolute magnitude of the TRGB as well as photometric errors. The distance PPDs obtained for each object both with, and without the aforementioned prior are made available to the reader in tabular form...

The Andromeda galaxy (M31) shows many tidal features in its halo, including the Giant Southern Stream (GSS) and a sharp ledge in surface density on its western side (the W Shelf). Using DEIMOS on the Keck telescope, we obtain radial velocities of M31's giant stars along its NW minor axis, in a radial range covering the W Shelf and extending beyond its edge. In the space of velocity versus radius, the sample shows the wedge pattern expected from a radial shell, which is detected clearly here for the first time. This confirms predictions from an earlier model of formation of the GSS, which proposed that the W Shelf is a shell from the third orbital wrap of the same tidal debris stream that produces the GSS, with the main body of the progenitor lying in the second wrap. We calculate the distortions in the shelf wedge pattern expected from its outward expansion and angular momentum, and show that these effects are echoed in the data. In addition, a hot, relatively smooth spheroid population is clearly present. We construct a bulge-disk-halo N-body model that agrees with surface brightness and kinematic constraints, and combine it with a simulation of the GSS. From the contrasting kinematic signatures of the hot spheroid and shelf components, we decompose the observed stellar metallicity distribution into contributions from each component using a non-parametric mixture model. The shelf component's metallicity distribution matches previous observations of the GSS superbly, further strengthening the evidence they are connected and bolstering the case for a massive progenitor of this stream.

We determine the velocity vector of M31 with respect to the Milky Way and use this to constrain the mass of the Local Group, based on HST proper-motion measurements presented in Paper I. We construct N-body models for M31 to correct the measurements for the contributions from stellar motions internal to M31. We also estimate the center-of-mass motion independently, using the kinematics of satellite galaxies of M31 and the Local Group. All estimates are mutually consistent, and imply a weighted average M31 heliocentric transverse velocity of (v_W,v_N) = (-125.2+/-30.8, -73.8+/-28.4) km/s. We correct for the reflex motion of the Sun using the most recent insights into the solar motion within the Milky Way. This implies a radial velocity of M31 with respect to the Milky Way of V_rad = -109.3+/-4.4 km/s, and a tangential velocity V_tan = 17.0 km/s (<34.3 km/s at 1-sigma confidence). Hence, the velocity vector of M31 is statistically consistent with a radial (head-on collision) orbit towards the Milky Way. We revise prior estimates for the Local Group timing mass, including corrections for cosmic bias and scatter. Bayesian combination with other mass estimates yields M_LG = M_MW(vir) + M_M31(vir) = (3.17 +/- 0.57) x 10^12 solar masses. The velocity and mass results imply at 95% confidence that M33 is bound to M31, consistent with expectation from observed tidal deformations. (Abridged)

Subhalo abundance matching (SHAM) is a technique for populating simulated dark matter distributions with galaxies, assuming a monotonic relation between a galaxy's stellar mass or luminosity and the mass of its parent dark matter halo or subhalo. We examine the accuracy of SHAM in two cosmological SPH simulations, one of which includes momentum-driven winds. The SPH simulations indeed show a nearly monotonic relation between stellar mass and halo mass provided that, for satellite galaxies, we use the mass of the subhalo at the epoch when it became a satellite. In each simulation, the median relation for central and satellite galaxies is nearly identical, though a somewhat larger fraction of satellites are outliers. SHAM-assigned masses (at z=0-2), luminosities (R-band at z=0), or star formation rates (at z=2) have a 68% scatter of 0.09-0.15 dex relative to the true simulation values. When we apply SHAM to the subhalo population of collisionless N-body simulation with the same initial conditions as the SPH runs, we find generally good agreement for the halo occupation distributions and halo radial profiles of galaxy samples defined by thresholds in stellar mass. However, because a small fraction of SPH galaxies suffer severe stellar mass loss after becoming satellites, SHAM slightly overpopulates high mass halos; this effect is more significant for the wind simulation, which produces galaxies that are less massive and more fragile. SHAM recovers the two-point correlation function of the SPH galaxies in the no-wind simulation to better than 10% at scales 0.1 < r < 10 Mpc/h. For the wind simulation, agreement is better than 15% at r > 2 Mpc/h, but overpopulation of massive halos increases the correlation function by a factor of ~2.5 on small scales.

The combination of large size, high stellar density, high metallicity, and Sersic surface brightness profile of the spheroidal component of the Andromeda galaxy (M31) within R_proj ~ 20 kpc suggest that it is unlike any subcomponent of the Milky Way. In this work we capitalize on our proximity to and external view of M31 to probe the kinematical properties of this "inner spheroid." We employ a Markov chain Monte Carlo (MCMC) analysis of resolved stellar kinematics from Keck/DEIMOS spectra of 5651 red giant branch stars to disentangle M31's inner spheroid from its stellar disk. We measure the mean velocity and dispersion of the spheroid in each of five spatial bins after accounting for a locally cold stellar disk as well as the Giant Southern Stream and associated tidal debris. For the first time, we detect significant spheroid rotation (v_rot ~ 50 km/s) beyond R_proj ~ 5 kpc. The velocity dispersion decreases from about 140 km/s at R_proj = 7 kpc to 120 km/s at R_proj = 14 kpc, consistent to 2 sigma with existing measurements and models. We calculate the probability that a given star is a member of the spheroid and find that the spheroid has a significant presence throughout the spatial extent of our sample. Lastly, we show that the flattening of the spheroid is due to velocity anisotropy in addition to rotation. Though this suggests that the inner spheroid of M31 more closely resembles an elliptical galaxy than a typical spiral galaxy bulge, it should be cautioned that our measurements are much farther out (2 - 14 r_eff) than for the comparison samples.

We examine the star-forming history (SFH) of the M31 disk during the past few hundred Myr. The luminosity functions (LFs) of main sequence stars at distances R_GC > 21 kpc (i.e. > 4 disk scale lengths) are matched by models that assume a constant star formation rate (SFR). However, at smaller R_GC the LFs suggest that during the past ~10 Myr the SFR was 2 - 3 times higher than during the preceding ~100 Myr. The rings of cool gas that harbor a significant fraction of the current star-forming activity are traced by stars with ages ~100 Myr, indicating that (1) these structures have ages of at least 100 Myr, and (2) stars in these structures do not follow the same relation between age and random velocity as their counterparts throughout the disks of other spiral galaxies, probably due to the inherently narrow orbital angular momentum distribution of the giant molecular clouds in these structures. The distribution of evolved red stars is not azimuthally symmetric, in the sense that their projected density along the north east segment of the major axis is roughly twice that on the opposite side of the galaxy. The north east arm of the major axis thus appears to be a fossil star-forming area that dates to intermediate epochs. Such a structure may be the consequence of interactions with a companion galaxy.

We study the behaviour of multiple radiative cooling algorithms implemented in seven Semi-Analytic Models (SAMs) of galaxy formation, including a new model we propose in this paper. We use versions of the models without feedback and apply them to dark matter haloes growing in a cosmological context, which have final masses that range from 10^{11}Msun to 10^{14}Msun. First, using simplified smoothly-growing halo models, we demonstrate that the different algorithms predict cooling rates and final cold gas masses that differ by a factor of ~5 for massive haloes (>10^{12}Msun). The algorithms are in better agreement for less massive haloes because they cool efficiently and, therefore, their cooling rates are largely limited by the halo accretion rate. However, for less massive haloes, all the SAMs predict less cooling than corresponding 1D hydrodynamic models. Second, we study the gas accretion history of the central galaxies of dark matter haloes using merger trees. The inclusion of mergers alters the cooling history of haloes by locking up gas in galaxies within small haloes at early times. For realistic halo models, the dispersion in the cold gas mass predicted by the algorithms is 0.5 dex for high mass haloes and 0.1 dex for low mass haloes, while the dispersion in the accretion rate is about two times larger. Comparing to cosmological SPH simulations, we find that most SAMs systematically under-predict the gas accretion rates for low-mass haloes but over-predict the gas accretion rates for massive haloes. Although the models all include both "rapid" and "slow" mode accretion, the transition between the two accretion modes varies between models and also differs from the simulations. Finally, we construct a new model that explicitly incorporates cold halo gas to illustrate that such a class of models can better match the results from cosmological hydrodynamic simulations.

We present a new approach for identifying the Tip of the Red Giant Branch (TRGB) which, as we show, works robustly even on sparsely populated targets. Moreover, the approach is highly adaptable to the available data for the stellar population under study, with prior information readily incorporable into the algorithm. The uncertainty in the derived distances is also made tangible and easily calculable from posterior probability distributions. We provide an outline of the development of the algorithm and present the results of tests designed to characterize its capabilities and limitations. We then apply the new algorithm to three M31 satellites: Andromeda I, Andromeda II and the fainter Andromeda XXIII, using data from the Pan-Andromeda Archaeological Survey (PAndAS), and derive their distances as $731^{(+ 5) + 18}_{(- 4) - 17}$ kpc, $634^{(+ 2) + 15}_{(- 2) - 14}$ kpc and $733^{(+ 13)+ 23}_{(- 11) - 22}$ kpc respectively, where the errors appearing in parentheses are the components intrinsic to the method, while the larger values give the errors after accounting for additional sources of error. These results agree well with the best distance determinations in the literature and provide the smallest uncertainties to date. This paper is an introduction to the workings and capabilities of our new approach in its basic form, while a follow-up paper shall make full use of the method's ability to incorporate priors and use the resulting algorithm to systematically obtain distances to all of M31's satellites identifiable in the PAndAS survey area.

The Pan Andromeda Archeological Survey (PAndAS) CFHT Megaprime survey of the M31-M33 system has found a star stream which extends about 120 kpc NW from the center of M31. The great length of the stream, and the likelihood that it does not significantly intersect the disk of M31, means that it is unusually well suited for a measurement of stream gaps and clumps along its length as a test for the predicted thousands of dark matter sub-halos. The main result of this paper is that the density of the stream varies between zero and about three times the mean along its length on scales of 2 to 20 kpc. The probability that the variations are random fluctuations in the star density is less than 10^-5. As a control sample we search for density variations at precisely the same location in stars with metallicity higher than the stream, [Fe/H]=[0, -0.5] and find no variations above the expected shot noise. The lumpiness of the stream is not compatible with a low mass star stream in a smooth galactic potential, nor is it readily compatible with the disturbance caused by the visible M31 satellite galaxies. The stream's density variations appear to be consistent with the effects of a large population of steep mass function dark matter sub-halos, such as found in LCDM simulations, acting on an approximately 10Gyr old star stream. The effects of a single set of halo substructure realizations are shown for illustration, reserving a statistical comparison for another study.

Recently Menard et al. detected a subtle but systematic change in the mean color of quasars as a function of their projected separation from foreground galaxies, extending to comoving separations of ~10Mpc/h, which they interpret as a signature of reddening by intergalactic dust. We present theoretical models of this remarkable observation, using SPH cosmological simulations of a (50Mpc/h)^3 volume. Our primary model uses a simulation with galactic winds and assumes that dust traces the intergalactic metals. The predicted galaxy-dust correlation function is similar in form to the galaxy-mass correlation function, and reproducing the MSFR data requires a dust-to-metal mass ratio of 0.24, about half the value in the Galactic ISM. Roughly half of the reddening arises in dust that is more than 100Kpc/h from the nearest massive galaxy. We also examine a simulation with no galactic winds, which predicts a much smaller fraction of intergalactic metals (3% vs. 35%) and therefore requires an unphysical dust-to-metal ratio of 2.18 to reproduce the MSFR data. In both models, the signal is dominated by sightlines with E(g-i)=0.001-0.1. The no-wind simulation can be reconciled with the data if we also allow reddening to arise in galaxies up to several x 10^10 Msun. The wind model predicts a mean visual extinction of A_V ~0.0133 mag out to z=0.5, with a sightline-to-sightline dispersion similar to the mean, which could be significant for future supernova cosmology studies. Reproducing the MSFR results in these simulations requires that a large fraction of ISM dust survive its expulsion from galaxies and its residence in the intergalactic medium. Future observational studies that provide higher precision and measure the dependence on galaxy type and environment will allow detailed tests for models of enriched galactic outflows and the survival of IG dust.

The kinematics of damped Lyman alpha absorbers (DLAs) are difficult to reproduce in hierarchical galaxy formation models, particularly the preponderance of wide systems. We investigate DLA kinematics at z=3 using high-resolution cosmological hydrodynamical simulations that include a heuristic model for galactic outflows. Without outflows, our simulations fail to yield enough wide DLAs, as in previous studies. With outflows, predicted DLA kinematics are in much better agreement with observations. Comparing two outflow models, we find that a model based on momentum-driven wind scalings provides the best match to the observed DLA kinematic statistics of Prochaska & Wolfe. In this model, DLAs typically arise a few kpc away from galaxies that would be identified in emission. Narrow DLAs can arise from any halo and galaxy mass, but wide ones only arise in halos with mass >10^11 Mo, from either large central or small satellite galaxies. This implies that the success of this outflow model originates from being most efficient at pushing gas out from small satellite galaxies living in larger halos. This increases the cross-section for large halos relative to smaller ones, thereby yielding wider kinematics. Our simulations do not include radiative transfer effects or detailed metal tracking, and outflows are modeled heuristically, but they strongly suggest that galactic outflows are central to understanding DLA kinematics. An interesting consequence is that DLA kinematics may place constraints on the nature and efficiency of gas ejection from high-z galaxies.

We study the nature of rapidly star-forming galaxies at z=2 in cosmological hydrodynamic simulations, and compare their properties to observations of sub-millimetre galaxies (SMGs). We identify simulated SMGs as the most rapidly star-forming systems that match the observed number density of SMGs. In our models, SMGs are massive galaxies sitting at the centres of large potential wells, being fed by smooth infall and gas-rich satellites at rates comparable to their star formation rates (SFR). They are not typically undergoing major mergers that significantly boost their quiescent SFR, but they still often show complex gas morphologies and kinematics. Our simulated SMGs have stellar masses of log M*/Mo~11-11.7, SFRs of ~180-500 Mo/yr, a clustering length of 10 Mpc/h, and solar metallicities. The SFRs are lower than those inferred from far-IR data by a factor of 3, which we suggest may owe to one or more systematic effects in the SFR calibrations. SMGs at z=2 live in ~10^13 Mo halos, and by z=0 they mostly end up as brightest group galaxies in ~10^14 Mo halos. We predict that higher-M* SMGs should have on average lower specific SFRs, less disturbed morphologies, and higher clustering. We also predict that deeper far-IR surveys will smoothly join SMGs onto the massive end of the SFR-M* relationship defined by lower-mass z=2 galaxies. Overall, our simulated rapid star-formers provide as good a match to available SMG data as merger-based scenarios, offering an alternative scenario that emerges naturally from cosmological simulations.

We analyse cosmological hydrodynamic simulations that include observationally-constrained prescriptions for galactic outflows. If these simulated winds accurately represent winds in the real Universe, then material previously ejected in winds provides the dominant source of gas infall for new star formation at redshifts z<1. This recycled wind accretion, or wind mode, provides a third physically distinct accretion channel in addition to the "hot" and "cold" modes emphasised in recent theoretical studies. Because of the interaction between outflows and gas in and around halos, the recycling timescale of wind material (t_rec) is shorter in higher-mass systems, which reside in denser gaseous environments. In these simulations, this differential recycling plays a central role in shaping the present-day galaxy stellar mass function (GSMF). If we remove all particles that were ever ejected in a wind, then the predicted GSMFs are much steeper than observed; galaxy masses are suppressed both by the direct removal of gas and by the hydrodynamic heating of their surroundings, which reduces subsequent infall. With wind recycling included, the simulation that incorporates our favoured momentum-driven wind scalings reproduces the observed GSMF for stellar masses 10^9 < M < 5x10^10 Msolar. At higher masses, wind recycling leads to excessive galaxy masses and excessive star formation rates relative to observations. In these massive systems, some quenching mechanism must suppress the re-accretion of gas ejected from star-forming galaxies. In short, as has long been anticipated, the form of the GSMF is governed by outflows; the unexpected twist here for our simulated winds is that it is not primarily the ejection of material but how the ejected material is re-accreted that governs the GSMF.

We examine the impact of gas pressure on the transverse coherence of high-redshift (2 <= z <= 4) Lyman-alpha forest absorption along neighboring lines of sight that probe the gas Jeans scale (projected separation Delta r <= 500 kpc/h comoving; angular separation Delta theta <= 30"). We compare predictions from two smoothed particle hydrodynamics (SPH) simulations that have different photoionization heating rates and thus different temperature-density relations in the intergalactic medium (IGM). We also compare spectra computed from the gas distributions to those computed from the pressureless dark matter. The coherence along neighboring sightlines is markedly higher for the hotter, higher pressure simulation, and lower for the dark matter spectra. We quantify this coherence using the flux cross-correlation function and the conditional distribution of flux decrements as a function of transverse and line-of-sight (velocity) separation. Sightlines separated by Delta theta <= 15" are ideal for probing this transverse coherence. Higher pressure decreases the redshift-space anisotropy of the flux correlation function, while higher thermal broadening increases the anisotropy. In contrast to the longitudinal (line-of-sight) structure of the Lya forest, the transverse structure on these scales is dominated by pressure effects rather than thermal broadening. With the rapid recent growth in the number of known close quasar pairs, paired line-of-sight observations offer a promising new route to probe the IGM temperature-density relation and test the unexpectedly high temperatures that have been inferred from single sightline analyses.

In the low density intergalactic medium (IGM) that gives rise to the Lyman-alpha forest, gas temperature and density are tightly correlated. The velocity scale of thermal broadening and the Hubble flow across the gas Jeans scale are of similar magnitude (Hlambda_J ~ sigma_th). To separate the effects of gas pressure support and thermal broadening on the Lya forest, we compare spectra extracted from two smoothed particle hydrodynamics (SPH) simulations evolved with different photoionization heating rates (and thus different Jeans scales), imposing different temperature-density relations on the evolved particle distributions. The turnover scales in the flux power spectrum and flux autocorrelation function are determined mainly by thermal broadening rather than pressure. However, the insensitivity to pressure arises partly from a cancellation effect with a sloped temperature-density relation (T ~ rho^{0.6} in our simulations): the high density peaks in the colder, lower pressure simulation are less smoothed by pressure support than in the hotter simulation, and it is this higher density gas that experiences the strongest thermal broadening. Changes in thermal broadening and pressure support have comparably important effects on the flux probability distribution (PDF), which responds directly to the gas overdensity distribution rather than the scale on which it is smooth. Tests on a lower resolution simulation show that our statistical results are converged even at this lower resolution. While thermal broadening generally dominates the longitudinal structure in the Lya forest, we show in Paper II that pressure support determines the transverse coherence of the forest observed towards close quasar pairs. [ABRIDGED]

Photometric surveys of M31's halo vividly illustrate the wreckage caused by hierarchical galaxy formation. Several of M31's satellites are being disrupted by M31's tidal field, among them M33 and And I, while other tidal structures are the corpses of satellites already destroyed. The extent to which M31's satellites have left battle scars upon it is unknown; to answer this we need accurate orbits and masses of the perturbers. I focus here on M31's 150-kpc-long Giant Southern Stream (GSS) as an example of how these can be determined even in the absence of a visible progenitor. Comparing N-body models to photometric and spectroscopic data, I find this stream resulted from the disruption of a large satellite galaxy by a close passage about 750 Myr ago. The GSS is connected to several other debris structures in M31's halo. Bayesian sampling of the simulations estimates the progenitor's initial mass as log(Mstar/Msun) = 9.5 +- 0.2, showing it was one of the most massive Local Group galaxies until quite recently. The stream model constrains M31's halo mass to be (1.8+-0.5)x10^{12} Msun. While these small uncertainties neglect several important degrees of freedom, they are likely to remain good even with a more complete model. Future work on M31's satellites and streams will provide independent constraints on M31's mass, and reveal the shared history of M31 and its halo components.

We present the discovery of two new dwarf galaxies, Andromeda XXI and Andromeda XXII, located in the surroundings of the Andromeda and Triangulum galaxies (M31 and M33). These discoveries stem from the first year data of the Pan-Andromeda Archaeological Survey (PAndAS), a photometric survey of the M31/M33 group conducted with the Megaprime/MegaCam wide-field camera mounted on the Canada-France-Hawaii Telescope. Both satellites appear as spatial overdensities of stars which, when plotted in a color-magnitude diagram, follow metal-poor, [Fe/H]=-1.8, red giant branches at the distance of M31/M33. Andromeda XXI is a moderately bright dwarf galaxy (M_V=-9.9+/-0.6), albeit with low surface brightness, emphasizing again that many relatively luminous M31 satellites still remain to be discovered. It is also a large satellite, with a half-light radius close to 1 kpc, making it the fourth largest Local Group dwarf spheroidal galaxy after the recently discovered Andromeda XIX, Andromeda II and Sagittarius around the Milky Way, and supports the trend that M31 satellites are larger than their Milky Way counterparts. Andromeda XXII is much fainter (M_V=-6.5+/-0.8) and lies a lot closer in projection to M33 than it does to M31 (42 vs. 224 kpc), suggesting that it could be the first Triangulum satellite to be discovered. Although this is a very exciting possibility in the context of a past interaction of M33 with M31 and the fate of its satellite system, a confirmation will have to await a good distance estimate to confirm its physical proximity to M33. Along with the dwarf galaxies found in previous surveys of the M31 surroundings, these two new satellites bring the number of dwarf spheroidal galaxies in this region to 20.

In hierarchical cosmological models, galaxies grow in mass through the continual accretion of smaller ones. The tidal disruption of these systems is expected to result in loosely bound stars surrounding the galaxy, at distances that reach $10 - 100$ times the radius of the central disk. The number, luminosity and morphology of the relics of this process provide significant clues to galaxy formation history, but obtaining a comprehensive survey of these components is difficult because of their intrinsic faintness and vast extent. Here we report a panoramic survey of the Andromeda galaxy (M31). We detect stars and coherent structures that are almost certainly remnants of dwarf galaxies destroyed by the tidal field of M31. An improved census of their surviving counterparts implies that three-quarters of M31's satellites brighter than $M_V < -6$ await discovery. The brightest companion, Triangulum (M33), is surrounded by a stellar structure that provides persuasive evidence for a recent encounter with M31. This panorama of galaxy structure directly confirms the basic tenets of the hierarchical galaxy formation model and reveals the shared history of M31 and M33 in the unceasing build-up of galaxies.

We compare the properties of galaxies that form in a cosmological simulation without strong feedback to observations at z=0. We confirm previous findings that models without strong feedback overproduce the observed galaxy baryonic mass function, especially at the low and high mass extremes. Through post-processing we investigate what kinds of feedback would be required to reproduce observed galaxy masses and star formation rates. To mimic an extreme form of "preventive" feedback (e.g., AGN radio mode) we remove all baryonic mass that was originally accreted via "hot mode" from shock-heated gas. This does not bring the high mass end of the galaxy mass function into agreement with observations because much of the stellar mass in these systems formed at high redshift from baryons that originally accreted via "cold mode" onto lower mass progenitors. An efficient "ejective" feedback mechanism, such as supernova driven winds, must reduce the masses of these progenitors. Feedback must also reduce the masses of lower mass z=0 galaxies, which assemble at lower redshifts and have much lower star formation rates. If we monotonically re-map galaxy masses to reproduce the observed mass function, but retain the simulation's predicted star formation rates, we obtain fairly good agreement with the observed sequence of star-forming galaxies but fail to recover the observed population of passive, low star formation rate galaxies. Suppressing all hot mode accretion improves agreement for high mass galaxies but worsens the agreement at intermediate masses. Reproducing these z=0 observations requires a feedback mechanism that dramatically suppresses star formation in a fraction of galaxies, increasing with mass, while leaving star formation rates of other galaxies essentially unchanged.

We study the formation of galaxies in a (50 Mpc/h)^3 cosmological simulation (2x288^3 particles), evolved using the entropy conserving SPH code Gadget-2. Most of the baryonic mass in galaxies of all masses is originally acquired through filamentary "cold mode" accretion of gas that was never shock heated to its halo virial temperature, confirming the key feature of our earlier results obtained with a different SPH code (Keres et al. 2005). Atmospheres of hot, virialized gas develop in halos above ~2.5e11 Msun, a transition mass that is nearly constant from z=3 to z=0. Cold accretion persists in halos above the transition mass, especially at z>=2. It dominates the growth of galaxies in low mass halos at all times, and it is the main driver of the cosmic star formation history. Satellite galaxies have accretion rates similar to central galaxies of the same baryonic mass at high redshifts, but they have less accretion than comparable central galaxies at low redshift. Relative to our earlier results, the Gadget-2 simulations predict much lower rates of "hot mode" accretion from the virialized gas component of massive halos. At z<=1, typical hot accretion rates in halos above 5e12 Msun are below 1 Msun/yr, even though our simulation does not include AGN heating or other forms of "preventive" feedback. The inner density profiles of hot gas in these halos are shallow, with long associated cooling times. The cooling recipes typically used in semi-analytic models can overestimate the accretion rates in these halos by orders of magnitude, so such models may overemphasize the role of preventive feedback in producing observed galaxy masses and colors. A fraction of the massive halos develop cuspy profiles and significant cooling rates between z=1 and z=0, a redshift trend similar to the observed trend in the frequency of cooling flow clusters.

We use N-body/gasdynamical LambdaCDM cosmological simulations to examine the effect of the assembly of a central galaxy on the shape and mass profile of its dark halo. Two series of simulations are compared; one that follows only the evolution of the dark matter component and a second one where a baryonic component is added. These simulations include radiative cooling but neglect star formation and feedback, leading most baryons to collect at the halo center in a disk which is too small and too massive when compared with typical spiral. This unrealistic model allows us, nevertheless, to gauge the maximum effect that galaxies may have in transforming their dark halos. We find that the shape of the halo becomes more axisymmetric: halos are transformed from triaxial into essentially oblate systems, with well-aligned isopotential contours of roughly constant flattening (c/a ~ 0.85). Halos always contract as a result of galaxy assembly, but the effect is substantially less pronounced than predicted by the "adiabatic contraction" hypothesis. The reduced contraction helps to reconcile LambdaCDM halos with constraints on the dark matter content inside the solar circle and should alleviate the long-standing difficulty of matching simultaneously the scaling properties of galaxy disks and the luminosity function. The halo contraction is also less pronounced than found in earlier simulations, a disagreement that suggests that halo contraction is not solely a function of the initial and final distribution of baryons. Not only how much baryonic mass has been deposited at the center of a halo matters, but also the mode of its deposition. It might prove impossible to predict the halo response without a detailed understanding of a galaxy's assembly history. (Abriged)

(Abridged) We use the simulations presented in Poole et al. 2006 to examine the effects of mergers on compact cool cores in X-ray clusters. We propose a scheme for classifying the morphology of clusters based on their surface brightness and entropy profiles. Three dominant morphologies emerge: two hosting compact cores and central temperatures which are cool (CCC systems) or warm (CWC systems) and one hosting extended cores which are warm (EWC systems). We find that CCC states are disrupted only after direct collisions between cluster cores in head-on collisions or during second pericentric passage in off-axis mergers. By the time they relax, our remnant cores have generally been heated to warm core (CWC or EWC) states but subsequently recover CCC states. The only case resulting in a long-lived EWC state is a slightly off-axis 3:1 merger for which the majority of shock heating occurs during the accretion of a low-entropy stream formed from the disruption of the secondary's core. Compression prevents core temperatures from falling until after relaxation thus explaining the observed population of relaxed CWC systems with no need to invoke AGN feedback. The morphological segregation observed in the L_x-T_x and beta-r_c scaling relations is reflected in our simulations as well. However, none of the cases we have studied produce sufficiently high remnant central entropies to account for the most under-luminous EWC systems observed. Lastly, systems which initially host central metallicity gradients do not yield merger remnants with flat metallicity profiles. Taken together, these results suggest that once formed, compact core systems are remarkably stable against disruption from mergers. It remains to be demonstrated exactly how the sizable observed population of extended core systems was formed.

The halo region of M31 exhibits a startling level of stellar inhomogeneities, the most prominent of which is the "giant southern stream". Our previous analysis indicates that this stream, as well as several other observed features, are products of the tidal disruption of a single satellite galaxy with stellar mass ~10^9 solar masses less than 1 Gyr ago. Here we show that the specific observed morphology of the stream and halo debris favors a cold, rotating, disk-like progenitor over a dynamically hot, non-rotating one. These observed characteristics include the asymmetric distribution of stars along the stream cross-section and its metal-rich core/metal-poor sheath structure. We find that a disk-like progenitor can also give rise to arc-like features on the minor axis at certain orbital phases that resemble the recently discovered minor-axis "streams", even reproducing the lower observed metallicity of these streams. Though interpreted by the discoverers as new, independent tidal streams, our analysis suggests that these minor-axis streams may alternatively arise from the progenitor of the giant southern stream. Overall, our study points the way to a more complete reconstruction of the stream progenitor and its merger with M31, based on the emerging picture that most of the major inhomogeneities observed in the M31 halo share a common origin with the giant stream.

We present the discovery of a kinematically-cold stellar population along the SE minor axis of the Andromeda galaxy (M31) that is likely the forward continuation of M31's giant southern stream. This discovery was made in the course of an on-going spectroscopic survey of red giant branch (RGB) stars in M31 using the DEIMOS instrument on the Keck II 10-m telescope. Stellar kinematics are investigated in eight fields located 9-30 kpc from M31's center (in projection). A likelihood method based on photometric and spectroscopic diagnostics is used to isolate confirmed M31 RGB stars from foreground Milky Way dwarf stars: for the first time, this is done without using radial velocity as a selection criterion, allowing an unbiased study of M31's stellar kinematics. The radial velocity distribution of the 1013 M31 RGB stars shows evidence for the presence of two components. The broad (hot) component has a velocity dispersion of 129 km/s and presumably represents M31's virialized spheroid. A significant fraction (19%) of the population is in a narrow (cold) component centered near M31's systemic velocity with a velocity dispersion that decreases with increasing radial distance, from 55.5 km/s at R_proj=12 kpc to 10.6 km/s at R_proj=18 kpc. The spatial and velocity distribution of the cold component matches that of the "Southeast shelf" predicted by the Fardal et al. (2007) orbital model of the progenitor of the giant southern stream. The metallicity distribution of the cold component matches that of the giant southern stream, but is about 0.2 dex more metal rich on average than that of the hot spheroidal component. We discuss the implications of our discovery on the interpretation of the intermediate-age spheroid population found in this region in recent ultra-deep HST imaging studies.

We use the suite of simulations presented in Poole et al (2006) to examine global X-ray and Sunyaev-Zel'dovich (SZ) observables for systems of merging relaxed X-ray clusters. In all of the scaling relations generated from the properties we have studied, we observe a generic evolution (in good qualitative agreement with previous authors): a rapid transient roughly along the mass scaling relations, a subsequent slow drift across the scatter until virialization, followed by a slow evolution along and up the mass scaling relations as cooling recovers in the cluster cores. However, this drift is not sufficient to account for the observed scatter in the scaling relations. We also study the effects of mergers on several theoretical temperature measures of the intracluster medium: emission weighted measures (T_ew), the spectroscopic-like measure proposed by Mazzotta et al (2004; T_sl) and plasma model fits to the integrated spectrum of the system (T_spec). We find that T_sl tracks T_spec for the entire duration of our mergers, illustrating that it remains a good tool for observational comparison even for highly disturbed systems. Furthermore, the transient temperature increases produced during first and second pericentric passage are 15-40% larger for T_ew than for T_sl or T_spec. This suggests that the effects of transient temperature increases on sigma_8 and Omega_M derived by Randall et al (2002) are over estimated. Lastly, we examine the X-ray SZ proxy proposed by Kravtsov et al (2006) and find that the tight mass scaling relation they predict remains secure through the entire duration of a merger event, independent of projection effects.

Published maps of red giant stars in the halo region of M31 exhibit a giant stellar stream to the south of this galaxy, as well as a giant "shelf" to the northeast of M31's center. Using these maps, we find that there is a fainter shelf of comparable size on the western side as well. By choosing appropriate structural and orbital parameters for an accreting dwarf satellite within the accurate M31 potential model of Geehan et al. (2006), we produce a very similar structure in an N-body simulation. In this scenario, the tidal stream produced at pericenter of the satellite's orbit matches the observed southern stream, while the forward continuation of this tidal stream makes up two orbital loops, broadened into fan-like structures by successive pericentric passages; these loops correspond to the NE and W shelves. The tidal debris from the satellite also reproduces a previously-observed "stream" of counterrotating PNe and a related stream seen in red giant stars. The debris pattern in our simulation resembles the shell systems detected around many elliptical galaxies, though this is the first identification of a shell system in a spiral galaxy and the first in any galaxy close enough to allow measurements of stellar velocities and relative distances. We discuss the physics of these partial shells, highlighting the role played by spatial and velocity caustics in the observations. We show that kinematic surveys of the tidal debris will provide a sensitive measurement of M31's halo potential, while quantifying the surface density of debris in the shelves will let us reconstruct the original mass and time of disruption of the progenitor satellite.

We report on the analysis of a suite of SPH simulations (incorporating cooling and star formation) of mergers involving idealised X-ray clusters whose initial conditions resemble relaxed clusters with cool compact cores observed by Chandra and XMM. The simulations sample the most interesting, theoretically plausible, range of impact parameters and progenitor mass ratios. We find that all mergers evolve via a common progression. We illustrate this progression in the projected gas density, X-ray surface brightness, SZ, temperature, and gas entropy maps. Several different classes of transient ``cold front''-like features can arise over the course of a merger. We find that all of these classes are present in Chandra and XMM observations of merging systems and propose a naming scheme for these features: ``comet-like'' tails, bridges, plumes, streams and edges. In none of the cases considered do the initial cool compact cores of the primary and the secondary get destroyed during the course of the mergers. We quantify the evolving morphology of our mergers using centroid variance, power ratios and offset between the X-ray and the projected mass maps. We find that the centroid variance best captures the dynamical state of the cluster. Placing the system at z=0.1, we find that all easily identified observable traces of the secondary disappear from a simulated 50 ks Chandra image following the second pericentric passage. The system, however, takes approximately 2 additional Gyrs to relax and virialize. Temperature fluctuations at the level of 20% can persist in the final systems well past the point of virialization, suggesting that that the existence of temperature fluctuations, in and of themselves, do not necessarily indicate a disturbed or unrelaxed system.

The total extragalactic background radiation can be an important test of the global star formation history (SFH). Using direct observational estimates of the SFH, along with standard assumptions about the initial mass function (IMF), we calculate the total extragalactic background radiation and the observed stellar density today. We show that plausible SFHs allow a significant range in each quantity, but that their ratio is very tightly constrained. Current estimates of the stellar mass and extragalactic background are difficult to reconcile, as long as the IMF is fixed to the Salpeter slope above 1 Msun. The joint confidence interval of these two quantities only agrees with that determined from the allowed range of SFH fits at the 3-sigma level, and for our best-fit values the discrepancy is about a factor of two. Alternative energy sources that contribute to the background, such as active galactic nuclei (AGN), Population III stars, or decaying particles, appear unlikely to resolve the discrepancy. However, changes to the IMF allow plausible solutions to the background problem. The simplest is an average IMF with an increased contribution from stars around 1.5--4 Msun. A ``paunchy'' IMF of this sort could emerge as a global average if low mass star formation is suppressed in galaxies experiencing rapid starbursts. Such an IMF is consistent with observations of star-forming regions, and would help to reconcile the fossil record of star formation with the directly observed SFH.

This paper is the first in a series which studies interactions between M31 and its satellites, including the origin of the giant southern stream. We construct accurate yet simple analytic models for the potential of the M31 galaxy to provide an easy basis for calculation of orbits in M31's halo. We use an NFW dark halo, an exponential disk, a Hernquist bulge, and a central black hole point mass to describe the galaxy potential. We constrain the parameters of these functions by comparing to existing surface brightness, velocity dispersion, and rotation curve measurements of M31. Our description provides a good fit to the observations, and agrees well with more sophisticated modeling of M31. While in many respects the parameter set is well constrained, there is substantial uncertainty in the outer halo potential and a near-degeneracy between the disk and halo components, producing a large, nearly two-dimensional allowed region in parameter space. We limit the allowed region using theoretical expectations for the halo concentration, baryonic content, and stellar $M/L$ ratio, finding a smaller region where the parameters are physically plausible. Our proposed mass model for M31 has $M_{bulge} = 3.2 \times 10^{10} \Msun$, $M_{disk} = 7.2 \times 10^{10} \Msun$, and $M_{200} = 7.1\times 10^{11} \Msun$, with uncorrected (for internal and foreground extinction) mass-to-light ratios of $M/L_R = 3.9$ and 3.3 for the bulge and disk, respectively. We present some illustrative test particle orbits for the progenitor of the stellar stream in our galaxy potential, highlighting the effects of the remaining uncertainty in the disk and halo masses.

We construct test-particle orbits and simple N-body models that match the properties of the giant stellar stream observed to the south of M31, using the model of M31's potential derived in the companion paper by Geehan et al. (2006). We introduce a simple approximation to account for the difference in position between the stream and the orbit of the progenitor; this significantly affects the best-fitting orbits. The progenitor orbits we derive have orbital apocenter $\sim 60 \kpc$ and pericenter $\sim 3 \kpc$, though these quantities vary somewhat with the current orbital phase of the progenitor which is as yet unknown. Our best combined fit to the stream and galaxy properties implies a mass within 125 kpc of M31 of $(7.4 \pm 1.2) \times 10^{11} \Msun$. Based on its length, width, luminosity, and velocity dispersion, we conclude that the stream originates from a progenitor satellite with mass $M_s \sim 10^9 \Msun$, and at most modest amounts of dark matter; the estimate of $M_s$ is again correlated with the phase of the progenitor. M31 displays a large number of faint features in its inner halo which may be progenitors or continuations of the stream. While the orbital fits are not constrained enough for us to conclusively identify the progenitor, we can identify several plausible candidates, of which a feature in the planetary nebula distribution found by Merrett et al. is the most plausible, and rule out several others. We make predictions for the kinematic properties of the successful candidates. These may aid in observational identification of the progenitor object, which would greatly constrain the allowed models of the stream.

High resolution X-ray observations indicate that the entropy profiles in the central regions of some massive cooling flow clusters are well approximated by powerlaws. McCarthy and coworkers recently accounted for this trend with an analytic model that includes the detailed effects of radiative cooling. Interestingly, these authors found that cooling (and subsequent inflow of the gas) naturally establishes approximate steady-state powerlaw entropy profiles in the cores of clusters. However, the origin of this behavior and its dependence on initial conditions have yet to be elucidated. In the present study, we explain this trend in the context of the self-similar cooling wave model developed previously by Bertschinger (1989). It is shown that the logarithmic slope of the entropy profile in the cores of relaxed cooling flow clusters is given by a simple analytic function that depends only on the logarithmic slopes of the local gravitational potential and the cooling function. We discuss a number of potentially interesting uses of the above result, including: (1) a way of measuring the shapes of gravitational potentials of cooling flow clusters (which may, for example, be compared with the results of the standard hydrostatic equilibrium method); (2) a simple method for constructing realistic analytic cluster models that include the effects of radiative cooling; and (3) a test of the reliability of cooling routines implemented in analytic models and hydrodynamic simulations.

We use a high resolution $\Lambda$CDM numerical simulation to calculate the mass function of dark matter haloes down to the scale of dwarf galaxies, back to a redshift of fifteen, in a 50 $h^{-1}$Mpc volume containing 80 million particles. Our low redshift results allow us to probe low $\sigma$ density fluctuations significantly beyond the range of previous cosmological simulations. The Sheth and Tormen mass function provides an excellent match to all of our data except for redshifts of ten and higher, where it overpredicts halo numbers increasingly with redshift, reaching roughly 50 percent for the $10^{10}-10^{11} \msun$ haloes sampled at redshift 15. Our results confirm previous findings that the simulated halo mass function can be described solely by the variance of the mass distribution, and thus has no explicit redshift dependence. We provide an empirical fit to our data that corrects for the overprediction of extremely rare objects by the Sheth and Tormen mass function. This overprediction has implications for studies that use the number densities of similarly rare objects as cosmological probes. For example, the number density of high redshift (z $\simeq$ 6) QSOs, which are thought to be hosted by haloes at 5$\sigma$ peaks in the fluctuation field, are likely to be overpredicted by at least a factor of 50%. We test the sensitivity of our results to force accuracy, starting redshift, and halo finding algorithm.

We study the predicted sub-mm emission from massive galaxies in a Lambda-CDM universe, using hydrodynamic cosmological simulations. Assuming that most of the emission from newly formed stars is absorbed and reradiated in the rest-frame far-IR, we calculate the number of galaxies that would be detected in sub-mm surveys conducted with SCUBA. The predicted number counts are strongly dependent on the assumed dust temperature and emissivity law. With plausible choices for SED parameters (e.g., T=35 K, beta=1.0), the simulation predictions reproduce the observed number counts above ~ 1 mJy. The sources have a broad redshift distribution with median z ~ 2, in reasonable agreement with observational constraints. However, the predicted count distribution may be too steep at the faint end, and the fraction of low redshift objects may be larger than observed. In this physical model of the sub-mm galaxy population, the objects detected in existing surveys consist mainly of massive galaxies (several M_*) forming stars fairly steadily over timescales ~ 10^8-10^9 years, at moderate rates ~100 Msun/yr. The typical descendants of these sub-mm sources are even more massive galaxies, with old stellar populations, found primarily in dense environments. While the resolution of our simulations is not sufficient to determine galaxy morphologies, these properties support the proposed identification of sub-mm sources with massive ellipticals in the process of formation. The most robust and distinctive prediction of this model, stemming directly from the long timescale and correspondingly moderate rate of star formation, is that the far-IR SEDs of SCUBA sources have a relative high 850 micron luminosity for a given bolometric luminosity. [Abridged]

(Abridged) We use a hydrodynamic simulation of a LambdaCDM model to predict the extragalactic X-ray background (XRB), focussing on emission from the intergalactic medium (IGM). We also include X-rays from point sources associated with galaxies in the simulation, and make maps of the angular distribution of the emission. We find that filaments in the maps are not evident, being diluted by projection. In the soft (0.5-2 keV) band, the mean intensity of radiation from intergalactic and cluster gas is 2.3*10^-12 ergdeg^-2cm^-2s^-1, 35% of the total soft band emission. This is compatible at the ~1 sigma level with estimates of the unresolved soft background from ROSAT and {\it Chandra}. Only 4% of the hard (2-10 keV) emission is associated with the IGM. Relative to AGN flux, the IGM component peaks at a lower redshift (median z~0.45) so its clustering makes an important contribution to that of the total XRB. The angular correlations on 0.1-10 arcmin scales are significant, with an amplitude roughly consistent with an extrapolation of recent ROSAT results to small scales. A cross-correlation of the XRB against nearby galaxies taken from a simulated redshift survey also yields a strong signal from the IGM. Although some recent papers have argued that the expected soft band intensity from gas in galaxy, group, and cluster halos would exceed XRB limits unless much of the gas is expelled by supernova feedback, we obtain reasonable compatibility with current observations in a simulation that incorporates cooling, star formation, and only modest feedback. A prediction of our model is that the unresolved portion of the soft XRB will remain mostly unresolved.

We examine the cooling radiation from forming galaxies in hydrodynamic simulations of the LCDM model (cold dark matter with a cosmological constant), focusing on the Ly-alpha line luminosities of high-redshift systems. Primordial composition gas condenses within dark matter potential wells, forming objects with masses and sizes comparable to the luminous regions of observed galaxies. As expected, the energy radiated in this process is comparable to the gravitational binding energy of the baryons, and the total cooling luminosity of the galaxy population peaks at z ~= 2. However, in contrast to the classical picture of gas cooling from the \sim 10^6 K virial temperature of a typical dark matter halo, we find that most of the cooling radiation is emitted by gas with T < 20,000 K. As a consequence, roughly 50% of this cooling radiation emerges in the Ly-alpha line. While a galaxy's cooling luminosity is usually smaller than the ionizing continuum luminosity of its young stars, the two are comparable in the most massive systems, and the cooling radiation is produced at larger radii, where the Ly-alpha photons are less likely to be extinguished by dust. We suggest, in particular, that cooling radiation could explain the two large (\sim 100 kpc), luminous (L_{Ly-alpha} \sim 10^{44} erg s^{-1}) ``blobs'' of Ly-alpha emission found in Steidel et al.'s (1999) narrow band survey of a z = 3 proto-cluster. Our simulations predict objects of the observed luminosity at about the right space density, and radiative transfer effects can account for the observed sizes and line widths. We discuss observable tests of this hypothesis for the nature of the Ly-alpha blobs, and we present predictions for the contribution of cooling radiation to the Ly-alpha luminosity function of galaxies as a function of redshift.

We compute the ionizing radiation field at low redshift, arising from Seyferts, QSOs, and starburst galaxies. This calculation combines recent Seyfert luminosity functions, extrapolated ultraviolet fluxes from our IUE-AGN database, and a new intergalactic opacity model based on Hubble Space Telescope and Keck Ly-alpha absorber surveys. At z = 0 for AGN only, our best estimate for the specific intensity at 1 Ryd is I_0 = 1.3 (+0.8/-0.5) x 10^-23 ergs/cm^2/s/Hz/sr, independent of H_0, Omega_0, and Lambda. The one-sided ionizing photon flux is Phi_ion = 3400 (+2100/-1300) photons/cm^2/s, and the H I photoionization rate is Gamma_HI = 3.2 (+2.0/-1.2) x 10^-14 s^-1 for alpha_s = 1.8. We also derive Gamma_ HI for z = 0 - 4. These error ranges reflect uncertainties in the spectral indexes for the ionizing EUV (alpha_s = 1.8 +/- 0.3) and the optical/UV (alpha_UV = 0.86 +/- 0.05), the IGM opacity model, the range of Seyfert luminosities (0.001 - 100 L*) and the completeness of the luminosity functions. Our estimate is a factor of three lower than the most stringent upper limits on the ionizing background (Phi_ion < 10^4 photons/cm^2/s) obtained from H-alpha observations in external clouds, and it lies within the range implied by other indirect measures. Starburst galaxies with a sufficiently large Lyman continuum escape fraction, f_ esc > 0.05, may provide a comparable background to AGN, I_0 (z=0) = 1.1 (+1.5/-0.7) x 10^{-23). An additional component of the ionizing background of this magnitude would violate neither upper limits from H-alpha observations nor the acceptable range from other measurements.

Absorption due to HeII Ly-alpha has been detected at low resolution in the spectra of four quasars between redshifts z = 2.74-3.29. We assess these observations, giving particular attention to the radiative transfer of the ionizing background radiation, cloud diffuse emission and ionization physics, and statistical fluctuations. We use high-resolution observations of HI absorption towards quasars to derive improved models for the opacity of the intergalactic medium (IGM) and use these models to calculate the HI and HeII photoionization history in both optically-thin and self-shielded clouds, and the average line-blanketing contribution to HeII absorption. The derived ionization rate at z = 2-4, Gamma(HI) = (1-3)x10^-12 s^-1, is consistent with the ionizing background intensity inferred from the proximity effect, but it remains larger than that inferred by N-body hydrodynamical simulations of the Lya absorber distribution. The HeII observations are con- sistent with line blanketing from clouds with N(HI) > 10^12 cm^-2, although a contribution from a more diffuse IGM would help to explain the observed opacity. We compute the expected HeII optical depth and examine the sizable fluc- tuations that arise from variations in the cloud numbers and ionizing radiation field. We assess how HeII absorption constrains the intensity and spectrum of the ionizing radiation and the fractional contributions of the dominant sources (quasars and starburst galaxies). Finally, we demonstrate how high-resolution UV observations can distinguish between absorption from the diffuse IGM and the Lya forest clouds and determine the source of the ionizing background.