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 .     

Large-scale galactic turbulence: can self-gravity drive the observed H i velocity dispersions?

Observations of turbulent velocity dispersions in the H  i component of galactic discs show a characteristic floor in galaxies with low star formation rates and within individual galaxies the dispersion profiles decline with radius. We carry out several high-resolution adaptive mesh simulations of gaseous discs embedded within dark matter haloes to explore the roles of cooling, star formation, feedback, shearing motions and baryon fraction in driving turbulent motions. In all simulations the disc slowly cools until gravitational and thermal instabilities give rise to a multiphase medium in which a large population of dense self-gravitating cold clouds are embedded within a warm gaseous phase that forms through shock heating. The diffuse gas is highly turbulent and is an outcome of large-scale driving of global non-axisymmetric modes as well as cloud–cloud tidal interactions and merging. At low star formation rates these processes alone can explain the observed H  i velocity dispersion profiles and the characteristic value of  ∼10 km s⁻¹  observed within a wide range of disc galaxies. Supernovae feedback creates a significant hot gaseous phase and is an important driver of turbulence in galaxies with a star formation rate per unit area  ≳10⁻³ M_⊙ yr⁻¹ kpc⁻²  .     

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.     

Simulations of the formation of stellar discs in the Galactic Centre via cloud–cloud collisions

Young massive stars in the central parsec of our Galaxy are best explained by star formation within at least one, and possibly two, massive self-gravitating gaseous discs. With help of numerical simulations, we here consider whether the observed population of young stars could have originated from a large angle collision of two massive gaseous clouds at   R ≃ 1 pc  from Sgr A*. In all the simulations performed, the post-collision gas flow forms an inner, nearly circular gaseous disc and one or two eccentric outer filaments, consistent with the observations. Furthermore, the radial stellar mass distribution is always very steep,  Σ_*∝ R ⁻²  , again consistent with the observations. All of our simulations produce discs that are warped by between 30° and 60°, in accordance with the most recent observations. The three-dimensional velocity structure of the stellar distribution is sensitive to initial conditions (e.g. the impact parameter of the clouds) and gas cooling details. For example, the runs in which the inner disc is fed intermittently with material possessing fluctuating angular momentum result in multiple stellar discs with different orbital orientations, contradicting the observed data. In all the cases the amount of gas accreted by our inner boundary condition is large, enough to allow Sgr A* to radiate near its Eddington limit over ∼10⁵ yr. This suggests that a refined model would have physically larger clouds (or a cloud and a disc such as the circumnuclear disc) colliding at a distance of a few parsecs rather than 1 pc as in our simulations.     

H i and dark matter in the windy starburst dwarf galaxy NGC 1705

We present 21-cm H  i line observations of the blue compact dwarf galaxy NGC 1705. Previous optical observations show a strong outflow powered by an ongoing starburst dominating the H  ii morphology and kinematics. In contrast, most of the H  i lies in a rotating disc. An extraplanar H  i spur accounts for ∼8 per cent of the total H  i mass, and is possibly associated with the H  ii outflow. The inferred mass loss rate out of the core of the galaxy is significant, ∼0.2 − 2 M_⊙ yr⁻¹, but does not dominate the H  i dynamics. Mass model fits to the rotation curve show that the dark matter (DM) halo is dominant at nearly all radii and has a central density ρ₀ ≈ 0.1 M_⊙ pc⁻³: ten times higher than typically found in dwarf irregular galaxies, but similar to the only other mass-modelled blue compact dwarf, NGC 2915. This large difference strongly indicates that there is little evolution between dwarf irregular and blue compact dwarf types. Instead, dominant DM haloes may regulate the morphology of dwarf galaxies by setting the critical surface density for disc star formation. Neither our data nor catalogue searches reveal any likely external trigger to the starburst in NGC 1705.     

3D spectroscopy of merger Seyfert galaxy Mrk 334: nuclear starburst, superwind and the circumnuclear cavern

We are presenting new results on kinematics and structure of the Mrk 334 Seyfert galaxy. Panoramic (3D) spectroscopy is performed at the 6-m telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences using the integral-field Multi-Pupil Fiber Spectrograph (MPFS) and scanning Fabry–Pérot interferometer. The deep images have revealed that Mrk 334 is observed during the final stage of its merging with a massive companion. A possible mass ratio ranges from 1/5 to 1/3. The merger has triggered mass redistribution in the disc resulting in an intensification of nuclear activity and in a burst of star formation in the inner region of the galaxy. The circumnuclear starburst is so intense that its contribution to the gas ionization exceeds that contribution of the active galactic nuclei (AGN). We interpret the nuclear gas outflow with velocities of  ∼200 km s⁻¹  as a galactic superwind that accompanies the violent star formation. This suggestion is consistent with the asymmetric X-ray brightness distribution in Mrk 334. The trajectory of the fragments of the disrupted satellite in the vicinity of the main galaxy nucleus can be traced. In the galaxy disc, a cavern is found that is filled with a low-density ionized gas. We consider this region to be the place where the remnants of the companion have recently penetrated through the gaseous disc of the main galaxy.     

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.     

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.     

