Hydrodynamics of photoionized columns in the Eagle Nebula, M 16

We present hydrodynamical simulations of the formation, structure and evolution of photoionized columns, with parameters based on those observed in the Eagle Nebula. On the basis of these simulations we argue that there is no unequivocal evidence that the dense neutral clumps at heads of the columns were cores in the pre-existing molecular cloud. In our simulations, a variety of initial conditions leads to the formation and maintenance of near-equilibrium columns. Therefore, it is likely that narrow columns will often occur in regions with large-scale inhomogeneities, but that observations of such columns can tell us little about the processes by which they formed. The manner in which the columns in our simulations develop suggests that their evolution may result in extended sequences of radiation-induced star formation.     

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.     

Galaxy formation: Warm dark matter, missing satellites, and the angular momentum problem

We present warm dark matter (WDM) as a possible solution to the missing satellites and angular momentum problem in galaxy formation and introduce improved initial conditions for numerical simulations of WDM models, which avoid the formation of unphysical haloes found in earlier simulations. There is a hint, that because of that the mass function of satellite haloes has been overestimated so far, pointing to higher values for the WDM particle mass. This revised version was published online in August 2006 with corrections to the Cover Date.     

Star formation and chemical evolution in smoothed particle hydrodynamics simulations: a statistical approach

In smoothed particle hydrodynamics (SPH) codes with a large number of particles, star formation as well as gas and metal restitution from dying stars can be treated statistically. This approach allows one to include detailed chemical evolution and gas re-ejection with minor computational effort. Here we report on a new statistical algorithm for star formation and chemical evolution, especially conceived for SPH simulations with large numbers of particles, and for parallel SPH codes.
For the sake of illustration, we also present two astrophysical simulations obtained with this algorithm, implemented into the Tree-SPH code by Lia & Carraro .
In the first simulation, we follow the formation of an individual disc-like galaxy, predict the final structure and metallicity evolution, and test resolution effects. In the second simulation we simulate the formation and evolution of a cluster of galaxies, to demonstrate the capabilities of the algorithm in investigating the chemo-dynamical evolution of galaxies and of the intergalactic medium in a cosmological context.     

Parameter tests within cosmological simulations of galaxy formation

Numerical simulations of galaxy formation require a number of parameters. Some of these are intrinsic to the numerical integration scheme (e.g., the time-step), while others describe the physical model (e.g., the gas metallicity). In this paper we present results of a systematic exploration of the effects of varying a subset of these parameters on simulations of galaxy formation. We use N -body and 'Smoothed Particle Hydrodynamics' techniques to follow the evolution of cold dark matter and gas in a small volume. We compare a fiducial model with 24 different simulations, in which one parameter at a time is varied, focusing on properties such as the relative fraction of hot and cold gas, and the abundance and masses of galaxies. We find that for reasonable choices of numerical values, many parameters have relatively little effect on the galaxies, with the notable exception of the parameters that control the resolution of the simulation and the efficiency with which gas cools.     

Disc-like objects in hierarchical hydrodynamical simulations: comparison with observations

We present results from a careful and detailed analysis of the structural and dynamical properties of a sample of 29 disc-like objects identified at z =0 in three AP3M–SPH fully consistent cosmological simulations. These simulations are realizations of a CDM hierarchical model, in which an inefficient Schmidt-law-like algorithm to model the stellar formation process has been implemented. We focus on properties that can be constrained with available data from observations of spiral galaxies, namely the bulge and disc structural parameters and the rotation curves. Comparison with data from Broeils, de Jong and Courteau gives satisfactory agreement, in contrast with previous findings using other codes. This suggests that the stellar formation implementation we have used has succeeded in forming compact bulges that stabilize disc-like structures in the violent phases of their assembly, while in the quiescent phases the gas has cooled and collapsed in accord with the Fall & Efstathiou standard model of disc formation.     

On the Reliability of Merger-Trees and the Mass-Growth Histories of Dark Matter Haloes

