EDP Sciences
Free Access
Issue
A&A
Volume 510, February 2010
Article Number A64
Number of page(s) 13
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/200913261
Published online 10 February 2010
A&A 510, A64 (2010)

Star formation in M 33: the radial and local relations with the gas

S. Verley1,2 - E. Corbelli1 - C. Giovanardi1 - L. K. Hunt1

1 - Osservatorio Astrofisico di Arcetri - INAF, Largo E. Fermi 5, 50125 Firenze, Italy
2 - Dept. de Física Teórica y del Cosmos, Facultad de Ciencias, Universidad de Granada, Spain

Received 7 September 2009 / Accepted 29 November 2009

Abstract
Aims. In the Local Group spiral galaxy M 33, we investigate the correlation between the star formation rate (SFR) surface density, $\Sigma _{\rm SFR}$, and the gas density $\Sigma_{\rm gas}$ (molecular, atomic, and total). We also explore whether there are other physical quantities, such as the hydrostatic pressure and dust optical depth, which establish a good correlation with $\Sigma _{\rm SFR}$.
Methods. We use the H${\alpha }$, far-ultraviolet (FUV), and bolometric emission maps to infer the SFR locally at different spatial scales, and in radial bins using azimuthally averaged values. Most of the local analysis is done using the highest spatial resolution allowed by gas surveys, 180 pc. The Kennicutt-Schmidt (KS) law, $\Sigma_{\rm SFR} \propto \Sigma_{\rm gas}^n$ is analyzed by three statistical methods.
Results. At all spatial scales, with H${\alpha }$ emission as a SFR tracer, the KS indices n are always steeper than those derived with the FUV and bolometric emissions. We attribute this to the lack of H${\alpha }$ emission in low luminosity regions where most stars form in small clusters with an incomplete initial mass function at their high mass end. For azimuthally averaged values the depletion timescale for the molecular gas is constant, and the KS index is $n_{\rm H_2}=1.1$ $\pm $ 0.1. Locally, at a spatial resolution of 180 pc, the correlation between $\Sigma _{\rm SFR}$ and $\Sigma_{\rm gas}$ is generally poor, even though it is tighter with the molecular and total gas than with the atomic gas alone. Considering only positions where the CO J=1-0 line is above the 2-$\sigma $ detection threshold and taking into account uncertainties in $\Sigma_{\rm H_2}$ and $\Sigma _{\rm SFR}$, we obtain a steeper KS index than obtained with radial averages: $n_{\rm H_2}=2.22 \pm 0.07$ (for FUV and bolometric SFR tracers), flatter than that relative to the total gas ( $n_{\rm H_{tot}}=2.59 \pm 0.05$). The gas depletion timescale is therefore larger in regions of lower $\Sigma _{\rm SFR}$. Lower KS indices ( $n_{\rm H_2}=1.46 \pm 0.34$ and $n_{\rm H_2}=1.12$) are found using different fitting techniques, which do not account for individual position uncertainties. At coarser spatial resolutions these indices get slightly steeper, and the correlation improves. We find an almost linear relation and a better correlation coefficient between the local $\Sigma _{\rm SFR}$ and the ISM hydrostatic pressure or the gas volume density. This suggests that the stellar disk, gravitationally dominant with respect to the gaseous disk in M 33, has a non-marginal role in driving the SFR. However, the tight local correlation that exists between the dust optical depth and the SFR sheds light on the alternative hypothesis that the dust column density is a good tracer of the gas that is prone to star formation.

Key words: galaxies: individual: M 33 - galaxies: ISM - Local Group - galaxies: spiral

1 Introduction

The gas-to-star conversion process is one of the most important ingredients for galaxy evolution. The rate at which stars form in a galaxy at a given epoch depends not only on the available gas reservoir but also on the ability of the gas to collapse and fragment. The first seminal papers (Schmidt 1963,1959) related the star formation rate (SFR) to the atomic gas densities using a power law. Comparing the gas density with the number of young stellar objects in the solar neighborhood Schmidt (1959) derived a power law index n=2, i.e.  $\Sigma_{\rm SFR} \propto \Sigma_{\rm gas}^2$. Values of n approximately 1.5 to 2.0 were further confirmed by Guibert et al. (1978), using more precise data on the radial and vertical distributions of the interstellar gas and a variety of young stellar objects.

Papers on external galaxies often use the term Schmidt law or Kennicutt-Schmidt (KS) law to relate the surface density of gas to the SFR per unit surface area, since these quantities are the observables for external galaxies. For half a century, numerous studies have been done on the KS law (see Kennicutt 1998a, for a review). For instance, Kennicutt (1989) studied how the globally averaged SFR in a galaxy traced by H${\alpha }$ emission correlates with the mean atomic, molecular, and total (H I+H2) gas surface densities. For a sample of 15 galaxies they found a good correlation for the atomic and total gas densities (KS index $n_{\rm H_{tot}} = 1.3$ $\pm $ 0.3), but not for the molecular gas density. One decade later, Kennicutt (1998b) studied the relation between the total gas surface density and the H${\alpha }$ SFR density averaged over the entire galaxy, using observations for 61 normal spiral and 36 starburst galaxies. Due to the rather large scatter in the KS relation, the KS index was highly dependent on the method used to fit the data. For the 61 normal galaxies, a least-squares fit on the SFR density yields $n_{\rm H_{tot}} = 1.29$ $\pm $ 0.18, while a bivariate least-squares regression leads to $n_{\rm H_{tot}} = 2.47$ $\pm $ 0.39. Likewise, these two methods to estimate the KS indices for the sample of 36 starburst galaxies lead to values of $n_{\rm H_{tot}} = 1.28$ $\pm $ 0.08 and $n_{\rm H_{tot}} = 1.40$ $\pm $ 0.13, respectively. The better agreement between the two fitting methods mainly reflects the higher gas and SFR dynamic ranges (three orders of magnitude) spanned by the sample of starburst galaxies with respect to the sample of normal galaxies (one order of magnitude). The derivation of the KS index is very sensitive to the method used to fit the data when there is not a wide dynamical range of the variables.

Today several questions regarding the KS law are still open. First, does the SFR surface density correlate better with the total (atomic plus molecular) gas surface density, as suggested by the globally averaged studies, or only with molecular surface density, since stars condense out of molecular gas? Second, which is the best tracer to characterize the SFR? Is H${\alpha }$ a good local current SFR tracer? Does the incompleteness of the IMF at the high mass end in regions where only small clusters form make H${\alpha }$ an unreliable SFR tracer (Corbelli et al. 2009)? The far-UV radiation traces recent SFR averaged over a longer period of time (up to 100 Myr): how well does it correlate with the gas density locally? On which scale and at which wavelength does the infrared radiation trace star formation? How much of the evolved stellar population is contributing to dust heating?

Third, the dependence of the KS law on the spatial scale considered is not yet clear. It is of interest to study at which scale the KS law might break down and how the power law index n changes as we sample smaller and smaller regions. A first step towards a spatially resolved KS law was performed using azimuthally averaged values of SFR and gas densities (e.g. Boissier et al. 2003; Wong & Blitz 2002; Boissier et al. 2007; Martin & Kennicutt 2001). For instance, Wong & Blitz (2002) studied a sample of molecule-rich spiral galaxies and found that the SFR density is more strongly correlated with the H2 surface density than with the total gas surface density. They derived $n_{\rm H_2}=1{-}2$ with an average value of 1.4 using the molecular gas surface density and a radially varying extinction correction. For the total gas surface density the average KS index found by the same authors is steeper, $n_{\rm H_{tot}}\sim1.7$.

As the resolution of the telescopes increased, the KS law has been examined locally at several spatial resolutions (e.g. Bigiel et al. 2008; Calzetti et al. 2005). However, feedback processes linked to SF activity cast doubt on the applicability of the KS law on very small scales. Photodissociation or photoionization radiation from massive stars or stellar winds can locally remove the molecular hydrogen and make H II regions not spatially coincident with peaks of CO emission. This effect is visible in M 33 for example when comparing the H${\alpha }$ emission map of this galaxy with CO J=1-0 line maps (Heyer et al. 2004; Engargiola et al. 2003). M 33 is a galaxy with a generally low extinction, and the recent Spitzer maps of this galaxy at 8 and 24 $\mu $m emission confirm that the displacement is not due to embedded H II regions (Corbelli et al. 2009; Verley et al. 2009).

In order to address these questions, recent space mission data in the UV (GALEX, Gil de Paz et al. 2007) and in the IR (Spitzer, Werner et al. 2004) have been used together with ground-based maps in the H${\alpha }$ recombination line. Such a multiwavelength database enables accurate diagnostics of the SFR at high angular resolution in nearby galaxies. The H I Nearby Galaxy Survey (THINGS, Walter et al. 2008) has provided H I and CO maps for many nearby galaxies at an angular resolution of 500-800 pc, thus offering the possibility to examine the relation between the gas surface density and SFR density locally, on smaller spatial scales. Bigiel et al. (2008) find that the best correlation is between the H2 surface density and the SFR (traced by a combination of the far-UV and IR surface brightness). The average KS index is: 1.0 $\pm $ 0.2. The value $n_{\rm H_2} = 1$ implies that the SFR scales with the mass surface density of molecular gas, i.e., that the timescale and the efficiency of SF (the fraction of gas mass converted into stars) is constant.

