Forming a large disc galaxy from a z < 1 major merger

Using high-resolution SPH simulations in a fully cosmological Λ cold dark matter context, we study the formation of a bright disc-dominated galaxy that originates from a 'wet' major merger at   z = 0.8  . The progenitors of the disc galaxy are themselves disc galaxies that formed from early major mergers between galaxies with blue colours. A substantial thin stellar disc grows rapidly following the last major merger and the present-day properties of the final remnant are typical of early-type spiral galaxies, with an i -band bulge-to-disc ratio ∼0.65, a disc scalelength of 7.2 kpc,   g − r = 0.5 mag  , an H  i linewidth ( W ₂₀/2) of 238 km s⁻¹ and total magnitude   i =−22.4  . The key ingredients for the formation of a dominant stellar disc component after a major merger are (i) substantial and rapid accretion of gas through cold flows followed at late times by cooling of gas from the hot phase, (ii) supernova feedback that is able to partially suppress star formation during mergers and (iii) relative fading of the spheroidal component. The gas fraction of the progenitors' discs does not exceed 25 per cent at   z < 3  , emphasizing that the continuous supply of gas from the local environment plays a major role in the regrowth of discs and in keeping the galaxies blue. The results of this simulation alleviate the problem posed for the existence of disc galaxies by the high likelihood of interactions and mergers for galaxy-sized haloes at relatively low z .     

Modelling star formation and feedback in simulations of interacting galaxies

We discuss a heuristic model to implement star formation and feedback in hydrodynamical simulations of galaxy formation and evolution. In this model, gas is allowed to cool radiatively and to form stars at a rate given by a simple Schmidt-type law. We assume that supernova feedback results in turbulent motions of gas below resolved scales, a process that can pressurize the diffuse gaseous medium effectively, even if it lacks substantial thermal support. Ignoring the complicated detailed physics of the feedback processes, we try to describe their net effect on the interstellar medium with a fiducial second reservoir of internal energy, which accounts for the kinetic energy content of the gas on unresolved scales. Applying the model to three-dimensional, fully self-consistent models of isolated disc galaxies, we show that the resulting feedback loop can be modelled with smoothed particle hydrodynamics such that converged results can be reached with moderate numerical resolution. With an appropriate choice of the free parameters, Kennicutt's phenomenological star formation law can be reproduced over many orders of magnitude in gas surface density. We also apply the model to mergers of equal-mass disc galaxies, typically resulting in strong nuclear starbursts. Confirming previous findings, the presence of a bulge can delay the onset of the starburst from the first encounter of the galaxies until their final coalescence. The final density profiles of the merger remnants are consistent with de Vaucouleurs profiles, except for the innermost region, where the newly created stars give rise to a luminous core with stellar densities that may be in excess of those observed in the cores of most elliptical galaxies. By comparing the isophotal shapes of collisionless and dissipative merger simulations we show that dissipation leads to isophotes that are more discy than those of corresponding collisionless simulations.     

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.     

Pre-enriched, not primordial ellipticals

We follow the chemical evolution of a galaxy through star formation and its feedback into the interstellar medium (ISM), starting from primordial gas and allowing for gas to inflow into the region being modelled. We attempt to reproduce observed spectral line strengths for early-type galaxies in order to constrain their star formation histories (SFH). The efficiencies and times of star formation are varied, as are the amount and duration of inflow. We evaluate the chemical enrichment and the mass of stars made with time. Single stellar population (SSP) data are then used to predict line strengths for composite stellar populations. The results are compared with observed line strengths in 10 ellipticals, including some features which help to break the problem of age–metallicity degeneracy in old stellar populations. We find that the elliptical galaxies modelled require high metallicity SSPs (> 3 Z⊙) at later times. In addition, the strong lines observed cannot be produced by an initial starburst in primordial gas, even if a large amount of inflow is allowed for during the first few × 10⁸ yr. This is because some pre-enrichment is required for lines in the bulk of the stars to approach the observed line strengths in ellipticals. These strong lines are better modelled by a system with a delayed burst of star formation, following an early SFH which can be a burst or more steady star formation. Such a model is representative of star formation in normal ellipticals or spirals, respectively, followed by a starburst and gas inflow during a merger or strong interaction with a gas-rich galaxy. Alternatively, a single initial burst of normal stars with a Salpeter initial mass function could produce the observed strong lines if it followed some pre-enrichment process which did not form long-lived stars (e.g. population III stars).     

