A simple model for the relationship between star formation and surface density

We investigate the relationship between the star formation rate per unit area and the surface density of the interstellar medium (ISM; the local Kennicutt–Schmitt law) using a simplified model of the ISM and a simple estimate of the star formation rate based on the mass of gas in bound clumps, the local dynamical time-scales of the clumps and an efficiency parameter of around  ε≈ 5  per cent. Despite the simplicity of the approach, we are able to reproduce the observed linear relation between star formation rate and surface density of dense (molecular) gas. We use a simple model for the dependence of H₂ fraction on total surface density to argue why neither total surface density nor the H  i surface density is a good local indicator of star formation rate. We also investigate the dependence of the star formation rate on the depth of the spiral potential. Our model indicates that the mean star formation rate does not depend significantly on the strength of the spiral potential, but that a stronger spiral potential, for a given mean surface density, does result in more of the star formation occurring close to the spiral arms. This agrees with the observation that grand design galaxies do not appear to show a larger degree of star formation compared to their flocculent counterparts.     

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.     

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.     

Application of the global modal approach to the spiral galaxies

We have tested the applicability of the global modal approach in the density wave theory of spiral structure for a sample of spiral galaxies with measured axisymmetric background properties. We report here the results of the simulations for four galaxies: NGC 488, NGC 628, NGC 1566, and NGC 3938. Using the observed radial distributions for the stellar velocity dispersions and the rotation velocities we have constructed the equilibrium models for the galactic disks in each galaxy and implemented two kinds of stability analyses - the linear global analysis and 2D-nonlinear simulations. In general, the global modal approach is able to reproduce the observed properties of the spiral arms in the galactic disks. The growth of spirals in the galactic disks can be physically understood in terms of amplification by over-reflection at the corotation resonance. Our results support the global modal approach as a theoretical explanation of spiral structure in galaxies. This revised version was published online in August 2006 with corrections to the Cover Date.     

Statistical properties of H II regions in the disc of M100

From a new mosaic image in the Hα line of the complete disc of the spiral galaxy M100, a catalogue is composed listing 1948 individual H II regions. I give details of the data collection and reduction procedure, and of the production of the H II region catalogue. For each H II region, the catalogue gives its position relative to the centre of the galaxy, its deprojected distance to the centre, its radius and its calibrated luminosity. An indication is included as to whether the H II region is located in the arms, between them, or in the circumnuclear star-forming region. I present the results of a statistical study of properties of the H II regions. The luminosity function of the complete ensemble of H II regions shows a characteristic shape well fitted by a power-law slope in the higher luminosity range, and complying with literature values for galaxies like M100. Luminosity function slopes for arm and interarm H II region populations separately are found to be equal within the errors of the fits, indicating that whereas the density wave accumulates material into the arm regions, and may trigger star formation there, it does not in fact change the mass distribution of the star-forming clouds, nor the statistical properties of the H II region population. Diameter distributions and the radial number density distribution are discussed. The latter indicates those areas where most star formation occurs: the circumnuclear region and the spiral arms. The huge number of H II regions allowed the construction of a number of independent luminosity functions at different distances to the nucleus. The slope of the luminosity function shows a marginal decrease with increasing distance from the centre, which could indicate a gradual change towards shallower IMF slopes with increasing galactocentric distance, or an evolutionary effect.     

A dynamical model for the extraplanar gas in spiral galaxies

Recent H  i observations reveal that the discs of spiral galaxies are surrounded by extended gaseous haloes. This extraplanar gas reaches large distances (several kiloparsecs) from the disc and shows peculiar kinematics (low rotation and inflow). We have modelled the extraplanar gas as a continuous flow of material from the disc of a spiral galaxy into its halo region. The output of our models is pseudo data cubes that can be directly compared to the H  i data. We have applied these models to two spiral galaxies (NGC 891 and NGC 2403) known to have a substantial amount of extraplanar gas. Our models are able to reproduce accurately the vertical distribution of extraplanar gas for an energy input corresponding to a small fraction (<4 per cent) of the energy released by supernovae. However, they fail in two important aspects: (1) they do not reproduce the right gradient in rotation velocity; (2) they predict a general outflow of the extraplanar gas, contrary to what is observed. We show that neither of these difficulties can be removed if clouds are ionized and invisible at 21 cm as they leave the disc but become visible at some point on their orbits. We speculate that these failures indicate the need for accreted material from the intergalactic medium that could provide the low angular momentum and inflow required.     

