We present Gemini-S and {\it Spitzer}-IRAC optical-through-near-IR observations in the field of the SPT2349-56 proto-cluster at $z=4.3$. We detect optical/IR counterparts for only nine of the 14 submillimetre galaxies (SMGs) previously identified by ALMA in the core of SPT2349-56. In addition, we detect four $z\sim4$ Lyman-break galaxies (LBGs) in the 30 arcsec diameter region surrounding this proto-cluster core. Three of the four LBGs are new systems, while one appears to be a counterpart of one of the nine observed SMGs. We identify a candidate brightest cluster galaxy (BCG) with a stellar mass of $(3.2^{+2.5}_{-1.4})\times10^{11}\,{\rm M}_{\odot}$. The stellar masses of the eight other SMGs place them on, above, and below the main sequence of star formation at $z\approx4.5$. The cumulative stellar mass for the SPT2349-56 core is at least $(11.5\pm2.9)\times10^{11}\,{\rm M}_{\odot}$, a sizeable fraction of the stellar mass in local BCGs, and close to the universal baryon fraction (0.16) relative to the virial mass of the core ($10^{13}\,{\rm M}_{\odot}$). As all 14 of these SMGs are destined to quickly merge, we conclude that the proto-cluster core has already developed a significant stellar mass at this early stage, comparable to $z=1$ BCGs. Importantly, we also find that the SPT2349-56 core structure would be difficult to uncover in optical surveys, with none of the ALMA sources being easily identifiable or constrained through $g,r,$ and $i$ colour-selection in deep optical surveys and only a modest overdensity of LBGs over the extended core structure. SPT2349-56 therefore represents a truly dust-obscured phase of a massive cluster core under formation.

We investigate the stellar and dynamical mass profiles in the centres of 25 brightest cluster galaxies (BCGs) at redshifts of 0.05 $\leq z \leq$ 0.30. Our spectroscopy enables us to robustly measure the Gauss-Hermite higher order velocity moments $h_{3}$ and $h_{4}$, which we compare to measurements for massive early-type galaxies, and central group galaxies. We measure positive central values for $h_{4}$ for all the BCGs. We derive the stellar mass-to-light ratio ($\Upsilon_{\star \rm DYN}$), and velocity anisotropy ($\beta$) based on a Multi-Gaussian Expansion (MGE) and axisymmetric Jeans Anisotropic Methods (JAM, cylindrically- and spherically-aligned). We explicitly include a dark matter halo mass component, which is constrained by weak gravitational lensing measurements for these clusters. We find a strong correlation between anisotropy and velocity dispersion profile slope, with rising velocity dispersion profiles corresponding to tangential anisotropy and decreasing velocity dispersion profiles corresponding to radial anisotropy. The rising velocity dispersion profiles can also indicate a significant contribution from the intracluster light (ICL) to the total light (in projection) in the centre of the galaxy. For a small number of BCGs with rising velocity dispersion profiles, a variable stellar mass-to-light ratio can also account for the profile shape, instead of tangential anisotropy or a significant ICL contribution. We note that, for some BCGs, a variable $\beta_{z}(r)$ (from radial to tangential anisotropy) can improve the model fit to the observed kinematic profiles. The observed diversity in these properties illustrates that BCGs are not the homogeneous class of objects they are often assumed to be.

Electromagnetic precursor waves generated by the synchrotron maser instability at relativistic magnetized shocks have been recently invoked to explain the coherent radio emission of Fast Radio Bursts. By means of two-dimensional particle-in-cell simulations, we explore the properties of the precursor waves in relativistic electron-positron perpendicular shocks as a function of the pre-shock magnetization $\sigma\gtrsim 1$ (i.e., the ratio of incoming Poynting flux to particle energy flux) and thermal spread $\Delta \gamma \equiv kT/mc^2=10^{-5}-10^{-1}$. We measure the fraction $f_\xi$ of total incoming energy that is converted into precursor waves, as computed in the post-shock frame. At fixed magnetization, we find that $f_{\xi}$ is nearly independent of temperature as long as $\Delta \gamma \lesssim 10^{-1.5}$ (with only a modest decrease of a factor of three from $\Delta\gamma =10^{-5}$ to $\Delta\gamma =10^{-1.5}$), but it drops by nearly two orders of magnitude for $\Delta \gamma \gtrsim 10^{-1}$. At fixed temperature, the scaling with magnetization $f_\xi\sim 10^{-3}\,\sigma^{-1}$ is consistent with our earlier one-dimensional results. For our reference $\sigma=1$, the power spectrum of precursor waves is relatively broad (fractional width $\sim 1-3$) for cold temperatures, whereas it shows pronounced line-like features with fractional width $\sim 0.2$ for $10^{-3} \lesssim \Delta \gamma \lesssim 10^{-1.5} $. For $\sigma\gtrsim 1$, the precursor waves are beamed within an angle $\simeq \sigma^{-1/2}$ from the shock normal (as measured in the post-shock frame), as required so they can outrun the shock. Our results can provide physically-grounded inputs for FRB emission models based on maser emission from relativistic shocks.

We analyze an extremely deep 450-$\mu$m image ($1\sigma=0.56$\,mJy\,beam$^{-1}$) of a $\simeq 300$\,arcmin$^{2}$ area in the CANDELS/COSMOS field as part of the SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). We select a robust (signal-to-noise ratio $\geqslant 4$) and flux-limited ($\geqslant 4$\,mJy) sample of 164 sub-millimeter galaxies (SMGs) at 450-$\mu$m that have $K$-band counterparts in the COSMOS2015 catalog identified from radio or mid-infrared imaging. Utilizing this SMG sample and the 4705 $K$-band-selected non-SMGs that reside within the noise level $\leqslant 1$\,mJy\,beam$^{-1}$ region of the 450-$\mu$m image as a training set, we develop a machine-learning classifier using $K$-band magnitude and color-color pairs based on the thirteen-band photometry available in this field. We apply the trained machine-learning classifier to the wider COSMOS field (1.6\,deg$^{2}$) using the same COSMOS2015 catalog and identify a sample of 6182 450-$\mu$m SMG candidates with similar colors. The number density, radio and/or mid-infrared detection rates, redshift and stellar mass distributions, and the stacked 450-$\mu$m fluxes of these SMG candidates, from the S2COSMOS observations of the wide field, agree with the measurements made in the much smaller CANDELS field, supporting the effectiveness of the classifier. Using this 450-$\mu$m SMG candidate sample, we measure the two-point autocorrelation functions from $z=3$ down to $z=0.5$. We find that the 450-$\mu$m SMG candidates reside in halos with masses of $\simeq (2.0\pm0.5) \times10^{13}\,h^{-1}\,\rm M_{\odot}$ across this redshift range. We do not find evidence of downsizing that has been suggested by other recent observational studies.

Recent observations show that the star formation rate (SFR) in the {\it Phoenix} cluster's central galaxy is $\sim 500$ M$_\odot$ yr$^{-1}$. Even though {\it Phoenix} is a massive cluster ($M_{200} \approx 2.0\times 10^{15}$ M$_\odot$; $z\approx 0.6$) such a high central SFR is not expected in a scenario in which feedback from an active galactic nucleus (AGN) maintains the intracluster medium (ICM) in a state of rough thermal balance. It has been argued that either AGN feedback saturates in very massive clusters or the central supermassive black hole (SMBH) is too small to produce enough kinetic feedback and hence is unable to quench the catastrophic cooling. In this work, we present an alternate scenario wherein intense short-lived cooling and star formation phases followed by strong AGN outbursts are part of the AGN feedback loop. Using results from a 3D hydrodynamic simulation of a standard cool-core cluster ($M_{200}\sim 7\times10^{14}$ M$_\odot$; $z=0$), scaled to account for differences in mass and redshift, we argue that {\it Phoenix} is at the end of a cooling phase in which an AGN outburst has begun but has not yet arrested core cooling. This state of high cooling rate and star formation is expected to last for $\lesssim$ 100 Myr in {\it Phoenix}.

We construct a SCUBA-2 450-$\mu$m map in the COSMOS field that covers an area of 300 arcmin$^{2}$ and reaches a 1$\sigma$ noise level of 0.65 mJy in the deepest region. We extract 256 sources detected at 450 $\mu$m with signal-to-noise ratio $>$ 4.0 and analyze the physical properties of their multi-wavelength counterparts. We find that most of the sources are at $z\lesssim3$, with a median of $z = 1.79^{+0.03}_{-0.15}$. About $35^{+32}_{-25}$% of our sources are classified as starburst galaxies based on their total star-formation rates (SFRs) and stellar masses ($M_{\ast}$). By fitting the far-infrared spectral energy distributions, we find that our 450-$\mu$m-selected sample has a wide range of dust temperatures (20 K $ \lesssim T_{\rm d} \lesssim$ 60 K), with a median of ${T}_{\rm d} = 38.3^{+0.4}_{-0.9}$ K. We do not find a redshift evolution in dust temperature for sources with $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ at $z<3$. However, we find a moderate correlation where dust temperature increases with the deviation from the SFR-$M_{\ast}$ relation. The increase in dust temperature also correlates with optical morphology, which is consistent with merger-triggered starbursts in sub-millimeter galaxies. Our galaxies do not show the tight IRX-$\beta_{\rm UV}$ correlation that has been observed in the local Universe. We construct the infrared luminosity functions of our 450-$\mu$m sources and measure their comoving SFR densities. The contribution of the $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ population to the SFR density rises dramatically from $z$ = 0 to 2 ($\propto$ ($1+z$)$^{3.9\pm1.1}$) and dominates the total SFR density at $z \gtrsim 2$.

We investigate the evolution of the gas mass fraction for galaxies in the COSMOS field using submillimetre emission from dust at 850$\mu$m. We use stacking methodologies on the 850$\mu$m S2COSMOS map to derive the gas mass fraction of galaxies out to high redshifts, 0 <= $z$ <= 5, for galaxies with stellar masses of $10^{9.5} < M_* (\rm M_{\odot}) < 10^{11.75}$. In comparison to previous literature studies we extend to higher redshifts, include more normal star-forming galaxies (on the main sequence), and also investigate the evolution of the gas mass fraction split by star-forming and passive galaxy populations. We find our stacking results broadly agree with scaling relations in the literature. We find tentative evidence for a peak in the gas mass fraction of galaxies at around $z$ ~ 2.5-3, just before the peak of the star formation history of the Universe. We find that passive galaxies are particularly devoid of gas, compared to the star-forming population. We find that even at high redshifts, high stellar mass galaxies still contain significant amounts of gas.

The current consensus on the formation and evolution of the brightest cluster galaxies is that their stellar mass forms early ($z \gtrsim 4$) in separate galaxies that then eventually assemble the main structure at late times ($z \lesssim 1$). However, advances in observational techniques have led to the discovery of protoclusters out to $z \sim 7$, suggesting that the late-assembly picture may not be fully complete. If these protoclusters assemble rapidly in the early universe, they should form the brightest cluster galaxies much earlier than suspected by the late-assembly picture. Using a combination of observationally constrained hydrodynamical and dark-matter-only simulations, we show that the stellar assembly time of a sub-set of the brightest cluster galaxies occurs at high redshifts ($z > 3$) rather than at low redshifts ($z < 1$), as is commonly thought. We find, using isolated non-cosmological hydrodynamical simulations, that highly overdense protoclusters assemble their stellar mass into brightest cluster galaxies within $\sim 1$ $\mathrm{Gyr}$ of evolution -- producing massive blue elliptical galaxies at high redshifts ($z \gtrsim 1.5$). We argue that there is a downsizing effect on the cluster scale wherein some of the brightest cluster galaxies in the cores of the most-massive clusters assemble earlier than those in lower-mass clusters. In those clusters with $z = 0$ virial mass $\geqslant 5\times10^{14}$ M$_\mathrm{\odot}$, we find that $9.8$% have their cores assembly early, and a higher fraction of $16.4$% in those clusters above $10^{15}$ M$_\mathrm{\odot}$. The James Webb Space Telescope will be able to detect and confirm our prediction in the near future, and we discuss the implications to constraining the value of $\sigma_\mathrm{8}$.

The intracluster medium (ICM) is a multi-phase environment, dynamically regulated by Active Galactic Nuclei (AGN), the motions of galaxies through it, and mergers with other clusters. AGN as a central heating source are key to preventing runaway cooling flows, but their role in heating cores in a cosmological context is still poorly understood. The activity of the AGN is strongly linked to star formation, especially in the Brightest Cluster Galaxy (BCG), likely because both rely on cold phase gas. A self-consistent model for AGN and star formation in galaxy clusters thus requires cosmological context, higher resolution, and a careful modeling of cooling and heating balance. In this paper, we use the high-resolution hydrodynamical cosmological simulation of the RomulusC galaxy cluster to study in detail the role of AGN and a major, head-on merger in shaping the cluster core. The unprecedented resolution of the RomulusC simulation captures the multiphase structure of the ICM. The realistic large-scale outflows launched by very small-scale thermal injections, the improved modeling of turbulent diffusion and mixing, and the particle nature of the simulation allow us to carefully separate different heating channels. We show that AGN activity, while efficient at regulating star formation, is incapable of destroying a CC. Instead, that process is facilitated by a head-on, 1:8 mass ratio merger. The merger generates bulk and turbulent motions, which in turn mix high entropy gas generated by AGN and merger driven shocks, turbulent dissipation and sloshing of the ICM by infalling substructures. While central cooling times remain shorter than the Hubble time, restoring a CC is made more difficult by the reduced precipitation rates at larger radii, emphasizing that the AGN-ICM connection is truly a multi-scale problem.

(Abridged) Palomar Gattini-IR is a new wide-field, near-infrared robotic time domain survey operating at Palomar Observatory. Using a 30 cm telescope mounted with a H2RG detector, Gattini-IR achieves a field of view of 25 sq. deg. with a pixel scale of 8.7" in J-band. Here, we describe the system design, survey operations, data processing system and on-sky performance of Palomar Gattini-IR. As a part of the nominal survey, Gattini-IR scans $\approx 7500$ square degrees of the sky every night to a median 5$\sigma$ depth of $15.7$ AB mag outside the Galactic plane. The survey covers $\approx 15000$ square degrees of the sky visible from Palomar with a median cadence of 2 days. A real-time data processing system produces stacked science images from dithered raw images taken on sky, together with PSF-fit source catalogs and transient candidates identified from subtractions within a median delay of $\approx 4$ hours from the time of observation. The calibrated data products achieve an astrometric accuracy (RMS) of $\approx 0.7$" with respect to Gaia DR2 for sources with S/N $> 10$, and better than $\approx 0.35$" for sources brighter than $\approx 12$ Vega mag. The photometric accuracy (RMS) achieved in the PSF-fit source catalogs is better than $\approx 3$% for sources brighter than $\approx 12$ Vega mag, as calibrated against the 2MASS catalog. With a field of view $\approx 40\times$ larger than any other existing near infrared imaging instrument, Gattini-IR is probing the reddest and dustiest transients in the local universe such as dust obscured supernovae in nearby galaxies, novae behind large columns of extinction within the galaxy, reddened micro-lensing events in the Galactic plane and variability from cool and dust obscured stars. We present results from transients and variables identified since the start of the commissioning period.

Using new 610 MHz and 235 MHz observations from the Giant Metrewave Radio Telescope (GMRT) in combination with archival GMRT and Very Large Array (VLA) survey data we present the radio properties of the dominant early-type galaxies in the low$-$richness sub-sample of the Complete Local-volume Groups Sample (CLoGS; 27 galaxy groups) and provide results for the radio properties of the full CLoGS sample for the first time. We find a high radio detection rate in the dominant galaxies of the low-richness sub-sample of 82% (22/27); for the full CLoGS sample, the detection rate is 87% (46/53). The group-dominant galaxies exhibit a wide range of radio power, 10$^{20}$ $-$ 10$^{25}$ W Hz$^{-1}$ in the 235 and 610 MHz bands, with the majority (53%) presenting point-like radio emission, 19% hosting currently active radio jets, 6% having remnant jets, 9% being diffuse and 13% having no detected radio emission. The mean spectral index of the detected radio sources in the 235$-$610 MHz frequency range is found to be $\alpha_{235}^{610}\sim$0.68, and $\alpha_{235}^{1400}\sim$0.59 in the 235$-$1400 MHz one. In agreement with earlier studies, we find that the fraction of ultra-steep spectrum sources ($\alpha>$1.3) is $\sim$4%, mostly dependent on the detection limit at 235 MHz. The majority of point-like systems are found to reside in dynamically young groups, whereas jet systems show no preference between spiral-rich and spiral-poor group environments. The mechanical power of the jet sources in the low$-$richness sample groups is estimated to be $\sim$10$^{42}$ $-$ 10$^{44}$ erg s$^{-1}$ with their black hole masses ranging between 2$\times$10$^{8}$ $-$ 5$\times$10$^{9}$ M$_{\odot}$. We confirm previous findings that, while radio jet sources tend to be associated with more massive black holes, black hole mass is not the decisive factor in determining jet activity or power.