Baryonic conversion tree: the global assembly of stars and dark matter in galaxies from the Sloan Digital Sky Survey

Using the spectroscopic sample of the Sloan Digital Sky Survey Data Release 1 (SDSS DR1), we measure how gas was transformed into stars as a function of time and stellar mass: the baryonic conversion tree (BCT). There is a clear correlation between early star formation activity and present-day stellar mass: the more massive galaxies have formed approximately 80 per cent of their stars at   z > 1  , while for the less massive ones the value is only approximately 20 per cent. By comparing the BCT with the dark matter merger tree, we find indications that star formation efficiency at   z > 1  had to be approximately a factor of two higher than today (∼10 per cent) in galaxies with present-day stellar mass larger than  2 × 10¹¹ M_⊙  , if this early star formation occurred in the main progenitor. Therefore, the λ cold dark matter (LCDM) paradigm can accommodate a large number of red objects. On the other hand, in galaxies with present-day stellar mass less than  10¹¹ M_⊙  , efficient star formation seems to have been triggered at   z ∼ 0.2  . We show that there is a characteristic mass  ( M _*∼ 10¹⁰ M_⊙)  for feedback efficiency (or lack of star formation). For galaxies with masses lower than this, feedback (or star formation suppression) is very efficient while for higher masses it is not. The BCT, determined here for the first time, should be an important observable with which to confront theoretical models of galaxy formation.     

The growth and assembly of a massive galaxy at z∼ 2

We study the stellar mass assembly of the Spiderweb galaxy  (MRC 1138−262)  , a massive   z = 2.2  radio galaxy in a protocluster and the probable progenitor of a brightest cluster galaxy. Nearby protocluster galaxies are identified and their properties are determined by fitting stellar population models to their rest-frame ultraviolet to optical spectral energy distributions. We find that within 150 kpc of the radio galaxy the stellar mass is centrally concentrated in the radio galaxy, yet most of the dust-uncorrected, instantaneous star formation occurs in the surrounding low-mass satellite galaxies. We predict that most of the galaxies within 150 kpc of the radio galaxy will merge with the central radio galaxy by   z = 0  , increasing its stellar mass by up to a factor of ≃2. However, it will take several hundred Myr for the first mergers to occur, by which time the large star formation rates are likely to have exhausted the gas reservoirs in the satellite galaxies. The tidal radii of the satellite galaxies are small, suggesting that stars and gas are being stripped and deposited at distances of tens of kpc from the central radio galaxy. These stripped stars may become intracluster stars or form an extended stellar halo around the radio galaxy, such as those observed around cD galaxies in cluster cores.     

Dissipationless mergers of elliptical galaxies and the evolution of the fundamental plane

We carry out numerical simulations of dissipationless major mergers of elliptical galaxies using initial galaxy models that consist of a dark matter haloes and a stellar bulge with properties consistent with the observed fundamental plane. By varying the density profile of the dark matter haloes [standard Navarro, Frenk & White (NFW) profile versus adiabatically contracted NFW profile], the global stellar to dark matter mass ratio and the orbit of the merging galaxies, we are able to assess the impact of each of these factors on the structure of the merger remnant. Our results indicate that the properties of the remnant bulge depend primarily on the angular momentum and energy of the orbit; for a cosmologically motivated orbit, the effective radius and velocity dispersion of the remnant bulge remain approximately on the fundamental plane. This indicates that the observed properties of elliptical galaxies are consistent with significant growth via late dissipationless mergers. We also find that the dark matter fraction within the effective radius of our remnants increases after the merger, consistent with the hypothesis that the tilt of the fundamental plane from the virial theorem is due to a varying dark matter fraction as a function of galaxy 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.     

