Sustaining star formation rates in spiral galaxies - Supernova-driven turbulent accretion disk models applied to THINGS galaxies

Gas disks of spiral galaxies can be described as clumpy accretion disks without a coupling of viscosity to the actual thermal state of the gas. The model description of a turbulent disk consisting of emerging and spreading clumps (Vollmer & Beckert 2003) contains free parameters, which can be constrained by observations of molecular gas, atomic gas and the star formation rate for individual galaxies. Radial profiles of 18 nearby spiral galaxies from THINGS, HERACLES, SINGS, and GALEX data are used to compare the observed star formation efficiency, molecular fraction, and velocity dispersion to the model. The observed radially decreasing velocity dispersion can be reproduced by the model. In the framework of this model the decrease in the inner disk is due to the stellar mass distribution which dominates the gravitational potential. Introducing a radial break in the star formation efficiency into the model improves the fits significantly. This change in star formation regime is realized by replacing the free fall time in the prescription of the star formation rate with the molecule formation timescale. Depending on the star formation prescription, the break radius is located near the transition region between the molecular-gas-dominated and atomic-gas-dominated parts of the galactic disk or closer to the optical radius. It is found that only less massive galaxies (log (M (M_solar)) <~ 10) can balance gas loss via star formation by radial gas accretion within the disk. These galaxies can thus access their gas reservoirs with large angular momentum. On the other hand, the star formation of massive galaxies is determined by the external gas mass accretion rate from a putative spherical halo of ionized gas or from satellite accretion.Comment: 26 pages, 5 figures, full figures 1 and 2 as ancillary pdf files available; accepted by A

Arm & Interarm Star Formation in Spiral Galaxies

We investigate the relationship between spiral arms and star formation in the grand-design spirals NGC 5194 and NGC 628 and in the flocculent spiral NGC 6946. Filtered maps of near-IR (3.6 micron) emission allow us to identify "arm regions" that should correspond to regions of stellar mass density enhancements. The two grand-design spirals show a clear two-armed structure, while NGC 6946 is more complex. We examine these arm and interarm regions, looking at maps that trace recent star formation - far-ultraviolet (GALEX NGS) and 24 micron emission (Spitzer, SINGS) - and cold gas - CO (Heracles) and HI (Things). We find the star formation tracers and CO more concentrated in the spiral arms than the stellar 3.6 micron flux. If we define the spiral arms as the 25% highest pixels in the filtered 3.6 micron images, we find that the majority (60%) of star formation tracers occurs in the interarm regions; this result persists qualitatively even when considering the potential impact of finite data resolution and diffuse interarm 24 micron emission. Even with a generous definition of the arms (45% highest pixels), interarm regions still contribute at least 30% to the integrated star formation rate tracers. We look for evidence that spiral arms trigger star or cloud formation using the ratios of star formation rate (SFR, traced by a combination of FUV and 24 micron emission) to H_2 (traced by CO) and H_2 to HI. Any enhancement of SFR / M(H_2) in the arm region is very small (less than 10%) and the grand design spirals show no enhancement compared to the flocculent target. Arm regions do show a weak enhancement in H_2/HI compared to the interarm regions, but at a fixed gas surface density there is little clear enhancement in the H_2/HI ratio in the arm regions. Thus, it seems that spiral arms may only act to concentrate the gas to higher densities in the arms.Comment: 11 pages, 9 Figures, accepted by Ap

Star formation efficiency in the Barred Spiral Galaxy NGC 4303

We present new $^{12}$CO(J=1-0) observations of the barred galaxy NGC 4303 using the Nobeyama 45m telescope (NRO45) and the Combined Array for Research in Millimeter-wave Astronomy (CARMA). The H$\alpha$ images of barred spiral galaxies often show active star formation in spiral arms, but less so in bars. We quantify the difference by measuring star formation rate and efficiency at a scale where local star formation is spatially resolved. Our CO map covers the central 2\farcm3 region of the galaxy; the combination of NRO45 and CARMA provides a high fidelity image, enabling accurate measurements of molecular gas surface density. We find that star formation rate and efficiency are twice as high in the spiral arms as in the bar. We discuss this difference in the context of the Kennicutt-Schimidt (KS) law, which indicates a constant star formation rate at a given gas surface density. The KS law breaks down at our native resolution ($\sim$ 250 pc), and substantial smoothing (to 500 pc) is necessary to reproduce the KS law, although with greater scatter.Comment: 17 pages, 10 figures, published by ApJ; http://adsabs.harvard.edu/abs/2010ApJ...721..383

