The collision between the Milky Way and Andromeda

We use an N -body/hydrodynamic simulation to forecast the future encounter between the Milky Way and the Andromeda galaxies, given present observational constraints on their relative distance, relative velocity, and masses. Allowing for a comparable amount of diffuse mass to fill the volume of the Local Group, we find that the two galaxies are likely to collide in a few billion years – within the Sun's lifetime. During the interaction, there is a chance that the Sun will be pulled away from its present orbital radius and reside in an extended tidal tail. The likelihood for this outcome increases as the merger progresses, and there is a remote possibility that our Sun will be more tightly bound to Andromeda than to the Milky Way before the final merger. Eventually, after the merger has completed, the Sun is most likely to be scattered to the outer halo and reside at much larger radii (>30 kpc). The density profiles of the stars, gas and dark matter in the merger product resemble those of elliptical galaxies. Our Local Group model therefore provides a prototype progenitor of late-forming elliptical galaxies.     

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.     

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.     

Milky Way type galaxies in a ΛCDM cosmology

We analyse a sample of 52 000 Milky Way (MW) type galaxies drawn from the publicly available galaxy catalogue of the Millennium Simulation with the aim of studying statistically the differences and similarities of their properties in comparison to our Galaxy. Model galaxies are chosen to lie in haloes with maximum circular velocities in the range 200–250 km s⁻¹ and to have bulge-to-disc ratios similar to that of the MW. We find that model MW galaxies formed 'quietly' through the accretion of cold gas and small satellite systems. Only ≈12 per cent of our model galaxies experienced a major merger during their lifetime. Most of the stars formed ' in situ ', with only about 15 per cent of the final mass gathered through accretion. Supernovae (SNe) and active galactic nuclei (AGN) feedback play an important role in the evolution of these systems. At high redshifts, when the potential wells of the MW progenitors are shallower, winds driven by SNe explosions blow out a large fraction of the gas and metals. As the systems grow in mass, SNe feedback effects decrease and AGN feedback takes over, playing a more important role in the regulation of the star formation activity at lower redshifts. Although model MW galaxies have been selected to lie in a narrow range of maximum circular velocities, they nevertheless exhibit a significant dispersion in the final stellar masses and metallicities. Our analysis suggests that this dispersion results from the different accretion histories of the parent dark matter haloes. Statistically, we also find evidences to support the MW as a typical Sb/Sc galaxy in the same mass range, providing a suitable benchmark to constrain numerical models of galaxy formation.     

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.     

A dimensional study of disc galaxies

We present a highly simplified model of the dynamical structure of a disc galaxy where only two parameters fully determine the solution, mass and angular momentum. We show through simple physical scalings that once the mass has been fixed, the angular momentum parameter λ is expected to regulate such critical galactic disc properties as colour, thickness of the disc and bulge-to-disc ratio. It is, hence, expected to be the determinant physical ingredient resulting in a given Hubble type. A simple analytic estimate of λ for an observed system is provided. An explicit comparison of the distribution of several galactic parameters against both Hubble type and λ is performed using observed galaxies. Both such distributions exhibit highly similar characteristics for all galactic properties studied, suggesting λ as a physically motivated classification parameter for disc galaxies.     

Chemical evolution of galaxies – I. A composition-dependent SPH model for chemical evolution and cooling

We describe an smooth particle hydrodynamics (SPH) model for chemical enrichment and radiative cooling in cosmological simulations of structure formation. This model includes: (i) the delayed gas restitution from stars by means of a probabilistic approach designed to reduce the statistical noise and, hence, to allow for the study of the inner chemical structure of objects with moderately high numbers of particles; (ii) the full dependence of metal production on the detailed chemical composition of stellar particles by using, for the first time in SPH codes, the   Q _ij   matrix formalism that relates each nucleosynthetic product to its sources and (iii) the full dependence of radiative cooling on the detailed chemical composition of gas particles, achieved through a fast algorithm using a new metallicity parameter ζ( T ) that gives the weight of each element on the total cooling function. The resolution effects and the results obtained from this SPH chemical model have been tested by comparing its predictions in different problems with known theoretical solutions. We also present some preliminary results on the chemical properties of elliptical galaxies found in self-consistent cosmological simulations. Such simulations show that the above ζ-cooling method is important to prevent an overestimation of the metallicity-dependent cooling rate, whereas the   Q _ij   formalism is important to prevent a significant underestimation of the [α/Fe] ratio in simulated galaxy-like objects.     

266 E+A galaxies selected from the Sloan Digital Sky Survey Data Release 2: the origin of E+A galaxies