The stellar populations of spiral galaxies

We have used a large sample of low-inclination spiral galaxies with radially resolved optical and near-infrared photometry to investigate trends in star formation history with radius as a function of galaxy structural parameters. A maximum-likelihood method was used to match all the available photometry of our sample to the colours predicted by stellar population synthesis models. The use of simplistic star formation histories, uncertainties in the stellar population models and considering the importance of dust all compromise the absolute ages and metallicities derived in this work; however, our conclusions are robust in a relative sense. We find that most spiral galaxies have stellar population gradients, in the sense that their inner regions are older and more metal rich than their outer regions. Our main conclusion is that the surface density of a galaxy drives its star formation history, perhaps through a local density dependence in the star formation law. The mass of a galaxy is a less important parameter; the age of a galaxy is relatively unaffected by its mass; however, the metallicity of galaxies depends on both surface density and mass. This suggests that galaxy‐mass-dependent feedback is an important process in the chemical evolution of galaxies. In addition, there is significant cosmic scatter suggesting that mass and density may not be the only parameters affecting the star formation history of a galaxy.     

Predicting the properties of the remnants of dissipative galaxy mergers

We construct a physically motivated model for predicting the properties of the remnants of gaseous galaxy mergers, given the properties of the progenitors and the orbit. The model is calibrated using a large suite of smoothed particle hydrodynamics (SPH) merger simulations. It implements generalized energy conservation while accounting for dissipative energy losses and star formation. The dissipative effects are evaluated from the initial gas fractions and from the orbital parameters via an 'impulse' parameter, which characterizes the strength of the encounter. Given the progenitor properties, the model predicts the remnant stellar mass, half-mass radius and velocity dispersion to an accuracy of 25 per cent. The model is valid for both major and minor mergers. We provide an explicit recipe for semi-analytic models of galaxy formation.     

Galactic winds and circulation of the interstellar medium in dwarf galaxies

We study, through 2D hydrodynamical simulations, the feedback of a starburst on the ISM of typical gas-rich dwarf galaxies. The main goal is to address the circulation of the ISM and metals following the starburst. We assume a single-phase rotating ISM in equilibrium in the galactic potential generated by a stellar disc and a spherical dark halo. The starburst is assumed to occur in a small volume in the centre of the galaxy, and it generates a mechanical power of 3.8×10³⁹ or 3.8×10⁴⁰ erg s⁻¹ for 30 Myr. We find, in accordance with previous investigations, that the galactic wind is not very effective in removing the ISM. The metal-rich stellar ejecta, however, can be efficiently expelled from the galaxy and dispersed in the intergalactic medium.
Moreover, we find that the central region of the galaxy is always replenished with cold and dense gas a few 100 million years after the starburst, achieving the requisite for a new star formation event in ≈0.5–1 Gyr. The hydrodynamical evolution of galactic winds is thus consistent with the episodic star formation regime suggested by many chemical evolution studies.
We also discuss the X-ray emission of these galaxies and find that the observable (emission-averaged) abundance of the hot gas underestimates the real one if thermal conduction is effective. This could explain the very low hot-gas metallicities estimated in starburst galaxies.     

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.     

Feedback and the formation of dwarf galaxy stellar haloes

Stellar population studies show that low-mass galaxies in all environments exhibit stellar haloes that are older and more spherically distributed than the main body of the galaxy. In some cases, there is a significant intermediate age component that extends beyond the young disc. We examine a suite of Smoothed Particle Hydrodynamic simulations and find that elevated early star formation activity combined with supernova feedback can produce an extended stellar distribution that resembles these haloes for model galaxies ranging from   v ₂₀₀= 15  to 35 km s⁻¹, without the need for accretion of subhaloes.     

Numerical simulations of hot halo gas in galaxy mergers