The THINGS result of unit slope differs from what has been found in other nearby galaxies by other groups using multiwavelength data. Thilker et al. (2007) derive a KS index of 1.64 (1.87) using $\Sigma_{\rm H_2}$ ( $\Sigma_{\rm H_I+H_2}$) in molecular dominated regions of NGC 7331. In M 51 a KS index of 1.56 has been found by Kennicutt et al. (2007) for the spatially resolved SFR density on 500 pc scales. The index refers to the total gas density and has been derived by taking into account only bright H II regions in the center or along the arms of M 51. The KS index for the H2 surface density in M 51 is slightly lower, 1.37, very similar to that found for azimuthally-averaged quantities in the BIMA sample or previously in M 33 (Heyer et al. 2004; Wong & Blitz 2002). In M 31 (Braun et al. 2009) the best correlation at 113 pc scales is established between the SFR density and the total gas density, with a KS index similar to that found for M 51. The different power law indices found in the literature, $1\le n \le 3.5$, imply that either the SF law differs from galaxy to galaxy (or from region to region such as arms versus whole disk) or that the KS law index is very sensitive to the method used for deriving it (spatial scale, SFR indicator, extinction corrections, CO to H2 conversion, background subtraction, etc.).

Due to its proximity, large angular size, and rather low inclination, the Local Group spiral galaxy M 33 is a unique target to investigate the physics underlying the KS law in a late-type galaxy. M 33 has about 1/3 of its baryons in gaseous form, a small molecular fraction and a low dust abundance (Verley et al. 2009; Corbelli 2003). Early attempts to test the KS law in M 33 considered only the atomic gas because of its low molecular content, given the limited spatial resolution and sensitivity of molecular gas observations. In the early 1970's, using stars and H II region counts together with neutral hydrogen gas, Madore et al. (1974), found a KS index $n_{\rm H_I} \sim 2.35$ $\pm $ 0.26. Newton (1980) noticed that the density of H II regions on scales of ${\sim}300$ pc had a weaker dependence on the H I surface density in the inner regions of the galaxy than in the outer ones. More recently, using the first unbiased census of 12CO J=1-0 line emission in M 33 and the far-infrared emission map provided by IRAS, Heyer et al. (2004) found a strong correlation between the azimuthally averaged SFR density and the average molecular gas surface computed for annular regions 250 pc wide. For the molecular gas surface density, the KS index is $n_{\rm H_2}=1.36$ $\pm $ 0.08, while for the total gas surface density, it is much steeper, $n_{\rm H_{tot}}=3.3$ $\pm $ 0.07. This steeper index is consistent with the molecular index, given the low molecular gas fraction in this galaxy which seems regulated by the balance between the gas pressure (acting on the H2 formation rate) and the dissociation radiation (Heyer et al. 2004; Elmegreen 1993a; Blitz & Rosolowsky 2006; Wong & Blitz 2002).

Another feature of M 33 that makes it interesting for such a study is the inconsistency of the disk stability with the ongoing SF (Kennicutt 1989; Martin & Kennicutt 2001). For many years most of the star-forming disk of M 33 has been known to be stable according to the simple Toomre gravitational stability criterion (Toomre 1964) if only the gas surface density is considered (Elmegreen 1993a). Corbelli (2003) has shown however that the Toomre criterion predicts correctly the size of the unstable star forming region of the M 33 disk when the stellar gravity is considered in addition to that of the gaseous disk. Apparently the stellar disk plays an important role in driving the disk instabilities which trigger SF. Its gravity compresses the gas and it can affect the SFR as well. If, for example, the density, or the free fall time, of a cloud depends on the disk gravity perpendicular to the plane, then the SFR density might not correlate with the gas density alone. Also if SF happens only in self gravitating clouds, then the presence of diffuse molecular material in CO all disk surveys might weaken the expected correlation.

The wealth and quality of newly available data for M 33 has drastically increased during the past years thanks to recent space missions in the UV (GALEX, Gil de Paz et al. 2007) and in the IR (Spitzer, Werner et al. 2004). The high resolution and sensitivity of the recent multiwavelength database for M 33 (Verley et al. 2009,2007) makes it an ideal target for investigating the radial and local relations between the various SFR tracers and gas densities. The spatial resolution of the CO J=1-0 survey (Heyer et al. 2004; Corbelli 2003) of M 33 is 45 arcsec, similar to the resolution achieved by Spitzer at 160 $\mu $m ($\sim$40 arcsec). We can therefore test the KS law at a spatial resolution of 180 pc. To shed light on the physical basis of the scaling relation between SFR and gas density in late type spiral galaxies, we test whether a simple relation between SFR and gas surface density exists in M 33. We do this using several SFR tracers and spatial resolutions applying different fitting methods. Galaxies such as M 33, which do not have a high SFR per unit surface area and have a nearly constant gas surface density, might not establish a tight KS relation; hence the difference between various fitting methods can be large.

This article is the fourth in a series dedicated to the SF in M 33, after Verley et al. (2007, hereafter Paper I), Verley et al. (2009, hereafter Paper II) and Corbelli et al. (2009, hereafter Paper III). It is organized as follows: Sect. 2 presents the data and the methodology used to derive SFRs and gas surface densities. Section 3 compares the neutral gas distribution, atomic and molecular, with the H${\alpha }$ emission line map. In Sect. 4 we study the KS law considering azimuthal averages of different SFR and gas surface density tracers, and in Sect. 5 we examine the local KS law at various spatial resolutions, using two different methods. In Sect. 6 we investigate whether the SFR per unit area establishes a better correlation with the gas volume density, and we outline possible biases when using H${\alpha }$ emission as SFR tracer on a local scale and the CO line luminosities as molecular gas tracer. Our conclusions are summarized in Sect. 7.

2 The data sets and methodology

Here we describe the multiwavelength data set that has been compiled and also the photometric methods and SFR diagnostics. These diagnostics will be used both for the azimuthally-averaged KS law and for the local KS law on a series of increasing spatial scales.

2.1 Ultraviolet and H${\alpha }$ line images

To investigate the continuum ultraviolet (UV) emission of M 33, we use Galaxy Evolution Explorer (GALEX) mission (Martin et al. 2005) data, in particular the data distributed by Gil de Paz et al. (2007). A description of GALEX observations in the far-UV (FUV, 1350-1750 Å) and near-UV (NUV, 1750-2750 Å) relative to M 33 and of the data reduction and calibration procedure can be found in Thilker et al. (2005).

To trace ionized gas, we adopt the narrow-line H${\alpha }$ image of M 33 obtained by Greenawalt (1998). The reduction process, using standard IRAF[*] procedures to subtract the continuum emission, is described in detail in Hoopes & Walterbos (2000). The total field of view of the image is 1.75 $\times$ 1.75 deg2 (2048 $\times$ 2048 pixels of 2 $.\!\!^{\prime\prime}$028 on a side).

2.2 Infrared images

Dust emission can be investigated through the mid- and FIR data of M 33 obtained with the Multiband Imaging Photometer for Spitzer (MIPS) instrument (Rieke et al. 2004; Werner et al. 2004). The complete set of MIPS (24, 70, and 160 $\mu $m) images of M 33 is described in Paper I: the Mopex software (Makovoz & Marleau 2005) was used to gather and reduce the Basic Calibrated Data (BCD). We chose a common pixel size equal to 1 $.\!\!^{\prime\prime}$2 for all images. The images were background subtracted, as explained in Paper I. The spatial resolutions measured on the images are 6'', 16'', and 40'' for MIPS 24, 70, and 160 $\mu $m, respectively. The complete field-of-view observed by Spitzer is very large and allows us to achieve high redundancy and a complete picture of the star forming disk of M 33, despite its relatively large extension on the sky.

2.3 21-cm and millimeter data