Galaxy clusters observed through the thermal Sunyaev-Zeldovich (tSZ) effect are a recent cosmological probe. The precision on the cosmological constraints is affected mainly by the current knowledge of cluster physics, which enters the analysis through the scaling relations. Here we aim to study one of the most important sources of systematic uncertainties, the mass bias, $b$. We have analysed the effects of a mass-redshift dependence, adopting a power-law parametrisation. We applied this parametrisation to the combination of tSZ number counts and power spectrum, finding a hint of redshift dependence that leads to a decreasing value of the mass bias for higher redshift. We tested the robustness of our results for different mass bias calibrations and a discrete redshift dependence. We find our results to be dependent on the clusters sample that we are considering, in particular obtaining an inverse (decreasing) redshift dependence when neglecting $z<0.2$ clusters. We analysed the effects of this parametrisation on the combination of cosmic microwave background (CMB) primary anisotropies and tSZ galaxy clusters. We find a preferred constant value of mass bias, having $(1-b) =0.62 \pm 0.05$. The corresponding value of $b$ is too high with respect to weak lensing and numerical simulations estimations. Therefore we conclude that this mass-redshift parametrisation does not help in solving the remaining discrepancy between CMB and tSZ clusters observations.

During the last decade, numerous and varied observations, along with increasingly sophisticated numerical simulations, have awakened astronomers to the central role the circumgalactic medium (CGM) plays in regulating galaxy evolution. It contains the majority of the baryonic matter associated with a galaxy, along with most of the metals, and must continually replenish the star forming gas in galaxies that continue to sustain star formation. And while the CGM is complex, containing gas ranging over orders of magnitude in temperature and density, a simple emergent property may be governing its structure and role. Observations increasingly suggest that the ambient CGM pressure cannot exceed the limit at which cold clouds start to condense out and precipitate toward the center of the potential well. If feedback fueled by those clouds then heats the CGM and causes it to expand, the pressure will drop and the "rain" will diminish. Such a feedback loop tends to suspend the CGM at the threshold pressure for precipitation. The coming decade will offer many opportunities to test this potentially fundamental principle of galaxy evolution.

The discovery of GeV gamma-rays from classical novae has led to a reassessment of these garden-variety explosions, and highlighted their importance for understanding radiative shocks, particle acceleration, and dust formation in more exotic, distant transients. Recent collaboration between observers and theorists has revealed that shocks in novae are energetically important, and can even dominate their bolometric luminosity. Shocks may also explain long-standing mysteries in novae such as dust production, super-Eddington luminosities, and `flares' in optical light curves. Here, we highlight the multi-wavelength facilities of the next decade that will further test our nova shock model and fulfill the promise of novae as powerful astrophysical laboratories.

A primary goal of the SCUBA-2 Web survey is to perform tomography of the early inter-galactic medium by studying systems containing some of the brightest quasi-stellar objects (QSOs; 2.5<z<3.0) and nearby submillimetre galaxies. As a first step, this paper aims to characterize the galaxies that host the QSOs. To achieve this, a sample of 13 hyper-luminous (L_AGN>10^14 L_odot) QSOs with previous submillimetre continuum detections were followed up with CO(3-2) observations using the NOEMA interferometer. All but two of the QSOs are detected in CO(3-2); for one non-detection, our observations show a tentative 2sigma line at the expected position and redshift, and for the other non-detection we find only continuum flux density an order of magnitude brighter than the other sources. In three of the fields, a companion potentially suitable for tomography is detected in CO line emission within 25 arcsec of the QSO. We derive gas masses, dynamical masses and far-infrared luminosities, and show that the QSOs in our sample have similar properties as compared to less luminous QSOs and SMGs in the literature, despite the fact that their black-hole masses (which are proportional to L_AGN) are 1-2 orders of magnitude larger. We discuss two interpretations of these observations: this is due to selection effects, such as preferential face-on viewing angles and picking out objects in the tail ends of the scatter in host-galaxy mass and black-hole mass relationships; or the black hole masses have been overestimated because the accretion rates are super-Eddington.

Modelling the turbulent diffusion of thermal energy, momentum, and metals is required in all galaxy evolution simulations due to the ubiquity of turbulence in galactic environments. The most commonly employed diffusion model, the Smagorinsky model, is known to be over-diffusive due to its strong dependence on the fluid velocity shear. We present a method for dynamically calculating a more accurate, locally appropriate, turbulent diffusivity: the dynamic localised Smagorinsky model. We investigate a set of standard astrophysically-relevant hydrodynamical tests, and demonstrate that the dynamic model curbs over-diffusion in non-turbulent shear flows and improves the density contrast in our driven turbulence experiments. In galactic discs, we find that the dynamic model maintains the stability of the disc by preventing excessive angular momentum transport, and increases the metal-mixing timescale in the interstellar medium. In both our isolated Milky Way-like galaxies and cosmological simulations, we find that the interstellar and circumgalactic media are particularly sensitive to the treatment of turbulent diffusion. We also examined the global gas enrichment fractions in our cosmological simulations, to gauge the potential effect on the formation sites and population statistics of Population III stars and supermassive black holes, since they are theorised to be sensitive to the metallicity of the gas out of which they form. The dynamic model is, however, not for galaxy evolution studies only. It can be applied to all astrophysical hydrodynamics simulations, including those modelling stellar interiors, planetary formation, and star formation.

We present the first results from RomulusC, the highest resolution cosmological hydrodynamic simulation of a galaxy cluster run to date. RomulusC, a zoom-in simulation of a halo with $z=0$ mass $10^{14}$ M$_{\odot}$, is run with the same sub-grid physics and resolution as Romulus25 (Tremmel et al. 2017). With unprecedented mass and spatial resolution, RomulusC represents a unique opportunity to study the evolution of galaxies in dense environments down to dwarf masses. We demonstrate that RomulusC results in an intracluster medium (ICM) consistent with observations. The star formation history and stellar mass of the brightest cluster galaxy (BCG) is consistent with observations and abundance matching results, indicating that our sub-grid models, optimized only to reproduce observations of field dwarf and Milky Way mass galaxies, are able to produce reasonable galaxy masses and star formation histories in much higher mass systems. Feedback from supermassive black holes (SMBHs) regulates star formation by driving large-scale, collimated outflows that coexist with a low entropy core. We find that non-BCG cluster member galaxies are substantially quenched compared to the field down to dwarf galaxy masses and, at low masses, quenching is seen to have no dependence on mass or distance from the cluster center. This enhanced quenched population extends beyond $R_{200}$ and is in place at high redshift. Similarly, we predict that SMBH activity is significantly suppressed within clusters outside of the BCG, but show how the effect could be lost when only focusing on the brightest AGN in the most massive galaxies.

The Pan-Andromeda Archaeological Survey is a survey of $>400$ square degrees centered on the Andromeda (M31) and Triangulum (M33) galaxies that has provided the most extensive panorama of a $L_\star$ galaxy group to large projected galactocentric radii. Here, we collate and summarise the current status of our knowledge of the substructures in the stellar halo of M31, and discuss connections between these features. We estimate that the 13 most distinctive substructures were produced by at least 5 different accretion events, all in the last 3 or 4 Gyrs. We suggest that a few of the substructures furthest from M31 may be shells from a single accretion event. We calculate the luminosities of some prominent substructures for which previous estimates were not available, and we estimate the stellar mass budget of the outer halo of M31. We revisit the problem of quantifying the properties of a highly structured dataset; specifically, we use the OPTICS clustering algorithm to quantify the hierarchical structure of M31's stellar halo, and identify three new faint structures. M31's halo, in projection, appears to be dominated by two `mega-structures', that can be considered as the two most significant branches of a merger tree produced by breaking M31's stellar halo into smaller and smaller structures based on the stellar spatial clustering. We conclude that OPTICS is a powerful algorithm that could be used in any astronomical application involving the hierarchical clustering of points. The publication of this article coincides with the public release of all PAndAS data products.

We use high-resolution cosmological hydrodynamical simulations of Milky Way-sized galaxies with varying supernovae feedback strengths and merger histories to investigate the formation of their gaseous halos and especially their hot ($>10^6$~K) X-ray luminous coronae. Our simulations predict the presence of significant hot gas in the halos as early as $z=3-4$, well before the halos ought to be able to sustain hot mode accretion in the conventional picture. The nascent coronae grow inside-out and initially do so primarily as a result of outflows from the central galaxies powered by merger-induced shock heating and strong supernovae feedback, both of which are elemental features of today's successful galaxy formation models. Furthermore, the outflows and the forming coronae also accelerate the transition from cold to hot mode accretion by contributing to the conditions for sustaining stable accretion shocks. They also disrupt the filamentary streams funneling cold gas onto the central galaxies by causing their mouths to fray into a broad delta, detach from the galaxies, and be pushed away to larger radii. And even though at early times the filaments repeatedly re-form, the hot gas and the outflows act to weaken the filaments and accelerate their ultimate disruption. Although galactic outflows are generally thought of as ejective feedback, we find that their action on the filaments suggests a preventive role as well.

We present IRAM 30m and APEX telescope observations of CO(1-0) and CO(2-1) lines in 36 group-dominant early-type galaxies, completing our molecular gas survey of dominant galaxies in the Complete Local-volume Groups Sample. We detect CO emission in 12 of the galaxies at >4sigma significance, with molecular gas masses in the range 0.01-6x10^8 Msol, as well as CO in absorption in the non-dominant group member galaxy NGC 5354. In total 21 of the 53 CLoGS dominant galaxies are detected in CO and we confirm our previous findings that they have low star formation rates (0.01-1 Msol/yr) but short depletion times (<1Gyr) implying rapid replenishment of their gas reservoirs. Comparing molecular gas mass with radio luminosity, we find that a much higher fraction of our group-dominant galaxies (60+-16%) are AGN-dominated than is the case for the general population of ellipticals, but that there is no clear connection between radio luminosity and the molecular gas mass. Using data from the literature, we find that at least 27 of the 53 CLoGS dominant galaxies contain HI, comparable to the fraction of nearby non-cluster early type galaxies detected in HI and significantly higher that the fraction in the Virgo cluster. We see no correlation between the presence of an X-ray detected intra-group medium and molecular gas in the dominant galaxy, but find that the HI-richest galaxies are located in X-ray faint groups. Morphological data from the literature suggests the cold gas component most commonly takes the form of a disk, but many systems show evidence of galaxy-galaxy interactions, indicating that they may have acquired their gas through stripping or mergers. We provide improved molecular gas mass estimates for two galaxies previously identified as being in the centres of cooling flows, NGC 4636 and NGC 5846, and find that they are relatively molecular gas poor compared to our other detected systems.

[abridged] Aims: We test the effects of re-orienting jets from an active galactic nucleus (AGN) on the intracluster medium in a galaxy cluster environment with short central cooling time. We investigate appearance and properties of the resulting cavities, and the efficiency of jets in providing near-isotropic heating to the cooling cluster core. Methods: We use numerical simulations to explore four models of jets over several active/inactive cycles. We keep the jet power and duration fixed, varying only the jet angle prescription. We track the total energy of the intracluster medium (ICM) in the cluster core over time, and the fraction of the jet energy transferred to the ICM, paying attention to where the energy is deposited. We also compare synthetic X-ray images of the simulated cluster to actual observations. Results: Jets whose re-orientation is minimal ($\lesssim 20^{\circ}$) typically produce conical structures of interconnected cavities, with the opening angle of the cones being $\sim 15-20^{\circ}$, extending to $\sim 300$ kpc from the cluster centre. Such jets transfer about $60\%$ of their energy to the ICM, yet they are not very efficient at heating the cluster core, as the jet energy is deposited further out. Jets that re-orient by $\gtrsim 20^{\circ}$ generally produce multiple pairs of detached cavities. Although smaller, these cavities are inflated within the central 50~kpc and are more isotropically distributed, resulting in more effective heating of the core. Such jets, over few hundreds Myr, can deposit up to $80\%$ of their energy where it is required. Consequently, these models come the closest to an heating/cooling balance and to mitigating runaway cooling of the core, even though all models have identical power/duration profiles. Additionally, the corresponding synthetic X-ray images exhibit structures closely resembling those seen in real cool-core clusters.

We present a study of the radio properties of the dominant early-type galaxies in 26 galaxy groups, the high-richness sub-sample of the Complete Local-volume Groups Sample (CLoGS). Combining new 610 MHz and 235 MHz observations of 21 groups from the Giant Metrewave Radio Telescope (GMRT) with archival GMRT and Very Large Array (VLA) survey data, we find a high detection rate, with 92% of the dominant galaxies hosting radio sources. The sources have a wide range of luminosities, 10$^{20}$ $-$ 10$^{24}$ W/Hz in the 235 and 610 MHz bands. The majority (54%) are point-like, but 23% have radio jets, and another 15% are diffuse radio sources with no clear jet/lobe morphology. The spectral index of the detected radio sources ranges from very flat values of ~0.2 to typical radio synchrotron spectra of ~0.9 with only two presenting steep radio spectra with $\alpha_{235}^{610}$ > 1. We find that jet sources are more common in X-ray bright groups, with radio non-detections found only in X-ray faint systems. Radio point sources appear in all group environments irrespective of their X-ray properties or spiral fraction. We estimate the mechanical power (Pcav) of the jet sources in the X-ray bright groups to be 10$^{41}$ $-$ 10$^{43}$ erg/s, with the two large-scale jet systems (NGC 193 and NGC 4261) showing jet powers two orders of magnitude greater than the radiative losses from the cool cores of their groups. This suggests that central AGN are not always in balance with cooling, but may instead produce powerful periodical bursts of feedback heating.

The recent EDGES collaboration detection of an absorption signal at a central frequency of $\nu = 78 \pm 1$ MHz points to the presence of a significant Lyman-$\alpha$ background by a redshift of $z=18$. The timing of this signal constrains the dark matter particle mass ($m_\chi$) in the warm dark matter (WDM) cosmological model. WDM delays the formation of small-scale structures, and therefore a stringent lower limit can be placed on $m_\chi,$ based on the presence of a sufficiently strong Ly-$\alpha$ background due to star formation at $z=18$. Our results show that the coupling the spin temperature to the gas through Ly-$\alpha$ pumping requires a minimum mass of $m_\chi>3$ keV if atomic cooling halos dominate the star formation rate at $z=18,$ and $m_\chi>2$ keV if ${\rm\,H_2}$ cooling halos also form stars efficiently at this redshift. These limits match or exceed the most stringent limits cited to date in the literature, even in the face of the many uncertainties regarding star-formation at high redshift.

Massive galaxy clusters are now found as early as 3 billion years after the Big Bang, containing stars that formed at even earlier epochs. The high-redshift progenitors of these galaxy clusters, termed 'protoclusters', are identified in cosmological simulations with the highest dark matter overdensities. While their observational signatures are less well defined compared to virialized clusters with a substantial hot intra-cluster medium (ICM), protoclusters are expected to contain extremely massive galaxies that can be observed as luminous starbursts. Recent claimed detections of protoclusters hosting such starbursts do not support the kind of rapid cluster core formation expected in simulations because these structures contain only a handful of starbursting galaxies spread throughout a broad structure, with poor evidence for eventual collapse into a protocluster. Here we report that the source SPT2349-56 consists of at least 14 gas-rich galaxies all lying at z = 4.31 based on sensitive observations of carbon monoxide and ionized carbon. We demonstrate that each of these galaxies is forming stars between 50 and 1000 times faster than our own Milky Way, and all are located within a projected region only $\sim$ 130 kiloparsecs in diameter. This galaxy surface density is more than 10 times the average blank field value (integrated over all redshifts) and $>$1000 times the average field volume density. The velocity dispersion ($\sim$ 410 km s$^{-1}$) of these galaxies and enormous gas and star formation densities suggest that this system represents a galaxy cluster core at an advanced stage of formation when the Universe was only 1.4 billion years old. A comparison with other known protoclusters at high redshifts shows that SPT2349-56 is a uniquely massive and dense system that could be building one of the most massive structures in the Universe today.

Recent analysis shows that it is important to explicitly include the gravitational potential of the central brightest central galaxy (BCG) to infer the acceleration due to gravity ($g$) and the free-fall time ($t_{\rm ff} \equiv [2r/g]^{1/2}$) in cool cluster cores. Accurately measuring $t_{\rm ff}$ is crucial because according to numerical simulations cold gas condensation and strong feedback occur in cluster cores with min($t_{\rm cool}/t_{\rm ff}$) below a threshold value close to 10. Recent observations which include the BCG gravity show that the observed threshold in min($t_{\rm cool}/t_{\rm ff}$) lies at a somewhat higher value, close to 10-30; there are only a few clusters in which this ratio falls much below 10. In this paper we compare numerical simulations of feedback AGN (Active Galactic Nuclei) jets interacting with the intracluster medium (ICM), with and without a BCG potential. We find that, for a fixed feedback efficiency, the presence of a BCG does not significantly affect the temperature but increases (decreases) the core density (entropy) on average. Most importantly, min($t_{\rm cool}/t_{\rm ff}$) is only affected slightly by the inclusion of the BCG gravity. Also notable is that the lowest value of min($t_{\rm cool}/t_{\rm ff}$) in the NFW+BCG runs are about twice larger than in the NFW runs. We also look at the role of depletion of cold gas due to star formation and show that it only affects the rotationally dominant component (torus), while the radially dominant component (which regulates the feedback cycle) remains largely unaffected. The distribution of metals due to AGN jets in our simulations is predominantly along the jet direction and the radial spread of metals is less. We also show that the turbulence in cool core clusters is weak, consistent with recent Hitomi results on Perseus cluster.

