Andreea S. Font

Ram pressure stripping the hot gaseous haloes of galaxies in groups and clusters

We use a large suite of carefully controlled full hydrodynamic simulations to study the ram pressure stripping of the hot gaseous haloes of galaxies as they fall into massive groups and clusters. The sensitivity of the results to the orbit, total galaxy mass, and galaxy structural properties is explored. For typical structural and orbital parameters, we find that ∼30 per cent of the initial hot galactic halo gas can remain in place after 10 Gyr. We propose a physically simple analytic model that describes the stripping seen in the simulations remarkably well. The model is analogous to the original formulation of Gunn & Gott, except that it is appropriate for the case of a spherical (hot) gas distribution (as opposed to a face-on cold disc) and takes into account that stripping is not instantaneous but occurs on a characteristic time-scale. The model reproduces the results of the simulations to within ≈10 per cent at almost all times for all the orbits, mass ratios, and galaxy structural properties we have explored. The one exception involves unlikely systems where the orbit of the galaxy is highly non-radial and its mass exceeds about 10 per cent of the group or cluster into which it is falling (in which case the model underpredicts the stripping following pericentric passage). The proposed model has several interesting applications, including modelling the ram pressure stripping of both observed and cosmologically simulated galaxies and as a way to improve present semi-analytic models of galaxy formation. One immediate consequence is that the colours and morphologies of satellite galaxies in groups and clusters will differ significantly from those predicted with the standard assumption of complete stripping of the hot coronae.

TRACING GALAXY FORMATION WITH STELLAR HALOS II: RELATING SUBSTRUCTURE IN PHASE- AND ABUNDANCE-SPACE TO ACCRETION HISTORIES

This paper explores the mapping between the observable properties of a stellar halo in phase- and abundance-space and the parent galaxy’s accretion history in terms of the characteristic epoch of accretion and mass and orbits of progenitor objects. The study utilizes a suite of eleven stellar halo models constructed within the context of a standard CDM cosmology. The results demonstrate that coordinate-space studies are sensitive to the recent (0-8 Gyears ago) merger histories of galaxies (this timescale corresponds to the last few to tens of percent of mass accretion for a Milky-Way-type galaxy). Specically, the frequency, sky coverage and fraction of stars in substructures in the stellar halo as a function of surface brightness are indicators of the importance of recent merging and of the luminosity function of infalling dwarfs. The morphology of features serves as a guide to the orbital distribution of those dwarfs. Constraints on the earlier merger history (> 8 Gyears ago) can be gleaned from the abundance patterns in halo stars: within our models, dramatic dierences in the dominant epoch of accretion or luminosity function of progenitor objects leave clear signatures in the [ /Fe] and [Fe/H] distributions of the stellar halo | halos dominated by very early accretion have higher average [ /Fe], while those dominated by high luminosity satellites have higher [Fe/H]. This intuition can be applied to reconstruct much about the merger histories of nearby galaxies from current and future data sets. Subject headings: Galaxy: evolution | Galaxy: formation | Galaxy:halo | Galaxy: kinematics and dynamics | galaxies: dwarf | galaxies: evolution | galaxies: formation | galaxies: halos | galaxies: kinematics and dynamics | Local Group | dark matter

The accretion of galaxies into groups and clusters