Several data sets are available to examine the atomic and molecular gas distributions: these include the Westerbork Radio Synthesis Telescope (WRST) array data (Deul & van der Hulst 1987, 24'' $\times$ 48'' spatial resolution) and Arecibo single dish survey (Corbelli & Schneider 1997, 4' sp. res.). For the molecular gas emission as traced by the CO J=1-0 rotational line, we can use the Berkeley Illinois Maryland Association (BIMA) array data (Engargiola et al. 2003, 13'' sp. res.) and the Five College Radio Astronomy Observatory (FCRAO) single dish data (Heyer et al. 2004; Corbelli 2003, 45'' sp. res.) or the map obtained by combining the two surveys as described by Rosolowsky et al. (2007). The interferometers in general recover less flux than single-dish data, since they tend to filter the diffuse emission from structures much more extended than the primary beam resolution. In fact, BIMA observations of M 33 contain roughly half of the single-dish flux (Engargiola et al. 2003); this implies that single-dish data are potentially more reliable to measure the total gas column density in M 33, even though it is not clear which one establishes a better correlation with the SFR. We shall use the WRST and the FCRAO data for the atomic and molecular gas distributions respectively, which have a comparable spatial resolution. Following Corbelli (2003), the CO measurements were converted to H2 mass column densities ($M_\odot $ pc-2), using the standard conversion factor, X=2.8 $\times$ 1020 cm-2 (K km s-1)-1.

2.4 Star formation tracers

We can use several SF tracers to test the KS law. Different tracers are sensitive to different timescales of SF, and their accuracy to trace SF episodes strongly depends on the SF history and dust content of the galaxy under scrutiny. H${\alpha }$ emission traces gas ionized by massive stars in recent bursts of star formation over timescales of 10 Myr or so. The FUV luminosity corresponds to relatively young stellar populations (${\leq}100$ Myr), so can be considered as complementary to H${\alpha }$ in terms of sensitivity to short timescales.

To convert H${\alpha }$ and FUV emission into SFR, we first correct both H${\alpha }$ and FUV for extinction using the formalism developed by Calzetti (2001). This empirical approach relates the extinction to the TIR and FUV luminosities. The TIR flux (in W m-2 pc-2) is the total IR flux from 3 to 1000 $\mu $m (Dale & Helou 2002) defined as:

\begin{displaymath}%
F{\rm (TIR)} = 10^{-14} \times [19.5\ F_\nu(24) \!+\! 3.3\ F_\nu(70) \!+\! 2.6\ F_\nu(160)],
\end{displaymath} (1)

where $F_\nu(24)$, $F_\nu(70)$, and $F_\nu(160)$ are the MIPS flux densities in Jy pc-2 (see also Paper II). The extinction correction can then be written as:

\begin{displaymath}%
A_{\rm FUV} = 2.5 \times C \times \log\left(\frac{1}{1.68} \times \frac{L(\rm TIR)}{L(\rm FUV)} + 1\right).
\end{displaymath} (2)

The value of C is unity if considering star forming regions. Since we sample star forming regions and the ISM with older populations, we adopt an average value C=0.7 (Paper II). For H${\alpha }$ extinction we shall use $A_{\rm H\alpha} = 0.3$ $\times$ $A_{\rm FUV}$. Following the calibrations given in Paper II, the SFR for the H${\alpha }$ and FUV wavelengths are given by the following formulae:

\begin{displaymath}%
{\rm SFR}({\rm H}\alpha)\ [{M}_\odot\ {\rm yr}^{-1}] = 8.3 \times 10^{-42}\ L({\rm H}\alpha)\ [{\rm erg~s}^{-1}],
\end{displaymath} (3)

\begin{displaymath}%
{\rm SFR}({\rm FUV})\ [M_\odot\ {\rm yr}^{-1}] = 8.8 \times 10^{-44} L(\rm {FUV)\ [erg~s}^{-1}].
\end{displaymath} (4)

Furthermore, Paper II has shown that in a galaxy of low dust content such as M 33, multifrequency diagnostics can help to determine the SFR more accurately. We adopt one ``hybrid'' SFR tracer in addition to H${\alpha }$ and FUV tracers: the bolometric luminosity, approximated by a linear combination of the FUV and TIR (Thilker et al. 2007). This is potentially a reliable SF tracer in regions which are more extended than individual H II regions. It reflects both the unobscured and obscured regions of the galaxy (see Paper III) and is sensitive to a variety of star-formation timescales. To estimate the bolometric luminosity, we used the following equation, given by Thilker et al. (2007):

\begin{displaymath}%
L{\rm (bol)} = \nu_{\rm FUV} L_{\nu, \rm obs} {\rm (FUV)} + (1 - \eta) L{\rm (TIR)}.
\end{displaymath} (5)

We know that the interstellar radiation field which heats the dust has a contribution from the old stellar population. Some of the IR emission may also come from dust heated by evolved stars (AGB), and thus not be directly associated with recent SF episodes (see Paper II). This heating might dominate for dust outside H II regions. The factor $\eta$ accounts for this (Iglesias-Páramo et al. 2006; Bell 2003), and we shall use $\eta = 0.3$. The bolometric SFR was calculated according to Thilker et al. (2007):

\begin{displaymath}%
{\rm SFR}{\rm (bol)}\ [M_\odot\ {\rm yr}^{-1}] = 4.6 \times 10^{-44} L{\rm (bol)}\ [{\rm erg~s}^{-1}].
\end{displaymath} (6)

This equation is based on the calibration of Iglesias-Páramo et al. (2006) derived from Starburst 99 models (Leitherer et al. 1999), with a Salpeter IMF ( $0.1{-}100~M_\odot$), solar metallicity, and a continuous SF.

2.5 Radial profiles and aperture photometry

In Sect. 4 we shall investigate the azimuthally averaged KS law by correlating the azimuthally averaged values of the SFR to gas surface density. Following Heyer et al. (2004), we compute $\Sigma _{\rm SFR}$ and $\Sigma_{\rm gas}$ as the mean values within elliptical annuli spaced by 0.24 kpc and centered on the galaxy center (01$^{\rm h}$33$^{\rm m}$50 $.\!\!^{\rm s}$90, +30$^\circ$39$^\prime$35 $.\!\!^{\prime\prime}$8). The annuli are assumed to be circular rings viewed at an inclination of 54$^\circ$ and with the line of nodes at a position angle of 22.5$^\circ$, thus representing the spatial orientation of the M 33 disk with respect to the line of sight (McConnachie et al. 2006). We performed this analysis on molecular, atomic and total gas tracers (from FCRAO and WRST data) and for each of the SFR tracers as described in the previous section. These azimuthal averages give $\Sigma _{\rm SFR}$ and $\Sigma_{\rm gas}$ as a function of radius, the central radius of each ring. We then correct $\Sigma_{\rm gas}$ and the radiation emitted from the newly formed stars for the disk inclination to obtain face-on values.

In addition to analyzing azimuthal averages, we shall investigate the local KS law using the highest spatial resolution possible for our dataset (see Sect. 5). This is the spatial resolution of FCRAO CO J=1-0 map, similar however to that of 21-cm and 160 $\mu $m maps. In order to compare $\Sigma _{\rm SFR}$ and $\Sigma_{\rm gas}$ locally we first degraded the 24, 70, 160 $\mu $m, FUV, and H${\alpha }$ images to the lowest resolution, i.e. to 45'', corresponding to the spatial resolution of FCRAO molecular map. This corresponds to sizes of 180 pc at the distance of M 33. This is still large enough not to resolve out single molecular clouds: molecular clouds sizes in M 33 are smaller than 100 pc (Rosolowsky et al. 2003). Then, at every position observed in the FCRAO map, we performed aperture photometry on each of the images with a beam of 45'' (full width half maximum). This gives a total of 7664 positions over the disk of M 33. Finally, we converted the photometric results to the appropriate surface density units to examine the local KS law.

To investigate how the KS law varies as a function of the spatial scale we evaluate the KS law on coarser spatial scales by averaging the surface densities over adjacent positions.

3 Comparison of neutral and ionized gas distributions

Figure 1 shows the H I contours (at logarithm levels of 0.55, 0.85, 1.03, 1.15 $M_\odot $ pc-2) overplotted on the H${\alpha }$ image. The H I atomic gas distribution in M 33 is very filamentary, and the H II regions lie on the high surface density H I filaments. But the bright H${\alpha }$ emission knots in M 33 are generally not coincident with the location of the neutral gas peaks. The H I contours follow the H${\alpha }$ emission on large spatial scales, but not locally, on small scales. This agrees with Wright et al. (1972), who found no correlation between the positions of the ten brightest H II regions and the H I density peaks. In the northern spiral arm, south of NGC 604, there is relatively strong H I emission, but little H${\alpha }$ emission associated with it, as already noted by Newton (1980). The H I emission extends very far from the center, much beyond the optical disk, with a slower radial decline than that of the H${\alpha }$ surface brightness. The H I contours even trace the northern plume where the H${\alpha }$ emission is faint. The maxima of the H I emission are in the southern half of the galaxy, exactly on the southern arm, and not in the geometrical center of M 33, where the gas mass density is dominated by molecular gas, and several areas appear to be totally devoid of H I.

\begin{figure}
\par\includegraphics[width=8cm,clip]{13261f1.eps}
\end{figure} Figure 1:

Equal density H I contours (in logarithm: 0.55, 0.85, 1.03, 1.15 $M_\odot $ pc-2) superimposed over the H${\alpha }$ image.

Open with DEXTER

We use the FCRAO CO J=1-0 emission map to trace the H2 molecular gas, with the conversion factor as given in Sect. 2.3. In Fig. 2, the levels of the contours are, in logarithm, 0.90, 1.11, 1.25, 1.36 $M_\odot $ pc-2. Above the detection threshold of FCRAO measurements ( $3\sigma \sim 0.8$ K km s-1), molecular gas is visible only in the inner disk of the galaxy, up to about 3-4 kpc. CO contours closely follow the southern spiral arm of the galaxy. In the north is a displacement of the CO emission peaks relative to the arm as traced by H${\alpha }$. This displacement is also evident in the BIMA data of M 33 (Engargiola et al. 2003). For the GMCs of BIMA, Engargiola et al. (2003) showed that there is clustering of GMCs and H II regions out to a separation of 150 pc. But only 67% of the detected GMCs have their centroid position within 50 pc of an H II region. Along the minor axis, on the western side of the galaxy, CO is absent, while there are bright H II regions visible in H${\alpha }$ emission as well as in the IR. No CO peaks are visible close to IC 133, although CO emission is detected around the bright H II region NGC 604.

\begin{figure}
\par\includegraphics[width=8cm,clip]{13261f2.eps}
\end{figure} Figure 2:

Equal density CO contours (in logarithm: 0.90, 1.11, 1.25, 1.36 $M_\odot $ pc-2) superimposed over the H${\alpha }$ image.

Open with DEXTER

Unlike luminous spiral galaxies, the gas component of M 33 is dominated by atomic gas (Heyer et al. 2004). The atomic-to-molecular fraction is extremely high (1 to 10) compared to molecule-rich galaxies (e.g., 0.05 to 0.4, Wong & Blitz 2002). Small molecular clouds have been found beyond 3-4 kpc (Gardan et al. 2007). The paucity of molecular gas at large radii compared to the atomic gas, which has no radial decline, can be explained as the balance between disk hydrostatic pressure and dissociating radiation (Heyer et al. 2004). Smaller clouds form in the outer disk and in some interarm regions, which might have a higher SF efficiency and a shorter dispersion time than elsewhere. GMCs in M 33 comprise less than 30% of the molecular mass, while in our Galaxy GMCs contain 80% of the molecular gas, and the fraction of molecular gas mass in GMCs decreases with radius.

4 Radial Kennicutt-Schmidt law

We first discuss the radial variations of the KS law by analyzing azimuthally averaged values of the SFR and gas surface density in bins of 0.24 kpc. To trace the SFR, we consider each of the three diagnostics derived in the previous section: extinction corrected H${\alpha }$ (Eq. (3)), extinction corrected FUV (Eq. (4)), and bolometric (Eq. (5)) luminosities. In Table 1 we display the KS indices for all three SFR diagnostics, with regard to each gas component and to the total gas. The KS indices are rather similar if one considers different SFR tracers; only slightly higher KS indices are found when the H${\alpha }$ emission is used to trace the SFR. We shall discuss this effect in Sect. 6.

The KS indices (slopes) are steeper for the total gas surface densities (Arecibo plus FCRAO) ( $n_{\rm H_{tot}} \approx 3.1$). Since the atomic gas, the dominant contributor to the total gas surface density, has a shallow radial distribution, we find a very high value of the KS index when we correlate the SFR with the atomic gas surface density alone. A tighter correlation with a shallower slope is found between the azimuthally averaged SFR and the molecular gas surface density with $n_{\rm H_2} = 1.2$ $\pm $ 0.1, from the inner disk out to 6 kpc. The Robertson & Kravtsov (2008) simulation results are in remarkable agreement with the radially averaged values of the indices found for the atomic and molecular gas in M 33.

The KS law is graphically presented in Fig. 3, for the molecular, atomic, and total gas as a function of the SFR surface densities calculated from the three diagnostics described above. The gas depletion timescales for a constant SF efficiency are represented by dotted lines. One can see that the depletion timescale for the molecular gas is rather constant and between 0.5 and 1 Gyr. On the other hand, the depletion timescales for the atomic and total gas vary widely across the disk of M 33. The fastest depletion timescales (about 0.5 Gyr) occur near the center of the galaxy where the SF is more pronounced; conversely, the lowest depletion timescales are found in the outer parts of M 33 with values reaching roughly 10 Gyr.

Table 1:   Radially averaged KS indices (n) and Pearson correlation coefficients (r).

\begin{figure}
\par\includegraphics[width=17cm,clip]{13261f3.eps}
\end{figure} Figure 3:

Radial KS law, elliptically averaged over bins of 0.24 kpc. The SFR is calculated from three different tracers: extinction-corrected H${\alpha }$ ( left), extinction-corrected FUV ( center), and bolometric ( right). The KS law taking into account only the molecular gas surface density is displayed with plus symbols. The KS law relative to the atomic gas is shown by open triangles. Filled circles corresponds to the KS law when the total gas surface density is considered. The dotted lines represent the gas depletion timescales (cf. labels in Gyr) considering a constant SF efficiency.

Open with DEXTER

Slightly higher indices than what we find using the FUV and bolometric emission have been found by Heyer et al. (2004) by examining radially averaged quantities and a SFR traced by IR emission alone (e.g.,  $n_{\rm H_2} \approx 1.36$). Our KS indices are compatible with those found by Wong & Blitz (2002) if we consider the molecular gas alone. On the contrary, they are much higher than Wong & Blitz (2002) for the atomic and total gas, because their sample of galaxies is dominated by molecular gas, while M 33 has less than 10% of its gas in molecular form.

5 Local Kennicutt-Schmidt law

In this section, we first examine the KS law locally at a spatial scale of 180 pc (45 arcsec), which is the resolution of the FCRAO CO J=1-0 dataset, and then at larger spatial scales. As in the previous section, we calculate the SFR using the three different tracers and evaluate the correlations with the atomic, molecular, and total gas surface densities. Extinction corrections are applied to the H${\alpha }$ and FUV fluxes at each beam position, and the SFR(H${\alpha }$) and SFR(FUV) are computed as described in Sect. 2.4. At all spatial scales, there is much dispersion in the $\Sigma_{\rm SFR} - \Sigma_{\rm gas}$ relations, and the dynamical range of the gas surface density is rather small. As a consequence, the KS indices depend on the fitting method used. We shall give the results from two different fitting methods, which we shall describe in the remainder of this section.

5.1 Local KS law from 180 to 1440 pc: first fitting method

To estimate the indices of the KS law at a spatial resolution of 180 pc, without taking into account outlier positions with very high or low SFR in our map, we used a recursive fit. The first linear fit includes all the 6400 spatial positions defined by the central square (80 $\times$ 80 positions) of the CO map (Heyer et al. 2004). We estimate the dispersion, $\sigma $, of the data with respect to this first fit and remove all points that are lying outside $2\sigma $ before attempting a second fit. We repeat this action recursively until no points lie outside the $2\sigma $ dispersion and then quote the slope of this final linear fit. The value of the threshold has been set equal to $2\sigma $, and changing this threshold could lead to slightly different estimates of the final KS index. Too small a threshold will remove many points at each step and will end up by discarding all the data. On the other hand, if the threshold is too large, the final set of data will remain the same as the initial one, as no points will lie outside the boundaries. Ideally we want to select the value of the threshold which leaves as many points as possible, representing the bulk of the distribution in the final sample while discarding the outliers; a threshold value of $2\sigma $ seems appropriate for this study (see Figs. 4 and 5) and gives stable results. In order to study the effect of the resolution on the KS indices, we used the same method at coarser spatial resolutions, by averaging the surface densities over adjacent positions: 360, 540, 720, 900, 1080, 1260, and 1440 pc. In Table 2, we show the final slopes obtained at each resolution for the three SFR tracers, with respect to the molecular, atomic, and total gas. Also in Table 2, we show the Pearson coefficients obtained for the raw samples of points. Clear trends appear: the most significant is that the correlation between $\log \Sigma_{\rm SFR}$ and $\log \Sigma_{\rm gas}$ gets better and better when the resolution gets coarser. This was expected due to the stochastic effects which take place at small spatial scales. Although less systematic, the Pearson coefficients are in general higher when considering the molecular gas, with respect to the total and atomic gas. This latter shows the weakest correlation, independently of the SFR tracer used for the correlation.

\begin{figure}
\par\includegraphics[width=13cm,clip]{13261f4.eps}
\end{figure} Figure 4:

Local KS law involving the extinction-corrected H${\alpha }$ luminosity as SFR tracer and molecular gas surface density at four spatial resolutions. The black points (inbetween the $2\sigma $ boundaries, evidenced by the two dashed lines) are the ones considered for the last iteration of the first fitting method. The solid line is the resulting best fit considering these positions. The grey points were discarded during the various iterations of the fit. The Kennicutt (1998b)'s index is depicted by a dash-dot line (slope of 1.4). The dotted lines represent the times of gas depletion (cf. labels in Gyr in the lower right panel) considering a constant SF efficiency.

Open with DEXTER

In Figs. 4 (molecular gas) and 5 (total gas) we display the final slope obtained (heavy line) with the final samples of positions included in the fit, for two cases. Positions discarded (shown in grey) during the several steps of the iteration process lie outside the $2\sigma $ boundaries (depicted by dashed lines). One can see that the predominant influence on the KS indices is the nature of the gas that is taken into account: the molecular gas gives values of $n_{\rm H_2}$ between 1 and 2, while atomic and total gas exhibit higher values: $n \approx 2{-}4$. As in the azimuthally averaged relation, the H${\alpha }$ SFR tracer displays higher KS indices than the FUV and bolometric SFR tracers; this will be discussed in the next section.

Indices get steeper on average as the resolution gets more coarse, i.e. as we average quantities over larger areas. For the total gas density the KS indices obtained radially are even higher than the ones obtained for the coarser local resolution. Steeper indices for radially averaged densities are also found by Wong & Blitz (2002) and Thilker et al. (2007). For the atomic gas the local KS indices are significantly lower, because the H I has a shallow radial falloff, while its distribution in the disk is filamentary. Hence local variations of the atomic gas are more significant than radial variations. This is not the case for the molecular gas, which declines radially with a scalelength of 2 kpc, similar to that of the H${\alpha }$ and FUV emission (Paper II). The local KS indices relative to the molecular gas are similar to those found for azimuthal averages in radial bins.

The number of iterations needed to reach the final slope and the data dispersion generally decreases as a coarser spatial resolution is considered (in Table 2 the value of the linear Pearson coefficient r with the original distributions of points, i.e. before the first iteration, is given). This is because stochastic effects weaken the relation between SFR and gas densities at small spatial scales. In the limit of H II region sizes the correlation could disappear altogether, because massive stellar winds and supernovae explosions could remove the molecular gas and quench subsequent SF. Using the combined BIMA-FCRAO dataset (Rosolowsky et al. 2007) at a spatial resolution of 30 arcsec (120 pc) around selected H II regions, we find in fact an even weaker correlation between $\Sigma _{\rm SFR}$ and $\Sigma_{\rm H_2}$ than over the 180 pc spatial scale.

\begin{figure}
\par\includegraphics[width=13cm,clip]{13261f5.eps}
\end{figure} Figure 5:

Local KS law involving the extinction corrected FUV luminosity as SFR tracer and total gas surface density at four spatial resolution. The black points (in between the $2\sigma $ boundaries, evidenced by the two dashed lines) are the ones considered for the last iteration of the first fitting method. The solid line is the resulting best fit considering these positions. The grey points were discarded during the various iterations of the fit. The Kennicutt (1998b)'s index is depicted by a dash-dot line (slope of 1.4). The dotted lines represent the times of gas depletion (cf. labels in Gyr in the lower right panel) considering a constant SF efficiency.

Open with DEXTER

5.2 Relative errors: second fitting method

One caveat of the procedure described in the previous paragraph is that we have considered all positions in the map to allow smoothing and averages over larger areas. This implies that even positions where the molecular gas detection was not above the noise level have been considered. Thus, the reliability of the results on the smallest spatial scale relative to the large-scale ones could be compromised because of low signal-to-noise; as the resolution is degraded, the signal-to-noise ratio increases because we are averaging over larger areas.

Table 2:   Local KS law indices (n) obtained using the first fitting method for the three SFR tracers.

We now fit the $\log \Sigma_{\rm SFR} - \log \Sigma_{\rm gas}$ local relation at the lowest spatial resolution (180 pc), taking into account the errors in the determination of the gas surface density as well as in the photometry and extinction correction for the determination of the $\Sigma _{\rm SFR}$. We use the GALEX FUV map to determine the photometric errors[*]. In the smoothed GALEX FUV image the pixel-to-pixel 1-$\sigma $ noise varies between 2.9 and 8.7 $\times$ 10-6 counts s-1, and we shall consider its maximum value. The maximum large-scale sky 1-$\sigma $ variation of sky uncertainty is 3 $\times$ 10-4 counts s-1. The sky noise is much smaller than the large-scale sky variation and can be neglected. To the photometric errors we add the errors for extinction corrections. Following Calzetti (2001), we compute the extinction in the GALEX FUV band as in Eq. (2). The largest source of uncertainty in the formula are not the photometric errors on TIR and FUV luminosities, but the bolometric correction factor (1.68). Since the sampled areas contain a mix of young and old stellar populations, it is not clear how much of the TIR luminosity is due to the young stellar population powering the FUV emission. We estimate the uncertainties in $A_{\rm FUV}$ to be on average 20%. The individual errors on CO flux estimates can be computed from the rms noise in individual spectra. To these we add the error on the atomic hydrogen surface density, which is on the order of 0.7 $M_\odot $ pc-2 (as from map noise), if the total gas surface density is considered.

We now consider all areas in M 33 for which the CO flux measured by the FCRAO survey is above $2\sigma $ uncertainty. We analyze the possible linear correlation between the $\log$ of molecular gas surface density and the $\log$ of SFR per unit surface as traced by the FUV emission corrected for extinction. A linear least square method with errors both in the gas surface density and in the SFR (subroutine fitexy in Press et al. 1992) gives a slope of 2.22 $\pm $ 0.07. The Pearson linear correlation coefficient is 0.42. The error on the slope is determined considering extinction correction error variations: 20% on average ${\pm}5$%. The best-fit regression slopes and correlation coefficients relating SFR and gas content in M 33 for the 180 pc regions using this method are given in Table 3. For comparison, the slopes obtained with the ordinary-least-square (OLS) method advocated by Isobe et al. (1990) are also listed in Table 3. This method does not take into account either uncertainties or outliers, so is less sophisticated than either of the two methods described previously. Indeed, the OLS slopes are shallower than either of the other two estimates; in the case of poor correlations with many data points, not considering properly errors or outliers can give too much relative weight to statistically insignificant data. In Fig. 6 we show the $\log \Sigma_{\rm H_{tot}} - \log \Sigma_{\rm SFR}$ relation with the relative best fitting linear relation (KS index = 2.59 $\pm $ 0.05). The correlation with total gas is of similar quality to that with the molecular gas alone, since both Pearson linear correlation coefficients are $\sim$0.43.

We now consider the correlation between the total gas surface density and the SFR per unit surface, including points which have the CO brightness below the detection threshold. For these we shall consider only the H I surface density. The correlation coefficient does not change, but the regression is shallower with a slope of 1.7 $\pm $ 0.2. But caution is called for in interpreting such a shallow slope as real or claiming that there are areas in the galaxy where the surface density decreases but the SFR does not (bi-modal distribution). The non-detection of molecular emission generates artificial tails with zero slopes at low surface densities and flattens the average slope. It is more feasible to consider positions without a reliable CO detection only when the corresponding molecular surface density threshold is well below the H I surface density, but there are not many positions in our data which satisfy this condition.

Table 3:   Pearson coefficients (r) and slopes (n) derived from the ``relative errors'' method for correlations between the two quantities shown in the first column.

\begin{figure}
\par\includegraphics[width=8cm,clip]{13261f6.eps}
\end{figure} Figure 6:

Local KS law involving the extinction corrected FUV luminosity as SFR tracer and the molecular gas surface density in  a), the total gas surface density in  b). The spatial resolution is 180 pc. The best fitting linear relation leads to a KS index of 2.22 $\pm $ 0.07 in a) and of 2.59 $\pm $ 0.05 in  b) shown by the heavy lines.

