Brought to you by:

Articles

THE UNIVERSAL INITIAL MASS FUNCTION IN THE EXTENDED ULTRAVIOLET DISK OF M83

, , , , , , , and

Published 2012 March 19 © 2012. The American Astronomical Society. All rights reserved.
, , Citation Jin Koda et al 2012 ApJ 749 20 DOI 10.1088/0004-637X/749/1/20

0004-637X/749/1/20

ABSTRACT

We report deep Subaru Hα observations of the extended ultraviolet (XUV) disk of M83. These new observations enable the first complete census of very young stellar clusters over the entire XUV disk. Combining Subaru and Galaxy Evolution Explorer data with a stellar population synthesis model, we find that (1) the standard, but stochastically sampled, initial mass function (IMF) is preferred over the truncated IMF because there are low-mass stellar clusters (102–3M) that host massive O-type stars; (2) the standard Salpeter IMF and a simple aging effect explain the counts of far-UV (FUV)-bright and Hα-bright clusters with masses >103M; and (3) the Hα-to-FUV flux ratio over the XUV disk supports the standard IMF. To reach conclusion (2), we assumed instantaneous cluster formation and a constant cluster formation rate over the XUV disk. The Subaru Prime Focus Camera covers a large area even outside the XUV disk—far beyond the detection limit of the H i gas. This enables us to statistically separate the stellar clusters in the disk from background contamination. The new data, model, and previous spectroscopic studies provide overall consistent results with respect to the internal dust extinction (AV ∼ 0.1 mag) and low metallicity (∼0.2 Z) using the dust extinction curve of the Small Magellanic Cloud. The minimum cluster mass for avoiding the upper IMF incompleteness due to stochastic sampling and the spectral energy distributions of O, B, and A stars are discussed in the Appendices.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

The initial mass function (IMF) is of central importance in modeling galaxy formation and evolution. Its universality has been supported observationally (e.g., Kroupa 2002; Bastian et al. 2010) and adopted in many applications (e.g., Kennicutt 1998). Observations of resolved stellar populations and integrated colors of most galaxies are consistent with a universal IMF, but a non-universal IMF has been inferred recently, possibly depending on environmental parameters (Bastian et al. 2010). The IMF could be steeper (or truncated), compared to the standard one, in low-density environments, such as dwarfs and low surface brightness (LSB) galaxies, thus producing fewer high-mass stars (Hoversten & Glazebrook 2008; Meurer et al. 2009; Lee et al. 2009). On the other hand, a top-heavy IMF is suggested in dense starburst galaxies (Rieke et al. 1993; Parra et al. 2007).

Pflamm-Altenburg & Kroupa (2008) explained the apparent non-universality as a purely statistical effect; the intrinsic IMF shape is invariant, but the maximum mass of the stars in a stellar cluster is limited by the mass of the parental cluster, the mass of which is in turn controlled by star formation activity. Thereby, the IMF integrated over a galaxy, integrated galactic IMF, may apparently differ in low-density environments. Krumholz & McKee (2008) suggested another model: a low gas density could result in fragmentation into gas clumps too small to allow high-mass star formation. The non-universality of the IMF potentially alters our understanding of star formation and galaxy evolution. Its confirmation is urgent.

The Galaxy Evolution Explorer (GALEX) provides a means to investigate the high-mass end of the IMF. The combination of GALEX and Hα data can separate the populations of O and B stars; the two GALEX bands, far-UV (FUV; 1539 Å) and near-UV (NUV; 2316 Å), are sensitive to both stellar clusters with O stars and with B stars, while Hα data are responsive predominantly to clusters with O stars. Meurer et al. (2009) analyzed the Hα-to-FUV flux ratio over galaxies, found a deficit of ionizing O stars, and concluded that the IMF is steeper (or truncated) in dwarfs and LSBs. This analysis, however, is sensitive to the procedure of dust extinction correction (see alternative analysis; Boselli et al. 2009). Another approach, the comparison of Hα flux with cluster mass, points to a universal IMF (Calzetti et al. 2010).

GALEX finds UV emission far beyond the optical edge of galaxies, i.e., their extended ultraviolet disks (XUV disks; Thilker et al. 2005; Gil de Paz et al. 2005; Thilker et al. 2007). Star formation in such low-density outskirts is a prime target to study the IMF variation. In fact, the radial UV profile extends beyond the Hα truncation radius (Boissier et al. 2007), indicating that O stars are extremely rare while B stars are abundant. The absence may be an indication of a non-universal IMF or could simply be the effect of mixed stellar ages (Thilker et al. 2005)—as O stars are the primary sources of H ii regions and have much shorter lives (∼6 Myr) than B stars (∼100 Myr). Indeed, Goddard et al. (2010) averaged Hα and FUV fluxes over XUV disks and found the ratio consistent with the prediction from the standard IMF, although the sensitivity and area coverage of their Hα data were limited.

We present Hα observations of M83 obtained with the Subaru Prime Focus Camera (Suprime-Cam; Miyazaki et al. 2002) on the Subaru telescope (Iye et al. 2004). With its high sensitivity and wide field of view (34' × 27'), we can cover the entire XUV disk and surrounding area with two pointings (Figure 1) and can detect H ii regions around individual massive stars. M83 is one of the nearest galaxies with a prominent XUV disk (d = 4.5 Mpc or mM = 28.27 mag; Thim et al. 2003) and is best suited for the detection of individual H ii regions. This study provides the first complete census of the UV-bright cluster population with and without H ii regions. We count the number of blue clusters in the UV and bright clusters in Hα and translate the number ratio into the ratio of durations when a cluster is UV-blue and Hα-bright. We compare the ratio with the prediction from a stellar population synthesis model using the standard IMF.