Formation and evolution of elliptical galaxies and QSO activity

We present the results of a numerical code that combines multi-zone chemical evolution with 1D hydrodynamics to follow in detail the evolution and radial behaviour of gas and stars during the formation of elliptical galaxies. We use the model to explore the links between the evolution and formation of elliptical galaxies and QSO activity. The knowledge of the radial gas flows in the galaxy allows us to trace metallicity gradients, and, in particular, the formation of a high-metallicity core in ellipticals. The high-metallicity core is formed soon enough to explain the metal abundances inferred in high-redshift quasars. The star formation rate and the subsequent feedback regulate the episodes of wind, outflow and cooling flow, thus affecting the recycling of the gas and the chemical enrichment of the intergalactic medium. The evolution of the galaxy shows several stages, some of which are characterized by a complex flow pattern, with inflow in some regions and outflow in other regions. All models, however, exhibit during their late evolution a galactic wind at the outer boundary and, during their early evolution, an inflow towards the galactic nucleus. The characteristics of the inner inflow could explain the bolometric luminosity of a quasar lodged at the galactic centre as well as the evolution of the optical luminosity of quasars.     

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 unusual distribution of molecular gas and star formation in Arp 140

We investigate the atomic and molecular interstellar medium and star formation of NGC 275, the late-type spiral galaxy in Arp 140, which is interacting with NGC 274, an early-type system. The atomic gas (H  i ) observations reveal a tidal tail from NGC 275 which extends many optical radii beyond the interacting pair. The H  i morphology implies a prograde encounter between the galaxy pair approximately ∼1.5 × 10⁸ yr ago. The Hα emission from NGC 275 indicates clumpy irregular star formation, clumpiness which is mirrored by the underlying mass distribution as traced by the K _s-band emission. The molecular gas distribution is striking in its anticorrelation with the H  ii regions. Despite the evolved nature of NGC 275's interaction and its barred potential, neither the molecular gas nor the star formation is centrally concentrated. We suggest that this structure results from stochastic star formation leading to preferential consumption of the gas in certain regions of the galaxy. In contrast to the often-assumed picture of interacting galaxies, NGC 275, which appears to be close to merger, does not display enhanced or centrally concentrated star formation. If the eventual merger is to lead to a significant burst of star formation it must be preceded by a significant conversion of atomic to molecular gas as at the current rate of star formation all the molecular gas will be exhausted by the time the merger is complete.     

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 SAURON project – VII. Integral-field absorption and emission-line kinematics of 24 spiral galaxy bulges

We present observations of the stellar and gas kinematics for a representative sample of 24 Sa galaxies obtained with our custom-built integral-field spectrograph SAURON operating on the William Herschel Telescope. The data have been homogeneously reduced and analysed by means of a dedicated pipeline. All resulting data cubes were spatially binned to a minimum mean signal-to-noise ratio of 60 per spatial and spectral resolution element. Our maps typically cover the bulge-dominated region. We find a significant fraction of kinematically decoupled components (12/24), many of them displaying central velocity dispersion minima. They are mostly aligned and co-rotating with the main body of the galaxies, and are usually associated with dust discs and rings detected in unsharp-masked images. Almost all the galaxies in the sample (22/24) contain significant amounts of ionized gas which, in general, is accompanied by the presence of dust. The kinematics of the ionized gas are consistent with circular rotation in a disc co-rotating with respect to the stars. The distribution of mean misalignments between the stellar and gaseous angular momenta in the sample suggests that the gas has an internal origin. The [O  iii ]/Hβ ratio is usually very low, indicative of current star formation, and shows various morphologies (ring-like structures, alignments with dust lanes or amorphous shapes). The star formation rates (SFRs) in the sample are comparable with that of normal disc galaxies. Low gas velocity dispersion values appear to be linked to regions of intense star formation activity. We interpret this result as stars being formed from dynamically cold gas in those regions. In the case of NGC 5953, the data suggest that we are witnessing the formation of a kinematically decoupled component from cold gas being acquired during the ongoing interaction with NGC 5954.     

Magnetic and optical spiral arms in the galaxy NGC 6946