We have used merger-trees realizations to study the formation of dark matter haloes. The construction of merger-trees is based on three different pictures about the formation of structures in the Universe. These pictures include the spherical collapse (SC), the ellipsoidal collapse (EC) and the non-radial collapse (NR). The reliability of merger-trees has been examined comparing their predictions related to the distribution of the number of progenitors, as well as the distribution of formation times, with the predictions of analytical relations. The comparison yields a very satisfactory agreement. Subsequently, the mass-growth histories (MGH) of haloes have been studied and their formation scale factors have been derived. This derivation has been based on two different definitions that are (a) the scale factor when the halo reaches half its present day mass and (b) the scale factor when the mass-growth rate falls below some specific value. Formation scale factors follow approximately power laws of mass. It has also been shown that MGHs are in good agreement with models proposed in the literature that are based on the results of N-body simulations. The agreement is found to be excellent for small haloes but, at the early epochs of the formation of large haloes, MGHs seem to be steeper than those predicted by the models based on N-body simulations. This rapid growth of mass of heavy haloes is likely to be related to a steeper central density profile indicated by the results of some N-body simulations.     

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.     

Populating a cluster of galaxies – I. Results at

We simulate the assembly of a massive rich cluster and the formation of its constituent galaxies in a flat, low-density universe. Our most accurate model follows the collapse, the star formation history and the orbital motion of all galaxies more luminous than the Fornax dwarf spheroidal, while dark halo structure is tracked consistently throughout the cluster for all galaxies more luminous than the SMC. Within its virial radius this model contains about     dark matter particles and almost 5000 distinct dynamically resolved galaxies. Simulations of this same cluster at a variety of resolutions allow us to check explicitly for numerical convergence both of the dark matter structures produced by our new parallel N -body and substructure identification codes, and of the galaxy populations produced by the phenomenological models we use to follow cooling, star formation, feedback and stellar aging. This baryonic modelling is tuned so that our simulations reproduce the observed properties of isolated spirals outside clusters. Without further parameter adjustment our simulations then produce a luminosity function, a mass-to-light ratio, luminosity, number and velocity dispersion profiles, and a morphology–radius relation which are similar to those observed in real clusters. In particular, since our simulations follow galaxy merging explicitly, we can demonstrate that it accounts quantitatively for the observed cluster population of bulges and elliptical galaxies.     

Dynamical and pressure structures in winds with multiple embedded evaporating clumps – I. Two-dimensional numerical simulations

Because of its key role in feedback in star formation and galaxy formation, we examine the nature of the interaction of a flow with discrete sources of mass injection. We show the results of two-dimensional numerical simulations in which we explore a range of configurations for the mass sources and study the effects of their proximity on the downstream flow. The mass sources act effectively as a single source of mass injection if they are so close together that the ratio of their combined mass injection rate is comparable to or exceeds the mass flux of the incident flow into the volume that they occupy. The simulations are relevant to many diffuse sources, such as planetary nebulae and starburst superwinds, in which a global flow interacts with material evaporating or being ablated from the surface of globules of cool, dense gas.     

H i and CO observations of Arp 104: a spiral–elliptical interacting pair

We present data probing the spatial and kinematical distribution of both the atomic (H  i ) and molecular (CO) gas in NGC 5218, the late-type barred spiral galaxy in the spiral–elliptical interacting pair, Arp 104. We consider these data in conjunction with far-infrared and radio-continuum data, and N -body simulations, to study the galaxies interactions, and the star formation properties of NGC 5218. We use these data to assess the importance of the bar and tidal interaction on the evolution of NGC 5218, and the extent to which the tidal interaction may have been important in triggering the bar. The molecular gas distribution of NGC 5218 appears to have been strongly affected by the bar; the distribution is centrally condensed with a very large surface density in the central region. The N -body simulations indicate a time-scale since perigalacticon of  ∼3 × 10⁸ yr  , which is consistent with the interaction having triggered or enhanced the bar potential in NGC 5218, leading to inflow and the large central molecular gas density observed. Whilst NGC 5218 appears to be undergoing active star formation, its star formation efficiency is comparable to a 'normal' SBb galaxy. We propose that this system may be on the brink of a more active phase of star formation.     

Merging time-scales of stellar subclumps in young star-forming regions

Recent observations and hydrodynamical simulations of star formation inside a giant molecular cloud have revealed that, within a star-forming region, stars do not form evenly distributed throughout this region, but rather in small subclumps. It is generally believed that these subclumps merge and form a young star cluster. The time-scale of this merging process is crucial for the evolution and the possible survival of the final star cluster. The key issue is whether this merging process happens faster than the time needed to remove the residual gas of the cloud. A merging time-scale shorter than the gas-removal time would enhance the survival chances of the resulting star cluster. In this paper, we show by means of numerical simulations that the time-scale of the merging is indeed very fast. Depending on the details of the initial subclump distribution, the merging may occur before the gas is expelled from the newly formed cluster via either supernovae or the winds from massive stars. Our simulations further show that the resulting merger objects have a higher effective star formation efficiency than the overall star-forming region and confirm the results that mass-segregated subclumps form mass-segregated merger objects.     