The first galaxies: assembly, cooling and the onset of turbulence

We investigate the properties of the first galaxies at   z ≳ 10  with highly resolved numerical simulations, starting from cosmological initial conditions and taking into account all relevant primordial chemistry and cooling. A first galaxy is characterized by the onset of atomic hydrogen cooling, once the virial temperature exceeds  ≃10⁴ K  , and its ability to retain photoheated gas. We follow the complex accretion and star formation history of a  ≃5 × 10⁷ M_⊙  system by means of a detailed merger tree and derive an upper limit on the number of Population III (Pop III) stars formed prior to its assembly. We investigate the thermal and chemical evolution of infalling gas and find that partial ionization at temperatures  ≳10⁴ K  catalyses the formation of  H₂  and hydrogen deuteride, allowing the gas to cool to the temperature of the cosmic microwave background. Depending on the strength of radiative and chemical feedback, primordial star formation might be dominated by intermediate-mass Pop III stars formed during the assembly of the first galaxies. Accretion on to the nascent galaxy begins with hot accretion, where gas is accreted directly from the intergalactic medium and shock heated to the virial temperature, but is quickly accompanied by a phase of cold accretion, where the gas cools in filaments before flowing into the parent halo with high velocities. The latter drives supersonic turbulence at the centre of the galaxy and could lead to very efficient chemical mixing. The onset of turbulence in the first galaxies thus likely marks the transition to Pop II star 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.     

Star formation and chemical evolution of damped Lyman α systems

In this paper, we investigate the star formation and chemical evolution of damped Lyman α systems (DLAs) based on the disc galaxy formation model developed by Mo, Mao & White. We propose that the DLAs are the central galaxies of less-massive dark haloes present at redshifts z ∼3, and they should inhabit haloes of moderately low circular velocity. The empirical Schmidt law of star formation rates, and closed box model of chemical evolution that an approximation known as instantaneous recycling is assumed, are adopted. In our models, when the predicted distribution of metallicity for DLAs is calculated, two cases are considered. One is that, using the closed-box model, empirical Schmidt law and star formation time, the distribution of metallicity can be directly calculated. The other is that, when the simple gravitational instability of a thin isothermal gas disc as first discussed by Toomre is considered, the star formation occurs only in the region where the surface density of gas satisfies the critical value, not everywhere of a gas disc. In this case, we first obtain the region where the star formation can occur by assuming that the disc has a flat rotation curve and rotational velocity is equal to the circular velocity of the surrounding dark matter halo, and then calculate the metallicity distribution as in case one. We assume that star formation in each DLA lasts for a period of 1 Gyr from redshifts z =3. There is only one output parameter in our models, i.e. the stellar yield, which relates to the time of star formation history and is obtained by normalizing the predicted distribution of metallicity to the mean value of 1/13 Z_⊙ as presented by Pettini et al.. The predicted metallicity distribution is consistent with the current (rather limited) observational data. A random distribution of galactic discs is taken into account.     

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.     

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.     

The extended rotation curve and the dark matter halo of M33

We present the 21-cm rotation curve of the nearby galaxy M33 out to a galactocentric distance of 16 kpc (13 disc scalelengths). The rotation curve keeps rising out to the last measured point and implies a dark halo mass ≳5×10¹⁰ M_⊙. The stellar and gaseous discs provide virtually equal contributions to the galaxy gravitational potential at large galactocentric radii, but no obvious correlation is found between the radial distribution of dark matter and the distribution of stars or gas.
Results of the best fit to the mass distribution in M33 picture a dark halo which controls the gravitational potential from 3 kpc outward, with a matter density which decreases radially as R ^−1.3. The density profile is consistent with the theoretical predictions for structure formation in hierarchical clustering cold dark matter (CDM) models, and favours lower mass concentrations than those expected in the standard cosmogony.     

Diffuse X-ray emission from late-type galaxy haloes

Current theories of galaxy formation predict that spiral galaxies are embedded in a reservoir of hot gas. This gas is able to cool on to the galaxy, replenishing cold gas that is consumed by star formation. Estimates of the X-ray luminosity emitted in the cooling region suggest a bolometric luminosity of the order of 10×10⁴¹ erg s⁻¹ in massive systems. We have used ROSAT PSPC data to search for extended X-ray emission from the haloes of three nearby, massive, late-type galaxies: NGC 2841, 4594 and 5529. We infer 95 per cent upper limits on the bolometric X-ray luminosities of the haloes of NGC 2841, 4594 and 5529 of 0.4, 1.2 and 3.8×10⁴¹ erg s⁻¹ respectively. Thus, the true luminosity lies well below the straightforward theoretical prediction. We discuss this discrepancy and suggest a number of ways in which the theoretical model might be brought into agreement with the observational results. A possible solution is that the gravitational potentials of the dark matter haloes of these galaxies are weaker than assumed in the current model. Alternatively, the present-day accretion may be substantially less than is required on average to build the disc over the Hubble time. Our results are, however, based on only three galaxies, none of which is ideal for this kind of study. A larger data set is required to explore this important problem further.     

