The formation of a disk galaxy within a slowly growing dark halo

The formation of a disk galaxy within a slowly growing dark halo is simulated with a new chemo-dynamical model. The model describes the evolution of the stellar populations, the multi-phase ISM and all important interaction. I find, that the galaxy forms radially from inside-out and vertically from top-to-bottom. The derived stellar age distributions show that the inner halo is the oldest component, followed by the outer halo, the triaxial bulge, the halo-disk transition region and the disk. Despite the still idealized model, the final galaxy resembles present-day disk galaxies in many aspects. In particular, the stellar metallicity distribution in the halo of the model resembles the one of M31. The bulge in the model shows, at least two stellar subpopulations, an early collapse population and a population that formed later out of accreted disk mass. In the stellar metallicity distribution of the disk, I find a pronounced ‘G-dwarf problem’ which is the result of a pre-enrichment of the disk ISM with metal-rich gas from the bulge. This revised version was published online in September 2006 with corrections to the Cover Date.     

The formation of galaxy discs

Galaxy disc formation must incorporate the multiphase nature of the interstellar medium. The resulting two-phase structure is generated and maintained by gravitational instability and supernova energy input, which yield a source of turbulent viscosity that is able to compete effectively in the protodisc phase with early angular momentum loss of the baryonic component via dynamical friction in the dark halo. Provided that star formation occurs on the viscous drag time-scale, this mechanism provides a means of accounting for disc sizes and radial profiles. The star formation feedback is self-regulated by turbulent gas pressure limited percolation of the supernova remnant heated hot phase, but can run away in gas-rich protodiscs to generate compact starbursts. A simple analytic model is derived for a Schmidt-like global star formation law in terms of the cold gas volume density.     

Galactic Archaeology: The dwarfs that survived and perished

From the archaeological point of view, the local dwarf galaxies are unique objects in which the imprint of the conditions that shaped the early structure formation can be studied today at high resolution. Over the last decade, this new window into the high redshift Universe has started to be exploited using deep wide-field imaging, high resolution spectroscopy and cutting edge N-body and hydro-dynamical simulations. We review the recent advances in the observational studies of the Milky Way dwarf galaxies, with the aim to understand the properties of the population as a whole and to assist an objective comparison between the models and the data.     

Oxygen abundance in local disk and bulge: chemical evolution with a strictly universal IMF

This paper has two parts: one about observational constraints related to the empirical differential oxygen abundance distribution (EDOD), and the other about inhomogeneous models of chemical evolution, in particular the theoretical differential oxygen abundance distribution (TDOD). In the first part, the EDOD is deduced from subsamples related to two different samples involving (i) N=532 solar neighbourhood (SN) stars within the range, −1.5<[Fe/H]<0.5, for which the oxygen abundance has been determined both in presence and in absence of the local thermodynamical equilibrium (LTE) approximation (Ramirez et al. in Astron. Astrophys. 465:271, 2007); and (ii) N=64 SN thick disk, SN thin disk, and bulge K-giant stars within the range, −1.7<[Fe/H]<0.5, for which the oxygen abundance has been determined (Melendez et al. in Astron. Astrophys. 484:L21, 2008). A comparison is made with previous results implying use of [O/H]–[Fe/H] empirical relations (Caimmi in Astron. Nachr. 322:241, 2001b; New Astron. 12:289, 2007) related to (iii) 372 SN halo subdwarfs (Ryan and Norris in Astron. J. 101:1865, 1991); and (iv) 268 K-giant bulge stars (Sadler et al. in Astron. J. 112:171, 1996). The EDOD of the SN thick + thin disk is determined by weighting the mass, for assumed SN thick to thin disk mass ratio within the range, 0.1–0.9. In the second part, inhomogeneous models of chemical evolution for the SN thick disk, the SN thin disk, the SN thick + thin disk, the SN halo, and the bulge, are computed assuming the instantaneous recycling approximation. The EDOD data are fitted, to an acceptable extent, by their TDOD counterparts with the exception of the thin or thick + thin disk, where two additional restrictions are needed: (i) still undetected, low-oxygen abundance thin disk stars exist, and (ii) a single oxygen overabundant star is removed from a thin disk subsample. In any case, the (assumed power-law) stellar initial mass function (IMF) is universal but gas can be inhibited from, or enhanced in, forming stars at different rates with respect to a selected reference case. Models involving a strictly universal IMF (i.e. gas neither inhibited from, nor enhanced in, forming stars with respect to a selected reference case) can also reproduce the data to an acceptable extent. Our main conclusions are (1) different models are necessary to fit the (incomplete) halo sample, which is consistent with the idea of two distinct halo components: an inner, flattened halo in slow prograde rotation, and an outer, spherical halo in net retrograde rotation (Carollo et al. in Nature 450:1020, 2007); (2) the oxygen enrichment within the inner SN halo, the SN thick disk, and the bulge, was similar and coeval within the same metallicity range, as inferred from observations (Prochaska et al. in Astron. J. 120:2513, 2000); (3) the fit to thin disk data implies an oxygen abundance range similar to its thick disk counterpart, with the extension of conclusion (2) to the thin disk, and the evolution of the thick + thin disk as a whole (Haywood in Mon. Not. R. Astron. Soc. 388:1175, 2008) cannot be excluded; (4) leaving outside the outer halo, a fit to the data related to different environments is provided by models with a strictly universal IMF but different probabilities of a region being active, which implies different global efficiencies of the star formation rate; (5) a special case of stellar migration across the disk can be described by models with enhanced star formation, where a fraction of currently observed SN stars were born in situ and a comparable fraction is due to the net effect of stellar migration, according to recent results based on high-resolution N-body + smooth particle hydrodynamics simulations (Roškar et al. in Astrophys. J. Lett. 684:L79, 2008).     

