A model of supernova feedback in galaxy formation

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

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⁻²  .     

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.     

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 .     

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.     

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.     

Numerical simulations of the formation and chemical evolution of galaxies

We present the results of a set of three-dimensional SPH-Treecode simulations which model the formation and early evolution of disc galaxies, including the generation and return of heavy elements to the interstellar medium by star formation. Starting from simple initial conditions which are given by a uniform density sphere of gas which is embedded in a dark matter halo and in solid-body rotation, we are able to form realistic disc galaxies, and find that an exponential gas disc is quickly formed. Star formation within this exponential disc naturally leads to the formation of abundance gradients which are in broad agreement with those observed, although they are slightly shallower than some observations.
We investigate the systematic effects of variation of mass, rotation and star formation parameters on the abundance gradients. We find that the abundance gradients are most sensitive to changes in the star formation parameters or rotation. Including a critical-density cut-off in the star formation law causes abundance gradients to be steepened.
Analysis of gas flows within the models shows radial flows which are a function of angle of azimuth around the galaxies, with alternating inward and outward flows. This motion is linked to the presence of a bar, whose strength is related to the amount of star formation in the models, and there is a gentle drift of mass inwards. The shallow abundance gradients may be linked to these radial flows.     

Storm fronts over galaxy discs: models of how waves generate extraplanar gas and its anomalous kinematics

The existence of partially ionized, diffuse gas and dust clouds at kiloparsec scale distances above the central planes of edge-on, galaxy discs was an unexpected discovery about 20 years ago. Subsequent observations showed that this extended or extraplanar diffuse interstellar gas (EDIG) has rotation velocities approximately 10–20 per cent lower than those in the central plane, and has been hard to account for. Here, we present results of hydrodynamic models, with radiative cooling and heating from star formation. We find that in models with star formation generated stochastically across the disc, an extraplanar gas layer is generated as long as the star formation is sufficiently strong. However, this gas rotates at nearly the same speed as the midplane gas. We then studied a range of models with imposed spiral or bar waves in the disc. EDIG layers were also generated in these models, but primarily over the wave regions, not over the entire disc. Because of this partial coverage, the EDIG clouds move radially, as well as vertically, with the result that observed kinematic anomalies are reproduced. The implication is that the kinematic anomalies are the result of three-dimensional motions when the cylindrical symmetry of the disc is broken. Thus, the kinematic anomalies are the result of bars or strong waves, and more face-on galaxies with such waves should have an asymmetric EDIG component. The models also indicate that the EDIG can contain a significant fraction of cool gas, and that some star formation can be triggered at considerable heights above the disc mid-plane. We expect all of these effects to be more prominent in young, forming discs, to play a role in rapidly smoothing disc asymmetries and in working to self-regulate disc structure.     

The effect of galaxy mass ratio on merger-driven starbursts

We employ numerical simulations of galaxy mergers to explore the effect of galaxy mass ratio on merger-driven starbursts. Our numerical simulations include radiative cooling of gas, star formation, and stellar feedback to follow the interaction and merger of four disc galaxies. The galaxy models span a factor of 23 in total mass and are designed to be representative of typical galaxies in the local universe. We find that the merger-driven star formation is a strong function of merger mass ratio, with very little, if any, induced star formation for large mass ratio mergers. We define a burst efficiency that is useful to characterize the merger-driven star formation and test that it is insensitive to uncertainties in the feedback parametrization. In accord with previous work we find that the burst efficiency depends on the structure of the primary galaxy. In particular, the presence of a massive stellar bulge stabilizes the disc and suppresses merger-driven star formation for large mass ratio mergers. Direct, coplanar merging orbits produce the largest tidal disturbance and yield the most intense burst of star formation. Contrary to naive expectations, a more compact distribution of gas or an increased gas fraction both decrease the burst efficiency. Owing to the efficient feedback model and the newer version of smoothed particle hydrodynamics employed here, the burst efficiencies of the mergers presented here are smaller than in previous studies.     

The star formation rate in disc galaxies: thresholds and dependence on gas amount

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

Star formation in southern Seyfert galaxies

We have produced radio maps, using the Australia Telescope Compact Array, of the central regions of six southern type 2 Seyfert galaxies (NGC 1365, 4945, 6221, 6810, 7582 and Circinus) with circumnuclear star formation, to estimate the relative contribution of star formation activity compared to activity from the active galactic nucleus (AGN). The radio morphologies range from extended diffuse structures to compact nuclear emission, with no evidence, even in the relatively compact sources, for synchrotron self-absorption. In each case the radio to far-infrared (FIR) ratio has a value consistent with star formation, and in all but one case the radio to [Fe  II ] ratio is also consistent with star formation. We derive supernova rates and conclude that, despite the presence of a Seyfert nucleus in these galaxies, the radio, FIR and [Fe  II ] line emissions are dominated by processes associated with the circumnuclear star formation (i.e. supernova remnants and H  II regions) rather than with the AGN.     

A high-resolution radio study of neutral gas in the starburst galaxy NGC 520