Figure 1.

Figure 1. Subaru coverage and definition of regions in M83. The background is a GALEX FUV image downloaded from Galexview. Solid overlapping rectangles are two pointings with Suprime-Cam (the field of view is 34' × 27' or 44.5 × 35.3 kpc2, in one pointing). The inner (yellow) contour is roughly at the traditional edge of the optical disk. The outer (red) counter is at an H i surface density of 1.5 × 1020 cm−2, using the H i data from Miller et al. (2009). We call FUV-bright objects inside the H i contour IN objects (excluding those in the inner disk—yellow contour) and those outside OUT objects. The Hα-bright clusters (NA659 − RC < −1 mag; Section 4) are also plotted: green for clusters with >103M and cyan for the ones with <103M.

Standard image High-resolution image

2. DATA

2.1. Subaru

Observations were made in the RC and NA659 (Hα) bands (Okamura et al. 2002; Yoshida et al. 2002; Hayashino et al. 2003) on 2010 September 7. Two Suprime-Cam fields cover the entire XUV disk of M83, i.e., the whole area with an H i surface density above 1.5 × 1020 cm−2 (roughly the detection limit in Miller et al. 2009) and a large surrounding area. Figure 1 shows the area coverage and the H i contour (red) on a GALEX FUV image. The two fields are marked with black overlapping solid boxes. The exposure time for each field is 5 × 12 minutes and 5 × 4 minutes in the NA659 and RC bands, respectively. The yellow contour indicates the edge of the optical disk (roughly, the 25 mag arcsec−2 isophote). The typical seeing was ∼1''.

Data were reduced in a standard way using the Suprime-Cam reduction software (Yagi et al. 2002; Ouchi et al. 2004). We performed overscan subtraction, flat fielding, distortion correction, background subtraction, and mosaicking. The flux calibration was performed against standard field stars (Oke 1990; Landolt 1992). RC and NA659 bands are practically at the same wavelength (Hayashino et al. 2003), and we do not expect any spatial variation in the RC/NA659 ratios of typical foreground field stars. However, residual errors remain at large spatial scales in the ratios (a few percent) which were likely introduced in flat fielding. We made a correction by making an NA659/RC ratio map using the average ratio of field stars in each 200'' × 200'' grid and by normalizing the NA659 image. The 1σ photometric limits with the ∼5'' aperture (i.e., GALEX resolution) are 24.90 and 25.56 AB mag in the NA659 (Hα) and RC, respectively.

2.2. GALEX

We use archival GALEX FUV and NUV images (Bigiel et al. 2010). The data are obtained from the GALEX data archive through Galexview in the multimission archive at STScI (MAST). The GALEX observations and data reduction were discussed in Bigiel et al. (2010). The FWHM sizes of the point-spread functions are measured with the IDL software ATV (Barth 2001) and are 4farcs7 for FUV and 5farcs2 for NUV. The 1σ photometric limits with the ∼5'' aperture are 26.34 and 26.20 AB mag in the FUV and NUV bands, respectively.

3. METHODS

3.1. Stellar Population Synthesis Model

Our analysis is guided by the stellar population synthesis model Starburst99 (Leitherer et al. 1999). The Padova stellar evolution tracks with asymptotic giant branch stars are adopted. Most UV-bright objects in the XUV disk of M83 are likely stellar clusters (Thilker et al. 2005; Gil de Paz et al. 2007). Therefore, we adopt an instantaneous burst approximation to model their photometric evolution.

We will show that the standard Salpeter IMF (Salpeter 1955) explains the observations well. Figure 2 demonstrates the differences between the standard and truncated IMFs. The cluster mass is set to 1000 M for the plot, which is roughly the minimum mass of a cluster to fully populate the IMF up to an O star (see Appendix A). Note that a factor of 10 increase in cluster mass shifts all the magnitude curves in Figure 2 (top) systematically upward by 2.5 mag, but it does not change the color curves. For these models, we assume the Salpeter IMF with a lower cutoff mass of ml = 0.1 M. The upper cutoff mass is set to mu = 100, 20, and 3 M, which corresponds to the masses of the most massive O-, B-, and A-type stars, respectively (Cox 2000). In other words, the last two models represent the truncated IMF (or, almost equivalently, the steeper IMF) without O or OB stars. We adopt 20% of the solar metallicity (0.2 Z) from the spectroscopic measurements (Gil de Paz et al. 2007; Bresolin et al. 2009). For reference, the spectra of individual O, B, and A stars are discussed in Appendix B.

Figure 2.

Figure 2. Photometric evolution of stellar clusters from Starburst99. A single-burst model with the Salpeter IMF is adopted. A low metallicity (0.2 Z) is adopted (red) for the analysis, but models of higher metalicities (0.4 Z and 1.0 Z) are also plotted (green and cyan, respectively). The lower mass boundary for the IMF is 0.1 M and the upper mass boundary is 100 M (including all types of stars), 20 M (no O stars), and 3 M (no OB stars). The adopted distance of M83 is d = 4.5 Mpc. Top: FUV magnitude for a cluster mass of 1000 M. Middle: FUV–NUV color. Bottom: NA659–RC color.