The spiral pattern in the nearby spiral galaxy NGC 6946 has been studied using the wavelet transformation technique, applied to galaxy images in polarized and total non-thermal radio emission at λλ 3.5 and 6.2 cm, in broadband red light, in the λ 21.1 cm H  i line and in the optical Hα line. Well-defined, continuous spiral arms are visible in polarized radio emission and red light, where we can isolate a multi-armed pattern in the range of galactocentric distances 1.5–12 kpc, consisting of four long arms and one short spiral segment. The 'magnetic arms' (visible in polarized radio emission) are localized almost precisely between the optical arms. Each magnetic arm is similar in length and pitch angle to the preceding optical arm (in the sense of galactic rotation) and can be regarded as its phase-shifted image. Even details like a bifurcation of an optical arm have their phase-shifted counterparts in the magnetic arms. The average relative amplitude of the optical spiral arms (the stellar density excess over the azimuthal average) grows with galactocentric radius up to 0.3–0.7 at r ≃5 kpc, decreases by a factor of two at r =5–6 kpc and remains low at 0.2–0.3 in the outer parts of the galaxy. By contrast, the magnetic arms have a constant average relative amplitude (the excess in the regular magnetic field strength over the azimuthal average) of 0.3–0.6 in a wide radial range r =1.5–12 kpc. We briefly discuss implications of our findings for theories of galactic magnetic fields.     

Chemical evolution with radial mixing

Models of the chemical evolution of our Galaxy are extended to include radial migration of stars and flow of gas through the disc. The models track the production of both iron and α-elements. A model is chosen that provides an excellent fit to the metallicity distribution of stars in the Geneva–Copenhagen survey (GCS) of the solar neighbourhood and a good fit to the local Hess diagram. The model provides a good fit to the distribution of GCS stars in the age–metallicity plane, although this plane was not used in the fitting process. Although this model's star formation rate is monotonically declining, its disc naturally splits into an α-enhanced thick disc and a normal thin disc. In particular, the model's distribution of stars in the ([O/Fe], [Fe/H]) plane resembles that of Galactic stars in displaying a ridge line for each disc. The thin-disc's ridge line is entirely due to stellar migration, and there is the characteristic variation of stellar angular momentum along it that has been noted by Haywood in survey data. Radial mixing of stellar populations with high  σ _z   from inner regions of the disc to the solar neighbourhood provides a natural explanation of why measurements yield a steeper increase of  σ _z   with age than predicted by theory. The metallicity gradient in the interstellar medium is predicted to be steeper than in earlier models, but appears to be in good agreement with data for both our Galaxy and external galaxies. The models are inconsistent with a cut-off in the star formation rate at low gas surface densities. The absolute magnitude of the disc is given as a function of time in several photometric bands, and radial colour profiles are plotted for representative times.     

Ram pressure stripping of spiral galaxies in clusters

We use three-dimensional SPH/ N -body simulations to study ram pressure stripping of gas from spiral galaxies orbiting in clusters. We find that the analytic expectation of Gunn & Gott, relating the gravitational restoring force provided by the disc to the ram pressure force, provides a good approximation to the radius at which gas will be stripped from a galaxy. However, at small radii it is also important to consider the potential provided by the bulge component. A spiral galaxy passing through the core of a rich cluster, such as Coma, will have its gaseous disc truncated to ∼4 kpc, thus losing ∼80 per cent of its diffuse gas mass. The time-scale for this to occur is a fraction of a crossing time ∼10⁷ yr. Galaxies orbiting within poorer clusters, or inclined to the direction of motion through the intracluster medium, will lose significantly less gas. We conclude that ram pressure alone is insufficient to account for the rapid and widespread truncation of star formation observed in cluster galaxies, or the morphological transformation of Sabs to S0s that is necessary to explain the Butcher–Oemler effect.     

An Analytic Model of Galactic Winds and Mass Outflows