We use the galaxy stellar mass and halo merger tree information from the semi-analyticmodel galaxy catalogue of Font et al. (2008) to examine the accretion of galaxies into a large sample of groups and clusters, covering a wide range in halo mass (1012.9 to 1015.3 h−1 M⊙), and selected from each of four redshift epochs (z=0, 0.5, 1.0 and 1.5). We find that clusters at all examined redshifts have accreted a significant fraction of their final galaxy populations through galaxy groups. A 1014.5 h−1 M⊙ mass cluster at z=0 has, on average, accreted_ 40% of its galaxies (Mstellar > 109 h−1 M⊙) from halos with masses greater than 1013 h−1 M⊙. Further, the galaxies which are accreted through groups are more massive, on average, than galaxies accreted through smaller halos or from the field population. We find that at a given epoch, the fraction of galaxies accreted from isolated environments is independent of the final cluster or group mass. In contrast, we find that observing a cluster of the same halo mass at each redshift epoch implies different accretion rates of isolated galaxies, from 5-6 % per Gyr at z=0 to 15% per Gyr at z=1.5. We find that combining the existence of a Butcher Oemler effect at z=0.5 and the observations that galaxies within groups display significant environmental effects with galaxy accretion histories justifies striking conclusions. Namely, that the dominant environmental process must begin to occur in halos of 1012 – 1013 h−1 M⊙, and act over timescales of > 2 Gyrs. This argues in favor of a mechanism like “strangulation”, in which the hot halo of a galaxy is stripped upon infalling into a more massive halo . This simple model predicts that by z=1.5 galaxy groups and clusters will display little to no environmental effects. This conclusion may limit the effectiveness of red sequence cluster finding methods at high redshift.

Cosmological simulations of the formation of the stellar haloes around disc galaxies

We use the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamical simulations to study the formation of stellar spheroids of Milky Way mass disc galaxies. The simulations contain accurate treatments of metal-dependent radiative cooling, star formation, supernova feedback and chemodynamics, and the large volumes that have been simulated yield an unprecedentedly large sample of ≈400 simulated ∼L∗ disc galaxies. The simulated galaxies are surrounded by low-mass, low surface brightness stellar haloes that extend out to ∼100 kpc and beyond. The diffuse stellar distributions bear a remarkable resemblance to those observed around the Milky Way, M31 and other nearby galaxies, in terms of mass density, surface brightness and metallicity profiles. We show that in situ star formation typically dominates the stellar spheroids by mass at radii of r 30 kpc, whereas accretion of stars dominates at larger radii and this change in origin induces a change in the slope of the surface brightness and metallicity profiles, which is also present in the observational data. The system-to-system scatter in the in situ mass fractions of the spheroid, however, is large and spans over a factor of 4. Consequently, there is a large degree of scatter in the shape and normalization of the spheroid density profile within r 30 kpc (e.g. when fitted by a spherical power-law profile, the indices range from −2.6 to −3.4). We show that the in situ mass fraction of the spheroid is linked to the formation epoch of the system. Dynamically, older systems have, on average, larger contributions from in situ star formation, although there is significant system-to-system scatter in this relationship. Thus, in situ star formation likely represents the solution to the long-standing failure of pure accretion-based models to reproduce the observed properties of the inner spheroid.

A COLD DARK MATTER, STELLAR FEEDBACK, AND THE GALACTIC HALO ABUNDANCE PATTERN

The hierarchical formation scenario for the stellar halo requires the accretion and disruption of dwarf galaxies, yet low-metallicity halo stars are enriched in � -elements compared to similar, low-metallicity stars in dwarf spheroidal (dSph) galaxies. We address this primary challenge for the hierarchical formation scenario for the stellar halo by combining chemical evolution modeling with cosmologically motivated mass accretion histories for the MilkyWaydarkhaloanditssatellites.Wedemonstratethatstellarhaloanddwarfgalaxyabundancepatternscanbe explained naturally within the CDM framework. Our solution relies fundamentally on the CDM model prediction that the majority of the stars in the stellar halo were formed within a few relatively massive, � 5 ;10 10 M� , dwarf irregular (dIrr) sized dark matter halos, which were accreted and destroyed � 10 Gyr in the past. These systems necessarily have short-lived, rapid star formation histories, are enriched primarily by Type II supernovae, and host stars with enhanced [� /Fe] abundances. In contrast, dwarf dSph galaxies exist within low-mass dark matter hostsof � 10 9 M� , wheresupernovaewindsareimportantin settingtheintermediate[� /Fe]ratios observed. Our model includes enrichment from Type Ia and Type II supernovae, as well as stellar winds, and includes a physicallymotivatedsupernovaefeedbackprescriptioncalibratedtoreproducethelocaldwarfgalaxystellarmass– metallicity relation. We use representative examples of the type of dark matter halos that we expect to host a destroyed ‘‘stellar halo progenitor’’ dwarf, a surviving dIrr, and a surviving dSph galaxy, and show that their derived abundance patterns, stellar masses, and gas masses are consistent with those observed for each type of system. Our model also self-consistently reproduces the observed stellar mass–vcirc relation for local group satellites and produces the correct cumulative mass for the Milky Way stellar halo. We predict that the lowest metallicity stars in intermediate-mass dIrr galaxies such as the SMC and LMC should follow abundance patterns similar to that observed in the stellar halo. Searches for accreted, disrupted, low-mass dwarfs may be enhanced by searching for unbound stars with dSph-like chemical abundance patterns.