Standard image High-resolution image

There are quantitative differences between clusters with and without OB stars, i.e., between the standard and truncated (steeper) IMFs. Clusters without OB stars are fainter by 6 mag in FUV for the first 10 Myr and never become as blue as FUV–NUV < 0.2–0.3 mag. Therefore, blue clusters (<0.2–0.3 mag) must have O and/or B stars. The FUV–NUV color is insensitive to the difference between clusters with only B stars (up to mu = 20 M) and with both O and B stars (100 M). This is primarily because the GALEX bands trace only the Rayleigh–Jeans side of stellar thermal radiation for both O and early-type B stars (Appendix B). The FUV–NUV color is therefore a useful indicator of the presence of O and/or B stars in a cluster.

The NA659 − RC color differentiates clusters with O stars from those without them, since O stars emit a significantly greater number of ionizing photons than B stars (Sternberg et al. 2003). This is demonstrated in Figure 2 (bottom). Clusters with NA659 − RC < −1 mag must have O stars, and those with NA659 − RC < −0.3 mag should have O and/or B stars. For this plot, we calculated the Hα line luminosity L using the Lyman photon flux QLyc from Starburst99 and the relation L[erg s−1] = 1.37 × 10−12QLyc[s−1] from case B recombination with an electron temperature of 104 K (Gavazzi et al. 2002). Case B recombination is perhaps a reasonable assumption even in the low-density environment, since the optical depths in the Lyman resonance lines are very likely quite large (Osterbrock & Ferland 2006) in any star-forming conditions. The NA659 magnitude is based on the stellar continuum luminosity plus L. We ignore the effect of [N ii] emission since the effect is very small (Gil de Paz et al. 2007; Kennicutt et al. 2008; Goddard et al. 2010; see Section 4.1).

3.2. Cluster Counts

Counting blue clusters in FUV–NUV and NA659 − RC constrains the IMF at its high-mass end. Clusters become bright in Hα ("blue" in the color NA659 − RC) predominantly due to O stars, while FUV–NUV is due to both O and B stars. Since O stars have much shorter lives than B stars, NA659 − RC changes faster than FUV–NUV (Figure 2). Under the assumptions of instantaneous cluster formation and a constant cluster formation rate, the numbers of blue clusters in FUV–NUV and NA659 − RC (NUV and N) should be proportional to the durations that clusters are blue (tUV and t, respectively). Thus, we expect N/NUV = t/tUV if the difference in the populations is due to a simple aging effect starting from the initial populations of the standard IMF. The assumption of a constant cluster formation rate should be approximately accurate when integrated over the large XUV disk and over a short timescale, i.e., the young ages of clusters which can be detected with our color selection criteria (below). In fact, the clusters that we count in Section 5.3 are distributed over the UV disk (Figure 1).

For the standard IMF (mu = 100 M and 0.2 Z), a cluster is blue (e.g., FUV–NUV < 0.0 mag) for tUV = 70.8 Myr and Hα-bright (e.g., NA659 − RC < −1.0 mag) for 5.8 Myr. Therefore, the number ratio should be N/NUV = 0.08 if each population is the consequence of the cluster aging effect. We will compare the observations with this theoretical prediction.

We identify and count stellar clusters at the GALEX resolution (∼5''; ∼ 109 pc at the distance of M83). It is possible that clusters may be blended at this resolution; however, even in that case, our analysis is valid if clusters formed in such proximity are typically coeval.

4. SAMPLE SELECTION

4.1. FUV-selected Objects

UV-bright objects are identified based on the FUV image using the SExtractor package (Bertin & Arnouts 1996). Saturated stars and image edges are masked before the identification. Objects within the traditional optical disk (yellow contour in Figure 1) are excluded. We confirmed the very low probability of false detection using the negative image technique. SExtractor is run against the sign-flipped FUV image and detected only one false target at 25.41 mag in FUV. Thus, false detections have a negligible impact on our number counts.

The aperture photometry is carried out in the four bands with the dual image mode of SExtractor using the FUV image as the reference. The 1σ photometric limits with the aperture of the GALEX resolution (∼5'') are 24.90, 25.56, 26.34, and 26.20 mag in the NA659 (Hα), RC, FUV, and NUV bands, respectively. The sensitivity in NA659 (FWHM = 120 Å) corresponds to a flux of 7.98 × 1034 erg s−1 at a distance d = 4.5 Mpc, which is sufficient to detect an H ii region around a single B star (e.g., 3.44 × 1035 erg s−1 in the case of B0; assuming case B recombination; see Appendix B). We note that the NA659 band includes stellar continuum emission as well as Hα emission.

The contamination from [N ii] λ6548, 6583 in the NA659 photometry is neglected. The line ratio [N ii]/Hα is small (≲ 0.1) in the outskirts of galaxies (Gil de Paz et al. 2007; Kennicutt et al. 2008; Goddard et al. 2010). Even in the extreme case that the stellar continuum has zero contribution to the NA659 flux, the error is only 0.1 mag (and likely less in more realistic conditions).

Figure 1 shows that the Subaru data cover both the inside and the outside of the H i contour (roughly the H i detection limit). Recently formed FUV-bright objects are more likely to reside inside the H i contour as star formation occurs in gas. For the background subtraction in Section 4.4, the objects inside the H i contour (hereafter, IN objects) and outside (OUT objects) are separately cataloged. Note that we intensively used the IMCAT software package (N. Kaiser 2004, private communication)8 for the following data analysis.