Galactic winds and mass outflows are observed both in nearby starburst galaxies and in high-redshift star-forming galaxies. We develop a simple analytic model to understand the observed superwind phenomenon with a discussion of the model uncertainties. Our model is built upon the model of McKee & Ostriker for the interstellar medium. It allows one to predict how properties of a superwind, such as wind velocity and mass outflow rate, are related to properties of its starforming host galaxy, such as size, gas density and star formation rate. The model predicts a threshold of star formation rate density for the generation of observable galactic winds. Galaxies with more concentrated star formation activities produce superwinds with higher velocities. The predicted mass outflow rates are comparable to (or slightly larger than) the corresponding star formation rates. We apply our model to both local starburst galaxies and high-redshift Lyman break galaxies, and find its predictions to be in good agreement with current observations. Our model is simple and so can be easily incorporated into numerical simulations and semi-analytical models of galaxy formation.     

Stellar and gas dynamics of late-type barred-spiral galaxies: NGC 3359, a test case

We study the dynamics of a model for the late-type barred-spiral galaxy NGC 3359 by using both observational and numerical techniques. The results of our modelling are compared with photometric and kinematical data. The potential used is estimated directly from observations of the galaxy. It describes with a single potential function, a barred-spiral system with an extended spiral structure. Thus, the study of the dynamics in this potential has an interest by itself. We apply orbital theory and response models for the study of the stellar component, and smoothed particle hydrodynamics for modelling the gas. In particular, we examine the pattern speed of the system and the orbital character (chaotic or ordered) of the spiral arms. We conclude that the spiral pattern rotates slowly, in the sense that its corotation is close to or even beyond the end of the arms. Although a single, slow pattern speed could, under certain assumptions, characterize the whole disc, the comparison with the observational data indicates that probably the bar and the spirals have different angular velocities. In our two pattern speeds model, the best fit is obtained with a bar ending close to its 4:1 resonance and a more slowly rotating spiral. Assuming an 11 Mpc distance to the galaxy, a match of our models with the observed data indicates a pattern speed of about  39 km s⁻¹ kpc⁻¹  for the bar and about  15 km s⁻¹ kpc⁻¹  for the spiral. We do not find any indication for a chaotic character of the arms in this barred-spiral system. The flow in the region of the spirals can best be described as a regular 'precessing-ellipses flow'.     

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.     

The origin of the light distribution in spiral galaxies

We analyse a high-resolution, fully cosmological, hydrodynamical disc galaxy simulation, to study the source of the double-exponential light profiles seen in many stellar discs, and the effects of stellar radial migration upon the spatiotemporal evolution of both the disc age and metallicity distributions. We find a 'break' in the pure exponential stellar surface brightness profile, and trace its origin to a sharp decrease in the star formation per unit surface area, itself produced by a decrease in the gas volume density due to a warping of the gas disc. Star formation in the disc continues well beyond the break. We find that the break is more pronounced in bluer wavebands. By contrast, we find little or no break in the mass density profile. This is, in part, due to the net radial migration of stars towards the external parts of the disc. Beyond the break radius, we find that ∼60 per cent of the resident stars migrated from the inner disc, while ∼25 per cent formed in situ . Our simulated galaxy also has a minimum in the age profile at the break radius but, in disagreement with some previous studies, migration is not the main mechanism producing this shape. In our simulation, the disc metallicity gradient flattens with time, consistent with an 'inside-out' formation scenario. We do not find any difference in the intensity or the position of the break with inclination, suggesting that perhaps the differences found in empirical studies are driven by dust extinction.     

Modelling gaseous and stellar kinematics in the disc galaxies NGC 772, 3898 and 7782

We present V -band surface photometry and major-axis kinematics of stars and ionized gas of three early-type spiral galaxies, namely NGC 772, 3898 and 7782. For each galaxy we present a self-consistent Jeans model for the stellar kinematics, adopting the light distribution of bulge and disc derived by means of a two-dimensional parametric photometric decomposition. This allows us to investigate the presence of non-circular gas motions, and derive the mass distribution of luminous and dark matter in these objects.
NGC 772 and 7782 have apparently normal kinematics with the ionized gas tracing the gravitational equilibrium circular speed. This is not true in the innermost region (| r |≲8 arcsec) of NGC 3898, where the ionized gas is rotating more slowly than the circular velocity predicted by dynamical modelling. This phenomenon is common in the bulge-dominated galaxies for which dynamical modelling enables us to make the direct comparison between the gas velocity and the circular speed, and it poses questions about the reliability of galaxy mass distributions derived by the direct decomposition of the observed ionized-gas rotation curve into the contributions of luminous and dark matter.     

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.