Galaxy merger simulations have explored the behaviour of gas within the galactic disc, yet the dynamics of hot gas within the galaxy halo have been neglected. We report on the results of high-resolution hydrodynamic simulations of colliding galaxies with metal-free hot halo gas. To isolate the effect of the halo gas, we simulate only the dark matter halo and the hot halo gas over a range of mass ratios, gas fractions and orbital configurations to constrain the shocks and gas dynamics within the progenitor haloes. We find that (i) a strong shock is produced in the galaxy haloes before the first passage, increasing the temperature of the gas by almost an order of magnitude to   T ∼ 10^6.3 K  . (ii) The X-ray luminosity of the shock is strongly dependent on the gas fraction; it is  ≳10³⁹ erg s⁻¹  for halo gas fractions larger than 10 per cent. (iii) The hot diffuse gas in the simulation produces X-ray luminosities as large as  10⁴² erg s⁻¹  . This contributes to the total X-ray background in the Universe. (iv) We find an analytic fit to the maximum X-ray luminosity of the shock as a function of merger parameters. This fit can be used in semi-analytic recipes of galaxy formation to estimate the total X-ray emission from shocks in merging galaxies. (v) ∼10–20 per cent of the initial gas mass is unbound from the galaxies for equal-mass mergers, while 3–5 per cent of the gas mass is released for the 3:1 and 10:1 mergers. This unbound gas ends up far from the galaxy and can be a feasible mechanism to enrich the intergalactic medium with metals.     

How galaxies lose their angular momentum

The processes are investigated by which gas loses its angular momentum during the protogalactic collapse phase, leading to disc galaxies that are too compact with respect to the observations. High-resolution N -body/SPH simulations in a cosmological context are presented including cold gas and dark matter (DM). A halo with quiet merging activity since redshift   z ∼ 3.8  and with a high-spin parameter is analysed that should be an ideal candidate for the formation of an extended galactic disc. We show that the gas and the DM have similar specific angular momenta until a merger event occurs at   z ∼ 2  with a mass ratio of 5:1. All the gas involved in the merger loses a substantial fraction of its specific angular momentum due to tidal torques and dynamical friction processes falls quickly into the centre. In contrast, gas infall through small subclumps or accretion does not lead to catastrophic angular momentum loss. In fact, a new extended disc begins to form from gas that was not involved in the 5:1 merger event and that falls in subsequently. We argue that the angular momentum problem of disc galaxy formation is a merger problem: in cold dark matter cosmology substantial mergers with mass ratios of 1:1 to 6:1 are expected to occur in almost all galaxies. We suggest that energetic feedback processes could in principle solve this problem, however only if the heating occurs at the time or shortly before the last substantial merger event. Good candidates for such a coordinated feedback would be a merger-triggered starburst or central black hole heating. If a large fraction of the low angular momentum gas would be ejected, late-type galaxies could form with a dominant extended disc component, resulting from late infall, a small bulge-to-disc ratio and a low baryon fraction, in agreement with observations.     

Simulating galactic outflows with kinetic supernova feedback

Feedback from star formation is thought to play a key role in the formation and evolution of galaxies, but its implementation in cosmological simulations is currently hampered by a lack of numerical resolution. We present and test a subgrid recipe to model feedback from massive stars in cosmological smoothed particle hydrodynamics simulations. The energy is distributed in kinetic form among the gas particles surrounding recently formed stars. The impact of the feedback is studied using a suite of high-resolution simulations of isolated disc galaxies embedded in dark haloes with total mass 10¹⁰ and  10¹²  h ⁻¹ M_⊙  . We focus, in particular, on the effect of pressure forces on wind particles within the disc, which we turn off temporarily in some of our runs to mimic a recipe that has been widely used in the literature. We find that this popular recipe gives dramatically different results because (ram) pressure forces on expanding superbubbles determine both the structure of the disc and the development of large-scale outflows. Pressure forces exerted by expanding superbubbles puff up the disc, giving the dwarf galaxy an irregular morphology and creating a galactic fountain in the massive galaxy. Hydrodynamic drag within the disc results in a strong increase in the effective mass loading of the wind for the dwarf galaxy, but quenches much of the outflow in the case of the high-mass galaxy.     