Chemical Abundance Distributions of Galactic Halos and Their Satellite Systems in a ΛCDM Universe

We present a cosmologically motivated model for the hierarchical formation of the stellar halo that includes a semianalytic treatment of galactic chemical enrichment coupled to numerical simulations that track the orbital evolution and tidal disruption of satellites. A major motivating factor in this investigation is the observed systematic difference between the chemical abundances of stars in satellite galaxies and those in the Milky Way halo. Specifically, for the same [Fe/H] values, stars in neighboring satellite galaxies display significantly lower [α/Fe] ratios than stars in the halo. We find that the observed chemical abundance patterns are a natural outcome of the process of hierarchical assembly of the Galaxy. This result follows because the stellar halo in this context is built up from satellite galaxies accreted early on (more than 8-9 Gyr ago) and enriched in α-elements produced in Type II supernovae. In contrast, satellites that still survive today are typically accreted late (within the last 4-5 Gyr) with nearly solar [α/Fe] values as a result of contributions from both Type II and Type Ia supernovae. We use our model to investigate the abundance distribution functions (using both [Fe/H] and [α/Fe] ratios) for stars in the halo and in surviving satellites. Our results suggest that the shapes and peaks in the abundance distribution functions provide a direct probe of the accretion histories of galaxies.

Photoevaporation of Circumstellar Disks around Young Stars

We examine the ability of photoevaporative disk winds to explain the low-velocity components observed in the forbidden line spectra of low-mass T Tauri stars. Using the analytic model of Shu and coworkers and Hollenbach and coworkers as a basis, we examine the characteristics of photoevaporative outflows with hydrodynamic simulations. General results from the simulations agree well with the analytic predictions, although some small differences are present. Most importantly, the flow of material from the disk surface develops at smaller radii than in the analytic approximations, and the flow velocity from the disk surface is only one-third the sound speed. A detailed presentation of observational consequences of the model is given, including predicted line widths, blueshifts, and integrated luminosities of observable sulfur and nitrogen emission lines. We demonstrate that these predictions are in agreement with current observational data on the low-velocity forbidden line emission of ionized species from T Tauri stars. This is in contrast to magnetic wind models, which systematically underpredict these forbidden line luminosities. However, the present model cannot easily account for the luminosities of neutral oxygen lines in T Tauri stars.

Mismatch and misalignment: dark haloes and satellites of disc galaxies