Open with DEXTER

6 Different tracers of the SFR and of the likelihood of the gas to form stars

In this last section we analyze whether the SFR correlates better with other physical quantities of the ISM, such as the hydrostatic pressure or gas volume density, rather than with the total gas surface density. We also discuss the uncertainties related to the use of the H${\alpha }$ as SFR tracer, which at first seems more appropriate than the FUV or bolometric surface brightness in tracing the most recent star formation. And in the last paragraph we show the tight correlation between the SFR and the dust optical depth. We briefly discuss whether the dust optical depth can trace the gas prone to form stars in the M 33 disk better than the CO J=1-0 line intensity.

6.1 The role of the stellar disk and the dependence of the SFR on the gas volume density

The widely used KS law relates the SFR per unit surface with the gas surface density. The SFR depends upon the amount of molecular gas able to collapse, fragment and form stars on a timescale  $t_{\rm SF}$. Often the relevant timescale has been referred to as the cloud collapse or free-fall timescale, which is inversely proportional to the square root of the gas volume density in the cloud: $t_{\rm SF}\propto t_{\rm ff} \propto \rho_{\rm H_2}^{-0.5}$. If  $\rho_{\rm H_2}$ scales with some power of the gas surface density or is constant, then it is conceivable that the SFR correlates with the surface density of the self-gravitating gas,  $\Sigma_{\rm H_2}$, to some power. However, since the free-fall time is much shorter than the molecular cloud formation timescale (e.g., Engargiola et al. 2003; Ward-Thompson et al. 1994; Krumholz et al. 2009; McKee & Tan 2003), the relevant time for star formation could be linked to the interstellar medium out of which the clouds form. For example Koyama & Ostriker (2009) have recently shown using numerical simulations of disks that the ISM structure plays a primary role in determining the actual SFR in galactic disks. To test this hypothesis in M 33, we relate the star formation timescale not to the molecular cloud free-fall time, but to their formation timescale. If  $t_{\rm SF}$ is the time needed to form self-gravitating molecular clouds out of the diffuse interstellar gas, $t_{\rm SF}$ will depend on the density of the ISM  $\rho_{\rm ISM}$, i.e.  $t_{\rm sf} \propto \rho_{\rm ISM}^{-0.5}$ (Koyama & Ostriker 2009). The value of $\rho_{\rm ISM}$ in the disk scales with the hydrostatic pressure, which in galaxies with a low gaseous content such as M 33 depends on the stellar disk surface density (Corbelli 2003). An important finding in M 33 is that both the outer SF threshold radius and the abundance of molecules can be explained if one includes the contribution of the stellar disk in the hydrostatic pressure equation (Corbelli 2003; Elmegreen 1993b).

