Stellar archaeology: Exploring the Universe with metal‐poor stars

The abundance patterns of the most metal‐poor stars in the Galactic halo and small dwarf galaxies provide us with a wealth of information about the early Universe. In particular, these old survivors allow us to study the nature of the first stars and supernovae, the relevant nucleosynthesis processes responsible for the formation and evolution of the elements, early star‐ and galaxy formation processes, as well as the assembly process of the stellar halo from dwarf galaxies a long time ago. This review presents the current state of the field of “stellar archaeology” – the diverse use of metal‐poor stars to explore the high‐redshift Universe and its constituents. In particular, the conditions for early star formation are discussed, how these ultimately led to a chemical evolution, and what the role of the most iron‐poor stars is for learning about Population III supernovae yields. Rapid neutron‐capture signatures found in metal‐poor stars can be used to obtain stellar ages, but also to constrain this complex nucleosynthesis process with observational measurements. Moreover, chemical abundances of extremely metal‐poor stars in different types of dwarf galaxies can be used to infer details on the formation scenario of the halo and the role of dwarf galaxies as Galactic building blocks. I conclude with an outlook as to where this field may be heading within the next decade. A table of ~ 1000 metal‐poor stars and their abundances as collected from the literature is provided in electronic format (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The X—ray Background from the Warm Gas of the Galactic Halo

LETTERS1 INTRODUCTIONIn the hierarchical scenario of structure formation, massive dark ha1os fOrm by gravitationalaggregation of individual low-mass objects, whi1e the stel1ar disks of spiral galaxies like theMilky Way form by accretion of gas which cools and falls onto the galaxies from an extendedsurrounding reservoir. FOr a massive galaxy of M ~ 10"MO, the surrounding gas can be heatedto temperature of T ~ 106 K by gravitational1y-driven shocks, the dominant cooling is thus dueto …     

The environment of formation as a second parameter for globular cluster classification

We perform an evolutionary multivariate analysis of a sample of 54 Galactic globular clusters with high-quality colour–magnitude diagrams and well-determined ages. The four parameters adopted for the analysis are: metallicity, age, maximum temperature on the horizontal branch and absolute V magnitude. Our cladistic analysis breaks the sample into three novel groups. An a posteriori kinematical analysis puts groups 1 and 2 in the halo, and group 3 in the thick disc. The halo and disc clusters separately follow a luminosity–metallicity relation of much weaker slope than galaxies. This property is used to propose a new criterion for distinguishing halo and disc clusters. A comparison of the distinct properties of the two halo groups with those of Galactic halo field stars indicates that the clusters of group 1 originated in the inner halo, while those of group 2 formed in the outer halo of the Galaxy. The inner halo clusters were presumably initially the most massive one, which allowed the formation of more strongly helium-enriched second generation stars, thus explaining the presence of Cepheids and of very hot horizontal-branch stars exclusively in this group. We thus conclude that the 'second parameter' is linked to the environment in which globular clusters form, the inner halo favouring the formation of the most massive clusters which subsequently become more strongly self-enriched than their counterparts of the galactic outer halo and disc.     

Mapping the substructure in the Galactic halo with the next generation of astrometric satellites

We run numerical simulations of the disruption of satellite galaxies in a Galactic potential to build up the entire stellar halo, in order to investigate what the next generation of astrometric satellites will reveal by observing the halo of the Milky Way. We generate artificial DIVA , FAME and GAIA halo catalogues, in which we look for the signatures left by the accreted satellites. We develop a method based on the standard Friends-of-Friends algorithm applied to the space of integrals of motion. We find this simple method can recover about 50 per cent of the different accretion events, when the observational uncertainties expected for GAIA are taken into account, even when the exact form of the Galactic potential is unknown. The recovery rate for DIVA and FAME is much smaller, but these missions, like GAIA , should be able to test the hierarchical formation paradigm on our Galaxy by measuring the amount of halo substructure in the form of nearby kinematically cold streams with, for example, a two-point correlation function in velocity space.     

The kinematics of lopsided galaxies

Lopsidedness is a common feature in galaxies, both in the distribution of light and in the kinematics. We investigate the kinematics of a model for lopsided galaxies that consists of a disc lying off-centre in a dark halo, and circling around the halo centre. We search for families of stable, closed, non-crossing orbits, and assume that gas in our galaxies moves on these orbits. Several of our models show strong lopsided gas kinematics, especially those in which the disc spins around its axis in a retrograde sense compared with its motion around the halo centre. We are able to reproduce the H  i velocity map of the kinematically lopsided galaxy NGC 4395.
The lopsidedness in our models is most pronounced in the models where the halo provides a relatively large fraction of the total mass at small radii. This may explain why the gas shows lopsidedness more frequently in late-type galaxies, which are dominated by dark matter. Surfaces of section show large regions of irregular orbits in the models where the halo density is low. This may indicate that these models are unstable.     

The Structure of the Galactic Halo

We used the star counts in 21 BATC fields obtained with the National Astronomical Observatories (NAOC) 60/90 cm Schmidt Telescope to study the structure of the Galactic halo.Adopting a de Vaucouleurs γ/4 law halo,we found that the halo is somewhat flatter (c/a～0.4) towards the Galactic center than in the anticentre and antirotation direction (c/a ＞ 0.4).We also notice that the axial ratios are smaller (flatter) towards the low latitude fields than the high latitude fields,except for a few fields.We provide robust limits on the large-scale flattening of the halo.Our analysis shows that the axial ratio of the halo may vary with distance and the observation direction.At large Galactocentric radii,the halo may not have a smooth density distribution,but rather,it may be largely composed of overlapping streams or substructures,which provides a support for the hybrid formation model.     

Dynamical evolution of galactic disks driven by interaction with a satellite

Dynamical evolution of galactic disks driven by interaction with satellite galaxies, particularly the problem of the disk warping and thickening is studied numerically. One of the main purpose of the study is to resolve the long standing problem of the origin of the disk warping. A possible cause of the warp is interaction with a satellite galaxy. In the case of the Milky Way, the LMC has been considered as the candidate. Some linear analysis have already given a positive result, but one had to wait for a fully self-consistent simulation as a proof. I have accomplished the numerical simulations with a million particles, by introducing a hybrid algorithm, SCF-TREE. Those simulations give us quantitative estimates for the Milky Way system. We have found an example in which large warp amplitudes are developed. We also found that the warp amplitudes depend on the halo distribution. Among our three models, the most massive and spherical halo is preferable for the observable warp excitation. This revised version was published online in August 2006 with corrections to the Cover Date.     

The Milky Way Galaxy as a strong gravitational lens

We study the gravitational lensing effects of spiral galaxies by taking a model of the Milky Way and computing its lensing properties. The model is composed of a spherical Hernquist bulge, a Miyamoto–Nagai disc and an isothermal halo. As a strong lens, a spiral galaxy like the Milky Way can give rise to four different imaging geometries. They are (i) three images on one side of the galaxy centre ('disc triplets'), (ii) three images with one close to the centre ('core triplets'), (iii) five images and (iv) seven images. Neglecting magnification bias, we show that the core triplets, disc triplets and fivefold imaging are roughly equally likely. Even though our models contain edge-on discs, their image multiplicities are not dominated by disc triplets. The halo is included for completeness, but it has a small effect on the caustic structure, the time delays and brightnesses of the images.
The Milky Way model has a maximum disc (i.e. the halo is not dynamically important in the inner parts). Strong lensing by nearly edge-on disc galaxies breaks the degeneracy between the relative contributions of the disc and halo to the overall rotation curve. If a spiral galaxy has a submaximum disc, then the astroid caustic shrinks dramatically in size, whilst the radial caustic shrinks more modestly. This causes changes in the relative likelihood of the image geometries, specifically (i) core triplets are now ∼9/2 times more likely than disc triplets, (ii) the cross-section for threefold imaging is reduced by a factor of ∼2/3, whilst (iii) the cross-section for fivefold imaging is reduced by ∼1/2. Although multiple imaging is less likely (the cross-sections are smaller), the average total magnification is greater. The time delays are smaller, as the total projected lensing mass is reduced.     

Tidal disruption of globular clusters in dwarf galaxies with triaxial dark matter haloes

We use N -body simulations to study the tidal evolution of globular clusters (GCs) in dwarf spheroidal (dSph) galaxies. Our models adopt a cosmologically motivated scenario in which the dSph is approximated by a static Navarro, Frenk & White halo with a triaxial shape. We apply our models to five GCs spanning three orders of magnitude in stellar density and two in mass, chosen to represent the properties exhibited by the five GCs of the Fornax dSph. We show that only the object representing Fornax's least dense GC (F1) can be fully disrupted by Fornax's internal tidal field – the four denser clusters survive even if their orbits decay to the centre of Fornax. For a large set of orbits and projection angles, we examine the spatial and velocity distribution of stellar debris deposited during the complete disruption of an F1-like GC. Our simulations show that such debris appears as shells, isolated clumps and elongated overdensities at low surface brightness (≥26 mag arcsec⁻²), reminiscent of substructure observed in several Milky Way dSphs. Such features arise from the triaxiality of the galaxy potential and do not dissolve in time. The kinematics of the debris depends strongly on the progenitor's orbit. Debris associated with box and resonant orbits does not display stream motions and may appear 'colder'/'hotter' than the dSph's field population if the viewing angle is perpendicular/parallel to the progenitor's orbital plane. In contrast, debris associated with loop orbits shows a rotational velocity that may be detectable out to a few kpc from the galaxy centre. Chemical tagging that can distinguish GC debris from field stars may reveal whether the merger of GCs contributed to the formation of multiple stellar components observed in dSphs.     

The existence and detection of optically dark galaxies by 21-cm surveys

One explanation for the disparity between cold dark matter (CDM) predictions of galaxy numbers and observations could be that there are numerous dark galaxies in the Universe. These galaxies may still contain baryons, but no stars, and may be detectable in the 21-cm line of atomic hydrogen. The results of surveys for such objects, and simulations that do/do not predict their existence, are controversial. In this paper, we use an analytical model of galaxy formation, consistent with CDM, to first show that dark galaxies are certainly a prediction of the model. Secondly, we show that objects like VIRGOHI21, a dark galaxy candidate recently discovered by us, while rare are predicted by the model. Thirdly, we show that previous 'blind' H  i surveys have placed few constraints on the existence of dark galaxies. This is because they have either lacked the sensitivity and/or velocity resolution or have not had the required detailed optical follow up. We look forward to new 21-cm blind surveys [Arecibo Legacy Fast ALFA (ALFALFA) survey and Arecibo Galactic Environments Survey (AGES)] using the Arecibo multibeam instrument which should find large numbers of dark galaxies if they exist.     

Star clusters around recoiled black holes in the Milky Way halo

Gravitational wave emission by coalescing black holes (BHs) kicks the remnant BH with a typical velocity of hundreds of  km s⁻¹  . This velocity is sufficiently large to remove the remnant BH from a low-mass galaxy but is below the escape velocity from the Milky Way (MW) galaxy. If central BHs were common in the galactic building blocks that merged to make the MW, then numerous BHs that were kicked out of low-mass galaxies should be freely floating in the MW halo today. We use a large statistical sample of possible merger tree histories for the MW to estimate the expected number of recoiled BH remnants present in the MW halo today. We find that hundreds of BHs should remain bound to the MW halo after leaving their parent low-mass galaxies. Each BH carries a compact cluster of old stars that populated the core of its original host galaxy. Using the time-dependent Fokker–Planck equation, we find that the present-day clusters are  ≲1 pc  in size, and their central bright regions should be unresolved in most existing sky surveys. These compact systems are distinguishable from globular clusters by their internal (Keplerian) velocity dispersion greater than 100 km s⁻¹ and their high mass-to-light ratio owing to the central BH. An observational discovery of this relic population of star clusters in the MW halo would constrain the formation history of the MW and the dynamics of BH mergers in the early Universe. A similar population should exist around other galaxies and may potentially be detectable in M31 and M33.     

Mass models of the Milky Way

A parametrized model of the mass distribution within the Milky Way is fitted to the available observational constraints. The most important single parameter is the ratio of the scalelength R _d* of the stellar disc to R ₀. The disc and bulge dominate v _c( R ) at R ≲ R ₀ only for R _d,*/ R ₀≲0.3. Since the only knowledge we have of the halo derives from studies like the present one, we allow it to contribute to the density at all radii. When allowed this freedom, however, the halo causes changes in assumptions relating to R  ≪  R ₀ to affect profoundly the structure of the best-fitting model at R  ≫  R ₀. For example, changing the disc slightly from an exponential surface-density profile significantly changes the form of v _c( R ) at R  ≫  R ₀, where the disc makes a negligible contribution to v _c. Moreover, minor changes in the constraints can cause the halo to develop a deep hole at its centre that is not physically plausible. These problems call into question the proposition that flat rotation curves arise because galaxies have physically distinct haloes rather than outwards-increasing mass-to-light ratios.   The mass distribution of the Galaxy and the relative importance of its various components will remain very uncertain until more observational data can be used to constrain mass models. Data that constrain the Galactic force field at z ≳ R and at R  >  R ₀ are especially important.     

Contrasting the chemical evolution of the Milky Way and Andromeda

The chemical evolution history of a galaxy hides clues about how it formed and has been changing through time. We have studied the chemical evolution history of the Milky Way (MW) and Andromeda (M31) to find which are common features in the chemical evolution of disc galaxies as well as which are galaxy-dependent. We use a semi-analytic multizone chemical evolution model. Such models have succeeded in explaining the mean trends of the observed chemical properties in these two Local Group spiral galaxies with similar mass and morphology. Our results suggest that while the evolution of the MW and M31 shares general similarities, differences in the formation history are required to explain the observations in detail. In particular, we found that the observed higher metallicity in the M31 halo can be explained by either (i) a higher halo star formation efficiency (SFE), or (ii) a larger reservoir of infalling halo gas with a longer halo formation phase. These two different pictures would lead to (i) a higher [O/Fe] at low metallicities, or (ii) younger stellar populations in the M31 halo, respectively. Both pictures result in a more massive stellar halo in M31, which suggests a possible correlation between the halo metallicity and its stellar mass.     

The correlation of star formation quenching with internal galaxy properties and environment

We investigate the correlation of star formation quenching with internal galaxy properties and large-scale environment (halo mass) in empirical data and theoretical models. We make use of the halo-based group catalogue of Yang and collaborators, which is based on the Sloan Digital Sky Survey. Data from the Galaxy evolution explorer are also used to extract the recent star formation rate. In order to investigate the environmental effects, we examine the properties of 'central' and 'satellite' galaxies separately. For central galaxies, we are unable to conclude whether star formation quenching is primarily connected with halo mass or stellar mass, because these two quantities are themselves strongly correlated. For satellite galaxies, a nearly equally strong dependence on halo mass and stellar mass is seen. We make the same comparison for five different semi-analytic models based on three independently developed codes. We find that the models with active galactic nuclei feedback reproduce reasonably well the dependence of the fraction of central red and passive galaxies on halo mass and stellar mass. However, for satellite galaxies, the same models badly overproduce the fraction of red/passive galaxies and do not reproduce the empirical trends with stellar mass or halo mass. This satellite overquenching problem is caused by the too-rapid stripping of the satellites' hot gas haloes, which leads to rapid strangulation of star formation.     

Phase-space density constraints on the formation of bulges from discs

The colours of stellar bulges and of inner stellar discs are comparable, and consistent with rather similar mean metallicities and ages. Indeed, the mean chemical abundances of the Milky Way bulge and old disc are approximately equal. Further, the scalelengths of discs and bulges are correlated. These observations imply a close relationship between discs and bulges, and may support models in which stellar bulges form from stellar discs. The present paper discusses constraints on this scenario from the stellar phase-space density of bulges and of discs. Phase-space density cannot increase in the absence of collisional processes. We show here that the maximum phase-space density of galactic bulges is higher than that of inner discs, arguing that instabilities of purely stellar discs cannot form bulges. Rather, the high densities of bulges probably reflect gaseous dissipation. Gas inflow from the disc would complicate the interpretation of the similarities in stellar colours between discs and bulges. Gas inflow from the stellar halo, if one exists, may be favoured on angular momentum grounds, but this means of formation of the bulge would provide no explanation for the relationships between disc and bulge in any one galaxy. At least in the Milky Way, the metallicity distribution of the bulge is not consistent with the bulge being built up from the dense regions of accreted satellite galaxies and/or globular clusters.     

The APM survey for cool carbon stars in the Galactic halo – I

A by-product of the APM high-redshift quasar survey was the discovery of several distant (20–100 kpc) N-type carbon stars at high galactic latitude. Following on from this, we have started a systematic all-sky survey at galactic latitudes ⊢ b ⊢>30° to find further examples of these rare objects, and we report here on the results from the first season of follow-up spectroscopy. Faint, high-latitude carbon (FHLC) giants make excellent probes of the kinematic structure of the outer Galactic halo. Therefore, in addition to detailed spectrophotometry covering a wide wavelength range, we have obtained high-resolution (∼1 Å) spectra centred on the CN bands at ∼8000 Å, and have derived accurate (≲10 km s⁻¹) radial velocities for the known FHLC stars. From the initial phase of our survey covering ≈6500 deg², we find a surface density of faint N-type carbon stars in the halo of ≈1 per 200 deg², roughly a factor of 4 less than the surface density of CH-type carbon stars in the halo. Intermediate-age, N-type carbon stars seem unlikely to have formed in the halo in isolation from other star-forming regions, and one possibility that we are investigating is that they either arise from the disruption of tidally captured dwarf satellite galaxies or are a manifestation of the long-sought optical component of the Magellanic Stream.     

Point spread function tails and the measurements of diffuse stellar halo light around edge-on disc galaxies

Measuring the integrated stellar halo light around galaxies is very challenging. The surface brightness of these haloes is expected to be many magnitudes below dark sky and the central brightness of the galaxy. Here, I show that in some of the recent literature the effect of very extended Point Spread Function (PSF) tails on the measurements of halo light has been underestimated; especially in the case of edge-on disc galaxies. The detection of a halo along the minor axis of an edge-on galaxy in the Hubble Ultra Deep Field can largely be explained by scattered galaxy light. Similarly, depending on filter and the shape one assumes for the uncertain extended PSF, 20–80 per cent of the halo light found along the minor axis of scaled and stacked Sloan Digital Sky Survey (SDSS) edge-on galaxy images can be explained by scattered galaxy light. Scattered light also significantly contributes to the anomalous halo colours of stacked SDSS images. The scattered light fraction decreases when looking in the quadrants away from the minor axis. The remaining excess light is well modelled with a Sérsic profile halo with shape parameters based on star count halo detections of nearby galaxies. Even though, the contribution from PSF scattered light does not fully remove the need for extended components around these edge-on galaxies, it will be very challenging to make accurate halo light shape and colour measurements from integrated light without very careful PSF measurements and scattered light modelling.     

Galaxy star formation in different environments

We use a semi-analytic model of galaxy formation to study signatures of large-scale modulations in the star formation (SF) activity in galaxies. In order to do this, we carefully define local and global estimators of the density around galaxies. The former are computed using a voronoi tessellation technique and the latter are parametrized by the normalized distance to haloes and voids, in terms of the virial and void radii, respectively. As a function of local density, galaxies show a strong modulation in their SF, a result that is in agreement with those from several authors. When taking subsamples of equal local density at different large-scale environments, we find relevant global effects whereby the fraction of red galaxies diminishes for galaxies in equal local density environments farther away from clusters and closer to voids. In general, the semi-analytic simulation is in good agreement with the available observational results, and offers the possibility to disentangle many of the processes responsible for the variation of galaxy properties with the environment; we find that the changes found in samples of galaxies with equal local environment but different distances to haloes or voids come from the variations in the underlying mass function of dark matter (DM) haloes. There is an additional possible effect coming from the host DM halo ages, indicating that halo assembly also plays a small but significant role (1.14σ) in shaping the properties of galaxies, and in particular, hints at a possible spatial correlation in halo/stellar mass ages. An interesting result comes from the analysis of the coherence of flows in different large-scale environments of fixed local densities; the neighbourhoods of massive haloes are characterized by lower coherences than control samples, except for galaxies in filament-like regions, which show highly coherent motions.     

Where do 'red and dead' early-type void galaxies come from?

Void regions of the Universe offer a special environment for studying cosmology and galaxy formation, which may expose weaknesses in our understanding of these phenomena. Although galaxies in voids are observed to be predominately gas rich, star forming and blue, a subpopulation of bright-red void galaxies can also be found, whose star formation was shutdown long ago. Are the same processes that quench star formation in denser regions of the Universe also at work in voids?
We compare the luminosity function of void galaxies in the 2dF Galaxy Redshift Survey, to those from a galaxy formation model built on the Millennium simulation. We show that a global star formation suppression mechanism in the form of low-luminosity 'radio-mode' active galactic nuclei (AGN) heating is sufficient to reproduce the observed population of void early types. Radio-mode heating is environment independent other than its dependence on dark matter halo mass, where, above a critical mass threshold of approximately   M _vir∼ 10^12.5 M_⊙  , gas cooling on to the galaxy is suppressed and star formation subsequently fades. In the Millennium simulation, the void halo mass function is shifted with respect to denser environments, but still maintains a high-mass tail above this critical threshold. In such void haloes, radio-mode heating remains efficient and red galaxies are found; collectively these galaxies match the observed space density without any modification to the model. Consequently, galaxies living in vastly different large-scale environments but hosted by haloes of similar mass are predicted to have similar properties, consistent with observations.