We analyse spatially-resolved deep optical spectroscopy of Brightest Cluster Galaxies (BCGs) located in 32 massive clusters with redshifts of 0.05 $\leq z \leq$ 0.30, to investigate their velocity dispersion profiles. We compare these measurements to those of other massive early-type galaxies, as well as central group galaxies, where relevant. This unique, large sample extends to the most extreme of massive galaxies, spanning M$_{K}$ between -25.7 to -27.8 mag, and host cluster halo mass M$_{500}$ up to 1.7 $\times$ 10$^{15}$ M$_{\odot}$. To compare the kinematic properties between brightest group and cluster members, we analyse similar spatially-resolved long-slit spectroscopy for 23 nearby Brightest Group Galaxies (BGGs) from the Complete Local-Volume Groups Sample (CLoGS). We find a surprisingly large variety in velocity dispersion slopes for BCGs, with a significantly larger fraction of positive slopes, unique compared to other (non-central) early-type galaxies as well as the majority of the brightest members of the groups. We find that the velocity dispersion slopes of the BCGs and BGGs correlate with the luminosity of the galaxies, and we quantify this correlation. It is not clear whether the full diversity in velocity dispersion slopes that we see is reproduced in simulations.

Galaxy clusters are expected to form hierarchically in a LCDM universe, growing primarily through mergers with lower mass clusters and the continual accretion of group-mass halos. Galaxy clusters assemble late, doubling their masses since z~0.5, and so the outer regions of clusters should be replete with infalling group-mass systems. We present an XMM-Newton survey to search for X-ray groups in the infall regions of 23 massive galaxy clusters at z~0.2, identifying 39 X-ray groups that have been spectroscopically confirmed to lie at the cluster redshift. These groups have mass estimates in the range 2x10^13-7x10^14Msun, and group-to-cluster mass ratios as low as 0.02. The comoving number density of X-ray groups in the infall regions is ~25x higher than that seen for isolated X-ray groups from the XXL survey. The average mass per cluster contained within these X-ray groups is 2.2x10^14Msun, or 19% of the mass within the primary cluster itself. We estimate that ~10^15Msun clusters increase their masses by 16% between z=0.223 and the present day due to the accretion of groups with M200>10^13.2Msun. This represents about half of the expected mass growth rate of clusters at these late epochs. The other half is likely to come from smooth accretion of matter not bound in halos. The mass function of the infalling X-ray groups appears significantly top-heavy with respect to that of field X-ray systems, consistent with expectations from numerical simulations, and the basic consequences of collapsed massive dark matter halos being biased tracers of the underlying large-scale density distribution.

We present the Complete Local-Volume Groups Sample (CLoGS), a statistically complete optically-selected sample of 53 groups within 80 Mpc. Our goal is to combine X-ray, radio and optical data to investigate the relationship between member galaxies, their active nuclei, and the hot intra-group medium (IGM). We describe sample selection, define a 26-group high-richness subsample of groups containing at least 4 optically bright (log L_B>=10.2 LBsol) galaxies, and report the results of XMM-Newton and Chandra observations of these systems. We find that 14 of the 26 groups are X-ray bright, possessing a group-scale IGM extending at least 65kpc and with luminosity >10^41 erg/s, while a further 3 groups host smaller galaxy-scale gas halos. The X-ray bright groups have masses in the range M_500=0.5-5x10^13 Msol, based on system temperatures of 0.4-1.4 keV, and X-ray luminosities in the range 2-200x10^41 erg/s. We find that ~53-65% of the X-ray bright groups have cool cores, a somewhat lower fraction than found by previous archival surveys. Approximately 30% of the X-ray bright groups show evidence of recent dynamical interactions (mergers or sloshing), and ~35% of their dominant early-type galaxies host AGN with radio jets. We find no groups with unusually high central entropies, as predicted by some simulations, and confirm that CLoGS is in principle capable of detecting such systems. We identify three previously unrecognized groups, and find that they are either faint (L_X,R500<10^42 erg/s) with no concentrated cool core, or highly disturbed. This leads us to suggest that ~20% of X-ray bright groups in the local universe may still be unidentified.

Several arguments suggest that stochastic condensation of cold gas and its accretion onto the central supermassive black hole (SMBH) is essential for active galactic nuclei (AGN) feedback to work in the most massive galaxies that lie at the centres of galaxy clusters. Our 3-D hydrodynamic AGN jet-ICM (intracluster medium) simulations, looking at the detailed angular momentum distribution of cold gas and its time variability for the first time, show that the angular momentum of the cold gas crossing $\lesssim 1$ kpc is essentially isotropic. With almost equal mass in clockwise and counter-clockwise orientations, we expect a cancellation of angular momentum on roughly the dynamical time. This means that a compact accretion flow with a short viscous time ought to form, through which enough accretion power can be channeled into jet mechanical energy sufficiently quickly to prevent a cooling flow. The inherent stochasticity, expected in feedback cycles driven by cold gas condensation, gives rise to a large variation in the cold gas mass at the centres of galaxy clusters, for similar cluster and SMBH masses, in agreement with the observations. Such correlations are expected to be much tighter for the smoother hot/Bondi accretion. The weak correlation between cavity power and Bondi power obtained from our simulations also match observations.

Observations with the Chandra X-ray Observatory are used to examine the hot gas properties within a sample of 10 galaxy groups selected from the Galaxy And Mass Assembly survey's optical Friends-of-Friends group catalogue. Our groups have been screened to eliminate spurious and unrelaxed systems, and the effectiveness of this procedure is demonstrated by the detection of intergalactic hot gas in 80 per cent of our sample. However, we find that 9 of the 10 are X-ray underluminous by a mean factor of $\sim$4 compared to typical X-ray-selected samples. Consistent with this, the majority of our groups have gas fractions that are lower and gas entropies somewhat higher than those seen in typical X-ray-selected samples. Two groups, which have high 2{\sigma} lower limits on their gas entropy, are candidates for the population of high-entropy groups predicted by some active galactic nucleus feedback models.

Radio halos are diffuse synchrotron sources on scales of ~1 Mpc that are found in merging clusters of galaxies, and are believed to be powered by electrons re-accelerated by the merger-driven turbulence. We present measurements of extended radio emission on similarly large scales in two clusters of galaxies hosting cool cores: Abell 2390 and Abell 2261. The analysis is based on interferometric imaging with the JVLA, VLA and GMRT. We present detailed radio images of the targets, subtract the compact emission components, and measure the spectral indices for the diffuse components. The radio emission in A2390 extends beyond a known sloshing-like brightness discontinuity, and has a very steep in-band spectral slope at 1.5 GHz that is similar to some known ultra-steep spectrum radio halos. The diffuse signal in A2261 is more extended than in A2390 but has lower luminosity. X-ray morphological indicators, derived from XMM-Newton X-ray data, place these clusters in the category of relaxed or regular systems, although some asymmetric features that can indicate past minor mergers are seen in the X-ray brightness images. If these two Mpc-scale radio sources are categorized as giant radio halos, they question the common assumption of radio halos occurring exclusively in clusters undergoing violent merging activity, in addition to commonly used criteria in distinguishing between radio halos and mini-halos.

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.

We present the first results from an extended survey of the Andromeda galaxy (M31) using 41.1 h of observations by Spitzer-IRAC at 3.6 and 4.5$\,\mu$m. This survey extends previous observations to the outer disc and halo, covering total lengths of 4$.\!\!^\circ$4 and 6$.\!\!^\circ$6 along the minor and major axes, respectively. We have produced surface brightness profiles by combining the integrated light from background-corrected maps with stellar counts from a new catalogue of point sources. Using auxiliary catalogues we have carried out a statistical analysis in colour-magnitude space to discriminate M31 objects from foreground Milky Way stars and background galaxies. The catalogue includes 426,529 sources, of which 66 per cent have been assigned probability values to identify M31 objects with magnitude depths of [3.6]$\,=\,$19.0$\,\pm\,$0.2, [4.5]$\,=\,$18.7$\,\pm\,$0.2. We discuss applications of our data for constraining the stellar mass and characterising point sources in the outer radii.