E+A galaxies are characterized as galaxies with strong Balmer absorption lines but without any [O  ii ] or Hα emission lines. The existence of strong Balmer absorption lines indicates that E+A galaxies have experienced starburst within the past one gigayear. However, the lack of [O  ii ] and Hα emission lines indicates that E+A galaxies do not have any on-going star formation. Therefore, E+A galaxies are interpreted as post-starburst galaxies. For many years, however, it has been a mystery why E+A galaxies started starburst and why they quenched star formation abruptly. Using one of the largest samples of 266 E+A galaxies carefully selected from the Sloan Digital Sky Survey Data Release 2, we have investigated the environment of E+A galaxies from 50 kpc to 8 Mpc scale, i.e. from a typical distance to satellite galaxies to the scale of large-scale structures. We found that E+A galaxies have an excess of local galaxy density only at a scale of <100 kpc (with a 2σ significance), but not at the cluster scale (∼1.5 Mpc) nor at the scale of large-scale structure (∼8 Mpc). These results indicate that E+A galaxies are not created by the physical mechanisms associated with galaxy clusters or the large-scale structure, but are likely to be created by dynamical interaction with closely accompanying galaxies at a <100 kpc scale. The claim is also supported by the morphology of E+A galaxies. We have found that almost all E+A galaxies have a bright compact core, and that ∼30 per cent of E+A galaxies have dynamically disturbed signatures or tidal tails, which quite strongly suggest the morphological appearance of merger/interaction remnants.     

Forming stars on a viscous time-scale: the key to exponential stellar profiles in disc galaxies?

We argue for implementing star formation on a viscous time-scale in hydrodynamical simulations of disc galaxy formation and evolution. Modelling two-dimensional isolated disc galaxies with the Bhatnagar–Gross–Krook (BGK) hydrocode, we verify the analytic claim of various authors that if the characteristic time-scale for star formation is equal to the viscous time-scale in discs, the resulting stellar profile is exponential on several scalelengths whatever the initial gas and dark matter profile. This casts new light on both numerical and semi-analytical disc formation simulations that either (a) commence star formation in an already exponential gaseous disc, (b) begin a disc simulation with conditions known to lead to an exponential, i.e. the collapse of a spherically symmetric nearly uniform sphere of gas in solid-body rotation under the assumption of specific angular momentum conservation, or (c) in simulations performed in a hierarchical context, tune their feedback processes to delay disc formation until the dark matter haloes are slowly evolving and without much substructure so that the gas has the chance to collapse under conditions known to give exponentials. In such models, star formation follows a Schmidt-like law, which for lack of a suitable time-scale, resorts to an efficiency parameter. With star formation prescribed on a viscous time-scale, however, we find gas and star fractions after ∼12 Gyr that are consistent with observations without having to invoke a 'fudge factor' for star formation. Our results strongly suggest that despite our gap in understanding the exact link between star formation and viscosity, the viscous time-scale is indeed the natural time-scale for star formation.     

The evolution of 3CR radio galaxies from z1

We present the results of a comprehensive re-analysis of the images of a virtually complete sample of 28 powerful 3CR radio galaxies with redshifts 0.6< z <1.8 from the Hubble Space Telescope ( HST ) archive. Using a two-dimensional modelling technique we have derived scalelengths and absolute magnitudes for a total of 16 3CR galaxies with a median redshift of z =0.8. Our results confirm the basic conclusions of Best, Longair & Röttgering in that we also find z =1 3CR galaxies to be massive, well-evolved ellipticals, the infrared emission of which is dominated by starlight. However, we in fact find that the scalelength distribution of 3CR galaxies at z ≃1 is completely indistinguishable from that derived for their low-redshift counterparts from our own recently completed HST study of active galactic nuclei hosts at z ≃0.2. There is thus no evidence that 3CR radio galaxies at z ≃1 are dynamically different from 3CR galaxies at low redshift. Moreover, for a 10-object subsample we have determined the galaxy parameters with sufficient accuracy to demonstrate, for the first time, that the z ≃1 3CR galaxies follow a Kormendy relation that is indistinguishable from that displayed by low-redshift ellipticals if one allows for purely passive evolution. The implied rather modest level of passive evolution since z ≃1 is consistent with that predicted from spectrophotometric models provided one assumes a high formation redshift ( z ≥4) within a low-density universe. We conclude that there is no convincing evidence for significant dynamical evolution among 3CR galaxies in the redshift interval 0< z <1, and that simple passive evolution remains an acceptable interpretation of the K – z relation for powerful radio galaxies.