The formation of the Hubble sequence of disc galaxies: the effects of early viscous evolution

We investigate a model of disc galaxies whereby viscous evolution of the gaseous disc drives material inwards to form a protobulge. We start from the standard picture of disc formation through the settling of gas into a dark halo potential well, with the disc initially coming into centrifugal equilibrium with detailed conservation of angular momentum. We derive generic analytic solutions for the disc–halo system after adiabatic compression of the dark halo, with free choice of the input virialized dark halo density profile and of the specific angular momentum distribution. We derive limits on the final density profile of the halo in the central regions. Subsequent viscous evolution of the disc is modelled by a variation of the specific angular momentum distribution of the disc, providing analytic solutions to the final disc structure. The assumption that the viscous evolution time-scale and the star formation time-scale are similar leads to predictions of the properties of the stellar components. Focusing on small 'exponential' bulges, i.e., ones that may be formed through a disc instability, we investigate the relationship between the assumed initial conditions, such as halo 'formation', or assembly, redshift z _f, spin parameter λ , baryonic fraction F , and final disc properties such as global star formation time-scale, gas fraction, and bulge-to-disc ratio. We find that the present properties of discs, such as the scalelength, are compatible with a higher initial formation redshift if the redistribution by viscous evolution is included than if it is ignored. We also quantify the dependence of final disc properties on the ratio F λ , thus including the possibility that the baryonic fraction varies from galaxy to galaxy, as perhaps may be inferred from the observations.     

The nature of high-redshift galaxies

Using semi-analytic models of galaxy formation set within the cold dark matter (CDM) merging hierarchy, we investigate several scenarios for the nature of the high-redshift     ) Lyman-break galaxies (LBGs). We consider a 'collisional starburst' model in which bursts of star formation are triggered by galaxy–galaxy mergers, and find that a significant fraction of LBGs are predicted to be starbursts. This model reproduces the observed comoving number density of bright LBGs as a function of redshift and the observed luminosity function at     and     with a reasonable amount of dust extinction. Model galaxies at     have star formation rates, half-light radii,     colours and internal velocity dispersions that are in good agreement with the data. Global quantities such as the star formation rate density and cold gas and metal content of the Universe as a function of redshift also agree well. Two 'quiescent' models without starbursts are also investigated. In one, the star formation efficiency in galaxies remains constant with redshift, while in the other, it scales inversely with disc dynamical time, and thus increases rapidly with redshift. The first quiescent model is strongly ruled out, as it does not produce enough high-redshift galaxies once realistic dust extinction is accounted for. The second quiescent model fits marginally, but underproduces cold gas and very bright galaxies at high redshift. A general conclusion is that star formation at high redshift must be more efficient than locally. The collisional starburst model appears to accomplish this naturally without violating other observational constraints.     

Disc galaxy evolution along the Hubble sequence

Galaxy discs are characterized by star formation histories that vary systematically along the Hubble sequence. We study global star formation, incorporating supernova feedback, gas accretion and enriched outflows in discs modelled by a multiphase interstellar medium in a fixed gravitational potential. The star formation histories, gas distributions and chemical evolution can be explained in a simple sequence of models which are primarily regulated by the cold gas accretion history.     

Semi-analytic modelling of galaxy formation: the local Universe

Using semi-analytic models of galaxy formation, we investigate galaxy properties such as the Tully–Fisher relation, the B - and K -band LFs, cold gas contents, sizes, metallicities and colours, and compare our results with observations of local galaxies. We investigate several different recipes for star formation and supernova feedback, including choices that are similar to the treatment by Kauffmann, White & Guiderdoni and Cole et al., as well as some new recipes. We obtain good agreement with all of the key local observations mentioned above. In particular, in our best models, we simultaneously produce good agreement with both the observed B - and K -band LFs and the I -band Tully–Fisher relation. Improved cooling and supernova feedback modelling, inclusion of dust extinction and an improved Press–Schechter model all contribute to this success. We present results for several variants of the CDM family of cosmologies, and find that models with values of Ω₀≃0.3–0.5 give the best agreement with observations.     

The evolution of galactic discs

We use recent observations of high-redshift galaxies to study the evolution of galactic discs over the redshift range 0 <  z ≲1. The data are inconsistent with models in which discs were already assembled at z  = 1 and have evolved only in luminosity since that time. Assuming that disc properties change with redshift as powers of 1 +   z and analysing the observations assuming an Einstein–de Sitter universe, we find that for given rotation speed, disc scalelength decreases with z as ∼ (1 +  z )⁻¹, total B -band mass-to-light ratio decreases with z as ∼ (1 +  z )⁻¹, and disc luminosity (again in B ) depends only weakly on z . These scalings are consistent with current data on the evolution of disc galaxy abundance as a function of size and luminosity. Both the scalings and the abundance evolution are close to the predictions of hierarchical models for galaxy formation. If different cosmogonies are compared, the observed evolution in disc size and disc abundance favours a flat low-Ω₀ universe over an Einstein–de Sitter universe.