We investigate the properties of halo gas using three cosmological `zoom-in' simulations of realistic Milky Way-galaxy analogs with varying sub-grid physics. In all three cases, the mass of hot ($T > 10^6$ K) halo gas is $\sim 1\%$ of the host's virial mass. The X-ray luminosity of two of the runs is consistent with observations of the Milky Way, while the third simulation is X-ray bright and resembles more closely a very massive, star-forming spiral. Hot halos extend to 140 kpc from the galactic center and are surrounded by a bubble of warm-hot ($T = 10^5 - 10^6$K) gas that extends to the virial radius. Simulated halos agree well outside 20-30 kpc with the $\beta$-model of Miller & Bregman (2013) based on OVII absorption and OVIII emission measurements. Warm--hot and hot gas contribute up to $80\%$ of the total gas reservoir, and contain nearly the same amount of baryons as the stellar component. The mass of warm-hot and hot components falls into the range estimated for $L^*$ galaxies. With key observational constraints on the density of the Milky Way corona being satisfied, we show that concealing of the ubiquitous warm-hot baryons, along with the ejection of just $20-30 \%$ of the diffuse gas out of the potential wells by supernova-driven outflows, can solve the "missing baryon problem". The recovered baryon fraction within 3 virial radii is $90\%$ of the universal value. With a characteristic density of $\sim 10^{-4}$ cm$^{-3}$ at $50-80$ kpc, diffuse coronae meet the requirement for fast and complete ram--pressure stripping of the gas reservoirs in dwarf galaxy satellites.

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.

A fraction of brightest cluster galaxies (BCGs) shows bright emission in the UV and the blue part of the optical spectrum, which has been interpreted as evidence of recent star formation. Most of these results are based on the analysis of broadband photometric data. Here, we study the optical spectra of a sample of 19 BCGs hosted by X-ray luminous galaxy clusters at 0.15 < z < 0.3, a subset from the Canadian Cluster Comparison Project (CCCP) sample. We identify plausible star formation histories of the galaxies by fitting Simple Stellar Populations (SSPs) as well as composite populations, consisting of a young stellar component superimposed on an intermediate/old stellar component, to accurately constrain their star formation histories. We detect prominent young (~200 Myr) stellar populations in 4 of the 19 galaxies. Of the four, the BCG in Abell 1835 shows remarkable A-type stellar features indicating a relatively large population of young stars, which is extremely unusual even amongst star forming BCGs. We constrain the mass contribution of these young components to the total stellar mass to be typically between 1% to 3%, but rising to 7% in Abell 1835. We find that the four of the BCGs with strong evidence for recent star formation (and only these four galaxies) are found within a projected distance of 5 kpc of their host cluster's X-ray peak, and the diffuse, X-ray gas surrounding the BCGs exhibit a ratio of the radiative cooling-to-free-fall time ($t_{c}/t_{ff}$) of < 10. These are also some of the clusters with the lowest central entropy. Our results are consistent with the predictions of the precipitation-driven star formation and AGN feedback model, in which the radiatively cooling diffuse gas is subject to local thermal instabilities once the instability parameter $t_{c}/t_{ff}$ falls below ~10, leading to the condensation and precipitation of cold gas. [Abridged]

We test the assumption of hydrostatic equilibrium in an X-ray luminosity selected sample of 50 galaxy clusters at $0.15<z<0.3$ from the Local Cluster Substructure Survey (LoCuSS). Our weak-lensing measurements of $M_{500}$ control systematic biases to sub-4 per cent, and our hydrostatic measurements of the same achieve excellent agreement between XMM-Newton and Chandra. The mean ratio of X-ray to lensing mass for these 50 clusters is $\beta_{\rm X}=0.95\pm0.05$, and for the 44 clusters also detected by Planck, the mean ratio of Planck mass estimate to LoCuSS lensing mass is $\beta_{\rm P}=0.95\pm0.04$. Based on a careful like-for-like analysis, we find that LoCuSS, the Canadian Cluster Comparison Project (CCCP), and Weighing the Giants (WtG) agree on $\beta_{\rm P}\simeq0.9-0.95$ at $0.15<z<0.3$. This small level of hydrostatic bias disagrees at $\sim5\sigma$ with the level required to reconcile Planck cosmology results from the cosmic microwave background and galaxy cluster counts.

Using high-resolution 3-D and 2-D (axisymmetric) hydrodynamic simulations in spherical geometry, we study the evolution of cool cluster cores heated by feedback-driven bipolar active galactic nuclei (AGN) jets. Condensation of cold gas, and the consequent enhanced accretion, is required for AGN feedback to balance radiative cooling with reasonable efficiencies, and to match the observed cool core properties. A feedback efficiency (mechanical luminosity $\approx \epsilon \dot{M}_{\rm acc} c^2$; where $\dot{M}_{\rm acc}$ is the mass accretion rate at 1 kpc) as small as $5 \times 10^{-5}$ is sufficient to reduce the cooling/accretion rate by $\sim 10$ compared to a pure cooling flow. This value is smaller compared to the ones considered earlier, and is consistent with the jet efficiency and the fact that only a small fraction of gas at 1 kpc is accreted on to the supermassive black hole (SMBH). We find hysteresis cycles in all our simulations with cold mode feedback: {\em condensation} of cold gas when the ratio of the cooling-time to the free-fall time ($t_{\rm cool}/t_{\rm ff}$) is $\lesssim 10$ leads to a sudden enhancement in the accretion rate; a large accretion rate causes strong jets and {\em overheating} of the hot ICM such that $t_{\rm cool}/t_{\rm ff} > 10$; further condensation of cold gas is suppressed and the accretion rate falls, leading to slow cooling of the core and condensation of cold gas, restarting the cycle. Therefore, there is a spread in core properties, such as the jet power, accretion rate, for the same value of core entropy or $t_{\rm cool}/t_{\rm ff}$. A fewer number of cycles are observed for higher efficiencies and for lower mass halos because the core is overheated to a longer cooling time. The 3-D simulations show the formation of a few-kpc scale, rotationally-supported, massive ($\sim 10^{11} M_\odot$) cold gas torus. (abstract abridged)

We present a study of the spatial distribution and kinematics of star-forming galaxies in 30 massive clusters at 0.15<z<0.30, combining wide-field Spitzer 24um and GALEX NUV imaging with highly-complete spectroscopy of cluster members. The fraction (f_SF) of star-forming cluster galaxies rises steadily with cluster-centric radius, increasing fivefold by 2r200, but remains well below field values even at 3r200. This suppression of star formation at large radii cannot be reproduced by models in which star formation is quenched in infalling field galaxies only once they pass within r200 of the cluster, but is consistent with some of them being first pre-processed within galaxy groups. Despite the increasing f_SF-radius trend, the surface density of star-forming galaxies actually declines steadily with radius, falling ~15x from the core to 2r200. This requires star-formation to survive within recently accreted spirals for 2--3Gyr to build up the apparent over-density of star-forming galaxies within clusters. The velocity dispersion profile of the star-forming galaxy population shows a sharp peak of 1.44-sigma_v at 0.3r500, and is 10--35% higher than that of the inactive cluster members at all cluster-centric radii, while their velocity distribution shows a flat, top-hat profile within r500. All of these results are consistent with star-forming cluster galaxies being an infalling population, but one that must also survive ~0.5--2Gyr beyond passing within r200. By comparing the observed distribution of star-forming galaxies in the stacked caustic diagram with predictions from the Millennium simulation, we obtain a best-fit model in which SFRs decline exponentially on quenching time-scales of 1.73\pm0.25 Gyr upon accretion into the cluster.

Masses of clusters of galaxies from weak gravitational lensing analyses of ever larger samples are increasingly used as the reference to which baryonic scaling relations are compared. In this paper we revisit the analysis of a sample of 50 clusters studied as part of the Canadian Cluster Comparison Project. We examine the key sources of systematic error in cluster masses. We quantify the robustness of our shape measurements and calibrate our algorithm empirically using extensive image simulations. The source redshift distribution is revised using the latest state-of-the-art photometric redshift catalogs that include new deep near-infrared observations. Nonetheless we find that the uncertainty in the determination of photometric redshifts is the largest source of systematic error for our mass estimates. We use our updated masses to determine b, the bias in the hydrostatic mass, for the clusters detected by Planck. Our results suggest 1-b=0.76+-0.05(stat)}+-0.06(syst)}, which does not resolve the tension with the measurements from the primary cosmic microwave background.

We present IRAM 30m telescope observations of the CO(1-0) and (2-1) lines in a sample of 11 group-dominant elliptical galaxies selected from the CLoGS nearby groups sample. Our observations confirm the presence of molecular gas in 4 of the 11 galaxies at >4 sigma significance, and combining these with data from the literature we find a detection rate of 43+-14%, comparable to the detection rate for nearby radio galaxies, suggesting that group-dominant ellipticals may be more likely to contain molecular gas than their non-central counterparts. Those group-dominant galaxies which are detected typically contain ~2x10^8 Msol of molecular gas, and although most have low star formation rates (<1 Msol/yr) they have short depletion times, indicating that the gas must be replenished on timescales ~100 Myr. Almost all of the galaxies contain active nuclei, and we note while the data suggest that CO may be more common in the most radio-loud galaxies, the mass of molecular gas required to power the active nuclei through accretion is small compared to the masses observed. We consider possible origin mechanisms for the gas, through cooling of stellar ejecta within the galaxies, group-scale cooling flows, and gas-rich mergers, and find probable examples of each type within our sample, confirming that a variety of processes act to drive the build up of molecular gas in group-dominant ellipticals.

We present a deep Suzaku observation of H1821+643, an extremely rare example of a powerful quasar hosted by the central massive galaxy of a rich cooling-core cluster of galaxies. Informed by previous Chandra studies of the cluster, we achieve a spectral separation of emission from the active galactic nucleus (AGN) and the intracluster medium (ICM). With a high degree of confidence, we identify the signatures of X-ray reflection/reprocessing by cold and slowly moving material in the AGN's immediate environment. The iron abundance of this matter is found to be significantly sub-solar (Z~0.4Zsun), an unusual finding for powerful AGN but in line with the idea that this quasar is feeding from the ICM via a Compton-induced cooling flow. We also find a subtle soft excess that can be described phenomenologically (with an additional black body component) or as ionized X-ray reflection from the inner regions of a high inclination (i=57 degrees) accretion disk around a spinning (a>0.4) black hole. We describe how the ionization state of the accretion disk can be used to constrain the Eddington fraction of the source. Applying these arguments to our spectrum implies an Eddington fraction of 0.25-0.5, with an associated black hole mass of 3-6x10^9Msun.

We present a study of multiwavelength X-ray and weak lensing scaling relations for a sample of 50 clusters of galaxies. Our analysis combines Chandra and XMM-Newton data using an energy-dependent cross-calibration. After considering a number of scaling relations, we find that gas mass is the most robust estimator of weak lensing mass, yielding 15 +/- 6% intrinsic scatter at r500 (the pseudo-pressure YX has a consistent scatter of 22%+/-5%). The scatter does not change when measured within a fixed physical radius of 1 Mpc. Clusters with small BCG to X-ray peak offsets constitute a very regular population whose members have the same gas mass fractions and whose even smaller <10% deviations from regularity can be ascribed to line of sight geometrical effects alone. Cool-core clusters, while a somewhat different population, also show the same (<10%) scatter in the gas mass-lensing mass relation. There is a good correlation and a hint of bimodality in the plane defined by BCG offset and central entropy (or central cooling time). The pseudo-pressure YX does not discriminate between the more relaxed and less relaxed populations, making it perhaps the more even-handed mass proxy for surveys. Overall, hydrostatic masses underestimate weak lensing masses by 10% on the average at r500; but cool-core clusters are consistent with no bias, while non-cool-core clusters have a large and constant 15-20% bias between r2500 and r500, in agreement with N-body simulations incorporating unthermalized gas. For non-cool-core clusters, the bias correlates well with BCG ellipticity. We also examine centroid shift variance and and power ratios to quantify substructure; these quantities do not correlate with residuals in the scaling relations. Individual clusters have for the most part forgotten the source of their departures from self-similarity.

We present the first scaling relation between weak-lensing galaxy cluster mass, $M_{WL}$, and near-infrared luminosity, $L_K$. Our results are based on 17 clusters observed with wide-field instruments on Subaru, the United Kingdom Infrared Telescope, the Mayall Telescope, and the MMT. We concentrate on the relation between projected 2D weak-lensing mass and spectroscopically confirmed luminosity within 1Mpc, modelled as $M_{WL} \propto L_{K}^b$, obtaining a power law slope of $b=0.83^{+0.27}_{-0.24}$ and an intrinsic scatter of $\sigma_{lnM_{WL}|L_{K}}=10^{+8}_{-5}\%$. Intrinsic scatter of ~10% is a consistent feature of our results regardless of how we modify our approach to measuring the relationship between mass and light. For example, deprojecting the mass and measuring both quantities within $r_{500}$, that is itself obtained from the lensing analysis, yields $\sigma_{lnM_{WL}|L_{K}}=10^{+7}_{-5}\%$ and $b=0.97^{+0.17}_{-0.17}$. We also find that selecting members based on their (J-K) colours instead of spectroscopic redshifts neither increases the scatter nor modifies the slope. Overall our results indicate that near-infrared luminosity measured on scales comparable with $r_{500}$ (typically 1Mpc for our sample) is a low scatter and relatively inexpensive proxy for weak-lensing mass. Near-infrared luminosity may therefore be a useful mass proxy for cluster cosmology experiments.

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 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 the results from a weak gravitational lensing study of the merging cluster A520 based on the analysis of Hubble Space Telescope/Advanced Camera for Surveys (ACS) data. The excellent data quality allows us to reach a mean number density of source galaxies of ~109 per sq. arcmin, which improves both resolution and significance of the mass reconstruction compared to a previous study based on Wide Field Planetary Camera 2 (WFPC2) images. We take care in removing instrumental effects such as the trailing of charge due to radiation damage of the ACS detector and the position-dependent point spread function (PSF). This new ACS analysis confirms the previous claims that a substantial amount of dark mass is present between two luminous subclusters. We examine the distribution of cluster galaxies and observe very little light at this location. We find that the centroid of the dark peak in the current ACS analysis is offset to the southwest by ~1 arcmin with respect to the centroid from the WFPC2 analysis. Interestingly, this new centroid is in better spatial agreement with the location where the X-ray emission is strongest, and the mass-to-light ratio estimated with this centroid is much higher 813+-78 M_sun/L_Rsun than the previous value; the aperture mass based on the WFPC2 centroid provides a slightly lower, but consistent mass. Although we cannot provide a definite explanation for the presence of the dark peak, we discuss a revised scenario, wherein dark matter with a more conventional range sigma_DM/m_DM < 1 cm^2/g of self-interacting cross-section can lead to the detection of this dark substructure. If supported by detailed numerical simulations, this hypothesis opens up the possibility that the A520 system can be used to establish a lower limit of the self-interacting cross-section of dark matter.

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.

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 an analysis of the levels and evolution of star formation activity in a representative sample of 30 massive galaxy clusters at 0.15<z<0.30 from the Local Cluster Substructure Survey (LoCuSS), combining wide-field Spitzer 24um data with extensive spectroscopy of cluster members. The specific-SFRs of massive (M>10^10 M_sun) star-forming cluster galaxies within r200 are found to be systematically 28% lower than their counterparts in the field at fixed stellar mass and redshift, a difference significant at the 8.7-sigma level. This is the unambiguous signature of star formation in most (and possibly all) massive star-forming galaxies being slowly quenched upon accretion into massive clusters, their SFRs declining exponentially on quenching time-scales in the range 0.7-2.0 Gyr. We measure the mid-infrared Butcher-Oemler effect over the redshift range 0.0-0.4, finding rapid evolution in the fraction (f_SF) of massive (M_K<-23.1) cluster galaxies within r200 with SFRs>3M_sun/yr, of the form f_SF (1+z)^7.6. We dissect the origins of the Butcher-Oemler effect, revealing it to be due to the combination of a ~3x decline in the mean specific-SFRs of star-forming cluster galaxies since z~0.3 with a ~1.5x decrease in number density. Two-thirds of this reduction in the specific-SFRs of star-forming cluster galaxies is due to the steady cosmic decline in the specific-SFRs among those field galaxies accreted into the clusters. The remaining one-third reflects an accelerated decline in the star formation activity of galaxies within clusters. The slow quenching of star-formation in cluster galaxies is consistent with a gradual shut down of star formation in infalling spiral galaxies as they interact with the intra-cluster medium via ram-pressure stripping or starvation mechanisms. We find no evidence for the build-up of cluster S0 bulges via major nuclear star-burst episodes.

The advent of high-resolution imaging of galaxy clusters using the Sunyaev-Zel'dovich Effect (SZE) provides a unique probe of the astrophysics of the intracluster medium (ICM) out to high redshifts. To investigate the effects of cluster mergers on resolved SZE images, we present a high-resolution cosmological simulation of a 1.5E15 M_sun adiabatic cluster using the TreeSPH code ChaNGa. This massive cluster undergoes a 10:3:1 ratio triple merger accompanied by a dramatic rise in its integrated Compton-Y, peaking at z = 0.05. By modeling the thermal SZE (tSZ) and kinetic SZE (kSZ) spectral distortions of the Cosmic Microwave Background (CMB) at this redshift with relativistic corrections, we produce various mock images of the cluster at frequencies and resolutions achievable with current high-resolution SZE instruments. The two gravitationally-bound merging subclusters account for 10% and 1% of the main cluster's integrated Compton-Y, and have extended merger shock features in the background ICM visible in our mock images. We show that along certain projections and at specific frequencies, the kSZ CMB intensity distortion can dominate over the tSZ due to the large line of sight velocities of the subcluster gas and the unique frequency-dependence of these effects. We estimate that a one-velocity assumption in estimation of line of sight velocities of the merging subclusters from the kSZ induces a bias of ~10%. This velocity bias is small relative to other sources of uncertainty in observations, partially due to helpful bulk motions in the background ICM induced by the merger. Our results show that high-resolution SZE observations, which have recently detected strong kSZ signals in subclusters of merging systems, can robustly probe the dynamical as well as the thermal state of the ICM.

AGN jets carry more than sufficient energy to stave off catastrophic cooling of the intracluster medium (ICM) in the cores of cool-core clusters. However, in order to prevent catastrophic cooling, the ICM must be heated in a near-isotropic fashion and narrow bipolar jets with $P_{\rm jet}=10^{44-45}$ ergs/s, typical of radio AGNs at cluster centres, are inefficient at heating the gas in the transverse direction to the jets. We argue that due to existent conditions in cluster cores, the SMBHs will, in addition to accreting gas via radiatively inefficient flows, experience short stochastic episodes of enhanced accretion via thin discs. In general, the orientation of these accretion discs will be misaligned with the spin axis of the black holes and the ensuing torques will cause the black hole's spin axis (and therefore, the jet axis) to slew and rapidly change direction. This model not only explains recent observations showing successive generations of jet-lobes-bubbles in individual cool-core clusters that are offset from each other in the angular direction with respect to the cluster center, but also shows that AGN jets {\it can} heat the cluster core nearly isotropically on the gas cooling timescale. Our model {\it does} require that the SMBHs at the centers of cool-core clusters be spinning relatively slowly. Torques from individual misaligned discs are ineffective at tilting rapidly spinning black holes by more than a few degrees. Additionally, since SMBHs that host thin accretion discs will manifest as quasars, we predict that roughly 1--2 rich clusters within $z<0.5$ should have quasars at their centers.

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.

The Canadian Cluster Comparison Project is a comprehensive multi-wavelength survey targeting 50 massive X-ray selected clusters of galaxies to examine baryonic tracers of cluster mass and to probe the cluster-to-cluster variation in the thermal properties of the hot intracluster medium. In this paper we present the weak lensing masses, based on the analysis of deep wide-field imaging data obtained using the Canada-France-Hawaii-Telescope. The final sample includes two additional clusters that were located in the field-of-view. We take these masses as our reference for the comparison of cluster properties at other wavelengths. In this paper we limit the comparison to published measurements of the Sunyaev-Zel'dovich effect. We find that this signal correlates well with the projected lensing mass, with an intrinsic scatter of 12\pm5% at ~r_2500, demonstrating it is an excellent proxy for cluster mass.

We analyze deep g' and r' band data of 97 galaxy clusters imaged with MegaCam on the Canada-France-Hawaii telescope. We compute the number of luminous (giant) and faint (dwarf) galaxies using criteria based on the definitions of de Lucia et al. (2007). Due to excellent image quality and uniformity of the data and analysis, we probe the giant-to-dwarf ratio (GDR) out to z ~ 0.55. With X-ray temperature (Tx) information for the majority of our clusters, we constrain, for the first time, the Tx-corrected giant and dwarf evolution separately. Our measurements support an evolving GDR over the redshift range 0.05 < z < 0.55. We show that modifying the (g'-r'), m_r' and K-correction used to define dwarf and giant selection do not alter the conclusion regarding the presence of evolution. We parameterize the GDR evolution using a linear function of redshift (GDR = alpha * z + beta) with a best fit slope of alpha = 0.88 +/- 0.15 and normalization beta = 0.44 +/- 0.03. Contrary to claims of a large intrinsic scatter, we find that the GDR data can be fully accounted for using observational errors alone. Consistently, we find no evidence for a correlation between GDR and cluster mass (via Tx or weak lensing). Lastly, the data suggest that the evolution of the GDR at z < 0.2 is driven primarily by dry merging of the massive giant galaxies, which when considered with previous results at higher redshift, suggests a change in the dominant mechanism that mediates the GDR.

IRAS 09104+4109 is a rare example of a dust enshrouded type 2 QSO in the centre of a cool-core galaxy cluster. Previous observations of this z=0.44 system showed that as well as powering the hyper-luminous infrared emission of the cluster-central galaxy, the QSO is associated with a double-lobed radio source. However, the steep radio spectral index and misalignment between the jets and ionised optical emission suggested that the orientation of the QSO had recently changed. We use a combination of new, multi-band Giant Metrewave Radio Telescope observations and archival radio data to confirm that the jets are no longer powered by the QSO, and estimate their age to be 120-160 Myr. This is in agreement with the ~70-200 Myr age previously estimated for star-formation in the galaxy. Previously unpublished Very Long Baseline Array data reveal a 200 pc scale double radio source in the galaxy core which is more closely aligned with the current QSO axis and may represent a more recent period of jet activity. These results suggest that the realignment of the QSO, the cessation of jet activity, and the onset of rapid star-formation may have been caused by a gas-rich galaxy merger. A Chandra X-ray observation confirms the presence of cavities associated with the radio jets, and we estimate the energy required to inflate them to be ~7.7x10^60 erg. The mechanical power of the jets is sufficient to balance radiative cooling in the cluster, provided they are efficiently coupled to the intra-cluster medium (ICM). We find no evidence of direct radiative heating and conclude that the QSO either lacks the radiative luminosity to heat the ICM, or that it requires longer than 100-200 Myr to significantly impact its environment. [Abridged]

We present a study of the distribution of X-ray AGN in a representative sample of 26 massive clusters at 0.15<z<0.30, combining Chandra observations with highly complete spectroscopy of cluster members down to M_K*+2. In total we identify 48 X-ray AGN among the cluster members, with luminosities 2x10^41-1x10^44erg/s. In the stacked caustic diagram, the X-ray AGN appear to preferentially lie along the caustics, suggestive of an infalling population. They also appear to avoid the region with lowest cluster-centric radii and relative velocities (r_proj<0.4 r_500; |v-<v>|/sigma_v<0.8), which is dominated by the virialized population of galaxies accreted earliest into the clusters. Moreover the velocity dispersion of the 48 X-ray AGN is 1.51x that of the overall cluster population, which is consistent with the sqrt(2) ratio expected by simple energetic arguments when comparing infalling versus virialized populations. This kinematic segregation is significant at the 4.66-sigma level. When splitting the X-ray AGN sample into two according to X-ray or infrared (IR) luminosity, both X-ray bright and IR-bright sub-samples show higher velocity dispersions than their X-ray dim and IR-dim counterparts at >2sigma significance. This is consistent with the nuclear activity responsible for the X-ray and IR emission being slowly shut down as the host galaxies are accreted into the cluster. Overall our results provide the strongest observational evidence to date that X-ray AGN found in massive clusters are an infalling population, and that the cluster environment very effectively suppresses radiatively-efficient nuclear activity in its member galaxies. These results are consistent with the view that for galaxies to host an X-ray AGN they should be the central galaxy within their dark matter halo and have a ready supply of cold gas.

We present a Hubble Space Telescope/Wide Field Planetary Camera 2 weak-lensing study of A520, where a previous analysis of ground-based data suggested the presence of a dark mass concentration. We map the complex mass structure in much greater detail leveraging more than a factor of three increase in the number density of source galaxies available for lensing analysis. The "dark core" that is coincident with the X-ray gas peak, but not with any stellar luminosity peak is now detected with more than 10 sigma significance. The ~1.5 Mpc filamentary structure elongated in the NE-SW direction is also clearly visible. Taken at face value, the comparison among the centroids of dark matter, intracluster medium, and galaxy luminosity is at odds with what has been observed in other merging clusters with a similar geometric configuration. To date, the most remarkable counter-example might be the Bullet Cluster, which shows a distinct bow-shock feature as in A520, but no significant weak-lensing mass concentration around the X-ray gas. With the most up-to-date data, we consider several possible explanations that might lead to the detection of this peculiar feature in A520. However, we conclude that none of these scenarios can be singled out yet as the definite explanation for this puzzle.

We present the results of Suzaku X-ray follow-up observation of weak-lensing-detected halos in the field around galaxy cluster ZwCl0823.2+0425. We clearly detected X-ray emission associated with most of these halos and determined their detailed physical parameters such as X-ray luminosity, temperature, and metal abundance, for the first time. We find that the X-ray luminosity - temperature relation for these halos agrees with former typical results. With mass determined from the weak gravitational lensing data, the mass-temperature relation for them is also investigated and found to be consistent with the prediction from a simple self-similar model and results of the previous studies with both lensing and X-ray data. We would like to emphasize that the self-similar scaling relation of mass and temperature is shown here for the first time using a weak-lensing selected sample, whereas previous studies of the mass scaling relation used X-ray-selected samples of clusters. Therefore, our study demonstrates importance of X-ray follow-up observations of shear-selected clusters, and shows that a joint X-ray and lensing analysis will be crucial for clusters discovered by the forthcoming weak-lensing surveys, such as the one planned with Subaru/Hyper-Suprime-Cam.

We present an early broad-brush analysis of Herschel/PACS observations of star-forming galaxies in 8 galaxy clusters drawn from our survey of 30 clusters at z~0.2. We define a complete sample of 192 spectroscopically confirmed cluster members down to L_TIR>3x10^10L_sun and L_K>0.25L_sun. The average K-band and bolometric infrared luminosities of these galaxies both fade by a factor of ~2 from clustercentric radii of 2r_200 to ~0.5r_200, indicating that as galaxies enter the clusters ongoing star-formation stops first in the most massive galaxies, and that the specific star-formation rate (SSFR) is conserved. On smaller scales the average SSFR jumps by 25%, suggesting that in cluster cores processes including ram pressure stripping may trigger a final episode of star-formation that presumably exhausts the remaining gas. This picture is consistent with our comparison of the Herschel-detected cluster members with the cluster mass distributions, as measured in our previous weak-lensing study of these clusters. For example, the spatial distribution of the Herschel sources is positively correlated with the structures in the weak-lensing mass maps at 5sigma significance, with the strongest signal seen at intermediate group-like densities. The strong dependence of the total cluster IR luminosity on cluster mass L_TIR propto M_virial^2 is also consistent with accretion of galaxies and groups of galaxies (i.e. the substructure mass function) driving the cluster IR luminosity. The most surprising result is that roughly half of the Herschel-detected cluster members have redder S_100/S_24 flux ratios than expected, based on the Rieke et al. models. On average cluster members are redder than non-members, and the fraction of red galaxies increases towards the cluster centers, both of which indicate that these colors are not attributable to systematic photometric errors. [Abridged]

We investigate the effect of three important processes by which AGN-blown bubbles transport material: drift, wake transport and entrainment. The first of these, drift, occurs because a buoyant bubble pushes aside the adjacent material, giving rise to a net upward displacement of the fluid behind the bubble. For a spherical bubble, the mass of upwardly displaced material is roughly equal to half the mass displaced by the bubble, and should be ~ 10^{7-9} solar masses depending on the local ICM and bubble parameters. We show that in classical cool core clusters, the upward displacement by drift may be a key process in explaining the presence of filaments behind bubbles. A bubble also carries a parcel of material in a region at its rear, known as the wake. The mass of the wake is comparable to the drift mass and increases the average density of the bubble, trapping it closer to the cluster centre and reducing the amount of heating it can do during its ascent. Moreover, material dropping out of the wake will also contribute to the trailing filaments. Mass transport by the bubble wake can effectively prevent the build-up of cool material in the central galaxy, even if AGN heating does not balance ICM cooling. Finally, we consider entrainment, the process by which ambient material is incorporated into the bubble. Abridged

We compare X-ray hydrostatic and weak-lensing mass estimates for a sample of 12 clusters that have been observed with both XMM-Newton and Subaru. At an over-density of \Delta=500, we obtain 1-M^{X}/M^{WL}=0.01+/-0.07 for the whole sample. We also divided the sample into undisturbed and disturbed sub-samples based on quantitative X-ray morphologies using asymmetry and fluctuation parameters, obtaining 1-M^{X}/M^{WL}=0.09+/-0.06 and -0.06+/-0.12 for the undisturbed and disturbed clusters, respectively. In addition to non-thermal pressure support, there may be a competing effect associated with adiabatic compression and/or shock heating which leads to overestimate of X-ray hydrostatic masses for disturbed clusters, for example, in the famous merging cluster A1914. Despite the modest statistical significance of the mass discrepancy, on average, in the undisturbed clusters, we detect a clear trend of improving agreement between M^{X} and M^{WL} as a function of increasing over-density, M^{X}/M^{WL}=(0.908+/-0.004)+(0.187+/-0.010) \log_{10}(\Delta /500). We also examine the gas mass fractions, f_{gas}=M^{gas}/M^{WL}, finding that they are an increasing function of cluster radius, with no dependence on dynamical state, in agreement with predictions from numerical simulations. Overall, our results demonstrate that XMM-Newton and Subaru are a powerful combination for calibrating systematic uncertainties in cluster mass measurements.

Unbounded outflows in the form of highly collimated jets and broad winds appear to be a ubiquitous feature of accreting black hole systems. The most powerful jets are thought to derive a significant fraction, if not the majority, of their power from the rotational energy of the black hole. Whatever the precise mechanism that causes them, these jets must therefore exert a braking torque on the black hole. We calculate the spin-up function for an accreting black hole, accounting for this braking torque. We find that the predicted black hole spin-up function depends only on the black hole spin and dimensionless parameters describing the accretion flow. Using recent relativistic magnetohydrodynamical numerical simulation results to calibrate the efficiency of angular momentum transfer in the flow, we find that an ADAF flow will spin a black hole up (or down) to an equilibrium value of about 96% of the maximal spin value in the absence of jets. Combining our ADAF system with a simple model for jet power, we demonstrate that an equilibrium is reached at approximately 93% of the maximal spin value, as found in the numerical simulation studies of the spin-up of accreting black holes, at which point the spin-up of the hole by accreted material is balanced by the braking torque arising from jet production. Our model also yields a relationship between jet efficiency and black hole spin that is in surprisingly good agreement with that seen in the simulation studies, indicating that our simple model is a useful and convenient description of ADAF inflow - jet outflow about a spinning black hole for incorporation in models of the formation and evolution of galaxies, groups and clusters of galaxies.

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.

In this paper, we undertake a study to determine what insight can be reliably gleaned from the comparison of the X-ray and the weak lensing mass maps of galaxy clusters. We do this by investigating the 2D substructure within three high-resolution cosmological simulations of galaxy clusters. Our main results focus on non-radiative gas dynamics, but we also consider the effects of radiative cooling at high redshift. For our analysis, we use a novel approach, based on unsharp-masking, to identify substructures in 2D surface mass density and X-ray surface brightness maps. At full resolution (~ 15 h^-1 kpc), this technique is capable of identifying almost all self-bound dark matter subhaloes with M>10^12 h^-1 M_sun. We also report a correlation between the mass of a subhalo and the area of its corresponding 2D detection; such a correlation, once calibrated, could provide a useful estimator for substructure mass. Comparing our 2D mass and X-ray substructures, we find a surprising number of cases where the matching fails: around one third of galaxy-sized substructures have no X-ray counterpart. Some interesting cases are also found at larger masses, in particular the cores of merging clusters where the situation can be complex. Finally, we degrade our mass maps to what is currently achievable with weak-lensing observations (~100 h^-1 kpc at z=0.2). While the completeness mass limit increases by around an order of magnitude, a mass-area correlation remains. Our paper clearly demonstrates that the next generation of lensing surveys should start to reveal a wealth of information on cluster substructure. (Abridged)

Virtually all massive galaxies, including our own, host central black holes ranging in mass from millions to billions of solar masses. The growth of these black holes releases vast amounts of energy that powers quasars and other weaker active galactic nuclei. A tiny fraction of this energy, if absorbed by the host galaxy, could halt star formation by heating and ejecting ambient gas. A central question in galaxy evolution is the degree to which this process has caused the decline of star formation in large elliptical galaxies, which typically have little cold gas and few young stars, unlike spiral galaxies.

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)

We present a joint analysis of near-ultraviolet (NUV) data from the GALEX mission and (optical) colour profiles for a sample of 7 Brightest Cluster Galaxies (BCGs) in the Canadian Cluster Comparison Project. We find that every BCG which has a blue rest-frame UV colour also shows a blue-core in its optical colour profile. Conversely, BCGs that lack blue cores and show monotonic colour gradients typical of old elliptical galaxies, are red in the UV. We interpret this as evidence that the NUV enhancement in the blue BCGs is driven by recent star formation and not from old evolved stellar populations such as horizontal branch stars. The recent star formation in the blue BCGs typically has an age less than 200 Myrs and contributes mass fractions of less than a percent. Although the sample studied here is small, we demonstrate, for the first time, a one-to-one correspondence between blue cores in elliptical galaxies (in particular BCGs) and a NUV-enhancement observed using GALEX. The combination of this one-to-one correspondence and the consistently young age of recent star formation, coupled with additional correlations with the host cluster's X-ray properties, strongly suggests that the star formation is fueled by gas cooling out of the intracluster medium. In turn, this implies that any AGN heating of the intracluster medium in massive clusters only acts to reduce the magnitude of the cooling flow and that once this flow starts, it is nearly always active. (abridged)

The number density of galaxy clusters provides tight statistical constraints on the matter fluctuation power spectrum normalization, traditionally phrased in terms of sigma_8, the root mean square mass fluctuation in spheres with radius 8 h^-1 Mpc. We present constraints on sigma_8 and the total matter density Omega_m0 from local cluster counts as a function of X-ray temperature, taking care to incorporate and minimize systematic errors that plagued previous work with this method. In particular, we present new determinations of the cluster luminosity - temperature and mass - temperature relations, including their intrinsic scatter, and a determination of the Jenkins mass function parameters for the same mass definition as the mass - temperature calibration. Marginalizing over the 12 uninteresting parameters associated with this method, we find that the local cluster temperature function implies sigma_8 (Omega_m0/0.32)^alpha = 0.86+/-0.04 with alpha = 0.30 (0.41) for Omega_m0 < 0.32 (Omega_mo > 0.32) (68% confidence for two parameters). This result agrees with a wide range of recent independent determinations, and we find no evidence of any additional sources of systematic error for the X-ray cluster temperature function determination of the matter power spectrum normalization. The joint WMAP5 + cluster constraints are: Omega_m0 = 0.30+0.03/-0.02 and sigma_8 = 0.85+0.04/-0.02 (68% confidence for two parameters).

We present measurements of surface brightness and colour profiles for the brightest cluster galaxies (BCGs) in a sample of 48 X-ray luminous galaxy clusters. These data were obtained as part of the Canadian Cluster Comparison Project (CCCP). The Kormendy relation of our BCGs is steeper than that of the local ellipticals, suggesting differences in the assembly history of these types of systems. We also find that while most BCGs show monotonic colour gradients consistent with a decrease in metallicity with radius, 25% of the BCGs show colour profiles that turn bluer towards the centre (blue-cores). We interpret this bluing trend as evidence for recent star formation. The excess blue light leads to a typical offset from the red sequence of 0.5 to 1.0 mag in (g'-r'), thus affecting optical cluster studies that may reject the BCG based on colour. All of the blue-core BCGs are located within ~10 kpc of the peak in the cluster X-ray emission. Furthermore, virtually all of the BCGs with recent star formation are in clusters that lie above the Lx-Tx relation. Based on photometry alone, these findings suggest that central star formation is a ubiquitous feature of BCGs in dynamically relaxed cool-core clusters. This implies that while AGNs and other heating mechanisms are effective at tempering cooling, they do not full compensate for the energy lost via radiation.

For a decade, N-body simulations have revealed a nearly universal dark matter density profile, which appears to be robust to changes in the overall density of the universe and the underlying power spectrum. Despite its universality, the physical origin of this profile has not yet been well understood. Semi--analytic models by Barnes et al. (2005) have suggested that the density structure of dark matter halos is determined by the onset of the radial orbit instability (ROI). We have tested this hypothesis using N-body simulations of collapsing dark matter halos with a variety of initial conditions. For dynamically cold initial conditions, the resulting halo structures are triaxial in shape, due to the mild aspect of the instability. We examine how variations in initial velocity dispersion affect the onset of the instability, and find that an isotropic velocity dispersion can suppress the ROI entirely, while a purely radial dispersion does not. The quantity sigma^2/vc^2 is a criterion for instability, where regions with sigma^2/vc^2 <~1 become triaxial due to the ROI or other perturbations. We also find that the radial orbit instability sets a scale length at which the velocity dispersion changes rapidly from isotropic to radially anisotropic. This scale length is proportional to the radius at which the density profile changes shape, as is the case in the semi--analytic models; however, the coefficient of proportionality is different by a factor of ~2.5. We conclude that the radial orbit instability is likely to be a key physical mechanism responsible for the nearly universal profiles of simulated dark matter halos.

(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.

(Abridged) X-ray clusters are conventionally divided into two classes: "cool core" (CC) clusters and "non-cool core" (NCC) clusters. Yet relatively little attention has been given to the origins of this dichotomy and, in particular, to the energetics and thermal histories of the two classes. We develop a model for the entropy profiles of clusters starting from the configuration established by gravitational shock heating and radiative cooling. At large radii, gravitational heating accounts for the observed profiles and their scalings well. However, at small and intermediate radii, radiative cooling and gravitational heating cannot be combined to explain the observed profiles of either type of cluster. The inferred entropy profiles of NCC clusters require that material is preheated prior to cluster collapse in order to explain the absence of low entropy (cool) material in these systems. We show that a similar modification is also required in CC clusters in order to match their properties at intermediate radii. In CC clusters, this modification is unstable, and an additional process is required to prevent cooling below a temperature of a few keV. We show that this can be achieved by adding a self-consistent AGN feedback loop in which the lowest-entropy, most rapidly cooling material is heated so that it rises buoyantly to mix with material at larger radii. The resulting model does not require fine tuning and is in excellent agreement with a wide variety of observational data. Some of the other implications of this model are briefly discussed.

We discuss the effects of energy injection into the intra-group medium of the group of galaxies associated with NGC 741. The X-ray emission reveals a large bubble, which in the absence of a currently bright central radio source, may have been inflated by a previous cycle of nuclear activity . If the bubble is filled with a light, relativistic fluid we calculate that in expanding, it has done more than sufficient work to counteract the energy lost from the intra-group medium via radiative cooling; the bubble can provide this energy as it expands and rises. Using upper limits on the flux density of the plasma filling the bubble at 330 MHz and 1.4 GHz, we derive constraints on its electron energy distribution and magnetic field strength. We show that the data require the high-energy cut-off of the electron spectrum to be very low compared to the cut-offs seen in more typical radio sources if the fluid filling the bubble is a conventional relativistic plasma. This suggests that the fluid filling the bubble may not have evolved by expansion or synchrotron losses consistent with a dead radio source, leaving a puzzle as to what the origin of the bubble may be.

Using a uniform analysis procedure, we measure spatially resolved weak gravitational lensing and hydrostatic X-ray masses for a sample of 18 clusters of galaxies. We find a radial trend in the X-ray to lensing mass ratio: at r2500 we obtain a ratio MX/ML=1.03+/-0.07 which decreases to MX/ML=0.78+/-0.09 at r500. This difference is significant at 3 sigma once we account for correlations between the measurements. We show that correcting the lensing mass for excess correlated structure outside the virial radius slightly reduces, but does not eliminate this trend. An X-ray mass underestimate, perhaps due to nonthermal pressure support, can explain the residual trend. The trend is not correlated with the presence or absence of a cool core. We also examine the cluster gas fraction and find no correlation with ML, an important result for techniques that aim to determine cosmological parameters using the gas fraction.

Various differences in galaxy cluster properties derived from X-ray and weak lensing observations have been highlighted in the literature. One such difference is the observation of mass concentrations in lensing maps which have no X-ray counterparts (e.g. Jee, White, Ford et al. 2005). We investigate this issue by identifying substructures in maps of projected total mass (analogous to weak lensing mass reconstructions) and maps of projected X-ray surface brightness for three simulated clusters. We then compare the 2D mass substructures with both 3D subhalo data and the 2D X-ray substructures. Here we present preliminary results from the first comparison, where we have assessed the impact of projecting the data on subhalo identification.

The rich cluster Abell 520 (z=0.201) exhibits truly extreme and puzzling multi-wavelength characteristics. It may best be described as a "cosmic train wreck." It is a major merger showing abundant evidence for ram pressure stripping, with a clear offset in the gas distribution compared to the galaxies (as in the bullet cluster 1E 0657-558). However, the most striking feature is a massive dark core (721 h_70 M_sun/L_sun) in our weak lensing mass reconstruction. The core coincides with the central X-ray emission peak, but is largely devoid of galaxies. An unusually low mass to light ratio region lies 500 kpc to the east, and coincides with a shock feature visible in radio observations of the cluster. Although a displacement between the X-ray gas and the galaxy/dark matter distributions may be expected in a merger, a mass peak without galaxies cannot be easily explained within the current collisionless dark matter paradigm. Interestingly, the integrated gas mass fraction (~0.15), mass-to-light ratio (220 h_70 M_sun/L_sun), and position on the X-ray luminosity-temperature and mass-temperature relations are unremarkable. Thus gross properties and scaling relations are not always useful indicators of the dynamical state of clusters.

The power of jets from black holes are expected to depend on both the spin of the black hole and the structure of the accretion disk in the region of the last stable orbit. We investigate these dependencies using two different physical models for the jet power: the classical Blandford-Znajek (BZ) model and a hybrid model developed by Meier. In the BZ case, the jets are powered by magnetic fields directly threading the spinning black hole while in the hybrid model, the jet energy is extracted from both the accretion disk as well as the black hole via magnetic fields anchored to the accretion flow inside and outside the hole's ergosphere. The hybrid model takes advantage of the strengths of both the Blandford-Payne and BZ mechanisms, while avoiding the more controversial features of the latter. We develop these models more fully to account for general relativistic effects and to focus on advection-dominated accretion flows (ADAF) for which the jet power is expected to be a significant fraction of the accreted rest mass energy. We apply the models to elliptical galaxies, in order to see if these models can explain the observed correlation between the Bondi accretion rates and the total jet powers. For typical values of the disk viscosity parameter alpha~0.04-0.3 and mass accretion rates consistent with ADAF model expectations, we find that the observed correlation requires j>0.9; i.e., it implies that the black holes are rapidly spinning. Our results suggest that the central black holes in the cores of clusters of galaxies must be rapidly rotating in order to drive jets powerful enough to heat the intracluster medium and quench cooling flows.

We provide a new framework for the joint analysis of cluster observations (JACO) using simultaneous fits to X-ray, Sunyaev-Zel'dovich (SZ), and weak lensing data. Our method fits the mass models simultaneously to all data, provides explicit separation of the gaseous, dark, and stellar components, and--for the first time--allows joint constraints on all measurable physical parameters. JACO includes additional improvements to previous X-ray techniques, such as the treatment of the cluster termination shock and explicit inclusion of the BCG's stellar mass profile. An application of JACO to the rich galaxy cluster Abell 478 shows excellent agreement among the X-ray, lensing, and SZ data. We find that Abell 478 is consistent with a cuspy dark matter profile with inner slope n=1. Accounting for the stellar mass profile of the BCG allows us to rule out inner dark matter slopes n>1.1 at the 99% confidence level. At large radii, an 1/r^3 asymptotic slope is preferred over an 1/r^4 behavior. All single power law dark matter models are ruled out at greater than the 99% confidence level. JACO shows that self-consistent modeling of multiwavelength data can provide powerful constraints on the shape of the dark profile.

The cosmic background due to the Sunyaev-Zeldovich (SZ) effect is expected to be the largest signal at mm and cm wavelengths at a resolution of a few arcminutes. We investigate some simple statistics of SZ maps and their scaling with the normalization of the matter power spectrum, sigma_8, as well as the effects of the unknown physics of the intracluster medium on these statistics. We show that the SZ background provides a significant background for SZ cluster searches, with the onset of confusion occurring around 10^{14} h^{-1} solar masses in a cosmology-dependent way, where confusion is defined as typical errors in recovered flux larger than 20%. The confusion limit, corresponds to the mass at which there are roughly ten clusters per square degree, with this number nearly independent of cosmology and cluster gas physics. Typical errors grow quickly as lower mass objects are included in the catalog. We also point out that there is nothing in particular about the rms of the filtered map that makes it especially well-suited for capturing aspects of the SZ effect, and other indicators of the one-point SZ probability distribution function are at least as well suited for the task. For example, the full width at half maximum of the one point probability distribution has a field-to-field scatter that is about 60% that of the rms. The simplest statistics of SZ maps are largely unaffected by cluster physics such aspreheating, although the impact of preheating is clear by eye in the maps.Studies aimed at learning about the physics of the intracluster medium will apparently require more specialized statistical indicators.

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.

(Abridged) The thermal history of the intracluster medium (ICM) is complex. Heat input from cluster mergers, AGN, and galaxy winds offsets and may even halt the cooling of the ICM. Consequently, the processes that set the properties of the ICM play a key role in determining how galaxies form. In this paper we focus on the shock heating of the ICM during cluster mergers, with the eventual aim of incorporating this mechanism into semi-analytic models of galaxy formation. We use a suite of hydrodynamic simulations to track the evolution of the ICM in idealised two-body mergers. We find the heating of the ICM can be understood relatively simply by considering the evolution of the gas entropy during the mergers. We examine the processes that generate the entropy in order to understand why previous analytic shock heating models failed. We find that: (1) The energy that is thermalised in the collision greatly exceeds the kinetic energy available when the systems first touch. The smaller system penetrates deep into the potential well before it is disrupted. (2) For unequal mass mergers, most of the energy is thermalised in the more massive component. The heating of the smaller system is minor and its gas sinks to the centre of the final system. (3) The bulk of the entropy generation occurs in two distinct episodes. The first episode occurs following core collision, when a shock wave is generated that propagates outwards from the centre. This causes the combined system to expand rapidly and overshoot hydrostatic equilibrium. The second episode occurs as this material is shock heated as it re-collapses. This revised model for entropy generation significantly improves our physical understanding of cosmological gas simulations.

Recently, Allen et al. measured a tight correlation between the Bondi accretion rates and jet powers of the nuclei of nearby X-ray luminous elliptical galaxies. We employ two models of jet powering to understand the above correlation and derive constraints on the spin and accretion rate of the central black holes. The first is the Blandford-Znajek model, in which the spin energy of the hole is extracted as jet power; the second model is an hybrid version of the Blandford-Payne and Blandford-Znajek processes, in which the outflow is generated in the inner parts of the accretion disk. We assume advection-dominated accretion flows (ADAF) and account for general relativistic effects. Our modelling implies that for typical values of the disk viscosity parameter alpha~0.01-1 the tight correlation implies the narrow range of spins j~0.7-1 and accretion rates Mdot(R_ms)~(0.01-1) Mdot_Bondi. Our results provide support for the "spin paradigm" scenario and suggest that the central black holes in the cores of clusters of galaxies must be rapidly rotating in order to drive radio jets powerful enough to quench the cooling flows.

(Abridged) Existing models invoking AGN activty to resolve the cooling flow conundrum in galaxy clusters focus exclusively on the role of the central galaxy. Such models require fine-tuning of highly uncertain microscopic transport properties to distribute the thermal thermal over the entire cluster cooling core. We propose that the ICM is instead heated by multiple, spatially distributed AGNs. There is mounting observational evidence for multiple AGNs in cluster environments. Active AGNs drive bubbles into the ICM. We identify three distinct interactions between the bubble and the ICM: (1) Upon injection, the bubbles expand rapidly in situ to reach pressure equilibrium with their surroundings, generating shocks and waves whose dissipation is the principal source of ICM heating. (2) Once inflated, the bubbles rise buoyantly at rate determined by a balance with the viscous drag force, which itself results in some additional heating. (3) Rising bubbles expand and compress their surroundings. This process is adiabatic and does not contribute to any additional heating; rather, the increased ICM density due to compression enhances cooling. Our model sidesteps the ``transport'' issue by relying on the spatially distributed galaxies to heat the cluster core. We include self regulation in our model by linking AGN activity in a galaxy to cooling characteristics of the surrounding ICM. We use a spherically symmetric one-dimensional hydrodynamical code to carry out a preliminary study illustrating the efficacy of the model. Our self-regulating scenario predicts that there should be enhanced AGN activity of galaxies inside the cooling regions compared to galaxies in the outer parts of the cluster. This prediction remains to be confirmed or refuted by observations.

There is no accepted mechanism that explains the equilibrium structures that form in collisionless cosmological N-body simulations. Recent work has identified nonextensive thermodynamics as an innovative approach to the problem. The distribution function that results from adopting this framework has the same form as for polytropes, but the polytropic index is now related to the degree of nonextensiveness. In particular, the nonextensive approach can mimic the equilibrium structure of dark matter density profiles found in simulations. We extend the investigation of this approach to the velocity structures expected from nonextensive thermodynamics. We find that the nonextensive and simulated N-body rms-velocity distributions do not match one another. The nonextensive rms-velocity profile is either monotonically decreasing or displays little radial variation, each of which disagrees with the rms-velocity distributions seen in simulations. We conclude that the currently discussed nonextensive models require further modifications in order to corroborate dark matter halo simulations. (adapted from TeX)

It has long been realized that dark matter halos formed in cosmological N-body simulations are characterized by density profiles rho(r) that, when suitably scaled, have similar shapes. Additionally, combining the density and velocity dispersion profiles sigma(r), each of which have decidedly nonpower-law shapes, leads to quantity rho/(sigma^3) that is a power-law in radius over 3 orders of magnitude in radius. Halos' velocity anisotropy profiles beta(r) vary from isotropic near the centers of halos to quite radially anisotropic near the virial radius. Finally, there appears to be a nearly linear correlation between beta and the logarithmic density slope gamma for a wide variety of halos. This work is part of a continuing investigation of the above interrelationships and their origins using analytical and semi-analytical techniques. Our finding suggest that the nearly linear beta--gamma relationship is not just another expression of scale-free rho/(sigma^3) behavior. We also note that simultaneously reproducing density and anisotropy profiles like those found in simulations requires beta(r) and gamma(r) to have similar shapes, leading to nearly linear beta--gamma correlations. This work suggests that the beta--gamma and power-law rho/(sigma^3) relations have distinct physical origins.

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.

We explore the possibility that extended disks, such as that recently discovered in M31, are the result of a single dwarf (10^9--10^10 Msun) satellite merger. We carry out N-body simulations of dwarf NFW halos with embedded spheriodal stellar components on co-planar, prograde orbits in a M31-like host galaxy. As the orbit decays due to dynamical friction and the system is disrupted, the stellar particles form an extended, exponential disk-like structure that spans the radial range 30--200 kpc once they have relaxed. The disk scale-length correlates with the initial extent of the stellar component within the satellite halo: the more embedded the stars, the smaller the resulting disk scale-length. If the progenitors start on circular orbits, the kinematics of the stars that make up the extended disk have an average rotational motion that is 30--50 km/s lower than the host's circular velocity. For dwarf galaxies moving on highly eccentric orbits (e=0.7), the stellar debris exhibits a much lower rotational velocity. Our results imply that extended galactic disks might be a generic feature of the hierarchical formation of spiral galaxies such as M31 and the Milky Way.

We investigate the connection between collisionless equilibria and the phase-space relation between density $\rho$ and velocity dispersion $\sigma$ found in simulations of dark matter halo formation, $F=\psd \propto r^{-\alpha}$. Understanding this relation will shed light on the physics relevant to collisionless collapse and on the subsequent structures formed. We show that empirical density profiles that provide good fits to N-body halos also happen to have nearly scale-free \psd distributions when in equilibrium. We have also done a preliminary investigation of variables other than $r$ that may match or supercede the correlation with $F$. In the same vein, we show that $\rho/\sigma^m$, where $m=3$ is the most appropriate combination to use in discussions of the power-law relationship. Since the mechanical equilibrium condition that characterizes the final systems does not by itself lead to power-law $F$ distributions, our findings prompt us to posit that dynamical collapse processes (such as violent relaxation) are responsible for the radial power-law nature of the \psd distributions of virialized systems.

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.

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.

(Abridged) We revisit the scaling relationships between the dark matter mass and observed X-ray luminosity and temperature of galaxy clusters and groups in the local Universe. Specifically, we compare recent observations with analytic models of the intracluster medium in which the gas entropy distribution has been shifted by a variable amount, K0, to investigate the origin of the scatter in these scaling relations, and its influence on the luminosity and temperature functions. We find that variations in halo concentration or formation epoch (which might determine the time available for low entropy gas to cool out) are insufficient to explain the amount of scatter in the mass-luminosity relation. Instead, a range of entropy floors at a fixed halo mass, spanning approximately 50 keV sq.cm to 700 keV sq.cm, is required to match the data. This range is likely related to the variance in heating and/or cooling efficiency from halo to halo. We demonstrate that these models are consistent with the observed temperature and luminosity functions of clusters, if scatter in the mass-luminosity relation is accounted for. Finally, we present predictions for the redshift evolution of these scaling relations and luminosity/temperature functions. Comparison with recent data at z<0.7 shows reasonable agreement with a model that assumes a median entropy floor of K0=200 keV sq.cm. When observations are extended to group scales (kT < 1 keV), this evolution will have the potential to discriminate between an entropy floor that is independent of redshift (for example, in a preheating scenario) and one that depends on the cooling time of the halo.

Starting with the well-known NFW dark matter halo distribution, we construct a simple polytropic model for the intracluster medium which is in good agreement with high resolution numerical hydrodynamical simulations, apply this model to a very large scale concordance dark matter simulation, and compare the resulting global properties with recent observations of X-ray clusters, including the mass-temperature and luminosity-temperature relations. We make allowances for a non-negligible surface pressure, removal of low entropy (short cooling time) gas, energy injection due to feedback, and for a relativistic (non-thermal) pressure component. A polytropic index n=5 (Gamma=1.2) provides a good approximation to the internal gas structure of massive clusters (except in the very central regions where cooling becomes important), and allows one to recover the observed M_500-T, L_x-T, and T/n_e^{2/3} ~ T^0.65 relations. Using these concepts and generalizing this method so that it can be applied to fully three-dimensional N-body simulations, one can predict the global X-ray and SZE trends for any specified cosmological model. We find a good fit to observations when assuming that twelve percent of the initial baryonic mass condenses into stars, the fraction of rest mass of this condensed component transferred back to the remaining gas (feedback) is 3.9E-5, and the fraction of total pressure from a nonthermal component is near ten percent.

Although N-body studies of dark matter halos show that the density profiles, rho(r), are not simple power-laws, the quantity rho/sigma^3, where sigma(r) is the velocity dispersion, is in fact a featureless power-law over ~3 decades in radius. In the first part of the paper we demonstrate, using the semi-analytic Extended Secondary Infall Model (ESIM), that the nearly scale-free nature of rho/sigma^3 is a robust feature of virialized halos in equilibrium. By examining the processes in common between numerical N-body and semi-analytic approaches, we argue that the scale-free nature of rho/sigma^3 cannot be the result of hierarchical merging, rather it must be an outcome of violent relaxation. The empirical results of the first part of the paper motivate the analytical work of the second part of the paper, where we use rho/sigma^3 proportional to r^{-alpha} as an additional constraint in the isotropic Jeans equation of hydrostatic equilibrium. Our analysis shows that the constrained Jeans equation has different types of solutions, and in particular, it admits a unique ``periodic'' solution with alpha=1.9444. We derive the analytic expression for this density profile, which asymptotes to inner and outer profiles of rho ~ r^{-0.78}, and rho ~ r^{-3.44}, respectively.

We investigate a hypothesis regarding the origin of the scalelength in halos formed in cosmological N-body simulations. This hypothesis can be viewed as an extension of an earlier idea put forth by Merritt and Aguilar. Our findings suggest that a phenomenon related to the radial orbit instability is present in such halos and is responsible for density profile shapes. This instability sets a scalelength at which the velocity dispersion distribution changes rapidly from isotropic to radially anisotropic. This scalelength is reflected in the density distribution as the radius at which the density profile changes slope. We have tested the idea that radially dependent velocity dispersion anisotropy leads to a break in density profile shape by manipulating the input of a semi-analytic model to imitate the velocity structure imposed by the radial orbit instability. Without such manipulation, halos formed are approximated by single power-law density profiles and isotropic velocity distributions. Halos formed with altered inputs display density distributions featuring scalelengths and anisotropy profiles similar to those seen in cosmological N-body simulations.

High resolution N-body simulations have all but converged on a common empirical form for the shape of the density profiles of halos, but the full understanding of the underlying physics of halo formation has eluded them so far. We investigate the formation and structure of dark matter halos using analytical and semi-analytical techniques. Our halos are formed via an extended secondary infall model (ESIM); they contain secondary perturbations and hence random tangential and radial motions which affect the halo's evolution at it undergoes shell-crossing and virialization. Even though the density profiles of NFW and ESIM halos are different their phase-space density distributions are the same: \rho/\sigma^3 ~ r^{-\alpha}, with \alpha=1.875 over ~3 decades in radius. We use two approaches to try to explain this ``universal'' slope: (1) The Jeans equation analysis yields many insights, however, does not answer why \alpha=1.875. (2) The secondary infall model of the 1960's and 1970's, augmented by ``thermal motions'' of particles does predict that halos should have \alpha=1.875. However, this relies on assumptions of spherical symmetry and slow accretion. While for ESIM halos these assumptions are justified, they most certainly break down for simulated halos which forms hierarchically. We speculate that our argument may apply to an ``on-average'' formation scenario of halos within merger-driven numerical simulations, and thereby explain why \alpha=1.875 for NFW halos. Thus, \rho/\sigma^3 ~ r^{-1.875} may be a generic feature of violent relaxation.

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.

In a previous paper, we described a new method for including detailed information about substructure in semi-analytic models of dark matter halo formation based on merger trees. In this paper, we present the basic predictions of our full model of halo formation. We first describe the overall properties of substructure in galaxy, group or cluster haloes at the present day. We then discuss the evolution of substructure, and the effect of the mass accretion history of an individual halo on the mass function and orbital grouping of its subhalo population. We show, in particular, that the shape of the subhalo mass function is strongly correlated with the formation epoch of the halo. In a third paper in this series we will compare the results of our semi-analytic method with the results of self-consistent numerical simulations of halo formation.

In a previous paper, we described a new method for including detailed information about substructure in semi-analytic models of halo formation based on merger trees. In this paper, we compare the predictions of our model with results from self-consistent numerical simulations. We find that in general the two methods agree extremely well, particularly once numerical effects and selection effects in the choice of haloes are taken into account. As expected from the original analyses of the simulations, we see some evidence for artificial overmerging in the innermost regions of the simulated haloes, either because substructure is being disrupted artificially or because the group-finding algorithms used to identify substructure are not detecting all the bound clumps in the highest-density regions. Our analytic results suggest that greater mass and force resolution may be required before numerical overmerging becomes negligible in all current applications. We discuss the implications of this result for observational and experimental tests of halo substructure, such as the analysis of discrepant magnification ratios in strongly lensed systems, terrestrial experiments to detect dark matter particles directly, or indirect detection experiments searching for positrons, gamma-rays, neutrinos or other dark matter decay products.

Observations of multiple-image gravitational lens systems suggest that the projected mass distributions of galaxy haloes may contain substantial inhomogeneities. The fraction of the halo mass in dense substructure is still highly uncertain, but could be as large as a few percent. While halo substructure is seen in numerical simulations of CDM haloes, little of this substructure survives in the innermost regions of haloes, and thus the observational claims for substructure at small projected radii are slightly surprising. There is evidence, however, that even the highest-resolution simulations published to-date are still limited by numerical effects that heat and disrupt substructure artificially in high-density regions. By comparing numerical and semi-analytic models of halo substructure, we show that current simulations probably underestimate the mass fraction in substructure at small projected radii, by a factor of at least 2-3. We discuss the prospects for using lensing observations as a fundamental test of the nature of dark matter.

(Abridged) Theoretical models that include only gravitationally-driven processes fail to match the observed mean X-ray properties of clusters. As a result, there has recently been increased interest in models in which either radiative cooling or entropy injection play a central role in mediating the properties of the intracluster medium. Both sets of models give reasonable fits to the mean properties of clusters, but cooling only models result in fractions of cold baryons in excess of observationally established limits and the simplest entropy injection models do not treat the "cooling core" structure present in many clusters and cannot account for entropy profiles revealed by recent X-ray observations. We consider models that marry radiative cooling with entropy injection, and confront model predictions for the global and structural properties of massive clusters with the latest X-ray data. The models successfully and simultaneously reproduce the observed L-T and L-M relations, yield detailed entropy, surface brightness, and temperature profiles in excellent agreement with observations, and predict a cooled gas fraction that is consistent with observational constraints. The model also provides a possible explanation for the significant intrinsic scatter present in the L-T and L-M relations and provides a natural way of distinguishing between clusters classically identified as "cooling flow" clusters and dynamically relaxed "non-cooling flow" clusters. The former correspond to systems that had only mild levels (< 300 keV cm^2) of entropy injection, while the latter are identified as systems that had much higher entropy injection. This is borne out by the entropy profiles derived from Chandra and XMM-Newton.

Models of structure formation based on cold dark matter predict that most of the small dark matter haloes that first formed at high redshift would have merged into larger systems by the present epoch. Substructure in present-day haloes preserves the remains of these ancient systems, providing the only direct information we may ever have about the low-mass end of the power spectrum. I describe some recent attempts to model halo substructure down to very small masses, using a semi-analytic model of halo formation. I make a preliminary comparison between the model predictions, observations of substructure in lensed systems, and the properties of local satellite galaxies.

We describe our recent attempts to model substructure in dark matter halos down to very small masses, using a semi-analytic model of halo formation. The results suggest that numerical simulations of halo formation may still be missing substructure in the central regions of halos due to purely numerical effects. If confirmed, this central 'overmerging' problem will have important consequences for the interpretation of lensing measurements of substructure. We also show that the spatial distribution of subhalos relative to the satellite companions of the Milky Way rules out at least one simple model of how dwarf galaxies form in low-mass halos.

Although high-resolution N-body simulations make robust empirical predictions for the density distribution within cold dark matter halos, these studies have yielded little physical insight into the origins of the distribution. We investigate the problem using analytic and semi-analytic approaches. Simple analytic considerations suggest that the inner slope of dark matter halos cannot be steeper than alpha=2 (rho ~ r^-alpha), with alpha=1.5-1.7 being a more realistic upper limit. Our analysis suggests that any number of effects, both real (angular momentum from tidal torques, secondary perturbations) and artificial (two-body interactions, the accuracy of the numerical integrator, round-off errors), will result in shallower slopes. We also find that the halos should exhibit a well-defined relation between r_peri/r_apo and j_theta/j_r. We derive this relation analytically and speculate that it may be "universal". Using a semi-analytic scheme based on Ryden & Gunn (1987), we further explore the relationship between the specific angular momentum distribution in a halo and its density profile. For now we restrict ourselves to halos that form primarily via nearly-smooth accretion of matter, and only consider the specific angular momentum generated by secondary perturbations associated with the cold dark matter spectrum of density fluctuations. Compared to those formed in N-body simulations, our ``semi-analytic'' halos are more extended, have flatter rotation curves and have higher specific angular momentum, even though we have not yet taken into account the effects of tidal torques. Whether the density profiles of numerical halos is indeed the result of loss in angular momentum outside the central region, and whether this loss is a feature of hierarchical merging and major mergers in particular, is under investigation.

Recent X-ray observations have highlighted clusters that lack entropy cores. At first glance, these results appear to invalidate the preheated ICM models. We show that a self-consistent preheating model, which factors in the effects of radiative cooling, is in excellent agreement with the observations. Moreover, the model naturally explains the intrinsic scatter in the L-T relation, with ``cooling flow'' and ``non-cooling flow'' systems corresponding to mildly and strongly preheated systems, respectively. We discuss why preheating ought to be favoured over merging as a mechanism for the origin of ``non-cooling flow'' clusters.

We describe the properties of a sample of galaxy groups with very unusual distributions of galaxy luminosities. The most extreme example has an X-ray luminosity similar to that of the Virgo cluster but has a very low richness, with only one galaxy brighter than L*, compared with six in Virgo. That one galaxy, however, is optically more luminous than any galaxy in Virgo and has an optical luminosity as bright as many of the central cD galaxies in rich Abell clusters. The characteristic feature of the fossil groups we study is that most of the light arises from one dominant, central galaxy. We define a fossil system and, based on this definition, construct a small X-ray selected, flux-limited sample of fossil groups with well known selection criteria. We confirm that these systems are indeed groups of galaxies, but dominated by one central luminous giant elliptical galaxy and with few, or no, L* galaxies. We find that fossil systems represent 8%-20% of all systems of the same X-ray luminosity. Fossil groups are at least as numerous as all poor and rich clusters combined, and are thus a possible site for the formation of luminous central cluster galaxies before infall into clusters occurs. The fossil systems in our sample have significantly higher X-ray luminosities than normal groups of similar total optical luminosities (or similar X-ray temperature, where the latter can be measured). These enhanced X-ray luminosities may be due to relatively cool gas in the innermost regions or due to a low central gas entropy. We interpret fossil groups as old, undisturbed systems which have avoided infall into clusters, but where galaxy merging of most of the L* galaxies has occurred. An early formation epoch, before that of most groups, could explain low central gas entropies and high X-ray luminosities.

Studies of cluster X-ray scaling relations have led to suggestions that non-gravitational processes, e.g., radiative cooling and/or "preheating", have significantly modified the entropy of the intracluster medium (ICM). For the first time, we test this hypothesis through a comparison of predicted thermal Sunyaev-Zeldovich (SZ) effect scaling relations with available data from the literature. One of the relations that we explore, in principle, depends solely on SZ effect observations, thus offering an X-ray independent probe of the ICM. A detailed comparison of the theoretical relations with the largest compilation of high redshift SZ effect data to date indicates that the presence of an entropy floor is favored by the data. Furthermore, the inferred level of that floor, 300 keV cm^2 or higher, is comparable to that found in studies of X-ray scaling relations of nearby massive clusters. Thus, we find no evidence for significant evolution of the entropy floor out to z = 0.7. We further demonstrate that the high quality data to be obtained from the upcoming Sunyaev-Zeldovich Array (SZA) and the (soon-to-be) upgraded Owens Valley Radio Observatory (OVRO) array will open powerful new windows into the properties of the ICM. Specifically, the new measurements will allow for accurate measurements of the ICM entropy for even the most distant galaxy clusters.

A common feature of the X-ray bubbles observed in Chandra images of some "cooling flow" clusters is that they appear to be surrounded by bright, cool shells. Temperature maps of a few nearby luminous clusters reveal that the shells consist of the coolest gas in the clusters -- much cooler than the surrounding medium. Using simple models, we study the effects of this cool emission on the inferred cooling flow properties of clusters. We find that the introduction of bubbles into model clusters that do NOT have cooling flows results in temperature and surface brightness profiles that resemble those seen in nearby "cooling flow" clusters. They also approximately reproduce the recent XMM-Newton and Chandra observations of a high minimum temperature of 1-3 keV. Hence, bubbles, if present, must be taken into account when inferring the physical properties of the ICM. In the case of some clusters, bubbles may account entirely for these observed features, calling into question their designation as clusters with cooling flows. However, since not all nearby "cooling flow" clusters show bubble-like features, we suggest that there may be a diverse range of physical phenomena that give rise to the same observed features.

X-ray observations of an "entropy floor" in nearby groups and clusters of galaxies offer evidence that important non-gravitational processes, such as radiative cooling and/or "preheating", have strongly influenced the evolution of the intracluster medium (ICM). We examine how the presence of an entropy floor modifies the thermal Sunyaev-Zeldovich (SZ) effect. A detailed analysis of scaling relations between X-ray and SZ effect observables and also between the two primary SZ effect observables is presented. We find that relationships between the central Compton parameter and the temperature or mass of a cluster are extremely sensitive to the presence of an entropy floor. The same is true for correlations between the integrated Compton parameter and the X-ray luminosity or the central Compton parameter. In fact, if the entropy floor is as high as inferred in recent analyses of X-ray data, a comparison of these correlations with both current and future SZ effect observations should show a clear signature of this excess entropy. Moreover, because the SZ effect is redshift-independent, the relations can potentially be used to track the evolution of the cluster gas and possibly discriminate between the possible sources of the excess entropy. To facilitate comparisons with observations, we provide analytic fits to these scaling relations.

The hierarchical mergers that form the haloes of dark matter surrounding galaxies, groups and clusters are not entirely efficient, leaving substantial amounts of dense substructure, in the form of stripped halo cores or `subhaloes', orbiting within these systems. Using a semi-analytic model of satellite dynamics, we study the evolution of haloes as they merge hierarchically, to determine how much substructure survives merging and how the properties of individual subhaloes change over time. We find that subhaloes evolve, due to mass loss, orbital decay, and tidal disruption, on a characteristic time-scale equal to the period of radial oscillations at the virial radius of the system. Subhaloes with realistic densities and density profiles lose 25-45 per cent of their mass per pericentric passage, depending on their concentration and on the circularity of their orbit. As the halo grows, the subhalo orbits also grow in size and become less bound. Based on these general patterns, we suggest a method for including realistic amounts of substructure in semi-analytic models based on merger trees. We show that the parameters in the resulting model can be fixed by requiring self-consistency between different levels of the merger hierarchy. In a companion paper, we will compare the results of our model with numerical simulations of halo formation.