In Section 5.3, we will discuss a selection of clusters based on their NA659 − RC color; but in advance, Figure 1 shows the locations of clusters with NA659 − RC < −1 mag (i.e., bright in Hα). They are distributed over the entire XUV disk, and all are within the H i contour (see Section 4.4).

4.2. Extinction Correction

The extinction curves of Pei (1992) are adopted to correct for Galactic and internal extinctions. The extinctions in FUV, NUV, RC, and NA659 are (AFUV, ANUV, $A_{R_{\rm C}}$, ANA659) = (2.58, 2.83, 0.81, 0.81)AV for the Galaxy, (3.37, 2.61, 0.80, 0.80)AV for the Large Magellanic Cloud (LMC), and (4.18, 2.59, 0.79, 0.79)AV for the Small Magellanic Cloud (SMC), where AV is the extinction in V-magnitude. We note that the FUV–NUV color becomes bluer under the Galactic extinction, but redder under the LMC and SMC extinction.

The Galactic extinction correction is made for all objects using AV = 0.218 mag (Schlegel et al. 1998). We do not apply an internal extinction correction until Section 5, but it will be corrected using AV = 0.1 mag and the SMC extinction curve for the following reason. The measured metallicity of most outer objects is low (∼0.2 Z; Gil de Paz et al. 2007), and therefore, the extinction curve of the SMC is perhaps the most appropriate. Using E(BV) values from Gil de Paz et al. (2007), we estimate the amount of internal extinction to be around AV ∼ 0.1 mag with outliers (e.g., 0.0–1.0 mag). We therefore adopt AV ∼ 0.1 mag for the statistical correction, and this degree of extinction will be supported by a comparison with a model color in Section 5.2.

4.3. Completeness

The mass detection limit is determined in color–magnitude space with the help of models. Figure 3 shows color–magnitude diagrams for IN objects. The overlaid solid lines are the 3σ and 5σ mag detection limits and model evolutionary sequences for clusters with masses of 102, 103, and 104M. The models are for clusters with the standard Salpeter IMF with (ml, mu) = (0.1, 100 M) and metallicity 0.2 Z. Models with lower mu do not reproduce clusters as blue as observed (Figure 2 and see discussion below). Clearly, the mass detection limit changes as a function of color in Figure 3; it decreases as clusters become redder. We note that the RC and NA659 bands are practically at the same wavelength, and therefore, NA659 − RC is zero without Hα emission. A negative NA659 − RC color indicates the presence of associated Hα emission.

Figure 3.

Figure 3. Color–magnitude diagrams of all the objects inside the H i contour (IN objects defined in Section 4.1). Top: NA659 − RC vs. NA659. 5σ and 3σ detection limits are illustrated with solid straight lines. The model evolutionary sequences with the standard Salpeter IMF are displayed with solid curves for clusters with 102, 103, 104M. Bottom: the same as the top panel, but for FUV–NUV vs. FUV. Single O3 and B0 star models are also plotted in both panels.

Standard image High-resolution image

The top panel shows that the Subaru observations detect clusters down to ∼102M. Clusters become redder and fainter as they age, and those around the detection boundary (102M) become too faint when they become redder than NA659 − RC ∼ −0.6 mag. In other words, the sample is complete down to 102M in the range NA659 − RC < −0.6 mag. It is complete down to 103M in NA659 − RC < 0.0 mag. The extinction is negligible, as shown by the extinction vector (for AV = 0.1 mag) in the plot.

The bottom panel shows the corresponding color–magnitude diagram in the UV. Again, the mass detection limit is a function of the color range (and thus, the cluster age). The detection limit is ∼103M for the color FUV–NUV < 0.0 mag and is even deeper (∼102M) for FUV–NUV < −0.2 mag. The extinction is not negligible in FUV and NUV and will be discussed later.

Note that the very blue objects in the bottom panel (FUV–NUV ≲ −0.3 mag and FUV ∼ 22–24 mag) are bluer than any model cluster with the standard IMF can reproduce and are located approximately in the region where single O stars appear (the models of O3 and B0 stars are also plotted). These objects are not included in our number count, since they appear below the mass selection criterion in Section 5.3. We will discuss these objects in Section 5.1.

4.4. Background Subtraction

The contamination of background (and foreground) objects has been an obstacle in identifying the objects within nearby galaxies (Werk et al. 2010; Goddard et al. 2010). The large field of view of Suprime-Cam is a clear advantage as it enables us to remove such contamination statistically. The objects outside the H i gas contour (OUT objects; see Section 4.1) are most likely in the background, whereas IN objects include stellar clusters within M83 as well as background objects.

The IN and OUT objects show different colors statistically. Figure 4 shows the plots of NA659–RC versus FUV–NUV for the IN and OUT objects above the 5σ mag detection limit and above the 103M threshold line in Figure 3 (bottom). Figure 5 is a histogram of FUV–NUV color for the IN (green) and OUT objects (blue). Clearly, the IN objects are bluer in FUV–NUV than the OUT objects and preferentially have FUV–NUV < 0.2–0.3 mag, while the OUT objects are more likely to have FUV–NUV > 0.2–0.3 mag. Figure 4 shows that virtually all objects with NA659 − RC < −0.3 mag (OB stars) are IN objects; almost no young stellar clusters exist outside of the H i gas. This clear distinction ensures that the objects outside of the H i gas are not the young clusters in M83 and represent the background population. Note that IN objects are young stellar clusters, not very old clusters with many evolved FUV-bright stars, such as extreme horizontal branch stars, blue stragglers, or planetary nebulae. The old clusters do not become as blue in FUV–NUV as our IN objects (Donovan et al. 2009).

