Origin and structure of the Galactic disc(s)

We examine the chemical and dynamical structure in the solar neighbourhood of a model Galaxy that is the endpoint of a simulation of the chemical evolution of the Milky Way in the presence of radial mixing of stars and gas. Although the simulation's star formation rate declines monotonically from its unique peak and no merger or tidal event ever takes place, the model replicates all known properties of a thick disc, as well as matching special features of the local stellar population such as a metal-poor extension of the thin disc that has high rotational velocity. We divide the disc by chemistry and relate this dissection to observationally more convenient kinematic selection criteria. We conclude that the observed chemistry of the Galactic disc does not provide convincing evidence for a violent origin of the thick disc, as has been widely claimed.     

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 star formation rate in disc galaxies: thresholds and dependence on gas amount

We reassess the applicability of the Toomre criterion in galactic discs and we study the local star formation law in 16 disc galaxies for which abundance gradients are published. The data we use consist of stellar light profiles, atomic and molecular gas (deduced from CO with a metallicity-dependent conversion factor), star formation rates (from Hα emissivities), metallicities, dispersion velocities and rotation curves. We show that the Toomre criterion applies successfully to the case of the Milky Way disc, but it has limited success with the data of our sample; depending on whether or not the stellar component is included in the stability analysis, we find average values for the threshold ratio of the gas surface density to the critical surface density in the range 0.5–0.7. We also test various star formation laws proposed in the literature, i.e. either the simple Schmidt law or modifications of it, that take into account dynamical factors. We find only small differences among them as far as the overall fit to our data is concerned; in particular, we find that all three star formation laws (with parameters derived from the fits to our data) match observations in the Milky Way disc particularly well. In all cases we find that the exponent n of our best-fitting star formation rate has slightly higher values than in other recent works and we suggest several reasons that may cause that discrepancy.     

A model of supernova feedback in galaxy formation

A model of supernova feedback during disc galaxy formation is developed. The model incorporates infall of cooling gas from a halo, and outflow of hot gas from a multiphase interstellar medium (ISM). The star formation rate is determined by balancing the energy dissipated in collisions between cold gas clouds with that supplied by supernovae in a disc marginally unstable to axisymmetric instabilities. Hot gas is created by thermal evaporation of cold gas clouds in supernova remnants, and criteria are derived to estimate the characteristic temperature and density of the hot component and hence the net mass outflow rate. A number of refinements of the model are investigated, including a simple model of a galactic fountain, the response of the cold component to the pressure of the hot gas, pressure-induced star formation and chemical evolution. The main conclusion of this paper is that low rates of star formation can expel a large fraction of the gas from a dwarf galaxy. For example, a galaxy with circular speed 50 km s¹ can expel 6080 per cent of its gas over a time-scale of 1 Gyr, with a star formation rate that never exceeds 0.1 M yr¹. Effective feedback can therefore take place in a quiescent mode and does not require strong bursts of star formation. Even a large galaxy, such as the Milky Way, might have lost as much as 20 per cent of its mass in a supernova-driven wind. The models developed here suggest that dwarf galaxies at high redshifts will have low average star formation rates and may contain extended gaseous discs of largely unprocessed gas. Such extended gaseous discs might explain the numbers, metallicities and metallicity dispersions of damped Lyman systems.     

The SAURON project – VII. Integral-field absorption and emission-line kinematics of 24 spiral galaxy bulges