We present subarcsec angular resolution observations of the neutral gas in the nearby starburst galaxy NGC 520. The central kpc region of NGC 520 contains an area of significantly enhanced star formation. The radio continuum structure of this region resolves into ∼10 continuum components. By comparing the flux densities of the brightest of these components at 1.4 GHz with published 15-GHz data we infer that these components detected at 1.4 and 1.6 GHz are related to the starburst and are most likely to be collections of several supernova remnants within the beam. None of these components is consistent with emission from an active galactic nuclei. Both neutral hydrogen (H  i ) and hydroxyl (OH) absorption lines are observed against the continuum emission, along with a weak OH maser feature probably related to the star formation activity in this galaxy. Strong H  i absorption  ( N _H∼ 10²² atoms cm⁻²)  traces a velocity gradient of 0.5 km s⁻¹ pc⁻¹ across the central kpc of NGC 520. The H  i absorption velocity structure is consistent with the velocity gradients observed in both the OH absorption and in CO emission observations. The neutral gas velocity structure observed within the central kpc of NGC 520 is attributed to a kpc-scale ring or disc. It is also noted that the velocity gradients observed for these neutral gas components appear to differ with the velocity gradients observed from optical ionized emission lines. This apparent disagreement is discussed and attributed to the extinction of the optical emission from the actual centre of this source hence implying that optical ionized emission lines are only detected from regions with significantly different radii to those sampled by the observations presented here.     

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.     

Accretion of gas on to nearby spiral galaxies

We present evidence for cosmological gas accretion on to spiral galaxies in the local universe. The accretion is seen through its effects on the dynamics of the extraplanar neutral gas. The accretion rates that we estimate for two nearby spiral galaxies are of the order of their star formation rates. Our model shows that most of the extraplanar gas is produced by supernova feedback (galactic fountain) and only 10–20 per cent comes from accretion. The accreting material must have low specific angular momentum about the disc's spin axis, although the magnitude of the specific angular momentum vector can be higher. We also explore the effects of a hot corona on the dynamics of the extraplanar gas and find that it is unlikely to be responsible for the observed kinematical pattern and the source of accreted gas. However, the interaction with the fountain flow should profoundly affect the hydrodynamics of the corona.     

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.     

Disc–halo interaction – I. Three-dimensional evolution of the Galactic disc

The results of a three-dimensional model for disc–halo interaction are presented here. The model considers explicitly the input of energy and mass by isolated and correlated supernovae in the disc. Once disrupted by the explosions, the disc never returns to its initial state. Instead it approaches a state where a thin H  i disc is formed in the Galactic plane, overlaid by thick H  i and H  ii gas discs with scaleheights of 500 pc and 1–1.5 kpc, respectively. The upper parts of the thick H  ii disc (the diffuse ionized medium) act as a disc–halo interface, and its formation and stability are directly correlated to the supernova rate per unit area in the simulated disc.     

A grid of chemical evolution models as a tool to interpret spiral and irregular galaxies data

We present a generalization of the multiphase chemical evolution model (CEM) applied to a wide set of theoretical galaxies with different masses and evolutionary rates. This generalized set of models has been computed using the so-called universal rotation curve from Persic, Salucci & Steel to calculate the radial mass distribution of 44 theoretical protogalaxies. This distribution is a fundamental input which, besides its own effect on the galaxy evolution, defines the characteristic collapse time-scale or gas infall rate on to the disc. We have adopted 10 sets of values, between 0 and 1, for the molecular cloud and star formation efficiencies, as corresponding to their probability nature, for each one of the radial distributions of total mass. Thus, we have constructed a biparametric grid of models, depending on those efficiency sets and on the rotation velocity, whose results are valid in principle for any spiral or irregular galaxy. The model results provide the time-evolution of different regions of the disc and the halo along galactocentric distance, measured by the gas (atomic and molecular) and stellar masses, the star formation rate (SFR) and chemical abundances of 14 elements, for a total of 440 models. This grid may be used to estimate the evolution of a given galaxy for which only present time information, such as radial distributions of elemental abundances, gas densities and/or star formation, which are the usual observational constraints of chemical evolution models (CEMs), is available.     

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.     

Maximum feedback and dark matter profiles of dwarf galaxies

The observed rotation curves of dark matter-dominated dwarf galaxies indicate low-density cores, contrary to the predictions of CDM models. A possible solution of this problem involves stellar feedback. A strong baryonic wind driven by vigorous star formation can remove a large fraction of the gas, causing the dark matter to expand. Using both numerical and analytical techniques, we explore the maximum effect of the feedback with an instantaneous removal of the gaseous disc. The energy input depends on the compactness of the disc, hence the specific angular momentum of the disc. For the plausible cosmological parameters and a wide range of the disc angular momenta, the feedback is insufficient to destroy the central halo cusp, while the inner density is lowered only by a modest factor of 2 to 6. Any realistic modelling of the feedback would have even lesser impact on dark matter. We find that no star formation effect can resolve the problems of CDM cusps.     

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.