Figure 4.

Figure 4. NA659 – RC vs. FUV–NUV color–color plot of the objects above 5σ detection and above 103M. Lines are model predictions with the Salpeter IMF with mu = 100 M for 0.2 Z. A reddening vector of AV = 0.1 mag from the SMC extinction curve is also plotted. We assume that stellar clusters (traced in RC) and the surrounding ionized gas (in NA659) suffer from extinction by the same amount. Top: FUV-bright objects within the H i contour (IN objects). Bottom: the FUV-bright objects outside (OUT objects).

Standard image High-resolution image
Figure 5.

Figure 5. Histograms of FUV–NUV color of the objects above 5σ detection and above 103M. Colors indicate the objects within the H i gas (green), the ones outside of the H i gas (background objects; blue), and the ones associated with M83 (red). The associated ones are calculated by subtracting OFF from ON with scaling for the area coverage difference.

Standard image High-resolution image

The background contamination in the IN objects can be removed statistically; we subtract the OUT objects from the IN. In Figure 5, the number of the OUT objects are scaled by the IN/OUT area ratio (55%) and subtracted from that of the IN objects. The derived histogram (red) represents the population of young clusters with mass >103M in the XUV disk.

Virtually all clusters after the subtraction show blue color in FUV–NUV (<0.2–0.3 mag; Figure 5), suggesting that all of these young clusters have O or, at least, B stars (see Figure 2). Thus, massive stars are populated even in the diffuse gas environment. It is also remarkable that almost no red objects (FUV–NUV > 0.2–0.3 mag) exist. This confirms that OUT objects are distributed everywhere across our field of view and represent the background objects; thus, our background subtraction is fairly successful.

There is a possibility that old stellar clusters (i.e., FUV–NUV > 0.2–0.3 mag before extinction correction) exist both in the IN and OUT regions and are removed from the IN population. However, it should not affect our cluster count analysis as we discuss only the young clusters. Such a population, if it exists, must be fairly uniformly distributed both in the IN and OUT regions, since the count after the subtraction is nearly zero in FUV–NUV > 0.2–0.3 mag. We also note that the apparent magnitudes of low-mass clusters (≲ 103M) fall below our detection limit after ∼100 Myr (Figure 2, top). Again, this is not a problem with the selection criteria that are adopted in our analysis (Section 5.3).

5. DISCUSSION

5.1. Truncated versus Stochastic IMFs

The presence of ionizing stars (O stars) is evident even in small clusters (<103M) in Figure 3, suggesting that the IMF is not truncated. Massive stars do not form under the truncated IMF at all. On the other hand, they should appear occasionally even in small clusters under the stochastic IMF, though they tend not to appear in a small sample due to their low probability of formation in small clusters. Such a stochastic effect becomes particularly apparent when the cluster mass is small (<1000 M). Our Subaru observations detect small clusters (down to ∼102M) and confirm that even those small clusters show Hα excess (i.e., NA659 − RC < 0.0 mag); and some even have NA659 − RC < −1.0 mag, a color which only O stars can produce. Thus, O stars are present in the 102–3M clusters. The IMF is not truncated even in the low-density and low-metallicity environment.

More precisely, it is necessary to account for the stochastic effects in discussions of very low mass clusters (<103M)—and even then, our result still supports the stochastic IMF over the truncated IMF, as we discuss below. We first need to consider the likelihood that a cluster with a given small mass has a massive star. In Appendix A, we estimate that a 103M cluster is likely to have as massive a star as 40 M; the cluster mass must be greater than ∼400 M to have a high probability of hosting at least one O star (>20 M); a 102M cluster likely hosts up to a ∼10 M star.

The NA659 − RC color is sensitive to the mass of the most massive star in a cluster when the cluster is young, whereas FUV–NUV does not distinguish O stars from B stars well (Section 3.1; Figure 2; see also Appendix B). If a cluster has no massive star, it should appear redder (NA659 − RC > −1 mag; Figure 2) than the observed clusters in the range 102–3M; yet some clusters appear in the range NA659 − RC = −1 ∼ −2 mag. Massive O stars in small clusters should be populated stochastically. In addition, a stochastically populated single O star can outshine a cluster with a small mass (102–3M; Figure 2). This may result in an overestimation of the cluster mass when we apply the standard IMF with mu = 100 M to the cluster (as done in Figure 3). Therefore, the masses of small clusters that appear in the range 102–3M in Figure 3 could be actually even less. This prefers the stochastic IMF, as even the less massive clusters have O stars.

The discussion so far is about clusters of stars. Another possible evidence for the presence of O stars is the population of very blue objects (FUV–NUV ≲ −0.3 mag) in Figure 3 (bottom). They have the apparent magnitude and color similar to those of a single O star. They are possibly isolated, single O stars or small clusters with an O star(s) but without many low-mass stars. Higher resolution images are necessary to confirm them; their presence, if confirmed, may be stronger evidence for O stars that are formed in the low-density environments.