We study the phase-space distribution of satellite galaxies associated with late-type galaxies in the GIMIC suite of simulations. GIMIC consists of resimulations of five cosmologically representative regions from the Millennium Simulation, which have higher resolution and incorporate baryonic physics. Whilst the disc of the galaxy is well aligned with the inner regions (r ∼ 0.1r200) of the dark matter halo, both in shape and angular momentum, there can be substantial misalignments at larger radii (r ∼r200). Misalignments of >45 ◦ are seen in ∼30 per cent of our sample. We find that the satellite population aligns with the shape (and angular momentum) of the outer dark matter halo. However, the alignment with the galaxy is weak owing to the mismatch between the disc and dark matter halo. Roughly 20 per cent of the satellite systems with 10 bright galaxies within r200 exhibit a polar spatial alignment with respect to the galaxy – an orientation reminiscent of the classical satellites of the Milky Way. We find that a small fraction (∼10 per cent) of satellite systems show evidence for rotational support which we attribute to group infall. There is a bias towards satellites on prograde orbits relative to the spin of the dark matter halo (and to a lesser extent with the angular momentum of the disc). This preference towards co-rotation is stronger in the inner regions of the halo where the most massive satellites accreted at relatively early times are located. We attribute the anisotropic spatial distribution and angular momentum bias of the satellites at z = 0 to their directional accretion along the major axes of the dark matter halo. The satellite galaxies have been accreted relatively recently compared to the dark matter mass and have experienced less phase-mixing and relaxation – the memory of their accretion history can remain intact to z = 0. Understanding the phase-space distribution of the z = 0 satellite population is key for studies that estimate the host halo mass from the line-of-sight velocities and projected positions of satellite galaxies. We quantify the effects of such systematics in estimates of the host halo mass from the satellite population.

The population of Milky Way satellites in the LambdaCDM cosmology

We present a model for the satellites of the Milky Way in which galaxy formation is followed using semi-analytic techniques applied to the six high-resolution N-body simulations of galactic halos of the Aquarius project. The model, calculated using the Galform code, incorporates improved treatments of the relevant physics in the LambdaCDM cosmogony, particularly a self-consistent calculation of reionization by UV photons emitted by the forming galaxy population, including the progenitors of the central galaxy. Along the merger tree of each halo, the model calculates gas cooling (by Compton scattering off cosmic microwave background photons, molecular hydrogen and atomic processes), gas heating (from hydrogen photoionization and supernova energy), star formation and evolution. The evolution of the intergalactic medium is followed simultaneously with that of the galaxies. Star formation in the more massive progenitor subhalos is suppressed primarily by supernova feedback, while for smaller subhalos it is suppressed primarily by photoionization due to external and internal sources. The model is constrained to match a wide range of properties of the present day galaxy population as a whole, but at high redshift it requires an escape fraction of UV photons near unity in order completely to reionize the universe by redshift z ~ 8. In the most successful model the local sources photoionize the pre-galactic region completely by z ~ 10. In addition to the luminosity function of Milky Way satellites, the model matches their observed luminosity-metallicity relation, their radial distribution and the inferred values of the mass within 300 pc, which in the models increase slowly but significantly with luminosity. There is a large variation in satellite properties from halo to halo, with the luminosity function, for example, varying by a factor of ~ 2 among the six simulations.

The population of Milky Way satellites in the Λ cold dark matter cosmology

We present a model for the satellites of the Milky Way in which galaxy formation is followed using semi-analytic techniques applied to the six high-resolution N-body simulations of galactic halos of the Aquarius project. The model, calculated using the Galform code, incorporates improved treatments of the relevant physics in the LambdaCDM cosmogony, particularly a self-consistent calculation of reionization by UV photons emitted by the forming galaxy population, including the progenitors of the central galaxy. Along the merger tree of each halo, the model calculates gas cooling (by Compton scattering off cosmic microwave background photons, molecular hydrogen and atomic processes), gas heating (from hydrogen photoionization and supernova energy), star formation and evolution. The evolution of the intergalactic medium is followed simultaneously with that of the galaxies. Star formation in the more massive progenitor subhalos is suppressed primarily by supernova feedback, while for smaller subhalos it is suppressed primarily by photoionization due to external and internal sources. The model is constrained to match a wide range of properties of the present day galaxy population as a whole, but at high redshift it requires an escape fraction of UV photons near unity in order completely to reionize the universe by redshift z ~ 8. In the most successful model the local sources photoionize the pre-galactic region completely by z ~ 10. In addition to the luminosity function of Milky Way satellites, the model matches their observed luminosity-metallicity relation, their radial distribution and the inferred values of the mass within 300 pc, which in the models increase slowly but significantly with luminosity. There is a large variation in satellite properties from halo to halo, with the luminosity function, for example, varying by a factor of ~ 2 among the six simulations.