In an attempt to improve the local correlation between the gas and the SFR in M 33, we consider the gravitational compression of the gas due to the stellar disk. If the hydrostatic pressure P sets the gas vertical scale height $h_{\rm g}$ of the disk, the volume density  $\rho_{\rm ISM}$ can be evaluated as  $P/c_{\rm g}^2$ or following Corbelli (2003) as:

\begin{displaymath}%
\rho_{\rm ISM} = {\Sigma_{\rm H_{tot}} \over 2 h_{\rm g}},
\end{displaymath} (7)

where

\begin{displaymath}%
h_{\rm g}={c_{\rm g}\over \pi G}\left({\Sigma_{\rm H_{tot}}\over c_{\rm g}} +
{\Sigma_{\rm s}\over c_{\rm s}}\right)^{-1}
\end{displaymath} (8)

and $\Sigma_{\rm s}$, $c_{\rm s}$and $c_{\rm g}$ are respectively the stellar surface density, the velocity dispersion of the stars and the gas perpendicular to the disk. If  $\rho_{\rm ISM}$ determines the SF timescale through the cloud formation timescale (or through $t_{\rm ff}$ if $\rho_{\rm H_2}$ scales with $\rho_{\rm ISM}$), we expect the following relation to hold:

\begin{displaymath}%
\Sigma_{\rm SFR} \propto {\Sigma_{\rm H_2} \over t_{\rm SF}} \propto
{\Sigma_{\rm H_2}\over \rho_{\rm ISM}^{-0.5}}\cdot
\end{displaymath} (9)