We present observations of the stellar and gas kinematics for a representative sample of 24 Sa galaxies obtained with our custom-built integral-field spectrograph SAURON operating on the William Herschel Telescope. The data have been homogeneously reduced and analysed by means of a dedicated pipeline. All resulting data cubes were spatially binned to a minimum mean signal-to-noise ratio of 60 per spatial and spectral resolution element. Our maps typically cover the bulge-dominated region. We find a significant fraction of kinematically decoupled components (12/24), many of them displaying central velocity dispersion minima. They are mostly aligned and co-rotating with the main body of the galaxies, and are usually associated with dust discs and rings detected in unsharp-masked images. Almost all the galaxies in the sample (22/24) contain significant amounts of ionized gas which, in general, is accompanied by the presence of dust. The kinematics of the ionized gas are consistent with circular rotation in a disc co-rotating with respect to the stars. The distribution of mean misalignments between the stellar and gaseous angular momenta in the sample suggests that the gas has an internal origin. The [O  iii ]/Hβ ratio is usually very low, indicative of current star formation, and shows various morphologies (ring-like structures, alignments with dust lanes or amorphous shapes). The star formation rates (SFRs) in the sample are comparable with that of normal disc galaxies. Low gas velocity dispersion values appear to be linked to regions of intense star formation activity. We interpret this result as stars being formed from dynamically cold gas in those regions. In the case of NGC 5953, the data suggest that we are witnessing the formation of a kinematically decoupled component from cold gas being acquired during the ongoing interaction with NGC 5954.     

Coupled galactic density-wave modes in a composite system of thin stellar and gaseous discs

A rotating disc galaxy is modelled as a composite system consisting of thin stellar and gaseous discs, which are described by a two-fluid modal formalism. The composite disc system is assumed to retain axisymmetry in the background equilibrium. General density-wave perturbations in the two discs are coupled through the mutual gravitational interaction. We study the basic properties of open and tight spiral density-wave modes in such a composite disc system. Within the Lindblad resonances, perturbation enhancements of surface mass density in stellar and gaseous discs are in phase; this is also true during the initial growth phase of density-wave perturbations. Outside the Lindblad resonances, there exists a possible spiral density-wave branch for which perturbation enhancements of surface mass density in stellar and gaseous discs are out of phase. We discuss implications of these results on the critical parameters for global star formation in barred and normal spiral galaxies and on magnetohydrodynamic density waves within the Lindblad resonances.     

The origin of the light distribution in spiral galaxies

We analyse a high-resolution, fully cosmological, hydrodynamical disc galaxy simulation, to study the source of the double-exponential light profiles seen in many stellar discs, and the effects of stellar radial migration upon the spatiotemporal evolution of both the disc age and metallicity distributions. We find a 'break' in the pure exponential stellar surface brightness profile, and trace its origin to a sharp decrease in the star formation per unit surface area, itself produced by a decrease in the gas volume density due to a warping of the gas disc. Star formation in the disc continues well beyond the break. We find that the break is more pronounced in bluer wavebands. By contrast, we find little or no break in the mass density profile. This is, in part, due to the net radial migration of stars towards the external parts of the disc. Beyond the break radius, we find that ∼60 per cent of the resident stars migrated from the inner disc, while ∼25 per cent formed in situ . Our simulated galaxy also has a minimum in the age profile at the break radius but, in disagreement with some previous studies, migration is not the main mechanism producing this shape. In our simulation, the disc metallicity gradient flattens with time, consistent with an 'inside-out' formation scenario. We do not find any difference in the intensity or the position of the break with inclination, suggesting that perhaps the differences found in empirical studies are driven by dust extinction.     

The formation and survival of discs in a ΛCDM universe

We study the formation of galaxies in a Λ cold dark matter (ΛCDM) universe using high-resolution hydrodynamical simulations with a multiphase treatment of gas, cooling and feedback, focusing on the formation of discs. Our simulations follow eight isolated haloes similar in mass to the Milky Way and extracted from a large cosmological simulation without restriction on spin parameter or merger history. This allows us to investigate how the final properties of the simulated galaxies correlate with the formation histories of their haloes. We find that, at   z = 0  , none of our galaxies contains a disc with more than 20 per cent of its total stellar mass. Four of the eight galaxies nevertheless have well-formed disc components, three have dominant spheroids and very small discs, and one is a spheroidal galaxy with no disc at all. The   z = 0  spheroids are made of old stars, while discs are younger and formed from the inside-out. Neither the existence of a disc at   z = 0  nor the final disc-to-total mass ratio seems to depend on the spin parameter of the halo. Discs are formed in haloes with spin parameters as low as 0.01 and as high as 0.05; galaxies with little or no disc component span the same range in spin parameter. Except for one of the simulated galaxies, all have significant discs at   z ≳ 2  , regardless of their   z = 0  morphologies. Major mergers and instabilities which arise when accreting cold gas is misaligned with the stellar disc trigger a transfer of mass from the discs to the spheroids. In some cases, discs are destroyed, while in others, they survive or reform. This suggests that the survival probability of discs depends on the particular formation history of each galaxy. A realistic ΛCDM model will clearly require weaker star formation at high redshift and later disc assembly than occurs in our models.     