Spherical galaxy models with power-law logarithmic slope

We present a new family of spherically symmetric models for the luminous components of elliptical and spiral galaxies and their dark matter haloes. Our starting point is a general expression for the logarithmic slope  α( r ) = d log ρ/d log  r   from which most of the cuspy models yet available in literature may be derived. We then dedicate our attention to a particular set of models where the logarithmic slope is a power-law function of the radius r investigating in detail their dynamics assuming isotropy in the velocity space. While the basic properties (such as the density profile and the gravitational potential) may be expressed analytically, both the distribution function and the observable quantities (the surface brightness and the line-of-sight velocity dispersion) have to be evaluated numerically. We also consider the extension to anisotropic models, trying two different parametrization. As the model recently proposed by Navarro et al. as the best fit to their sample of numerically simulated haloes belongs to the family we present here, analytical approximations are given for the most useful quantities.     

Milky Way Gas Dynamics

The Milky Way is made up of a central bar, a disk with embedded spiral arms, and a dark matter halo. Observational and theoretical constraints for the characteristic parameters of these components will be presented, with emphasis on the constraints from the dynamics of the Milky Way gas. In particular, the fraction of dark matter inside the solar radius, the location of the main resonances, and the evidence for multiple pattern speeds will be discussed.Invited talk at the AAS Division on Dynamical Astronomy meeting, Santa Barbara, April 2005     

Galaxy secular mass flow rate determination using the potential-density phase shift approach: Application to six nearby spiral galaxies

Using the potential-density phase shift approach developed by the present authors in earlier publications, we estimate the magnitude of radial mass accretion/excretion rates across the disks of six nearby spiral galaxies (NGC 628, NGC 3351, NGC 3627, NGC 4321, NGC 4736, and NGC 5194) having a range of Hubble types. Our goal is to examine these rates in the context of bulge building and secular morphological evolution along the Hubble sequence. Stellar surface density maps of the sample galaxies are derived from SINGS 3.6 μm and SDSS i-band images using colors as an indicator of mass-to-light ratios. Corresponding molecular and atomic gas surface densities are derived from published CO (1-0) and HI interferometric observations of the BIMA SONG, THINGS, and VIVA surveys. The mass flow rate calculations utilize a volume-type torque integral to calculate the angular momentum exchange rate between the basic state disk matter and what we assume to be density wave modes in the observed galaxies. This volume-type integral contains the contributions from both the gravitational surface torque couple and the advective surface torque couple at the nonlinear, quasi-steady state of the wave modes, in sharp contrast to its behavior in the linear regime, where it contains only the contribution from the gravitational surface torque couple used by Lynden-Bell & Kalnajs in 1972. The potential-density phase shift approach yields angular momentum transport rates several times higher than those estimated using the Lynden-Bell and Kalnajs approach. And unlike Lynden-Bell and Kalnajs, whose approach predicts zero mass redistribution across the majority of the disk surface (apart from the isolated locations of wave-particle resonances) for quasi-steady waves, the current approach leads to predictions of significant mass redistribution induced by the quasi-steady density wave modes, enough for the morphological types of disks to evolve substantially within its lifetime. This difference with the earlier conclusions of Lynden-Bell and Kalnajs reflects the dominant role played by collisionless shocks in the secular evolution of galaxies containing extremely non-linear, quasi-steady density wave modes, thus enabling significant morphological transformation along the Hubble sequence during a Hubble time. We show for the first time also, using observational data, that stellar mass accretion/excretion is just as important, and oftentimes much more important, than the corresponding accretion/excretion processes in the gaseous component, with the latter being what had been emphasized in most of the previous secular evolution studies.     

Galactic satellite systems: radial distribution and environment dependence of galaxy morphology