We shall consider a gas velocity dispersion of 6 km s-1 and a ratio of gas-to-stellar velocity dispersion equal to 0.3 (Corbelli & Walterbos 2007). By looking at the correlation between $\Sigma _{\rm SFR}$ and $\Sigma_{\rm H_2}/\rho_{\rm ISM}^{-0.5}$ we find that this is tighter than that between $\Sigma _{\rm SFR}$ and the gas surface density alone. The Pearson linear correlation coefficient is 0.62, and the slope of the linear fit is 1.16 $\pm $ 0.04 in the $\log - \log$ plane (see Fig. 7), close to linear as expected by the above equation. It seems clear that the gravity of the stellar disk, dominant in M 33, plays a major role in regulating the SFR. Variations of the gas surface density alone are insufficient to explain the scatter in  $\Sigma _{\rm SFR}$ and imply a steep KS index as well as a loose correlation when considering only the gas surface density.

Finally, it is noteworthy that the Pearson linear correlation coefficient is even higher, r=0.71, when considering the $\log \Sigma_{\rm SFR} - \log \rho_{\rm ISM}$ relation. The slope of the linear correlation, shown in Fig. 7, is 1.07 $\pm $ 0.02. The tighter relation between $\Sigma _{\rm SFR}$ and $\rho_{\rm ISM}$ (or the hydrostatic pressure since $\rho_{\rm ISM}\propto P$) is well known. It can be interpreted as follows: suppose that  $t_{\rm SF}$ is proportional to $\rho_{\rm ISM}^{-0.5}$ because the cloud formation timescale is proportional to  $\rho_{\rm ISM}^{-0.5}$ and regulates the SF. If the molecular gas surface density (or volume density if SF is confined into a layer of constant scale height), which is in the form of bound, self-gravitating units, is proportional to  $\rho_{\rm ISM}^{0.57}$, the expected relation is that we find $\Sigma_{\rm SFR} \propto \rho_{\rm ISM}^{1.07}$. This would be the case for a constant star formation efficiency.

The surface density of molecular gas in the form of bound, self-gravitating units can be different from $\Sigma_{\rm H_2}$ since $\Sigma_{\rm H_2}$ includes the non negligible contribution of diffuse molecular gas. Unfortunately there are no detailed surveys of the molecular gas in M 33 available as yet which would be able to constrain the fraction of gas in the form of bound, self-gravitating units. Only results of GMCs surveys are available, but the molecular mass spectrum in M 33 is steeper than -2 and hence dominated by molecular clouds of small mass (Blitz & Rosolowsky 2005).

\begin{figure}
\par\includegraphics[width=8cm,clip]{13261f7.eps}
\end{figure} Figure 7:

The left-hand panel shows the linear correlation between log  $\Sigma _{\rm SFR}$ and log  $\Sigma _{\rm H_2}/\rho ^{-0.5}_{\rm ISM}$ (in units of $M_\odot $/pc-2 cm-1.5). The Pearson linear correlation coefficient is 0.62, and the slope of the linear fit is 1.16 $\pm $ 0.04. The slope is derived taking into account uncertainties on both axes, but they are not shown in the graph for clarity. On the right-hand panel we show the linear correlation between  $\log \Sigma_{\rm SFR}$ and $\log \rho_{\rm ISM}$ with the Pearson coefficient 0.71.

Open with DEXTER

\begin{figure}
\par\includegraphics[width=8.28cm,clip]{13261f8a.eps}\hspace*{-1.3mm}
\includegraphics[width=7cm,clip]{13261f8b.eps}
\end{figure} Figure 8:

$\Sigma _{\rm SFR}$ derived from H${\alpha }$ plotted against $\Sigma _{\rm SFR}$(FUV) ( left panel) and $\Sigma _{\rm SFR}$(bol) ( right panel). Only the 180 pc regions are shown. The dotted line corresponds to equality in the two SFRs.

Open with DEXTER

6.2 Is H${\alpha }$ luminosity a good SFR tracer on a local scale?

We compare the three different SFR tracers in Fig. 8, which shows $\Sigma _{\rm SFR}$(H${\alpha }$) plotted against $\Sigma _{\rm SFR}$ derived from the FUV (left panel) and bolometric luminosity (right) (see Sect. 2). At high SFR densities, H${\alpha }$ and FUV give quite similar values, but H${\alpha }$ underpredicts SFR relative to FUV at low SFR. SFR as inferred from bolometric luminosity is almost always larger than that from H${\alpha }$, and at low SFR densities the effect is quite strong.

With the H${\alpha }$ emission as SFR indicator, the correlation of $\Sigma _{\rm SFR}$ with $\Sigma_{\rm gas}$ or  $\rho_{\rm ISM}$ becomes looser, and steeper indices are found. This is to be expected when one samples regions, as we do, with a total bolometric luminosity below 1039 erg s-1, and we explain why. The bolometric luminosity of a single cluster or of an ensemble of clusters measures the mass of the clusters only when the IMF is fully populated. Hence, when $L_{\rm bol}$ is comparable to or lower than the luminosity of the most massive star in the cluster, the IMF cannot be fully sampled up to its high mass end.

As Corbelli et al. (2009) pointed out, there are two types of models for populating the stellar clusters when the IMF is incompletely sampled. Either the IMF is truncated to a limiting mass which depends on the cluster mass, or the IMF is incomplete and stochastically sampled, but maintains its original shape and completeness up to its high mass limit over the whole galaxy. For both cases the $L_{\rm bol}$/ $L_{\rm H\alpha}$ ratio rises when the most massive stars are lacking because the recombination-line luminosity, $L_{\rm H\alpha}$ has a steeper dependence on the stellar mass than the bolometric luminosity does. For the first model, the truncation model, a deviation from a simple scaling law between log  $L_{\rm bol}$ and log  $L_{\rm H\alpha}$ is expected when  $L_{\rm bol}$ is on the order of the luminosity of the most massive star. For 120 $M_\odot $ this is about $7\times 10^{39}$ erg s-1. For the second model, the randomly sampled cluster model, the stochastic character of the IMF implies that low luminosity clusters can be made either by populating the IMF up to a certain stellar mass, or by just one single star, or by something in between. The presence of outliers makes the average  $L_{\rm bol}$/ $L_{\rm H\alpha}$ deviate less from its constant value at high luminosities (where the IMF is fully populated). In practice it is only for $L_{\rm bol}\le 1.5$ $\times$ 1039 erg s-1 that $L_{\rm bol}$/ $L_{\rm H\alpha}$ drops dramatically. In this regime  $L_{\rm H\alpha}$ is no longer an adequate measure of the cluster mass, and for this model a larger scatter is expected in the bolometric luminosity - cluster mass relation.