Tidal evolution of discy dwarf galaxies in the Milky Way potential: the formation of dwarf spheroidals

We conduct high-resolution collisionless N -body simulations to investigate the tidal evolution of dwarf galaxies on an eccentric orbit in the Milky Way (MW) potential. The dwarfs originally consist of a low surface brightness stellar disc embedded in a cosmologically motivated dark matter halo. During 10 Gyr of dynamical evolution and after five pericentre passages, the dwarfs suffer substantial mass loss and their stellar component undergoes a major morphological transformation from a disc to a bar and finally to a spheroid. The bar is preserved for most of the time as the angular momentum is transferred outside the galaxy. A dwarf spheroidal (dSph) galaxy is formed via gradual shortening of the bar. This work thus provides a comprehensive quantitative explanation of a potentially crucial morphological transformation mechanism for dwarf galaxies that operates in groups as well as in clusters. We compare three cases with different initial inclinations of the disc and find that the evolution is fastest when the disc is coplanar with the orbit. Despite the strong tidal perturbations and mass loss, the dwarfs remain dark matter dominated. For most of the time, the one-dimensional stellar velocity dispersion, σ, follows the maximum circular velocity, V _max, and they are both good tracers of the bound mass. Specifically, we find that   M _bound∝ V ^3.5_max  and     in agreement with earlier studies based on pure dark matter simulations. The latter relation is based on directly measuring the stellar kinematics of the simulated dwarf, and may thus be reliably used to map the observed stellar velocity dispersions of dSphs to halo circular velocities when addressing the missing satellites problem.     

Chemo-spectrophotometric evolution of spiral galaxies – V. Properties of galactic discs at high redshift

We explore the implications for the high-redshift universe of 'state-of-the-art' models for the chemical and spectrophotometric evolution of spiral galaxies. The models are based on simple 'scaling relations' for discs, obtained in the framework of cold dark matter models for galaxy formation, and were 'calibrated' so as to reproduce the properties of the Milky Way and of nearby discs (at redshift z ∼0) . In this paper, we compare the predictions of our 'hybrid' approach to galaxy evolution to observations at moderate and high redshift. We find that the models are in fairly good agreement with observations up to z ∼1 , while some problems appear at higher redshift (provided there is no selection bias in the data); these discrepancies may suggest that galaxy mergers (not considered in this work) played a non-negligible role at z >1 . We also predict the existence of a 'universal' correlation between abundance gradients and disc scalelengths, independent of redshift.     

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.     

Simulations of minor mergers – II. The phase-space structure of thick discs

We analyse the phase-space structure of simulated thick discs that are the result of a 5:1 mass-ratio merger between a disc galaxy and a satellite. Our main goal is to establish what would be the imprints of a merger origin for the Galactic thick disc. We find that the spatial distribution predicted for thick-disc stars is asymmetric, seemingly in agreement with recent observations of the Milky Way thick disc. Near the Sun, the accreted stars are expected to rotate more slowly, to have broad velocity distributions and to occupy preferentially the wings of the line-of-sight velocity distributions. The majority of the stars in our model thick discs have low eccentricity orbits (in clear reference to the pre-existing heated disc) which give rise to a characteristic (sinusoidal) pattern for their line-of-sight velocities as a function of galactic longitude. The z -component of the angular momentum of thick-disc stars provides a clear discriminant between stars from the pre-existing disc and those from the satellite, particularly at large radii. These results are robust against the particular choices of initial conditions made in our simulations.     

Chemo-spectrophotometric evolution of spiral galaxies – III. Abundance and colour gradients in discs