Dynamical friction and galaxy merging time-scales

The time-scale for galaxies within merging dark matter haloes to merge with each other is an important ingredient in galaxy formation models. Accurate estimates of merging time-scales are required for predictions of astrophysical quantities such as black hole binary merger rates, the build-up of stellar mass in central galaxies and the statistical properties of satellite galaxies within dark matter haloes. In this paper, we study the merging time-scales of extended dark matter haloes using N -body simulations. We compare these results to standard estimates based on the Chandrasekhar theory of dynamical friction. We find that these standard predictions for merging time-scales, which are often used in semi-analytic galaxy formation models, are systematically shorter than those found in simulations. The discrepancy is approximately a factor of 1.7 for M _sat/ M _host≈ 0.1 and becomes larger for more disparate satellite-to-host mass ratios, reaching a factor of ∼3.3 for M _sat/ M _host≈ 0.01. Based on our simulations, we propose a new, easily implementable fitting formula that accurately predicts the time-scale for an extended satellite to sink from the virial radius of a host halo down to the halo's centre for a wide range of M _sat/ M _host and orbits. Including a central bulge in each galaxy changes the merging time-scale by ≲10 per cent. To highlight one concrete application of our results, we show that merging time-scales often used in the literature overestimate the growth of stellar mass by satellite accretion by ≈40 per cent, with the extra mass gained in low mass ratio mergers.     

The flattening and the orbital structure of early-type galaxies and collisionless N-body binary disc mergers

We use oblate axisymmetric dynamical models including dark haloes to determine the orbital structure of intermediate mass to massive early-type galaxies in the Coma galaxy cluster. We find a large variety of orbital compositions. Averaged over all sample galaxies the unordered stellar kinetic energy in the azimuthal and the radial direction are of the same order, but they can differ by up to 40 per cent in individual systems. In contrast, both for rotating and non-rotating galaxies the vertical kinetic energy is on average smaller than in the other two directions. This implies that even most of the rotating ellipticals are flattened by an anisotropy in the stellar velocity dispersions. Using three-integral axisymmetric toy models, we show that flattening by stellar anisotropy maximizes the entropy for a given density distribution. Collisionless disc merger remnants are radially anisotropic. The apparent lack of strong radial anisotropy in observed early-type galaxies implies that they may not have formed from mergers of discs unless the influence of dissipational processes was significant.     

Dry mergers: a crucial test for galaxy formation

We investigate the role that dry mergers play in the build-up of massive galaxies within the cold dark matter paradigm. Implementing an empirical shut-off mass scale for star formation, we find a nearly constant dry merger rate of  ∼6 × 10⁻⁵ Mpc⁻³ Gyr⁻¹  at   z ≤ 1  and a steep decline at larger z . Less than half of these mergers are between two galaxies that are morphologically classified as early-types, and the other half is mostly between an early- and late-type galaxy. Latter are prime candidates for the origin of tidal features around red elliptical galaxies. The introduction of a transition mass scale for star formation has a strong impact on the evolution of galaxies, allowing them to grow above a characteristic mass scale of   M _{*, c} ∼ 6.3 × 10¹⁰ M_⊙  by mergers only. As a consequence of this transition, we find that around   M _{*, c}   , the fraction of 1:1 mergers is enhanced with respect to unequal mass major mergers. This suggests that it is possible to detect the existence of a transition mass scale by measuring the relative contribution of equal mass mergers to unequal mass mergers as a function of galaxy mass. The evolution of the high-mass end of the luminosity function is mainly driven by dry mergers at low z . We however find that only 10–20 per cent of galaxies more massive than   M _{*, c}   experience dry major mergers within their last Gyr at any given redshift   z ≤ 1  .