\begin{figure}
\par\includegraphics[width=8cm,clip]{13261f9.eps}
\end{figure} Figure 9:

The theoretical birthline (Paper III) which predicts a decreasing H${\alpha }$ flux for decreasing luminosity of young stellar clusters due to incompleteness of the IMF is shown with filled squares for the randomly sampled cluster model. The filled circles are the model prediction for clusters with a truncated IMF. Cross symbols refer to M 33 data, areas observed at 180 pc resolution.

Open with DEXTER

Figure 9 shows the theoretical $L_{\rm bol}-L_{\rm H\alpha}$ relation, referred to as the cluster birthline. The filled square symbols correspond to the randomly sampled cluster model, and filled circles indicate the prediction for clusters modelled with a truncated IMF (see above and Paper III). Cross symbols in Fig. 9 show the data for the 180 pc regions in M 33. The luminosities of these regions are clearly sufficiently low to enter the regime of an incomplete IMF. These data cannot be used to judge which of the two cluster models applies, because regions 180 pc wide may contain more than one cluster and not necessarily young (see Corbelli et al. 2009, to address this issue correctly).

The fast decrease of the H${\alpha }$ luminosity in regions of low luminosity, which do not contain massive clusters, implies that using a constant H${\alpha }$ luminosity - SFR conversion factor, the inferred SFR will be much lower than the effective one. If the inferred SFR in low density regions is lower, the index of the KS relation will be artificially higher. This explains the steeper indices we find when using H${\alpha }$ as the SFR indicator in the local study of M 33. We conclude that the lack of massive stars in low luminosity regions makes H${\alpha }$ an unreliable tracer of the SFR.

The birthline is the line where very young stellar clusters lie. Aging or leakage of ionizing photons bring the clusters above the birthline. Aging decreases $L_{\rm H\alpha}$ faster than  $L_{\rm bol}$ and moves $L_{\rm bol}$/ $L_{\rm H\alpha}$ to higher values more or less vertically above the birthline (see Paper III). Also leakage of ionizing photons from H II regions raises the $L_{\rm bol}$/ $L_{\rm H\alpha}$ values above the birthline. All our sampled regions (M 33 areas at 180 pc resolution) lie above the birthline (cross symbols) as they should, implying that we have not underestimated their bolometric luminosity. However, most of the SF regions we sample are well above the birthline, and this suggests that there is a mix of young and aged H II regions in each area and/or that there is leakage of ionizing photons. The aging and leakage processes and the approximate formula used to estimate $L_{\rm bol}$ and extinction correction to $L_{\rm H\alpha}$ explains the observed scatter in the $L_{\rm bol}-L_{\rm H\alpha}$ relation. M 33 is known to have a high diffuse fraction of H${\alpha }$ emission due to the leakage process from individual star-forming regions (Hoopes & Walterbos 2000). But also in the continuum UV radiation M 33 has high diffuse fractions (Thilker et al. 2005), and if radiation escapes from the disk from these diffuse regions or from individual star-forming sites, both non-ionizing and ionizing photons will contribute to this leakage. Hence although leakage can explain a possible steepening of the $\Sigma _{\rm SFR}$ (from all tracers) at low gas surface densities, it does not explain the non-linear scaling between $L_{\rm bol}$ or  $L_{\rm FUV}$ and  $L_{\rm H\alpha}$.

6.3 The correlation between the SFR and the dust optical depth

We now use the dust opacity at 160 $\mu $m, $\tau_{160}$ derived as in Paper II, as an unbiased indicator of the gas column density. In M 33 the metallicity gradient is very shallow (Magrini et al. 2007b; Rosolowsky & Simon 2008), and hence it is probable that the gas-to-dust ratio does not change much going radially outwards. The declining behavior of the dust-to-gas ratio found in Paper II when using the 21-cm emission line brightness to infer the H I column density and the CO J=1-0 emission line brightness to infer the H2 column density, can then be due to an underestimation of the gas surface density at large galactocentric radii. In particular the CO-to-H2 conversion factor can be different from the value assumed in the M 33 FCRAO survey (2.8 $\times$ 1020 cm-2 (K km s-1)-1), and can vary across the M 33 disk according to local variations of the metal abundance and molecular clouds properties. Many papers in the literature (e.g. Leroy et al. 2009,2007, and references therein) show that the CO-to-H2 conversion factor is effectively higher in low metallicity galaxies and can vary even in single molecular cloud complexes according to the cloud self-shielding conditions. As an alternative there could be opaque H I gas unaccounted for by the usual optically thin assumption for the 21-cm line emission, as found in M 31 by Braun et al. (2009).

\begin{figure}
\par\includegraphics[width=8cm,clip]{13261f10.eps}
\vspace{4mm}\end{figure} Figure 10:

Correlation between the SFR (FUV) surface density and the dust opacity at 160 $\mu $m. Cross symbols are data for areas in M 33 which have a 180 pc wide diameter; only those points are shown with the CO luminosity above 2$\sigma $ detection threshold.

Open with DEXTER

We show in Fig. 10 the correlation between the SFR surface density and the dust opacity. The Pearson linear correlation coefficient is much higher (r=0.81) when considering the correlation in the $\log \Sigma_{\rm SFR} - \log \tau_{160}$ plane than when using the brightness of the 21-cm and CO J=1-0 lines to derive $\log \Sigma_{\rm SFR} - \log \Sigma_{\rm H_{tot}}$ relation (r=0.43). As shown in Fig. 10 the distribution seems bimodal; there is an inflection in the slope around $\log \tau_{160}\approx-4$. The average slope (given in Table 3) and the correlation coefficient do not depend on whether we include positions where the CO luminosity is below the 2$\sigma $ limit. We have assumed a 20% error on  $\tau_{160}$ to determine n, but given the high Pearson coefficient, the slope is not much dependent on the relative errors.

Although the correlation $\log \Sigma_{\rm SFR} - \log \tau_{160}$ is quite tight, some more work is needed to check the assumptions made and better address the implications of this result. For example even though the average radial variation of the metallicity is small ( ${\sim}-0.4$ dex), metallicity rises in the center (U et al. 2009), and at a given radius the variation in the oxygen abundance can be almost one order of magnitude ( ${\sim}-0.8$ dex) (Magrini et al. 2007a). Right now the relative errors on metal abundances are large, but if future measurements confirm this scatter, the assumption of a constant dust-to-gas ratio, implicit in the correlation, may not be appropriate. Also, the actual estimate of the optical depth is based on the intensity of the emission at 160 $\mu $m and on the dust temperature derived using the 70 and 160 $\mu $m emission. If the emission at 70 $\mu $m is due to dust heated not only by the interstellar radiation field but also by radiation from star-forming regions, then we would expect a higher colour temperature. This would give a larger thermal intensity, and thus we would derive a smaller  $\tau_{160}$, anti-correlated with the SFR. High-resolution IR observations at longer wavelengths, which will be available in the near future with Herschel, will allow a more accurate estimate of the dust optical depth. The resulting correlation coefficient involving $\log \tau_{160}$ and  $\log \Sigma_{\rm SFR}$ could have a higher value than the one derived here.

7 Summary and conclusions