We study the relations between luminosity and chemical-abundance profiles of spiral galaxies, using detailed models for the chemical and spectrophotometric evolution of galactic discs. The models are 'calibrated' on the Milky Way disc and are successfully extended to other discs with the help of simple 'scaling' relations, obtained in the framework of semi-analytic models of galaxy formation. We find that our models exhibit oxygen abundance gradients that increase in absolute value with decreasing disc luminosity (when expressed in dex kpc⁻¹) and are independent of disc luminosity (when expressed in dex scalelength⁻¹), both in agreement with observations. We notice an important strong correlation between abundance gradient and disc scalelength. These results support the idea of 'homologous evolution' of galactic discs.     

Chemo-spectrophotometric evolution of spiral galaxies – II. Main properties of present-day disc galaxies

We study the chemical and spectrophotometric evolution of galactic discs with detailed models calibrated on the Milky Way and using simple scaling relations, based on currently popular semi-analytic models of galaxy formation. We compare our results with a large body of observational data on present-day galactic discs, including disc sizes and central surface brightness, Tully–Fisher relations in various wavelength bands, colour–colour and colour–magnitude relations, gas fractions versus magnitudes and colours and abundances versus local and integrated properties, as well as spectra for different galactic rotational velocities. Despite the extremely simple nature of our models, we find satisfactory agreement with all those observables, provided that the time-scale for star formation in low-mass discs is longer than for more massive ones. This assumption is apparently in contradiction with the standard picture of hierarchical cosmology. We find, however, that it is extremely successful in reproducing major features of present-day discs, like the change in the slope of the Tully–Fisher relation with wavelength, the fact that more massive galaxies are on average 'redder' than low-mass ones (a generic problem of standard hierarchical models) and the metallicity–luminosity relation for spirals. It is concluded that, on a purely empirical basis, this new picture is at least as successful as the standard one. Observations at high redshifts could help to distinguish between the two possibilities.     

NGC 3741: the dark halo profile from the most extended rotation curve

We present new H  i observations of the nearby dwarf galaxy NGC 3741. This galaxy has an extremely extended H  i disc, which allows us to trace the rotation curve out to unprecedented distances in terms of the optical disc: we reach 42 B -band exponential scalelengths or about 7 kpc. The H  i disc is strongly warped, but the warp is very symmetric. The distribution and kinematics are accurately derived by building model data cubes, which closely reproduce the observations. In order to account for the observed features in the data cube, radial motions of the order of 5–13 km s⁻¹ are needed. They are consistent with an inner bar of several hundreds of pc and accretion of material in the outer regions.
The observed rotation curve was decomposed into its stellar, gaseous and dark components. The Burkert dark halo (with a central constant density core) provides very good fits. The dark halo density distribution predicted by the Λ cold dark matter (CDM) theory fails to fit the data, unless NGC 3741 is a 2.5σ exception to the predicted relation between concentration parameter and virial mass and at the same time a high value of the virial mass (though poorly constrained) of  10¹¹ M_⊙  . Noticeably, modified Newtonian dynamics (MOND) seems to be consistent with the observed rotation curve. Scaling up the contribution of the gaseous disc also gives a good fit.     

Compact stellar systems around NGC 1399

We have obtained spectroscopic redshifts of colour-selected point sources in four wide area VLT-FLAMES (Very Large Telescope-Fibre Large Array Multi Element Spectrograph) fields around the Fornax cluster giant elliptical galaxy NGC 1399, identifying as cluster members 27 previously unknown faint     compact stellar systems (CSS), and improving redshift accuracy for 23 previously catalogued CSS.
By amalgamating our results with CSS from previous 2dF observations and excluding CSS dynamically associated with prominent (non-dwarf) galaxies surrounding NGC 1399, we have isolated 80 'unbound' systems that are either part of NGC  1399's globular cluster (GC) system or intracluster GCs. For these unbound systems, we find (i) they are mostly located off the main stellar locus in colour–colour space; (ii) their projected distribution about NGC  1399 is anisotropic, following the Fornax cluster galaxy distribution, and there is weak evidence for group rotation about NGC  1399; (iii) their completeness-adjusted radial surface density profile has a slope similar to that of NGC  1399's inner GC system; (iv) their mean heliocentric recessional velocity is between that of NGC  1399's inner GCs and that of the surrounding dwarf galaxies, but their velocity dispersion is significantly lower; (v) bright CSS  ( M _V < −11)  are slightly redder than the fainter systems, suggesting they have higher metallicity; (vi) CSS show no significant trend in   g '− i '  colour index with radial distance from NGC  1399.     