5.2. Internal Extinction Correction

The internal extinction is assumed to be AV = 0.1 mag from spectroscopic measurements (Section 4.2; Gil de Paz et al. 2007). This degree of extinction on average is also supported by a comparison with the model curve in Figure 4 (top). The extinction vector, corresponding to AV = 0.1, is also plotted. An extinction correction by this amount creates agreement between the UV data and the model curve, providing additional support for the adopted AV.

We assume that a stellar cluster (traced in RC) and the surrounding ionized gas (in NA659) suffer from extinction by the same amount. The two components may have different spatial distributions and a slightly different attenuation (e.g., E(BV)continuum ∼ 0.44E(BV)gas in case of more active star-forming environments; Calzetti 2001). However, the difference is negligible for our analysis since the extinction is small (AV = 0.1).

5.3. Universal IMF

Detected clusters with a color of NA659 − RC < −1 mag are evidence that some clusters host O stars. Here we test whether their abundance, relative to the ones without O stars, can be explained by a simple aging effect with the standard IMF. Under the assumptions of instantaneous formation of star clusters and a constant cluster formation rate over the XUV disk, the numbers and durations of blue clusters should show the relation N/NUV = t/tUV, if the relative cluster populations are due to the aging effect.

Consistent sample selection between GALEX and Subaru data is crucial in this analysis. Each of the four photometric bands has its own detection limit (Figure 3), and a bias could be introduced if only magnitude detection limits are taken into account. We therefore use the cluster mass and color as thresholds for cluster counts. As discussed in Section 4.3, both GALEX and Subaru data are complete at ≳ 5σ, in colors NA659 − RC < −0.3 mag and FUV–NUV< 0.0 mag down to a cluster mass of 103M. This is true even after the internal extinction correction (Section 5.2) if a slightly lower detection limit (3.4σ) is adopted. The extinction correction is made for all four bands, which is equivalent to shifting the data points along the extinction vectors in Figure 3. We select the mass-limited sample after the extinction correction using the technique discussed in Section 4.

Almost all of the clusters are in the range FUV–NUV < 0 mag when the extinction is taken into account (Figure 5); thus, we adopt this color range for our analysis. This is the range where B stars definitely, and O stars possibly, exist in the clusters; we count how many of them have O stars below. Clusters are blue, FUV–NUV < 0 mag, for only a short duration. With the standard IMF of (ml, mu) = (0.1, 100 M) and 0.2 Z, this is tUV = 71 Myr (Figure 2). The short tUV is advantageous, since the assumption of a constant cluster formation rate is likely valid for a short period over the large area (i.e., the entire XUV disk)—no local events likely bias the statistics. On the other hand, clusters with O stars show NA659 − RC < −1 mag for only 5.8 Myr. The model therefore predicts $t_{\rm {\rm H}{\alpha }}/t_{\rm UV} \sim 0.08$.

The number of blue clusters with FUV–NUV < 0 mag is counted using Figure 5 (red)—after the internal extinction correction, it is NUV = 88 ± 9. From Figure 3, N = 9 ± 3 for NA659 − RC < −1 mag (see Figure 6). The errors are based on the Poisson distribution, and the counts are made after the internal extinction correction. Therefore, the observed ratio is N/NUV ∼ 0.10 ± 0.03. This is consistent with the prediction of the standard IMF (within a 1σ deviation). The counts of stellar clusters support the standard IMF with a simple aging effect.

Figure 6.

Figure 6. NA659 images of the Hα-bright clusters with NA659 − RC < −1 mag and >103M. O-type stars must be present to produce this cluster color. The size of each image is 24farcs2 × 24farcs2.

Standard image High-resolution image

It is worth noting that the slight deviation does not likely have a statistical significance. However, if anything, it is inclined toward the top-heavy IMF side, as opposed to the previous suggestion of a truncated IMF.

The combination of the parameters we adopt (i.e., low metallicity (0.2 Z), SMC extinction, and AV = 0.1 mag) is the best to consistently explain the new data, model, and previous spectroscopic studies. The result for the standard IMF does not change even if a higher metallicity (∼0.4 Z), and hence the LMC extinction curve, is adopted. If we ignore the spectroscopic studies, the comparison of the data and model (e.g., Figure 4) provides AV ∼ 0.2 mag. With this extinction, we count NUV = 86 ± 9 and N = 11 ± 3 and obtain N/NUV ∼ 0.13 ± 0.04. The model ratio is 0.10 with the metallicity 0.4 Z, where $t_{\rm {\rm H}{\alpha }} = 5.2$ Myr and tUV = 51.3 Myr. Therefore, the observations and model are consistent even if a different metallicity and extinction are assumed.

5.4. Hα-to-FUV Flux Ratio

The Hα-to-FUV flux ratio (F/fν FUV) over the entire XUV disk is another measure of the high-mass end of the IMF. We find that the observed ratio is consistent with the one predicted from the standard IMF. A lower ratio may suggest a truncated IMF (Meurer et al. 2009; Lee et al. 2009), though an extinction correction introduces an ambiguity (Boselli et al. 2009). We can estimate the ratio more accurately since we derived the extinction of AV = 0.1 mag statistically (Section 5.2).