We test the KS law in the Local Group spiral galaxy M 33, in azimuthally averaged areas (240 pc wide) and locally from $\sim$1.5 kpc spatial resolution down to a resolution of 180 pc. Starting from a multiwavelength set of observations, we used the H${\alpha }$, FUV, and bolometric (FUV+TIR) luminosities to estimate the SFR. For gas surface density we consider the molecular, atomic, and total gas phases. We use extinction-corrected SFR even though extinction marginally affects the observed properties of the optical or ultraviolet emissions in M 33 because of the rather low dust content of the galaxy. The most important results are summarized below:

  • At every spatial scale we find that the H${\alpha }$ KS indices are always higher than the FUV and bolometric ones, and we explain this as due to the lack of H${\alpha }$ emission in low luminosity regions where most of stars form in small clusters with an incomplete initial mass function at their high mass end. We use the cluster birthline to support this, which implies a non-linear relation between the H${\alpha }$ and bolometric luminosities. The birthline also shows that most regions, even at the highest spatial resolution, contain a mixture of ages which on average are smaller than 10 Myr.

  • For azimuthally averaged values, the depletion timescale for the molecular gas is radially constant and the KS index is 1.1 $\pm $ 0.1 for both FUV and bolometric SFR tracer, lower than that found by Heyer et al. (2004) using IR SFR tracers alone. The depletion time for the molecular gas is relatively constant, with a value of about 1 Gyr. The correlations with the molecular and total gas are tighter than with the H I gas, despite the fact that most of the gaseous mass of M 33 is in H I form. Denser filaments of neutral gas are found where star formation takes place, but their H I surface density does not correlate with the SFR. The KS radial index for the total gas surface density is 2.9 $\pm $ 0.2. The depletion time for the total gas density widely varies from $\sim$0.5 Gyr in the center of the galaxy to about 10 Gyr in the outer parts.

  • Locally the dispersion in the SFR-gas density relation per unit surface is high, and results are very sensitive to the statistical method used to fit the data. The scatter in the KS relation increases as the spatial resolution increases. At 180 pc resolution, fitting the $\log \Sigma_{\rm SFR} - \log \Sigma_{\rm H_2}$ data to a straight line using a recursive fitting method which removes outlier points, we obtain the same the KS index than in the radial average analysis (1.1). There is a slight general tendency for the local KS laws to show higher slopes for coarser spatial resolutions. As for radial averages, the local KS index for the molecular gas is much lower than that for the total gas.

  • For our finest resolution (180 pc) we also derive the best fitting straight line to the data in the $\log - \log$ plane, taking into account the errors in the determination of the gas surface density as well as the errors in the determination of the SFR as traced by FUV emission. Considering only the positions where the CO detection is above 2$\sigma $ noise, a bivariate regression gives the KS index for the molecular gas of 2.22 $\pm $ 0.07 with a Pearson linear correlation coefficient of 0.42. We believe that given the large scatter in the $\log \Sigma_{\rm SFR} - \log \Sigma_{\rm H_2}$ relation at the highest spatial resolution examined in this paper, this method gives more robust result. This KS index is higher than the one we find using the recursive fit and all positions in our maps. It is also higher than what the ordinary-least-square method gives for the same set of data (no uncertainties considered). The bivariate regression gives a slope for the $\log \Sigma_{\rm SFR} - \log \Sigma_{\rm H_{tot}}$ relation only slightly higher than for the molecular gas alone: 2.64 $\pm $ 0.07.

  • Given the rather poor correlations between the gas and the SFR per unit surface, we analyze whether the SFR correlates with other physical quantities. A good correlation is found with the hydrostatic pressure, i.e. with the interstellar medium volume density considering a constant sound speed. This implies that the stellar disk, gravitationally dominant with respect to the gaseous disk in M 33, plays a major role in driving the SFR. The slope of the correlation is close to unity, suggesting that the SFR per unit area (or per unit volume if the thickness of the SF disk does not vary) is proportional to the ISM volume gas density. Since the correlation is tighter, its slope is less dependent on the statistical method used.

  • There is a good correlation between the dust optical depth at 160 $\mu $m and the local SFR density. This can be interpreted in terms of a $\Sigma_{\rm SFR} - \Sigma_{\rm gas}$ relation if the dust opacity is used as an unbiased indicator of the surface density of gas prone to star formation. This might be the case for example if the gas-to-dust ratio were constant in M 33 while the CO-to-H2 conversion factor were not, or if there were opaque H I gas unaccounted for by the usual optically thin assumption for the 21-cm line emission. The tightness of the $\log \Sigma_{\rm SFR} - \log \tau_{160}$ correlation and the shallower slope compared to the $\log \Sigma_{\rm SFR} - \log \Sigma_{\rm gas}$ relation would then imply that the KS law in low luminosity galaxies still holds, but the surface density of the gas in the process of forming stars has a wider dynamical range than shown by the intensity of the 21-cm and CO J=1-0 lines. Still, some caution is called for in interpreting this result because of some of the assumptions made in deriving it. A more detailed analysis is required together with a more accurate determination of the dust optical depth, which will be possible in the near future thanks for example to the Herschel satellite data.

Acknowledgements
We would like to thank Rene Walterbos for providing us the H${\alpha }$ image of M 33, R. Kennicutt and G. Helou for interesting discussions, and the anonymous referee for comments that helped to improve our work. The work of S. V. was supported by a INAF - Osservatorio Astrofisico di Arcetri fellowship. The Spitzer Space Telescope is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This research has made use of the NASA/IPAC Extragalactic Database, which is operated by JPL/Caltech, under contract with NASA.

References

Footnotes

... IRAF[*]
IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation.
... errors[*]
It is not possible to determine the errors for the H${\alpha }$ map since the map we have is in emission measure units and it is not possible to recover the original errors in counts.

All Tables

Table 1:   Radially averaged KS indices (n) and Pearson correlation coefficients (r).

Table 2:   Local KS law indices (n) obtained using the first fitting method for the three SFR tracers.

Table 3:   Pearson coefficients (r) and slopes (n) derived from the ``relative errors'' method for correlations between the two quantities shown in the first column.

All Figures

  \begin{figure}
\par\includegraphics[width=8cm,clip]{13261f1.eps}
\end{figure} Figure 1:

Equal density H I contours (in logarithm: 0.55, 0.85, 1.03, 1.15 $M_\odot $ pc-2) superimposed over the H${\alpha }$ image.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8cm,clip]{13261f2.eps}
\end{figure} Figure 2:

Equal density CO contours (in logarithm: 0.90, 1.11, 1.25, 1.36 $M_\odot $ pc-2) superimposed over the H${\alpha }$ image.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=17cm,clip]{13261f3.eps}
\end{figure} Figure 3:

Radial KS law, elliptically averaged over bins of 0.24 kpc. The SFR is calculated from three different tracers: extinction-corrected H${\alpha }$ ( left), extinction-corrected FUV ( center), and bolometric ( right). The KS law taking into account only the molecular gas surface density is displayed with plus symbols. The KS law relative to the atomic gas is shown by open triangles. Filled circles corresponds to the KS law when the total gas surface density is considered. The dotted lines represent the gas depletion timescales (cf. labels in Gyr) considering a constant SF efficiency.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=13cm,clip]{13261f4.eps}
\end{figure} Figure 4:

Local KS law involving the extinction-corrected H${\alpha }$ luminosity as SFR tracer and molecular gas surface density at four spatial resolutions. The black points (inbetween the $2\sigma $ boundaries, evidenced by the two dashed lines) are the ones considered for the last iteration of the first fitting method. The solid line is the resulting best fit considering these positions. The grey points were discarded during the various iterations of the fit. The Kennicutt (1998b)'s index is depicted by a dash-dot line (slope of 1.4). The dotted lines represent the times of gas depletion (cf. labels in Gyr in the lower right panel) considering a constant SF efficiency.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=13cm,clip]{13261f5.eps}
\end{figure} Figure 5:

Local KS law involving the extinction corrected FUV luminosity as SFR tracer and total gas surface density at four spatial resolution. The black points (in between the $2\sigma $ boundaries, evidenced by the two dashed lines) are the ones considered for the last iteration of the first fitting method. The solid line is the resulting best fit considering these positions. The grey points were discarded during the various iterations of the fit. The Kennicutt (1998b)'s index is depicted by a dash-dot line (slope of 1.4). The dotted lines represent the times of gas depletion (cf. labels in Gyr in the lower right panel) considering a constant SF efficiency.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8cm,clip]{13261f6.eps}
\end{figure} Figure 6:

Local KS law involving the extinction corrected FUV luminosity as SFR tracer and the molecular gas surface density in  a), the total gas surface density in  b). The spatial resolution is 180 pc. The best fitting linear relation leads to a KS index of 2.22 $\pm $ 0.07 in a) and of 2.59 $\pm $ 0.05 in  b) shown by the heavy lines.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8cm,clip]{13261f7.eps}
\end{figure} Figure 7:

The left-hand panel shows the linear correlation between log  $\Sigma _{\rm SFR}$ and log  $\Sigma _{\rm H_2}/\rho ^{-0.5}_{\rm ISM}$ (in units of $M_\odot $/pc-2 cm-1.5). The Pearson linear correlation coefficient is 0.62, and the slope of the linear fit is 1.16 $\pm $ 0.04. The slope is derived taking into account uncertainties on both axes, but they are not shown in the graph for clarity. On the right-hand panel we show the linear correlation between  $\log \Sigma_{\rm SFR}$ and $\log \rho_{\rm ISM}$ with the Pearson coefficient 0.71.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8.28cm,clip]{13261f8a.eps}\hspace*{-1.3mm}
\includegraphics[width=7cm,clip]{13261f8b.eps}
\end{figure} Figure 8:

$\Sigma _{\rm SFR}$ derived from H${\alpha }$ plotted against $\Sigma _{\rm SFR}$(FUV) ( left panel) and $\Sigma _{\rm SFR}$(bol) ( right panel). Only the 180 pc regions are shown. The dotted line corresponds to equality in the two SFRs.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8cm,clip]{13261f9.eps}
\end{figure} Figure 9:

The theoretical birthline (Paper III) which predicts a decreasing H${\alpha }$ flux for decreasing luminosity of young stellar clusters due to incompleteness of the IMF is shown with filled squares for the randomly sampled cluster model. The filled circles are the model prediction for clusters with a truncated IMF. Cross symbols refer to M 33 data, areas observed at 180 pc resolution.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8cm,clip]{13261f10.eps}
\vspace{4mm}\end{figure} Figure 10:

Correlation between the SFR (FUV) surface density and the dust opacity at 160 $\mu $m. Cross symbols are data for areas in M 33 which have a 180 pc wide diameter; only those points are shown with the CO luminosity above 2$\sigma $ detection threshold.

Open with DEXTER
In the text


Copyright ESO 2010

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.