Bulges versus discs: the evolution of angular momentum in cosmological simulations of galaxy formation

We investigate the evolution of angular momentum in simulations of galaxy formation in a cold dark matter universe. We analyse two model galaxies generated in the N -body/hydrodynamic simulations of Okamoto et al. Starting from identical initial conditions, but using different assumptions for the baryonic physics, one of the simulations produced a bulge-dominated galaxy and the other one a disc-dominated galaxy. The main difference is the treatment of star formation and feedback, both of which were designed to be more efficient in the disc-dominated object. We find that the specific angular momentum of the disc-dominated galaxy tracks the evolution of the angular momentum of the dark matter halo very closely: the angular momentum grows as predicted by linear theory until the epoch of maximum expansion and remains constant thereafter. By contrast, the evolution of the angular momentum of the bulge-dominated galaxy resembles that of the central, most bound halo material: it also grows at first according to linear theory, but 90 per cent of it is rapidly lost as pre-galactic fragments, into which gas had cooled efficiently, merge, transferring their orbital angular momentum to the outer halo by tidal effects. The disc-dominated galaxy avoids this fate because the strong feedback reheats the gas, which accumulates in an extended hot reservoir and only begins to cool once the merging activity has subsided. Our analysis lends strong support to the classical theory of disc formation whereby tidally torqued gas is accreted into the centre of the halo conserving its angular momentum.     

Pattern speed of the stellar bar in NGC 7479

We study the pattern speed of the bar in NGC 7479 by comparing observations with numerical simulations of gas flow in a two-dimensional gravitational potential, derived from observations. The best agreement between the observations and the modelling is achieved for the fast bar pattern speed of 27 km s⁻¹ kpc⁻¹, when the corotation radius is at 50 arcsec, i.e. 1.1 times the radial length of the bar. This result is supported by the gas and dust lane morphologies, star formation distribution, projected velocity field and overall morphology. We find that star formation is most likely to be triggered close to the large-scale shocks and dust lanes in the bar. The net gas inflow rate in the simulations at 1-kpc radius is 4–6 M⊙ yr⁻¹ at intermediate times.     

Testing the modified Press–Schechter model against N-body simulations

A modified version of the extended Press–Schechter model for the growth of dark-matter haloes was introduced in two previous papers, with the aim of explaining the mass–density relation shown by haloes in high-resolution cosmological simulations. In this model, major mergers are well separated from accretion, thereby allowing a natural definition of halo formation and destruction. This makes it possible to derive analytic expressions for halo formation and destruction rates, the mass accretion rate and the probability distribution functions of halo formation times and progenitor masses. The stochastic merger histories of haloes can be readily derived and easily incorporated into semi-analytical models of galaxy formation, thus avoiding the usual problems encountered in the construction of Monte Carlo merger trees from the original extended Press–Schechter formalism. Here we show that the predictions of the modified Press–Schechter model are in good agreement with the results of N -body simulations for several scale-free cosmologies.     

Gas in simulations of high-redshift galaxies and minihaloes

We study the gas content of haloes in the early universe using high-resolution hydrodynamical simulations. We extract from the simulations and also predict, based on linear theory, the halo mass for which the enclosed baryon fraction equals half of the mean cosmic fraction. We find a rough agreement between the simulations and the predictions, which suggests that during the high-redshift era before stellar heating, the minimum mass needed for a minihalo to keep most of its baryons throughout its formation was  ∼3 × 10⁴ M_⊙  . We also carry out a detailed resolution analysis and show that in order to determine a halo's gas fraction even to 20 per cent accuracy, the halo must be resolved into at least 500 dark matter particles.     

Three-dimensional evolution of worms and chimneys in the Galactic disc

Three-dimensional modelling of the flow of gas and plasma in a section of the Galaxy has been carried out to study the evolution and formation of Galactic chimneys and worms. It is found that clustered supernovae located on either side of the Galactic plane are sources for the formation of well-collimated chimneys, having widths of ∼     . The thick gas disc may have a role in the collimation of chimneys. Channel maps of disc gas, obtained from the simulations, show the presence of sheet-like structures running perpendicular to the Galactic plane and resembling worms. Worms are believed to result from the break-up of the shells and supershells. However, the simulations show that although some worms correlate well with the debris of broken shells/supershells, others do not. They are cold gas that has been accelerated in the disc and rise on to the thick gas disc.