We use a semi-analytic model of halo formation to study the dynamical history of giant field galaxies like the Milky Way. We find that in a concordance LCDM cosmology, most isolated disk galaxies have remained undisturbed for 8-10 Gyr, such that the age of the Milky Way's thin disk is unremarkable. Many systems also have older disk components which have been thickened by minor mergers, consistent with recent observations of nearby field galaxies. We do have a considerable problem, however, reproducing the morphological mix of nearby galaxies. In our fiducial model, most systems have disk-to-bulge mass ratios of order 1, and look like S0s rather than spirals. This result depends mainly on merger statistics, and is unchanged for most reasonable choices of our model parameters. We discuss two possible solutions to this morphology problem in LCDM.

Recent efforts to account for the observed X-ray luminosity - temperature relation of galaxy clusters has led to suggestions that the ICM has an apparent ``entropy floor'' at or above the level of 300 keV cm^2. Here, we propose new tests based on the thermal Sunyaev-Zeldovich effect and on the cluster gas mass - temperature trend (from X-ray data) to probe the level of excess entropy in the ICM. We show that these new tests lend further support to the case for a high entropy floor in massive clusters.

Recent X-ray observations have been used to demonstrate that the cluster gas mass-temperature relation is steeper than theoretical self-similar predictions drawn from numerical simulations that consider the evolution of the cluster gas through the effects of gravity and shock heating alone. One possible explanation for this is that the gas mass fraction is not constant across clusters of different temperature, as is usually assumed. Observationally, however, there is no compelling evidence for gas mass fraction variation, especially in the case of hot clusters. Seeking an alternative physical explanation for the observed trends, we investigate the role in the cluster gas mass - temperature relation of the preheating of the intracluster medium by some arbitrary source for clusters with emission-weighted mean temperatures of greater than about 3 keV. Making use of the physically-motivated, analytic models developed in 2002 by Babul and coworkers, we find that preheating does, indeed, lead to a steeper relation. This is in agreement with previous theoretical studies on the relation. However, in apparent conflict with these studies, we argue that a ``high'' level of entropy injection is required to match observations. In particular, an entropy floor of greater than about 300 keV cm^2 is required. We also present a new test, namely, the study of the relation within different fixed radii. This allows one to indirectly probe the density profiles of clusters, since it samples different fractions of the virial radius for clusters of different temperature. This test also confirms that a high level of preheating is required to match observations.