The role of feedback in shaping the structure of the interstellar medium

We present an analysis of the role of feedback in shaping the neutral hydrogen (H I) content of simulated disc galaxies. For our analysis, we have used two realizations of two separate Milky Way-like (similar to L star) discs - one employing a conservative feedback scheme (McMaster Unbiased Galaxy Survey), the other significantly more energetic [Making Galaxies In a Cosmological Context (MaGICC)]. To quantify the impact of these schemes, we generate zeroth moment (surface density) maps of the inferred H I distribution; construct power spectra associated with the underlying structure of the simulated cold interstellar medium, in addition to their radial surface density and velocity dispersion profiles. Our results are compared with a parallel, self-consistent, analysis of empirical data from The H I Nearby Galaxy Survey (THINGS). Single power-law fits (P proportional to k(gamma)) to the power spectra of the stronger feedback (MaGICC) runs (over spatial scales corresponding to similar to 0.5 to similar to 20 kpc) result in slopes consistent with those seen in the THINGS sample (gamma similar to -2.5). The weaker feedback (MUGS) runs exhibit shallower power-law slopes (gamma similar to -1.2). The power spectra of the MaGICC simulations are more consistent though with a two-component fit, with a flatter distribution of power on larger scales (i.e. gamma similar to -1.4 for scales in excess of similar to 2 kpc) and a steeper slope on scales below similar to 1 kpc (gamma similar to -5), qualitatively consistent with empirical claims, as well as our earlier work on dwarf discs. The radial H I surface density profiles of the MaGICC discs show a clear exponential behaviour, while those of the MUGS suite are essentially flat; both behaviours are encountered in nature, although the THINGS sample is more consistent with our stronger (MaGICC) feedback runs

The dark matter halo shape of edge-on disk galaxies - I. HI observations

This is the first paper of a series in which we will attempt to put constraints on the flattening of dark halos in disk galaxies. We observe for this purpose the HI in edge-on galaxies, where it is in principle possible to measure the force field in the halo vertically and radially from gas layer flaring and rotation curve decomposition respectively. In this paper, we define a sample of 8 HI-rich late-type galaxies suitable for this purpose and present the HI observations.Comment: Accepted for publication by Astronomy & Astrophysics. For a higher resolution version see http://www.astro.rug.nl/~vdkruit/jea3/homepage/12565.pd

Regulation of Star Formation Rates in Multiphase Galactic Disks: a Thermal/Dynamical Equilibrium Model

We develop a model for regulation of galactic star formation rates Sigma_SFR in disk galaxies, in which ISM heating by stellar UV plays a key role. By requiring simultaneous thermal and (vertical) dynamical equilibrium in the diffuse gas, and star formation at a rate proportional to the mass of the self-gravitating component, we obtain a prediction for Sigma_SFR as a function of the total gaseous surface density Sigma and the density of stars + dark matter, rho_sd. The physical basis of this relationship is that thermal pressure in the diffuse ISM, which is proportional to the UV heating rate and therefore to Sigma_SFR, must adjust to match the midplane pressure set by the vertical gravitational field. Our model applies to regions where Sigma < 100 Msun/pc^2. In low-Sigma_SFR (outer-galaxy) regions where diffuse gas dominates, the theory predicts Sigma_SFR \propto Sigma (rho_sd)^1/2. The decrease of thermal equilibrium pressure when Sigma_SFR is low implies, consistent with observations, that star formation can extend (with declining efficiency) to large radii in galaxies, rather than having a sharp cutoff. The main parameters entering our model are the ratio of thermal pressure to total pressure in the diffuse ISM, the fraction of diffuse gas that is in the warm phase, and the star formation timescale in self-gravitating clouds; all of these are (in principle) direct observables. At low surface density, our model depends on the ratio of the mean midplane FUV intensity (or thermal pressure in the diffuse gas) to the star formation rate, which we set based on Solar neighborhood values. We compare our results to recent observations, showing good agreement overall for azimuthally-averaged data in a set of spiral galaxies. For the large flocculent spiral galaxies NGC 7331 and NGC 5055, the correspondence between theory and observation is remarkably close.Comment: 49 pages, 7 figures; accepted by the Ap.