We have studied the radial distribution of the early (E/S0) and late (S/Irr) types of satellites around bright host galaxies. We made a volume-limited sample of 4986 satellites brighter than   M _r=−18.0  associated with 2254 hosts brighter than   M _r=−19.0  from the Sloan Digital Sky Survey Data Release 5 sample. The morphology of satellites is determined by an automated morphology classifier, but the host galaxies are visually classified. We found segregation of satellite morphology as a function of the projected distance from the host galaxy. The amplitude and shape of the early-type satellite fraction profile are found to depend on the host luminosity. This is the morphology–radius/density relation at the galactic scale. There is a strong tendency for morphology conformity between the host galaxy and its satellites. The early-type fraction of satellites hosted by early-type galaxies is systematically larger than that of late-type hosts, and is a strong function of the distance from the host galaxies. Fainter satellites are more vulnerable to the morphology transformation effects of hosts. Dependence of satellite morphology on the large-scale background density was detected. The fraction of early-type satellites increases in high-density regions for both early- and late-type hosts. It is argued that the conformity in morphology of galactic satellite system is mainly originated by the hydrodynamical and radiative effects of hosts on satellites.     

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.     

The disc-dominated host galaxy of FR-I radio source B2 0722+30

We present new observational results that conclude that the nearby radio galaxy B2 0722+30 is one of the very few known disc galaxies in the low-redshift Universe that host a classical double-lobed radio source. In this paper, we use H  i observations, deep optical imaging, stellar population synthesis modelling and emission-line diagnostics to study the host galaxy, classify the active galactic nucleus (AGN) and investigate environmental properties under which a radio-loud AGN can occur in this system. Typical for spiral galaxies, B2 0722+30 has a regularly rotating gaseous disc throughout which star formation occurs. Dust heating by the ongoing star formation is likely responsible for the high infrared luminosity of the system. The optical emission-line properties of the central region identify a Low Ionization Nuclear Emission-line Region (LINER)-type nucleus with a relatively low [O  iii ] luminosity, in particular when compared with the total power of the Fanaroff & Riley type-I radio source that is present in this system. This classifies B2 0722+30 as a classical radio galaxy rather than a typical Seyfert galaxy. The environment of B2 0722+30 is extremely H  i -rich, with several nearby interacting galaxies. We argue that a gas-rich interaction involving B2 0722+30 is a likely cause for the triggering of the radio AGN and/or the fact that the radio source managed to escape the optical boundaries of the host galaxy.     

The evolution of the galaxy red sequence in simulated clusters and groups

N -body/hydrodynamical simulations of the formation and evolution of galaxy groups and clusters in a Λ cold dark matter (ΛCDM) cosmology are used in order to follow the building-up of the colour–magnitude relation in two clusters and in 12 groups. We have found that galaxies, starting from the more massive, move to the red sequence (RS) as they get aged over times and eventually set upon a 'dead sequence' (DS) once they have stopped their bulk star formation activity. Fainter galaxies keep having significant star formation out to very recent epochs and lie broader around the RS. Environment plays a role as galaxies in groups and cluster outskirts hold star formation activity longer than the central cluster regions. However, galaxies experiencing infall from the outskirts to the central parts keep star formation on until they settle on to the DS of the core galaxies. Merging contributes to mass assembly until z ∼ 1, after which major events only involve the brightest cluster galaxies.
The emerging scenario is that the evolution of the colour–magnitude properties of galaxies within the hierarchical framework is mainly driven by star formation activity during dark matter haloes assembly. Galaxies progressively quenching their star formation settle to a very sharp 'red and dead' sequence, which turns out to be universal, its slope and scatter being almost independent of the redshift (since at least z ∼ 1.5) and environment.
Differently from the DS, the operatively defined RS evolves more evidently with z , the epoch when it changes its slope being closely corresponding to that at which the passive galaxies population takes over the star-forming one: this goes from z ≃ 1 in clusters down to 0.4 in normal groups.     

The mass of the Andromeda galaxy

This paper argues that the Milky Way galaxy is probably the largest member of the Local Group. The evidence comes from estimates of the total mass of the Andromeda galaxy (M31) derived from the three-dimensional positions and radial velocities of its satellite galaxies, as well as the projected positions and radial velocities of its distant globular clusters and planetary nebulae. The available data set comprises 10 satellite galaxies, 17 distant globular clusters and nine halo planetary nebulae with radial velocities. We find that the halo of Andromeda has a mass of together with a scalelength of 90 kpc and a predominantly isotropic velocity distribution. For comparison, our earlier estimate for the Milky Way halo is Although the error bars are admittedly large, this suggests that the total mass of M31 is probably less than that of the Milky Way . We verify the robustness of our results to changes in the modelling assumptions and to errors caused by the small size and incompleteness of the data set.
Our surprising claim can be checked in several ways in the near future. The numbers of satellite galaxies, planetary nebulae and globular clusters with radial velocities can be increased by ground-based spectroscopy, while the proper motions of the companion galaxies and the unresolved cores of the globular clusters can be measured using the astrometric satellites Space Interferometry Mission ( SIM ) and Global Astrometric Interferometer for Astrophysics ( GAIA ). Using 100 globular clusters at projected radii 20 R 50 kpc with both radial velocities and proper motions, it will be possible to estimate the mass within 50 kpc to an accuracy of 20 per cent. Measuring the proper motions of the companion galaxies with SIM and GAIA will reduce the uncertainty in the total mass caused by the small size of the data set to 22 per cent.