In a Starburst99 simulation with a constant star formation rate, a metallicity of 0.2 Z, and the Salpeter IMF with (ml, mu) = (0.1, 100) M, the ratio approaches an approximate constant and is, e.g., log (F/fν FUV) ∼ 13.10 in unit of log (Hz) after 1 Gyr assuming a star formation rate of 1 M yr−1 (though the constant does not depend on these values). Fumagalli et al. (2011) calculate the ratio with their own model using the Kroupa IMF and derive ∼13.22 (from their Figure 1). The scatter between the models is about 30%. Compared to these values, the observed ratios previously reported in LSB galaxies by Meurer et al. (2009)—as low as one-third of the values calculated in these models—deviate at a significant level with respect to the model uncertainties.

We calculate the Hα-to-FUV flux ratio over the XUV disk of M83 by summing up the fluxes of all >3σ detections in the FUV image (Section 4.1). Assuming that OUT objects represent the background population (Section 4.4) and subtracting their total fluxes form those of IN objects (accounting for the IN/OUT area ratio), we obtain F = 8.36 × 1028 erg s−1 and $f_{\nu \,\rm FUV}=7.23\times 10^{25} \,\rm erg\, s^{-1} \, Hz^{-1}$. Adopting AV = 0.1 mag and the SMC extinction curve (Section 4.2), the intrinsic ratio after the extinction correction is log (F/fν FUV)int = 12.93. This value is only ∼30% lower than our model prediction, which is not as much as what the previous study found (Meurer et al. 2009). As the difference between this value and our model calculation is similar to the scatter between the models, we consider a truncated IMF to be unwarranted in the case of M83.

6. SUMMARY

Using new deep Subaru Hα data and archival GALEX data of the XUV disk of M83, we show that the standard Salpeter IMF explains the number counts of stellar clusters in the low-density, low-metallicity environment. Even in clusters with low-mass (102–3M) O stars are populated, at least occasionally, suggesting that the stochastic IMF is preferable to the truncated IMF. The Hα-to-FUV flux ratio over the XUV disk also supports the standard IMF.

The analyses were guided by the stellar population synthesis model Starburst99. The new data, together with the population synthesis model and previous spectroscopic studies, provide overall very consistent results at the metallicity (∼0.2 Z), internal extinction (AV ∼ 0.1 mag), and colors of the stellar clusters in the XUV disk. We counted the numbers of FUV-bright clusters (which should have O and/or B stars) and Hα-bright clusters (having O stars), whose ratio constrains the high-mass end of the IMF. The ratio of number counts can be translated to the ratio of their lifetimes with Starburst99, under the two assumptions of instantaneous cluster formation and a constant cluster formation rate in the XUV disk. As a result, the number counts are consistent with the Standard IMF and a simple aging effect.

This study benefits from the combination of the unparalleled wide field of view and sensitivity of Suprime-Cam. The new image covers an area far beyond the detection limit of the H i gas (thus, beyond the boundary of the XUV disk), enabling us to remove the background contaminations statistically, which had been a difficulty in previous studies. We found that virtually all blue clusters (FUV–NUV ≲ 0.0 mag after the extinction correction or ≲0.16 mag before the correction) are within the H i gas.

We thank Sadanori Okamura and Tomoki Hayashino for providing the redshifted Hα filter for Suprime-Cam. J.K. thanks Peter Capak for his help with the IMCAT software and Deane Peterson and Fred Walter for discussion on stellar spectra. A.G.d.P is funded by the "Ramoń y Cajal" program of the Spanish MICINN and partly supported by the AYA2009-10368 and Consolider-Ingenio 2010 CSD2006-70 projects. M.I. is supported by Grants-in-Aid for Scientific Research (No. 22012006, 23540273). Data analyses were in part carried out on the common use data analysis computer system at the Astronomical Data Center (ADC) of the National Astronomical Observatory of Japan. This work is based in part on data collected at Subaru Telescope, which is operated by the National Astronomical Observatory of Japan. Some of the data presented in this paper were obtained from the Multimission Archive at the Space Telescope Science Institute (MAST). STScI is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. Support for MAST for non-HST data is provided by the NASA Office of Space Science via grant NNX09AF08G and by other grants and contracts.

Facilities: Subaru - Subaru Telescope, GALEX - Galaxy Evolution Explorer satellite

APPENDIX A: MINIMUM STELLAR CLUSTER MASS TO AVOID STOCHASTICITY

We consider the threshold mass of a stellar cluster where the incomplete sampling (filling) of the IMF (i.e., stochastic effect) becomes important. It depends on the definition of "stochastic sampling," because in reality an O star with the highest mass ∼100 M is rarely found (at least in the Milky Way), and even a massive cluster (e.g., 106M) may not host any star with 100 M in a statistical sense. On the contrary, we could define the fully sampled IMF as the one in which at least one O star of any subclass (>20 M) is populated. Of course, a cluster is more likely to host a 20 M star than a 100 M star, and the threshold mass should be lower in the former case. One could also take a different approach; for example, Calzetti et al. (2010) found ∼103M as the threshold mass where the Hα luminosity over cluster mass starts to deviate from a constant.

Here, we take a simpler approach to illustrate a range of threshold cluster mass that may provide an important guideline. The threshold mass necessary to populate the IMF to some upper mass cutoff depends on that maximum stellar mass

The stellar IMF (Ψ(m)) with the upper and lower stellar mass cuts (mu and ml, respectively) is normalized as

Equation (A1)

where Mcl is the total mass of a stellar cluster. Assuming a power-law IMF (Ψ(m) = Am−α) the normalization coefficient (A) is