Mass models from high-resolution HI data of the dwarf galaxy NGC 1560

We present HI observations performed at the GMRT of the nearby dwarf galaxy NGC 1560. This Sd galaxy is well-known for a distinct "wiggle" in its rotation curve. Our new observations have twice the resolution of the previously published HI data. We derived the rotation curve by taking projection effects into account, and we verified the derived kinematics by creating model datacubes. This new rotation curve is similar to the previously published one: we confirm the presence of a clear wiggle. The main differences are in the innermost ~100 arcsec of the rotation curve, where we find slightly (<~ 5 km/s) higher velocities. Mass modelling of the rotation curve results in good fits using the core-dominated Burkert halo (which however does not reproduce the wiggle), bad fits using the a Navarro, Frenk & White halo, and good fits using MOND (Modified Newtonian Dynamics), which also reproduces the wiggle.Comment: Accepted for publication in MNRAS. 11 pages, 13 figures. High-resolution version available at http://users.ugent.be/~ggianfra/1560_final.pd

Observational Evidence Against Long-Lived Spiral Arms in Galaxies

We test whether the spiral patterns apparent in many large disk galaxies should be thought of as dynamical features that are stationary in a co-rotating frame for > t_{dyn}, as implied by the density wave approach for explaining spiral arms. If such spiral arms have enhanced star formation (SF), observational tracers for different stages of the SF sequence should show a spatial ordering, from up-stream to downstream in the corotating frame: dense HI, CO, tracing molecular hydrogen gas, 24 micron emission tracing enshrouded SF and UV emission tracing unobscured young stars. We argue that such a spatial ordering should be reflected in the angular cross-correlation (CC, in polar coordinates) using all azimuthal positions among pairs of these tracers; the peak of the CC should be offset from zero, in different directions inside and outside the corotation radius. Recent spiral SF simulations by Dobbs & Pringle, show explicitly that for the case of a stationary spiral arm potential such angular offsets between gas and young stars of differing ages should be observable as cross-correlation offsets. We calculate the angular cross-correlations for different observational SF sequence tracers in 12 nearby spiral galaxies, drawing on a data set with high quality maps of the neutral gas HI, THINGS), molecular gas (CO, HERACLES) along with 24 micron emission (Spitzer, SINGS); we include FUV images (GALEX) and 3.6 $\mu$m emission (Spitzer, IRAC) for some galaxies, tracing aging stars and longer timescales. In none of the resulting tracer cross-correlations for this sample do we find systematic angular offsets, which would be expected for a stationary dynamical spiral pattern of well-defined pattern speed. This result indicates that spiral density waves in their simplest form are not an important aspect of explaining spirals in large disk galaxies.Comment: 13 pages, 16 figure

The Fine-Scale Structure of the neutral Interstellar Medium in nearby Galaxies

We present an analysis of the properties of HI holes detected in 20 galaxies that are part of "The HI Nearby Galaxy Survey" (THINGS). We detected more than 1000 holes in total in the sampled galaxies. Where they can be measured, their sizes range from about 100 pc (our resolution limit) to about 2 kpc, their expansion velocities range from 4 to 36 km/s, and their ages are estimated to range between 3 and 150 Myr. The holes are found throughout the disks of the galaxies, out to the edge of the HI; 23% of the holes fall outside R25. We find that shear limits the age of holes in spirals (shear is less important in dwarf galaxies) which explains why HI holes in dwarfs are rounder, on average than in spirals. Shear, which is particularly strong in the inner part of spiral galaxies, also explains why we find that holes outside R25 are larger and older. We derive the scale height of the HI disk as a function of galactocentric radius and find that the disk flares up in all galaxies. We proceed to derive the surface and volume porosity (Q2D and Q3D) and find that this correlates with the type of the host galaxy: later Hubble types tend to be more porous. The size distribution of the holes in our sample follows a power law with a slope of a ~ -2.9. Assuming that the holes are the result of massive star formation, we derive values for the supernova rate (SNR) and star formation rate (SFR) which scales with the SFR derived based on other tracers. If we extrapolate the observed number of holes to include those that fall below our resolution limit, down to holes created by a single supernova, we find that our results are compatible with the hypothesis that HI holes result from star formation.Comment: 142 pages, 55 figures, accepted for publication in the Astronomical Journa