We have developed a new semi-analytic model for the formation and evolution of structure on galaxy, group and cluster scales. The model combines merger trees with a detailed, spatially resolved description of the dynamical evolution of halo substructure. It reproduces the results of numerical simulations of mergers and halo formation at a fraction of the computational expense. We use this model to study the formation of the dwarf galaxies in the Local Group and the structural components of the Milky Way. We find that reionization can explain the scarcity of Local Group dwarfs, although the reionization epoch is constrained to be high -- z_ri >= 12. We also find that disk disruption at recent times is rare, such that many galactic disks have an old, thin component.

Within the standard framework of structure formation, where clusters and groups of galaxies are built up from the merging of smaller systems, the physical properties of the intracluster medium, such as the gas temperature and the total X-ray luminosity, are predicted to possess well defined self-similar scaling relations. Observed clusters and groups, however, show strong deviations from these predicted relations. We argue that these deviations are unlikely to be entirely due to observational biasses; we assume they are physically based, due to the presence of excess entropy in the intracluster medium in addition to that generated by accretion shocks during the formation of the cluster. Several mechanisms have been suggested as a means of generating this entropy. Focussing on those mechanisms that preheat the gas before it becomes a constituent of the virialized cluster environment, we present a simple, intuitive, physically motivated, analytic model that successfully captures the important physics associated with the accretion of high entropy gas onto group and cluster-scale systems. We use the model to derive the new relationships between the observable properties of clusters and groups of galaxies, as well as the evolution of these relations. These include the luminosity-temperature and luminosity-$\sigma$ relations, as well as the temperature distribution function and X-ray luminosity function. These properties are found to be a more accurate description of the observations than those predicted from the standard framework. Future observations will that further test the efficacy of the preheated gas scenario are also discussed.