Equation (A2)

If we sum up the masses of all stars with masses greater than a threshold stellar mass mt, the sum is

Equation (A3)

In order to populate at least one star as massive as mt, this sum should be greater than mt [i.e., M(> mt) > mt]; and therefore,

Equation (A4)

in the case of a power-law IMF. This condition is necessary to unstochastically populate the IMF up to mt in the cluster, providing a threshold cluster mass to have a star up to mass mt.

A similar calculation also provides a guideline of the maximum stellar mass that a cluster with Mcl is likely to host—the IMF would be filled up to mt if it is populated smoothly from the lowest stellar mass ml toward the highest until the total mass reaches the cluster mass Mcl. Thus, the maximum stellar mass in the cluster is around mt.

Figure 7 shows the threshold cluster mass Mcl as a function of the maximum stellar mass mt that should be populated. We assumed the Salpeter IMF (α = 2.35), ml = 0.1 M, and mu = 100 M. In order to have a 40 M star the cluster mass should be at minimum 103M, while to have 80 M it needs to be 104M. All stars with >20 M are classified as O-type stars; and to have at least one O star the cluster mass should be greater than ∼400 M. To fill the IMF literally up to ∼100 M, the cluster mass should be on the order of 105M. If we adopt 50 M (arbitrarily) as a typical mass of O stars, the cluster mass should be >2000 M to fully populate the IMF up to this typical O star mass and avoid the stochastic effect.

Figure 7.

Figure 7. Minimum cluster mass to avoid the stochastic sampling effect of the IMF, as a function of the maximum stellar mass that one wants to populate unstochastically. The cluster mass needs to be 103, 104, and a few 105M to populate an O star up to 40, 80, and 100 M, respectively.

Standard image High-resolution image

These threshold cluster masses are meaningful only in a statistical sense. A stochastic effect could appear even above this threshold mass if a Monte Carlo simulation is performed. Even if a cluster mass is above the threshold and the cluster can host a star with mt, this mass could be split and allocated to multiple stars with smaller masses (whose total is mt). On the other hand, a cluster with mass below the threshold could occasionally have a star with mass >mt at the expense of many less massive stars. The threshold cluster mass should be taken as a minimum mass to avoid the stochastic effect, only in a statistical sense. Despite this limitation, our estimate provides a crude guideline for discussion of the stochastic sampling of the IMF.

APPENDIX B: UV SPECTRA OF O, B, AND A STARS

The spectra of O3-, B0-, and A0-type stars from the Kurucz's ATLAS9 grid (Kurucz 1992) are plotted in Figure 8. The masses of these stars are approximately 87.6, 17.5, and 2.9 M, respectively (Sternberg et al. 2003; Cox 2000). The effective temperature (Teff) and gravity (log g) are required to find a model from the grid, and the stellar radius R is necessary to calculate the total flux. From the ATLAS9 grid, we choose the models whose parameter sets are the closest to the stellar parameters presented in Sternberg et al. (2003) or Cox (2000): (Teff, log g, R) = (50,000 K, 5.0, 13.2 R); (30,000, 4.0, 7.4); and (10,000, 4.0, 2.4) for O3, B0, and A0, respectively. We note that the spectra of high-mass stars are not calibrated very well as only few of those exist in the solar neighborhood for calibration. In addition, an adjustment of the parameters is necessary to fit observed spectra of individual stars (Sternberg et al. 2003; Leitherer et al. 2010). The spectra that we present here are only a reference to understand our observations.

Figure 8.

Figure 8. Top: UV spectra of O3 (87.6 M), B0 (17.5 M), and A0 (2.9 M) stars. The y-axis is an apparent magnitude at the distance of M83 (d = 4.5 Mpc; mM = 28.27 mag). Bottom: transmission curves of GALEX FUV and NUV bands.

Standard image High-resolution image

The continuum slopes are more or less the same between the O and B stars in the wavelength range of GALEX bands. On the other hand, the A star spectrum is rising from FUV to NUV. The FUV absolute magnitude (AB) and FUV–NUV color of the individual stars for three metallicity [Fe/H] = +0.0 (1.0 Z), −0.5 (0.3 Z), −1.0 (0.1 Z) are in Table 1. We also list the Lyman continuum photon emission rate (QLyc), which can be used to calculate an Hα luminosity (L[erg s−1] = 1.37 × 10−12QLyc[s−1] in case B recombination). The QLyc is calculated by integrating the model spectra over the frequency above the ionizing frequency (i.e., the wavelength of 912 Å). These QLyc values are four to five times smaller than those in Sternberg et al. (2003), reflecting the uncertainty in model spectra of high-mass stars.

Table 1. UV Color of High-mass Star Models

    [Fe/H] = +0.0   −0.5   −1.0
Type log QLyc FUV FUV–NUV   FUV FUV–NUV   FUV FUV–NUV
  (s−1) (mag) (mag)   (mag) (mag)   (mag) (mag)
O3 49.16 −6.47 −0.54   −6.46 −0.57   −6.44 −0.58
B0 47.37 −4.44 −0.31   −4.47 −0.37   −4.48 −0.41
A0 37.24 4.30 0.58   4.16 0.45   4.10 0.40

Note. FUV magnitudes are the absolute magnitude.

Download table as:  ASCIITypeset image

Footnotes

Please wait… references are loading.
10.1088/0004-637X/749/1/20