Do high-velocity clouds form by thermal instability?

We examine the proposal that the H  i 'high-velocity' clouds (HVCs) surrounding the Milky Way and other disc galaxies form by condensation of the hot galactic corona via thermal instability. Under the assumption that the galactic corona is well represented by a non-rotating, stratified atmosphere, we find that for this formation mechanism to work the corona must have an almost perfectly flat entropy profile. In all other cases, the growth of thermal perturbations is suppressed by a combination of buoyancy and thermal conduction. Even if the entropy profile were nearly flat, cold clouds with sizes smaller than  10 kpc  could form in the corona of the Milky Way only at radii larger than  100 kpc  , in contradiction with the determined distances of the largest HVC complexes. Clouds with sizes of a few kpc can form in the inner halo only in low-mass systems. We conclude that unless even slow rotation qualitatively changes the dynamics of a corona, thermal instability is unlikely to be a viable mechanism for formation of cold clouds around disc galaxies.     

Inhomogeneous chemical evolution of dwarf spheroidal galaxies

Stellar abundance pattern of n-capture elements such as barium is used as a powerful tool to infer how the star formation proceeded in dwarf spheroidal (dSph) galaxies. It is found that the abundance correlation of barium with iron in stars belonging to dSph galaxies orbiting the Milky Way, i.e., Draco, Sextans, and Ursa Minor have a feature similar to that in Galactic metal-poor stars. The common feature of these two correlations can be realized by our in homogeneous chemical evolution model based on the supernova-driven star formation scenario if dSph stars formed from gas with a velocity dispersion of ∼ 26 km s^-1. This velocity dispersion together with the stellar luminosities strongly suggest that dark matter dominated dSph galaxies. The tidal force of the Milky Way links this velocity dispersion with the currently observed value ≲ 10 km s^-1 by stripping the dark matter in dSph galaxies. As a result, the total mass of each dSph galaxy is found to have been originally ∼ 25 times larger than at present. In this model, supernovae immediately after the end of the star formation can expel the remaining gas over the gravitational potential of the dSph galaxy. This revised version was published online in August 2006 with corrections to the Cover Date.     

The origin of the density distribution of disc galaxies: a new problem for the standard model of disc formation

We present new models for the formation of disc galaxies that improve upon previous models by following the detailed accretion and cooling of the baryonic mass, and by using realistic distributions of specific angular momentum. Under the assumption of detailed angular momentum conservation, the discs that form have density distributions that are more centrally concentrated than an exponential. We examine the influence of star formation, bulge formation, and feedback on the outcome of the surface brightness distributions of the stars. Low angular momentum haloes yield disc galaxies with a significant bulge component and with a stellar disc that is close to exponential, in good agreement with observations. High angular momentum haloes, on the other hand, produce stellar discs that are much more concentrated than an exponential, in clear conflict with observations. At large radii, the models reveal distinct truncation radii in both the stars and the cold gas. The stellar truncation radii result from our implementation of star formation threshold densities, and are in excellent agreement with observations. The truncation radii in the density distribution of the cold gas reflect the maximum specific angular momentum of the gas that has cooled. We find that these truncation radii occur at H  i surface densities of roughly 1 M_⊙ pc⁻², in conflict with observations. We examine various modifications to our models, including feedback, viscosity, and dark matter haloes with constant-density cores, but show that the models consistently fail to produce bulge less discs with exponential surface brightness profiles. This signals a new problem for the standard model of disc formation: if the baryonic component of the protogalaxies out of which disc galaxies form has the same angular momentum distribution as the dark matter, discs are too compact.