We have developed a simple yet surprisingly accurate analytic scheme for tracking the dynamical evolution of substructure within larger dark halos. The scheme incorporates the effects of dynamical friction, tidal mass loss and tidal heating via physically motivated approximations. Using our scheme, we can predict the orbital evolution and mass-loss history of individual subhalos in detail. We are also able to determine the impact and importance of the different physical processes on the dynamical evolution of the subhalos. To test and calibrate this model, we compare it with a set of recent high-resolution numerical simulations of mergers between galaxies and small companions. We find that we can reproduce the orbits and mass-loss rates seen in all of these simulations with considerable accuracy, using a single set of values for the three free parameters in our model. Computationally, our scheme is more than 1000 times faster than the simplest of the high-resolution numerical simulations. This means that we can carry out detailed and statistically meaningful investigations into the characteristics of the subhalo population in different cosmologies, the stripping and disruption of the subhalos, and the interactions of the subhalos with other dynamical structures such as a thin disk. This last point is of particular interest given the ubiquity of minor mergers in hierarchical models. In this regard, our method's simplicity and speed makes it particularly attractive for incorporation into semi-analytic models of galaxy formation.

We present the analysis of a suite of simulations of a Virgo mass galaxy cluster. Undertaken within the framework of standard cold dark matter cosmology, these simulations were performed at differing resolutions and with increasingly complex physical processes, with the goal of identifying the effects of each on the evolution of the cluster. We focus on the cluster at the present epoch and examine properties including the radial distributions of density, temperature, entropy and velocity. We also map `observable' projected properties such as the surface mass density, X-ray surface brightness and SZ signature. We identify significant differences between the simulations, which highlights the need for caution when comparing numerical simulations to observations of galaxy clusters. While resolution affects the inner density profile in dark matter simulations, the addition of a gaseous component, especially one that cools and forms stars, affects the entire cluster. We conclude that both resolution and included physical processes play an important role in simulating the formation and evolution of galaxy clusters. Therefore, physical inferences drawn from simulations that do not include a gaseous component that can cool and form stars present a poor representation of reality. (Abridged)

We have carried out 2D hydrodynamical simulations in order to study the interaction between supernova-powered gas outflows from low-mass galaxies and the local intergalactic medium (IGM). We are specifically interested in investigating whether a high pressure IGM, such as that in clusters of galaxies, can prevent the gas from escaping from the galaxy. The interface between the outflow and ambient IGM is demarcated by a dense expanding shell formed by the gas swept-up by the outflow. A sufficiently high IGM pressure can bring the shell to a halt well before it escapes the galaxy. Galaxies in such high pressure environments are, however, to be ploughing through the IGM at relatively high velocities. Hence, they will also be subject to ram pressure, which acts to strip the gas from the galaxy. We have carried out simulations that take into account the combined impact of ram pressure and thermal pressure. We find that ram pressure deforms the shell into a tail-like structure, fragments it into dense clouds and eventually drags the clouds away from the galaxy. The clouds are potential sites of star formation and if viewed during this transient phase, the galaxy will appear to have a low-surface brightness tail much like the galaxies with diffuse comet-like tail seen in z=1.15 cluster 3C324. In contrast, the relatively unhindered outflows in low density, low temperature environments can drive the shells of swept-up gas out to large distances from the galaxy. Such shells, if they intersect a quasar line-of-sight, would give rise to Ly $\alpha$ absorption lines of the kind seen in quasar spectra.

We confirm that the standard assumption of isothermal, shock-heated gas in cluster potentials is unable to reproduce the observed X-ray luminosity- temperature relation of groups of galaxies. As an alternative, we construct a physically motivated model for the adiabatic collapse of pre-heated gas into an isothermal potential that improves upon the original work of Kaiser (1991). The luminosity and temperature of the gas is calculated, assuming an appropriate distribution of halo formation times and radiation due to both bremsstrahlung and recombination processes. This model successfully reproduces the slope and dispersion of the luminosity-temperature relation of galaxy groups. We also present calculations of the temperature and luminosity functions for galaxy groups under the prescription of this model. This model makes two strong predictions for haloes with total masses M<10^13 M_sun, which are not yet testable with current data: (1) the gas mass fraction will increase in direct proportion to the halo mass; (2) the gas temperature will be larger than the virial temperature of the mass. The second effect is strong enough that group masses determined from gas temperatures will be overestimated by about an order of magnitude if it is assumed that the gas temperature is the virial temperature. The entropy required to match observations can be obtained by heating the gas at the turnaround time, for example, to about 3 X 10^6 K at z=1, which is too high to be generated by a normal rate of supernova explosions. This model breaks down on the scale of low mass clusters, but this is an acceptable limitation, as we expect accretion shocks to contribute significantly to the entropy of the gas in such objects.

We analyse parallel N-body simulations of three Cold Dark Matter (CDM) universes to study the abundance and clustering of galaxy clusters. The simulations cover a volume comparable to the forthcoming SDSS. We are able to make robust measurements of cluster properties to a redshift larger than unity. We extract halos using two independent, public domain group finders (FOF & HOP) and find consistent results. The correlation function of clusters is in very good agreement with a simple analytic prescription based upon a Lagrangian biasing scheme developed by Mo & White (1996) and the Press-Schechter (PS) formalism for the mass function. The R_0--D_c relation for the open CDM model is in good agreement with the results from the APM Cluster Survey. The SCDM universe shows a robust deviation in the shape and evolution of the mass function when compared with that predicted by the PS formalism. Critical models with a low sigma_8 normalization or small shape parameter Gamma show an excess of massive clusters compared with the PS prediction. When cluster normalized, the SCDM universe at z =1 contains 10 times more clusters with temperatures greater than 7keV, compared with the PS prediction. The agreement between the analytic and N-body mass functions of SCDM can be improved if the value of the delta_c (the extrapolated linear theory threshold for collapse) is revised to be $ delta_c(z) = 1.685[(0.7/sigma_8)(1+z)]^{-0.125}. Our best estimate for the amplitude of fluctuations inferred from the local cluster abundance for SCDM is sigma_{8} = 0.5 \pm 0.04. However, the discrepancy between the temperature function predicted in a critical density universe and that observed at z=0.33 (Henry et al. 1998) remains. (abridged)

The coherent image distortions induced by weak gravitational lensing can be used to measure the power spectrum of density inhomogeneities in the universe. We present our on-going effort to detect this effect with the FIRST radio survey, which currently contains about 400,000 sources over 4,200 square degrees, and thus provides a unique resource for this purpose. We discuss the sensitivity of our measurement in the context of various cosmological models. We then discuss the crucial issue of systematic effects, the most serious of which are source fragmentation, image-noise correlation, and VLA-beam anisotropy. After accounting for these effects, we expect our experiment to yield a detection, or at least a tight upper limit, for the weak lensing power spectrum on 0.2-20 degree scales.

Several independent lines of reasoning suggest that many of the very faint (B > 24) blue galaxies are low-mass objects that experienced a short epoch of star formation at redshifts 0.5 < z < 1 and have since faded into low luminosity, low surface brightness objects. Such a scenario, which arises naturally if star formation in dwarf galaxies is delayed by photoionisation due to the metagalactic UV radiation field, provides an attractive way to reconcile the Einstein-de Sitter cosmological model to the steeply rising galaxy counts observed at blue wavelengths. The Hubble Deep Field (HDF) images provide a stringent test of this model. We compare the Babul & Ferguson (1996) model to the data by constructing using simulated images and carrying out source detection and photometry for the simulations in the same way they were carried out for the real data. We compare the model predictions for the counts, sizes, and colours of galaxies observed in the HDF, and to the predictions from a low q_0 pure-luminosity-evolution (PLE) model. Both models fail to reproduce the observations. The low q_0 model predicts far more Lyman-break ``dropouts'' than are seen in the data. The fading dwarf model predicts too many remnants: faded dwarf galaxies in the redshift range 0.2 < z < 0.5 that should be detectable in the HDF as low-surface brightness red objects but are not seen. If fading dwarf galaxies are to reconcile the Einstein-de Sitter geometry to the counts, then the dwarf population must (a) form earlier than z ~ 1, with a higher initial luminosity; (b) have an initial-mass function more heavily weighted toward massive stars than the Salpeter IMF; or (c) expand much more than assumed during the supernova wind phase.

Weak lensing by large-scale mass inhomogeneities in the Universe induces correlations in the observed ellipticities of distant sources. We first review the harmonic analysis and statistics required of these correlations and discuss calculations for the predicted signal. We consider the ellipticity correlation function, the mean-square ellipticity, the ellipticity power spectrum, and a global maximum-likelihood analysis to isolate a weak-lensing signal from the data. Estimates for the sensitivity of a survey of a given area, surface density, and mean intrinsic source ellipticity are presented. We then apply our results to the FIRST radio-source survey. We predict an rms ellipticity of roughly 0.011 in $1^\circ \times 1^\circ$ pixels and 0.018 in $20' \times 20'$ pixels if the power spectrum is normalized to $\sigma_8 \Omega^{0.53} = 0.6$, as indicated by the cluster abundance. The signal is significantly larger in some models if the power spectrum is normalized instead to the COBE anisotropy. The uncertainty in the predictions from imprecise knowledge of the FIRST redshift distribution is about 25% in the rms ellipticity. We show that FIRST should be able to make a statistically significant detection of a weak-lensing signal for cluster-abundance-normalized power spectra.

One of the largest uncertainties in understanding the effect of a background UV field on galaxy formation is the intensity and evolution of the radiation field with redshift. This work attempts to shed light on this issue by computing the quasi-hydrostatic equilibrium states of gas in spherically symmetric dark matter halos (roughly corresponding to dwarf galaxies) as a function of the amplitude of the background UV field. We integrate the full equations of radiative transfer, heating, cooling and non-equilibrium chemistry for nine species: H, H^+, H^-,H_2, H_2^+, He, He^+, He^{++}, and e^-. As the amplitude of the UV background is decreased the gas in the core of the dwarf goes through three stages characterized by the predominance of ionized (H^+), neutral (H) and molecular (H_2) hydrogen. Characterizing the gas state of a dwarf galaxy with the radiation field allows us to estimate its behavior for a variety of models of the background UV flux. Our results indicate that a typical radiation field can easily delay the collapse of gas in halos corresponding to 1-$\sigma$ CDM perturbations with circular velocities less than 30 km/s.

(Abridged) In this paper we refine the theory of microlensing for a planar distribution of point masses. We derive the macroimage magnification distribution P(A) at high magnification (A-1 >> tau^2) for a low optical depth (tau << 1) lens distribution by modeling the illumination pattern as a superposition of the patterns due to individual ``point mass plus weak shear'' lenses. We show that a point mass plus weak shear lens produces an astroid- shaped caustic and that the magnification cross-section obeys a simple scaling property. By convolving this cross-section with the shear distribution, we obtain a caustic-induced feature in P(A) which also exhibits a simple scaling property. This feature results in a 20% enhancement in P(A) at A approx 2/tau. In the low magnification (A-1 << 1) limit, the macroimage consists of a bright primary image and a large number of faint secondary images formed close to each of the point masses. Taking into account the correlations between the primary and secondary images, we derive P(A) for low A. The low-A distribution has a peak of amplitude ~ 1/tau^2 at A-1 ~ tau^2 and matches smoothly to the high-A distribution. We combine the high- and low-A results and obtain a practical semi-analytic expression for P(A). This semi-analytic distribution is in qualitative agreement with previous numerical results, but the latter show stronger caustic-induced features at moderate A for tau as small as 0.1. We resolve this discrepancy by re-examining the criterion for low optical depth. A simple argument shows that the fraction of caustics of individual lenses that merge with those of their neighbors is approx 1-exp(-8 tau). For tau=0.1, the fraction is surprisingly high: approx 55%. For the purpose of computing P(A) in the manner we did, low optical depth corresponds to tau << 1/8.

In Paper I we studied the theory of gravitational microlensing for a planar distribution of point masses. In this second paper, we extend the analysis to a three-dimensional lens distribution. First we study the lensing properties of three-dimensional lens distributions by considering in detail the critical curves, the caustics, the illumination patterns, and the magnification cross-sections sigma(A) of multiplane configurations with two, three, and four point masses. For N_\ast point masses that are widely separated in Lagrangian space, we find that there are ~ 2^{N_\ast} - 1 critical curves in total, but that only ~ N_\ast of these produce prominent caustic-induced features at the high magnification end of sigma(A). In the case of a low optical depth random distribution of point masses, we show that the multiplane lens equation near a point mass can be reduced to the single plane equation of a point mass perturbed by weak shear. This allows us to calculate the caustic-induced feature in the macroimage magnification distribution P(A) as a weighted sum of the semi-analytic feature derived in Paper I for a planar lens distribution. The resulting semi-analytic caustic-induced feature is similar to the feature in the planar case, but it does not have any simple scaling properties, and it is shifted to higher magnification. The semi-analytic distribution is compared to the results of previous numerical simulations for optical depth tau ~ 0.1, and they are in better agreement than a similar comparison in the planar case. We explain this by estimating the fraction of caustics of individual lenses that merge with those of their neighbors. For tau=0.1, the fraction is approx 20%, much less than the approx 55% for the planar case.

We report on the detection of dark matter in the cluster of galaxies Abell~2163 using the weak gravitational distortion of background galaxies, and an analysis of the cluster X-ray emission. We find that while the qualitative distributions of the cluster light and the dark matter are similar -- shallow and extended, with significant substructure -- the X-ray morphology shows a more regular overall appearance. We interpret the joint lensing and X-ray observations as a signature of a merger event in the cluster. We present new ROSAT/HRI data and reanalyze ROSAT/PSPC data, accounting for the effect of a varying background to determine the best fit parameters in the $\beta$-model formalism. We combine the surface brightness fits with two determinations of the radial temperature profile to determine the total mass. Although there are slight variations in the total mass determinations introduced by the uncertainties in the $\beta$-fit, the main contributor to the error arises from the uncertainties in the temperature determinations. Even though the morphologies of the dark matter/light and X-ray gas are quite different, we find that the total mass determined from the X-ray and weak lensing estimates are consistent with each other within the $2\sigma$ error bars, with the X-ray inferred mass a factor of $\simeq 2$ larger. However, as the lensing mass estimates are differential (the surface density at any point is determined relative to the mean in a control annulus), the shallow, extended nature of the mass profile biases the lensing inferred mass downwards. We estimate the correction for this effect and find very good agreement between the corrected lensing and X-ray results. We determine the gas mass fraction and find $f_g \simeq 0.07h^{-3/2}$ at all radii and a constant mass-to-light ratio of $M/L_V

We report on the detection of dark matter in the cluster Abell 2390 using the weak gravitational distortion of background galaxies. We find that the cluster light and total mass distributions are quite similar over an angular scale of $\simeq 7^\prime \;(1 \Mpc$). The cluster galaxy and mass distributions are centered on the cluster cD galaxy and exhibit elliptical isocontours in the central $\simeq 2^\prime \; (280 \kpc)$. The major axis of the ellipticity is aligned with the direction defined by the cluster cD and a ``straight arc'' located $\simeq 38^{\prime\prime}$ to the northwest. We determined the radial mass-to-light profile for this cluster and found a constant value of $(320 \pm 90) h\; M_\odot/L_{\odot V}$, which is consistent with other published determinations. We also compared our weak lensing azimuthally averaged radial mass profile with a spherical mass model proposed by the CNOC group on the basis of their detailed dynamical study of the cluster. We find good agreement between the two profiles, although there are weak indications that the CNOC density profile may be falling more steeply for $\theta\geq 3^\prime$ $(420\kpc)$.

Both theoretical and observational lines of reasoning suggest that the very faint (B > 24) galaxies seen in deep images of the sky are small low-mass galaxies that experienced starburst at 0.5<z<1 and have since faded into low luminosity, low surface brightness objects. We scrutinize this hypothesis in order to determine whether such a model can account for the observed wavelength- dependent number counts as well as redshift, color and size distributions. Low-mass galaxies generically arise in large numbers in hierarchical clustering scenarios. Although expected to form at high redshifts, the metagalactic UV radiation at high redshifts will delay significant star formation (SF) activity in these galaxies until z < 1. We combine these two elements, along with simple assumptions about SF histories and efficiency, to construct our bursting dwarf model. We further augment the model with a phenomenological prescription for the formation and evolution of the locally observed galaxies, and use spectral synthesis/Monte-Carlo methods to generate realistic galaxy catalogs for comparison with observations. For reasonable choices of SF histories for the low-mass galaxies, the model results are in very good agreement with the results of the deep galaxy surveys. Moreover, the results of deep HST imaging studies, of gravitational lensing and galaxy clustering studies, and of regular galaxies at intermediate/high redshifts favor the dwarf-dominated model. We also discuss some generic predictions of the model. For example, the model predicts that the K-band number counts should begin to rise more steeply at magnitudes fainter than K_{AB} = 24-25 and the local field luminosity function should steepen at magnitudes fainter than M_B \approx -16. We also consider the possibility of detecting the low surface

We report on the detection of dark matter in the cluster Abell 2218 using the weak gravitational distortion of background galaxies. We find a highly significant, coherent detection of the distortion in the images of the background galaxies. The inferred 2D mass distribution has a peak that is coincident with the optical and X-ray centroid. The qualitative distributions of the cluster light, the X-ray emission and the dark matter are similar and the projected total mass, gas, and light surface densities are consistent with a $r^{-1}$ profile at distance of $r > 180^{\prime\prime}$ from the cluster cD galaxy. Using the weak lensing technique, we determine a lower bound for the total mass in A2218 of $(3.9 \pm 0.7) \times 10^{14}$~h$^{-1}$~M$_\odot$ within a fiducial aperture of radius 0.4~h$^{-1}$Mpc. The associated cluster mass-to-light ratio is $(440 \pm 80)$~h~$M_\odot/L_{\odot B}$. The mass estimated by the weak lensing method is consistent with that inferred from the X-ray data under the assumption of hydrostatic equilibrium and we derive an upper bound for the gas-to-total mass ratio at 400~h$^{-1}$kpc of $M_{gas}/M_{tot} = (0.04 \pm 0.02)$h$^{-3/2}$.

Long arcs in clusters of galaxies, produced by gravitational lensing, can be used to estimate the mass interior to the arcs and therefore, constrain the cluster mass distribution. The radial density distribution of the intracluster gas (ICM) can be extracted from the X-ray surface brightness observations. If the gas temperature is also known, it is then possible to probe the dynamical state of the gas and test whether the ICM is in hydro- static equilibrium within the gravitational potential of the cluster as a result of thermal pressure support. We analyze three clusters that exhibit large arcs, whose X-ray surface brightness profiles have been observed, and whose gas temperatures have been determined. In two of the clusters, A2218 and A1689, the central mass implied by lensing is a factor of $2$--$2.5$ too large for the gas at the observed temperature to be in hydrostatic equilibrium solely due to thermal pressure support. In other words, if we accept the mass estimate derived from the lensing analysis and demand that the X-ray surface brightness profile be consistent with the observations, the required gas temperature is a factor of $2$--$2.5$ higher than observed. The results for the third cluster, A2163 (the most luminous and the hottest cluster known), are more ambiguous. The discrepancy between the X-ray and the lensing mass estimates arise because the presence of arcs imply a highly concentrated cluster mass distribution whereas the observed X-ray profiles imply a more extended mass distribution. The large X-ray core radii are not the result of the limited resolution of the X-ray detectors. We consider various possibilities that could account for the discrepancy.

We present preliminary results of our investigation into the influence of shear fields on the evolution of galactic scale fluctuations in a primordial Gaussian random density field. Specifically, we study how the matter associated with a galaxy-scale peak evolves, to determine whether the shear can affect the peak's ability to form a virialized structure. We find that the evolution of the mass distribution in the immediate surroundings of $\geq 1\sigma$ initial density peaks is sensitive to the nature and magnitude of the shear. Its final fate and configuration spans a plethora of possibilities, ranging from accretion onto the peak to complete disruption. On the other hand, the mass defining the peak itself always tends to form a virialized object, though a given galaxy-scale peak need not necessarily form only one halo. Furthermore, galaxy-size halos need not necessarily be associated with initial galaxy-scale density peaks. Under certain conditions, the shear field is capable of breaking up a single primordial peak into two (and perhaps, more) distinct halos, or of promoting the growth of smaller-scale peaks into galaxy-size halos.

We challenge a widely accepted assumption of observational cosmology: that successful reconstruction of observed galaxy density fields from measured galaxy velocity fields (or vice versa), using the methods of gravitational instability theory, implies that the observed large-scale structures and large-scale flows were produced by the action of gravity. This assumption is false, in that there exist non-gravitational theories that pass the reconstruction tests and gravitational theories with certain forms of biased galaxy formation that fail them. Gravitational instability theory predicts specific correlations between large-scale velocity and mass density fields, but the same correlations arise in any model where (a) structures in the galaxy distribution grow from homogeneous initial conditions in a way that satisfies the continuity equation, and (b) the present-day velocity field is irrotational and proportional to the time-averaged velocity field. We demonstrate these assertions using analytical arguments and N-body simulations of gravitational and non-gravitational models. We also show examples of gravitational and non-gravitational models that {\it fail} reconstruction tests because galaxy formation is modulated (``biased'') in a way that violates the continuity equation. We discuss the relation between the value of $\Omega$ inferred from velocity-density comparisons and the true cosmological value.