<<

arXiv:astro-ph/0611586 v1 19 Nov 2006 ag rcino tr omi lses sdeter- as clusters, in form of fraction large loigo,I 47405 IN Bloomington, 3 E. 727 University Indiana Astronomy of Department S102 7 .Wlo Ave. 91125 Wilson CA Pasadena, S. 770 100-22 MS Technology of Center Institute Analysis California and Processing Infrared lcrncades [email protected]fl.edu address: Electronic [email protected] address: Electronic [email protected] address: Electronic 1 2 trfraini ihycutrdpoes A process. clustered highly a is formation rpittpstuigL using 2006 typeset 21, Preprint November version Draft H NTA LSE ASFNTO FSPRSA LSESI I IN CLUSTERS STAR SUPER OF FUNCTION MASS CLUSTER INITIAL THE rsn Address: Present rsn Address: Present rd g n n infiatrdcini einetnto after extinction median in reduction significant find and age asv trcutr odses rmv wyfo hi nat their from away identify move We or headings: disperse Subject . to irregular clusters and star spiral massive both in clusters ealct n aatcserrt.W xmn h evoluti the form examine We cluster star rate. of shear theories galactic for and result our of implications eeotie o h CFo h prlglx ape with sample, spiral the of ICMF the for obtained were rmtefit aarlaeo h la iia k uvy(SD Survey Sky Digital Sloan the of release data i gal ( nearby ICMFs fifth thirteen the on the based compare from is and galaxies measure irregular of we sample universality, its gauge To otn.Teetntos gs n asswr determined were masses and ages, extinctions, The routine. h aewyt eieterIM ae nasmlrnme fy of number similar and a irregular of on ICMFs based the ICMF that spir find their nearby We three derive ICMF. analyzed galaxy to We irregular way small. same be the to found inserte and analyzi cor by were galaxies assessed Completeness clusters was contamination artificial galactic models. which background in starburst simulations from Carlo generated those to r ttsial nitnusal.Frcutr oema more clusters For indistinguishable. statistically are aaisi esnbywl tb oe law power a by fit well reasonably is galaxies . Street h nta lse asfnto IM)i udmna pr fundamental a is (ICMF) function mass cluster initial The 0Mr,msie( massive Myr), 20 1. A INTRODUCTION T E tl mltajv 10/09/06 v. emulateapj style X aais reua,glxe:sia,glxe:sa clu star galaxies: spiral, galaxies: irregular, galaxies: & 3 × 10 4 M eateto hsc Astronomy & Physics of Department nt fAtooy et fPhysics of Dept. Astronomy, of Inst. ⊙ rf eso oebr2,2006 21, November version Draft lseswr eetduiga uoae oreextracti source automated an using selected were clusters ) 00E nvriyAe and Ave. University E. 1000 rn .Buckalew A. Brent nvriyo lrd and Florida of University 03Z¨urich, Switzerland 8093 n et fAstronomy of Dept. and ac .Dowell D. Jayce ansil,F 32611 FL Gainesville, nvriyo Wyoming of University oahnC Tan C. Jonathan aai,W 82071 WY Laramie, oso,T 77005 TX Houston, 30S anS and St Main S. 6320 ABSTRACT GALAXIES ieUniversity Rice dN T Z¨urichETH dM ( and M ) ∝ oie a on pe Tn20) h tla densi- stellar the The of excess 2005). in to (Tan are due speed that sound pressures speeds gas high their escape ionized very that and have gravity mean must self sizes clouds small new gas relatively present natal and their properties example, masses For large physical theories. formation extreme star for challenges Their al. sons. et Mengel 2000). 2005; Ho al. & et Beck, (Clark all Turner, in galaxy 2002; found of are types examples these and of (SSCs), clusters star super 03,i neatn aaislk h Antennae the 2001). clusters, al. star- star et (Tremonti dwarf like massive 5253 most in NGC The as and galaxies such 2005), galaxies burst interacting Whitmore & Lada Chandra, & in (Fall, (Lada neighborhood Galactic 2003), local our in mined M h td fSC sipratfranme frea- of number a for important is SSCs of study The sv hn10 than ssive − noec aay oerudselrand stellar Foreground galaxy. each into d xe eetdfo ouelmtdsample volume-limited a from selected axies to,wihapast eidpnetof independent be to appears which ation, gSS mgso ed rudtesample the around fields of images SDSS ng α 1 hstmsaea htrqie o young, for required that as timescale this M no iulextinction, visual of on 2 lgas. al ycmaigtheir comparing by with etoswr efre sn Monte using performed were rections lglxe ihSS aai exactly in data SDSS with galaxies al ∼ pryo trfraini galaxies. in formation star of operty reua n prlglxe.Our galaxies. spiral and irregular n α prlglxe o masses for galaxies spiral 5 S,fo hc bu 8 young 580 about which from SS), M α − 5 M 1 = M 0Mrb bu . a for mag 0.5 about by Myr 10 sters ug asv lsesa the as clusters massive oung, 1 = ⊙ . h CFo h irregular the of ICMF the , 74 . 62 ± RGLRADSPIRAL AND RREGULAR ± 0 . 0 7 edsusthe discuss We 07. u . 8 iia results Similar 08. ′ g A ′ & i V ′ z ihcluster with , 10 ′ magnitudes 5 > M ⊙ 10 r nw as known are , 5 M on ⊙ 2 ties of the final clusters are not likely to be that much SSCs (Mengel et al. 2002; McCrady, Gilbert, & Graham greater than typical rich young Galactic star clusters, like 2003). Most previous studies of cluster populations have the Orion Nebula Cluster, but if the initial mass func- concentrated on the luminosity function (LF), which also tion is similar, there are more massive OB stars present, often resembles a power law, potentially thousands within just a few of each dN(L) other. ∝ L−αL . (2) Young SSCs are likely to be similar to the progeni- dL tors of globular clusters (Gonz´ales-Delgado et al. 1997; From a sample of spiral galaxies, Larsen (2002) found Zhang & Fall 1999) although this requires an evolution 2.0 . αL . 2.4. In the , of the initial cluster mass function (ICMF) from a power- Whitmore et al. (1999) found αL = 2.6 ± 0.2 for MV < 5 law form to a log-normal form, likely to occur through −10.4 (equivalent to M & 10 M⊙), and αL =1.7 ± 0.2 preferential disruption of lower-mass clusters. Thus, for fainter, lower-mass clusters. Fall (2006) has discussed studying nearby populations of SSCs may shed light on the relations between luminosity and mass functions in processes similar to those occurring in the early phases more detail. of galaxy formation. For unresolved clusters, conversion from luminosity to From a practical point of view, young SSCs are the eas- mass requires knowledge of the age, metallicity, and red- iest clusters to identify in magnitude limited surveys of dening of the cluster together with a model for the clus- external galaxies, and so corrections for incompleteness ter’s stellar evolution and an assumed initial mass func- are relatively small, allowing for more accurate deter- tion (IMF) for the stars. Multi-color photometry is a mination of the ICMF. However, a disadvantage is that minimum observational requirement for performing this SSCs are relatively rare, so a large sample of galaxies conversion on a cluster by cluster basis and is the ap- needs to be surveyed to build up a good statistical sam- proach adopted in this paper using Sloan Digital Sky Sur- ple. vey (SDSS) data. Billett, Hunter, & Elmegreen (2002) While SSCs occur in many different types of galac- performed a similar analysis for nearby dwarf irregular tic environments, there are many prominent examples galaxies using a heterogeneous dataset from the HST in irregular and dwarf irregular galaxies. For example, archives. Their study has the advantage of high spa- 6 NGC 5253 hosts a 10 M⊙ SSC (Turner et al. 2000), tial resolution, but the disadvantage of a limited spec- Henize 2-10 hosts several forming SSCs as revealed in the tral coverage (and thus more of an uncertain conver- radio and mm continuum (Kobulnicky & Johnson 1999; sion from luminosity to mass). Hunter et al. (2003) used Johnson & Kobulnicky 2003), NGC 1569 hosts 3 SSCs ground-based data to study the cluster populations of (Gilbert & Graham 2003), SBS 0335-052 hosts a 2 × 106 the Large and Small Magellanic Clouds, concluding that M⊙ cluster (Plante & Sauvage 2002), and NGC 1705 2.0 . αM . 2.4. and IC 4662 host SSCs (Johnson, Indebetouw, & Pisano Alternatively, the conversion from luminosity to mass 2003). Several other dwarf starburst galaxies with SSCs can be made on a statistical basis, given assumptions are described by Beck, Turner & Kovo (2000). about the history. Larsen (2002) showed Billett, Hunter, & Elmegreen (2002) discussed how the that a cluster population with mass function with αM = 3 5 most massive cluster in a galaxy depends on the mass 2 from 10 M⊙

7 various rejection methods used in this paper can be found 3.0 × 10 M⊙ (Ma et al. 2006). Hasegan et al. (2005) in Table 3. studied massive clusters in M87, concluding that above 6 2 × 10 M⊙ there may be a transition to 4. DERIVATION OF THE ICMF AND COMPARISON BETWEEN IRREGULAR AND SPIRAL GALAXIES dominated, ultracompact dwarf galaxies. Other possible explanations for the ultra massive clus- 4.1. Completeness Corrections ters in our sample are that they are simply blends of The completeness and sample bias corrections de- lower-mass clusters that happen to be relatively close to- scribed by Larsen (1999) were applied to all galaxies with gether. The derived sizes of these objects are on the order twenty or more detected clusters. For these tests, 30 sets of hundreds of parsecs, which lends some credence to this of randomly generated artificial clusters with masses be- explanation. Another possibility is that some of these ul- 4.4 6.5 ′ tween ∼ 10 − 10 M⊙ were inserted into the g im- tra massive systems are foreground stellar or background age of each galaxy. The number of artificial clusters was galactic contaminants that have passed through the vari- set to one-quarter the number of clusters found. This ous color filters. This is supported by the unusually large dynamic factor helps to reduce overcrowding in smaller reddening values associated with the clusters. galaxies while better sampling the area of larger galax- The small number of ultra massive cluster candidates ies. The magnitudes used for these test clusters came permits us to examine them in more detail. The mas- from the appropriate metallicity Starburst99 models. We sive cluster in NGC 4449 is located about two optical limited the completeness and bias corrections to clusters radii from the galaxy center and is listed in the SDSS 4.4 more massive than ∼ 10 M⊙ to avoid contamination photometric database and the USNO catalog as a star. 8 contributed by supergiants and stellar blends. The arti- The 1.12 × 10 M⊙ cluster in NGC 5457 is associated ficial clusters used in this test were uniformly distributed with a nuclear region of diffuse emission rather than a in age between 2 and 20 Myr. Masses were sampled from single object. The other cluster in NGC 5457 is listed a power law distribution with an index of −1.7. The in the SDSS photometric database as a galaxy. Thus, spatial distribution of these clusters was determined by the three candidate ultra massive clusters in our sample fitting an ellipse to the observed cluster distribution and are contaminants rather than true clusters and have been altering the semi-major axis until & 90% of the observed excluded from the statistical analysis of the ICMF. clusters were contained in this ellipse. Once these param- 4.3. eters of the ellipse were determined, the artificial clusters Variation of Cluster Frequency within the Galaxy were randomly distributed throughout its interior. As a Samples final measure to insure that clusters that we had gen- It should be noted that the frequency at which mas- erated were as authentic as possible, random amounts sive clusters appear in our galaxy samples varies greatly. of reddening were added to the model magnitudes. The From Table 1 we see that three galaxies dominate the 5 reddening was uniformly distributed between an AV of irregular sample of massive (> 10 M⊙) clusters: NGC 0.02 mag and the maximum reddening magnitudes de- 4449, NGC 4485 and NGC 4656. These are all inter- termined for the particular galaxy. At this point, the acting galaxies (see §5). Using our derived Hubble flow artificial clusters were added to the images using IRAF’s distances and the apparent B-band magnitudes from the ‘mkobjects’ routine. RC3 yields absolute magnitudes of -18.8, -17.5, and - Once the artificial clusters were added to the images, 18.7, respectively. This range of magnitudes places these the SExtractor routine was executed on the images us- galaxies at the high end of the irregular galaxy luminos- ing the parameters found in §2.1. The returned object ity distribution. Only one of these galaxies, NGC 4656, list was then compared with the list of artificial clusters has a classification other than ImB. Its classification of to determine the number of clusters recovered. These SB/SM(p) is likely due to its optical appearance (see Fig- results were then averaged across the various samples to ure 1) being almost linear, with the possible remnant of arrive at the completeness corrections presented in Table a spiral arm towards the northeast. 4 for both the irregular and spiral samples. For the lowest In the small spiral sample, NGC 5457 contributes 4.40 4.64 mass bin that we consider, 10 −10 M⊙, we estimate about 90% of the massive clusters. Thus our compar- that our method of finding young clusters is 70% com- ison of the mass functions of the two galaxy samples is plete for the irregular galaxy sample and 71% complete in essence a comparison of the clusters in NGC 5457 with for the spiral galaxy sample. Higher angular resolution those in NGC 4449, NGC 4485 and NGC 4656. Larger optical searches by HST would help us to probe regions samples of galaxies are needed to improve upon this sit- of diffuse stellar light deeper and detect the less massive uation. clusters with a higher detection rate. 4.4. The Initial Cluster Mass Functions 4.2. Ultra Massive Clusters The young (. 20 Myr) clusters were binned by mass Our method of cluster identification finds a small num- using the bin sizing method described in Scott (1979). 7.5 ber of ultra massive (Mcluster > 10 M⊙) clusters. If This sizing procedure takes into consideration the stan- these objects are true star clusters, then their derived dard deviation of the data along with the total number masses are comparable to that of the ultra massive clus- of data points to reduce any bias introduced by binning. ter W3 in NGC 7252 (Maraston et al. 2004), which is For our data set, this resulted in a bin size of 0.24 dex. 7 estimated to have a mass of 8 × 10 M⊙ and a radius We applied completeness corrections to each bin and as- of about 18 pc. Fellhauer & Kroupa (2005) suggested sumed a combination of Poisson and uniform systematic that W3 may have formed from the merger of lower- uncertainties of 30% in the bin amplitudes. Histograms mass clusters. The most massive cluster in the Local for both samples are presented in Figure 3. In each case, Group is thought to be a in M31 with the distribution of clusters spans approximately the same 6 range in mass. Thus, despite being physically smaller higher frequency of SSCs, but larger samples of clusters and possessing different internal conditions, irregulars are needed to confirm this. are capable of producing clusters with the same mass We have also compared the spiral and irregular raw range as spirals. ICMFs (i.e. with no completeness corrections applied) dN M 5 ( ) −αM for masses above 10 M⊙ using a Kolmogorov-Smirnov We fit power laws, dM ∝ M , to the ICMF in 5 7.5 test. There are 243 such clusters in the irregular galax- the range 10 M⊙

Irregular galaxies, being typically of later Hubble type, galaxies in this subsample are observed with clusters with 6.5 smaller mass, and higher gas mass fraction compared to masses ≥ 10 M⊙. This is tentative evidence in favor normal spiral galaxies, tend to have much lower metal- of a more top heavy ICMF in interacting galaxies. licities. This does not appear to have any bearing on the ICMF. Note that the range in metallicity of the Star- 5.4. Continuity of the ICMF burst99 models considered, Z∼ 0.004 − 0.02, is not an accurate assessment of the range in metallicities of our Since the derived star cluster masses presented for the irregular galaxy sample cover a wide range of masses, we galaxy samples, and the metallicity recorded in Table 2 are able to examine the question of the continuity of the for each galaxy is a somewhat arbitrary choice intended to maximize the number of clusters found by the color3 ICMF discussed in Billett, Hunter, & Elmegreen (2002). The significant linear correlation p-value (< 0.01) for the method (see §3.2). Good measurements of the metallicity power law fits implies that the ICMF is well modeled by of the galaxies in our samples are not generally available, 5 7 a single power law in the ∼ 10 −10 M⊙ mass range. We so a more quantitative discussion of the range of metal- licity over which the ICMF is invariant is not possible at conclude that the ICMF in our sample of irregulars is well described by a continuous power law. Due to complete- this stage. ness and stellar contamination issues in the lower mass One mechanism by which metallicity might be ex- pected to affect star formation involves the degree of bins, the only conclusion that can be drawn about the is- sue of continuity is that if a break in the power law exists, ionization and heating inside molecular clouds due to × 4 Far UV photons. In a higher metallicity cloud with a it occurs for clusters with masses < 3 10 M⊙. There is tentative evidence for a cut-off at the high-mass end higher dust-to-gas ratio and higher mean extinction, the Far UV photons do not penetrate as easily. If magnetic of the irregular ICMF. In Figure 4, we see that this sam- ple has essentially no clusters more massive than about fields are important for preventing GMC collapse and 6 3 × 10 M⊙. The fact that the number of clusters in the thus regulating star formation (McKee 1989) and if the ≫ creation of magnetically super-critical regions occurs pri- mass bin just below this mass is about 20 1, suggests that this is a real, physical effect, and not due to small marily via ambipolar diffusion in regions of low ioniza- tion, then one expects such regions to require smaller number statistics (see discussion in Williams & McKee 1997). amounts of surrounding, shielding gas in a higher metal- licity galaxy. In other words, shielding a given mass of It should be noted that a power law index of αM < 2 gas in a low-metallicity galaxy requires it to be in the implies that the total star formation occurring in clus- ters (in a given mass range) is dominated by that in the center of a more massive GMC, i.e. with a larger total largest clusters. If the ICMF does not steepen signifi- mass surface density, Σ, than in a high-metallicity galaxy. 4 cantly below 5 × 10 M⊙, then our results imply that the The pressure due to the self-gravity of the gas will be of order ≃ GΣ2, and the combination of higher pressure star formation in irregular and spiral galaxies is domi- nated by that occurring in super star clusters or super and more surrounding material may lead to higher star formation efficiencies and higher-mass clusters. See Tan associations. & McKee (2004) for a more extensive discussion of feed- 5.5. back processes in forming super star clusters. Since we do Evolution of Extinction not observe such a difference in the ICMFs of (presum- We may expect that younger clusters should suffer ably lower-metallicity) irregulars and (presumably high- from a greater amount of local extinction relative to metallicity) spirals, we conclude that either magnetically older clusters. To test this hypothesis, we examined super-critical regions are created by turbulent diffusion how the local extinction estimate for our clusters var- of field strength rather than ambipolar diffusion or that ied with cluster age. In Figure 5, we see an apparent other processes, such as turbulent fragmentation, set the decrease in the cluster extinction with age. We also ICMF on scales smaller that the super-critical regions of notice that there are relatively few (1.2% for irregu- GMCs. lars; 1.3% for spirals) clusters younger than 5 Myr that 5.3. have an AV less than 0.1 mag while 30% of clusters Effect of Galaxy-Galaxy Interactions in irregulars and 42% of clusters in spirals older than It has been suggested that galaxy-galaxy interactions 5 Myr have an AV less than 0.1 mag. This implies may trigger enhanced star formation (e.g. Lisenfeld et al. that 5-10 Myr is the characteristic timescale in which 2004), and that this could lead to the presence of more a cluster destroys or moves away from its natal molec- massive star clusters. We examined the SDSS g′ images ular cloud, which is consistent with estimates of this of the eight irregular galaxies that had sufficient clusters timescale based on studies of young Galactic star clusters for completeness corrections to look for evidence of an (Leisawitz, Bash, & Thaddeus 1989). Performing a lin- interaction. This resulted in the galaxies NGC 4485 and ear fit between the extinctions and ages we derived that NGC 4656 being identified as interacting. NGC 4190 was an average cluster in an irregular galaxy has AV decrease also identified as interacting by Dahari (1985) and NGC by 0.091 mag Myr−1, while an average cluster in a spiral −1 4449 by Hunter et al. (1998). To determine the effect galaxy has AV decrease by 0.053 mag Myr . of galaxy-galaxy interactions on the ICMF, we ratio the We examined the mass dependence of these relations, 6.0 4.4 5 number of star clusters with masses ≥ 10 M⊙ to the defining a “low-mass” sample from 10 − 10 M⊙ and 5 6 7.5 number with masses ≥ 10 M⊙. The lower mass limit a “high-mass” sample from 10 − 10 M⊙. The high- was chosen to insure that each galaxy is complete to 90%. mass clusters in irregular and spiral galaxies start with From this, we see that the average non- high values of mean AV at zero age: 1.69 and 1.62 mag, has a ratio of 0.118 while the average interacting galaxy respectively. The corresponding values for the low-mass has a ratio of 0.143. We also see that only the interacting clusters are 0.83 and 0.41 mag. The rates of decrease 8

−1 of AV are -0.145 and -0.144 mag Myr for the high- adjacent smaller clusters. Nevertheless, this does not af- mass clusters in irregular and spiral galaxies, respec- fect our conclusion on the similarity of the irregular and tively, while they are -0.073 and -0.026 mag Myr−1 for spiral ICMFs since these data were analyzed in the same the low-mass clusters, respectively. Thus the 5-10 Myr way and share the same systematic errors. Furthermore, timescale for reduction in extinction is similar for all we are able to conclude that most stars in irregular and these sub-samples. spiral galaxies form in massive clusters or massive asso- ciations. 6. SUMMARY We find a characteristic timescale of about 5-10 Myr We have presented the age and masses determined for for clusters to move away from or destroy their natal gas 439 young (. 20 Myr) star clusters in thirteen irregular and dust clouds. galaxies and 611 young star clusters in three spiral galax- ies based on an automated analysis of SDSS data. From this, we have found the irregular ICMF to be well fit by dN(M) −αM a power law dM ∝ M with αM = 1.62 ± 0.08 We thank an anonymous referee for detailed comments 5 7.5 for clusters over the 10 − 10 M⊙ mass range. The that led to a significantly improved paper. This research equivalent index for the spiral ICMF was found to be has made use of the NASA/IPAC Extragalactic Database αM = 1.74 ± 0.07. The irregular and spiral ICMFs are (NED) as well as data from the fourth data release of the statistically indistinguishable in spite of expected differ- Sloan Digital Sky Survey. This work has been supported ences in galactic shear and metallicity. by the University of Wyoming Research Experience for Our derived ICMFs are slightly top heavy compared to Undergraduates program under NSF REU grant AST- some previous studies, perhaps because of the generally 0353760. JCT acknowledges support from a Zwicky fel- poorer spatial resolution of the SDSS data; it is likely lowship from the Inst. of Astronomy, ETH Z¨urich and that many of our more massive “clusters” are blends of CLAS, University of Florida..

REFERENCES Adelman-McCarthy, J., et al. 2007, ApJS, submitted. Leroy, A., Bolatto, A., Simon, J., & Blitz, L. 2005, ApJ, 625, 763. Arp, H. & Sandage, A. 1985, AJ, 90, 1163. Lisenfeld, U., Verdes-Montenegro, L., Espada, D., Garcia, E., Leon, Beck, S., Turner, J., & Kovo, O. 2000, AJ, 120, 244 S., Sabater, J., Sulentic, J., & Verley, S. 2004, in ?The Evolution Bertin, E. & Arnouts, S. 1996, A&AS, 117, 393. of Starbursts?, ed. S. Huettemeister et al., AIP Proceedings, Billett, O., Hunter, D., & Elmegreen, B. 2002, AJ, 123, 1454. (astro-ph/0411443). Clark, J., Negueruela, I., Crowther, P., & Goodwin, S. 2005, A&A, Lupton, R., Gunn, J., & Szalay, A. 1999, AJ, 118, 1406. 434, 949. Ma, J., de Grijs, R., Yang, Y., Zhou, X., Chen, J., Jiang, Z., Wu, Cresci, G., Vanzi, L., & Sauvage, M. 2004, A&A, in press Z., & Wu, J. 2006, MNRAS, 368, 1443 [astro-ph/0411486] Ma´ız-Alpell´aniz, J., Cieza, L., & MacKenty, J. 2002, AJ, 123, 1307. Dahari, O. 1985, ApJS, 57, 643. Maraston, C., Bastian, N., Saglia, R. P., Kissler-Patig, M., de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H., Buta, R., Schweizer, F., Goudfrooij, P. 2004, A&A, 416, 467 Paturel, G., & Fouqu´e, P. 1991, Third Reference Catalogue of McCrady, N., Gilbert, A., & Graham, J. 2003, ApJ, 596, 240. Bright Galaxies (New York:Springer) McKee, C. 1989, ApJ, 345, 782 Elmegreen, B. & Efremov, Y. 1997, ApJ, 480, 235. Mengel, S., Lehnert, M., Thatte, N., & Genzel, R. 2002, A&A, 383, Fall, S. 2006, ApJ, in press (astro-ph/0609201) 137. Fall, S., Chandar, R., & Whitmore, B. 2005, ApJ, 631, L133 Patureal G., Petit C., Prugniel P., Theureau G., Rousseau J., Fellhauer, M. & Kroupa, P. 2005, MNRAS, 359, 223. Brouty M., Dubois P., Cambresy L. 2003, A&A, 412, 45. Gilbert, A., & Graham, J. 2003, in: Extragalactic Globular Clusters Plante, S. & Sauvage, M. 2002, AJ, 124, 1995 and their Host Galaxies, 25th meeting of the IAU, Joint Schlegel, D., Finkbeiner, D., & Davis, M. 1998, ApJ, 500, 525. Discussion 6, 17 July 2003, Sydney, Australia. Scott, D. 1979, Biometrika, 66, 605. Gonz´alez-Delgado, R., Leitherer, C., Heckman, T., & Cervi˜o, M. Sufue, Y., Tutui, U., Honma, M., Tomita, A., Takamiya, T., Koda, 1997 ApJ, 483, 705. J., & Takeda, Y. 1999, ApJ, 523, 136. Hasegan, M., et al.. 2005, ApJ, 627, 203. Smith, J. et al. 2002, AJ, 123, 2121. Hunt, L., Thuan, T., & Izotov, Y. 2003, ApJ, 588, 281. Tan, J. 2000, ApJ, 536, 173 Hunter, D., Wilcots, E., van Woerden, H., Gallager, J., & Kohle, Tan, J. 2005, in Cores to Clusters: Star Formation with Next S. 1998, ApJ, 495, L47. Generation Telescopes, eds. M. S. N. Kumar, M. Tafalla, P. Hunter, D., O’Connell, R., Gallagher, J., & Smecher-Hane, T. 2000, Caselli, (Springer), p87, (astro-ph/0504256) AJ, 120, 2383. Tan, J., & McKee, C. 2004, in The Formation and Evolution of Hunter, D., Elmegreen, B., Dupy, T., & Mortonson, M. 2003, AJ, Massive Young Clusters in Giant Molecular Clouds, ASP Conf. 126, 1836. Ser. 322, eds. H.J.G.L.M. Lamers, L.J. Smith, and A. Nota. (San Johnson, K., Indebetouw, R., & Pisano, D. 2003, AJ, 126, 101. Francisco: ASP), 263 Johnson, K. & Kobulnicky, H. 2003, ApJ, 597, 923. Terry, J., Paturel, G., & Ekholm, T. 2002, A&A, 393, 57. Kennicutt, R. 1998, ApJ, 498, 541. Tremonti, C., Calzetti, D., Leitherer, C., & Heckman, T. 2001, Kobulnicky, H. & Johnson, K. 1999, ApJ, 527, 154. ApJ, 555, 322 Lada, C. & Lada, E. 2003, ARA&A, 41, 57. Turner, J. L., Beck, S., & Ho, P. 2000, ApJ, 532, 109 Larsen, S. 1999, A&AS, 139, 393. Whitmore, B., Zhang, Q., Leitherer, C., Fall, S., Schweizer, F., & Larsen, S. 2002, AJ, 124, 1393. Miller, B. 1999, AJ, 118, 1551. Leisawitz, D., Bash, F., & Thaddeus, P. 1989, ApJS, 70, 731. Williams, J., & McKee, C. F. 1997, ApJ, 476, 166 Leitherer, C., Schaerer, D., Goldader, J., Delgado, R., Robert, C., Zhang, Q. & Fall, S. 1999, ApJ, 527, L81. Kune, D., de Mello, D., Devost, D., & Heckman, T. 1999, ApJS, 123, 3. 9

Fig. 1.— g′ images of irregular galaxies NGC 2552 (upper left), NGC 3738 (upper right), NGC 4190 (center left), NGC 4214 (center right), NGC 4449 (lower left), and NGC 4485 (lower right) with candidate clusters circled. The circles indicate the location of young (< 20 Myr) clusters while the squares indicate the location of older clusters. In each case the size of the marker is proportional to the photometric aperture used. The black bar shown in each image is 30”. 10

Fig. 1.— Continued — g′ images of NGC 4656 (upper left), IC 3521 (upper right), UGC 6541 (center left), UGC 7408 (center right), DDO 165 (lower left), and DDO 167 (lower right). 11

Fig. 1.— Continued — g′ image of DDO 168 (upper left). Also shown are g′ images of spiral galaxies NGC 4571 (upper right), NGC 4713 (lower left), and NGC 5457 (lower right). 12

Fig. 2.— Color3 plots of g′ − i′ vs. u′ − g′ colors (top), g′ − i′ vs. i′ − z′ colors (center), and u′ − g′ vs. i′ − z′ colors (bottom). Left Panel – This series demonstrates how color3 derives an internal reddening of 0.48 mag and an age of 4 Myr for cluster 2 in DDO 165. The ellipses represent projections of the 1σ uncertainty ellipsoid. Ages along the solar metallicity model are indicated by crosses at 1, 5, 10, 20, 100, and 1,000 Myr. The arrow in each plot represents reddening of AV = 1.00 mag. Right Panel – The locations of dereddened clusters (plus signs) with 1σ uncertainties for all clusters in DDO 165. Also shown in this series are the initial (Galactic reddening corrected) cluster locations (crosses), clusters rejected by the color filters (diamonds), and clusters reject by color3 as being more than 3σ and 0.50 mag from the model (triangles). 13

Fig. 3.— Initial cluster mass function for irregular (left) and spiral galaxies (right), based on the mass histograms of young (< 20 Myr) clusters identified in these galaxy samples. The uncorrected data are shown with a dashed lines, while the completeness corrected data (§4.1) are shown with solid lines. 14

Fig. 4.— Plots of the best power law fits for the ICMFs in irregular (solid lines) and spiral galaxies (dashed lines). Assuming the uncertainties in bin amplitude is a combination of Poisson plus 30% fixed systematic added in quadrature, the best fit power laws over the 5 7.5 mass range 10 − 10 M⊙ are shown with a solid line with power law index −1.62 ± 0.08 for irregulars and a dashed line with power law index −1.74 ± 0.07 for spirals. Statistically, these two distributions are indistinguishable. 15

Fig. 5.— Visual extinction versus age of young massive clusters in irregular (left) and spiral (right) galaxies. We see an apparent evolution in the extinction of clusters as they age. Clusters younger than 5 Myr always have at least a few tenths of a magnitude of local V band extinction. Clusters with ages 5 Myr or older have a &50% chance of being relatively free (AV < 0.1 mag) of extinction. This 5-10 Myr timescale is thus the characteristic time for massive clusters to destroy or move away from their natal dust clouds. Note that some of the relatively high values of extinction at later ages are likely caused by our line-of-sight intercepting local, but unrelated, clouds in the host galaxies. The concentration of clusters at particular ages is due to extremities of the tracks of model cluster evolution in color space, to which the color3 method concentrates clusters with relatively large photometric uncertainties by adjusting the amount of extinction. Thus some of the scatter in the derived extinction is also likely to be due to a combination of photometric errors coupled with systematic uncertainties in the positions of the tracks of model clusters in color space. Nevertheless the evolution of the mean cluster extinction does appear to show a real decrease to very small values over about 5 to 10 Myr. 16

TABLE 1 Properties of the Galaxies Used in this Survey

Galaxy R.A. Dec. Hubble Stage MB v Distance Ncluster −1 5 (J2000.0) (Type) (km s ) (Mpc) M > 10 M⊙

Irregulars NGC 2552 08h19m20s.53 +50◦00′34′′. 7 9.0 (SA(s)m) -17.3 699 9.32 5 NGC 3738 11h35m48s.79 +54◦31′26′′. 0 10.0 (Im) -16.8 464 6.19 8 NGC 4190 12h13m44s.73 +36◦38′02′′. 9 10.0 (Im(p)) -14.8 415 5.53 7 NGC 4214 12h15m39s.17 +36◦19′36′′. 8 10.0 (IAB(s)m) -17.1 ··· 2.98 9 NGC 4449 12h28m11s.90 +44◦05′39′′. 6 10.0 (ImB) -18.8 422 5.63 55 NGC 4485 12h30m31s.13 +41◦42′04′′. 2 10.0 (ImB(sp)) -17.5 690 9.20 90 NGC 4656 12h43m57s.73 +32◦10′05′′. 3 9.0 (SB/Sm(p)) -18.7 646 8.61 65 IC 3521 12h34m39s.50 +07◦09′37′′. 0 10.0 (ImB) -16.0 663 8.84 1 UGC 6541 11h33m28s.90 +49◦14′14′′. 0 10.0 (Im) -13.0 249 3.32 3 UGC 7408 12h21m15s.01 +45◦48′40′′. 8 10.0 (ImA) -16.4 683 9.11 2 DDO 165 13h06m24s.85 +67◦42′25′′. 0 10.0 (Im) -15.3 312 4.16 1 DDO 167 13h13m22s.73 +46◦19′13′′. 2 10.0 (Im) ··· 400 5.33 0 DDO 168 13h14m27s.95 +45◦55′08′′. 9 10.0 (ImB) -16.1 426 5.68 1

Spirals NGC 4571 12h36m56s.37 +14◦13′02′′. 5 6.5 (SAd(r)) -17.0 435 5.80 2 NGC 4713 12h49m57s.87 +05◦18′41′′. 1 7.0 (SABd(rs)) -17.7 715 9.54 33 NGC 5457 14h03m12s.59 +54◦20′56′′. 7 6.0 (SABcd(rs)) -20.8 503 6.71 379

References. — Coordinates and radial velocities were taken from NASA/IPAC Extragalactic Database (NED), and morphologies are from RC3. Distances are from LEDA (Patureal et al. 2003) except for Ma´ız-Alpell´aniz, Cieza, & MacKenty (2002) (NGC 4214). MB s are derived from the apparent B-band magnitudes found in RC3 and our adpoted distances. 17 ) ⊙ 05 18 19 15 04 32 10 22 11 28 94 18 39 52 18 26 34 55 20 46 31 45 74 45 ...... 0 0 0 0 0 0 0 0 0 0 0 0 − +0 − +0 − +0 − − +0 +0 − − +0 +0 − +0 − +0 − +0 − +0 − +0 M/M 73 20 11 64 51 51 22 10 12 87 58 42 ( ...... 5 5 7 4 6 4 5 5 6 4 4 5 log 09 10 . 00 00 60 60 00 10 00 00 00 00 00 40 20 00 60 99 00 00 . 30 ...... 99 00 00 . . . 1 1 3 3 1 0 1 1 . 1 1 0 1 − +1 − +1 − +1 − − +1 +1 − +3 − +0 − +1 − +1 +18 − − +1 +177 − Age Mass 4 91 71 61 11 21 61 . 11 71 61 (Myr) ...... 91 . o 0.04 3 0.03 7 0.02 3 0.03 5 0.13 9 0.02 5 0.15 10 0.08 5 ) 0.51 3 0.12 13 0.22 0 0.14 7 ′ ± ± ± ± ± ± ± ± z ± ± ± ± ks. M ( 1 Data o 0.02 -11.47 0.02 -11.09 0.01 -13.19 0.03 -9.98 0.01 -12.96 0.03 -10.49 0.02 -12.27 0.02 -11.26 ) 0.07 -8.94 0.05 -9.20 0.04 -10.11 0.06 -8.68 ′ ± ± ± ± ± ± ± ± i ± ± ± ± M ( 008) 020) . . = 0 = 0 fit fit ,Z ,Z o 0.02 -11.87 0.02 -11.16 0.02 -13.20 0.02 -10.47 0.04 -9.45 0.01 -13.10 0.03 -10.43 0.02 -12.02 0.01 -11.23 1 1 ) 0.05 -8.86 0.03 -9.72 0.08 -8.43 ′ ± ± ± ± ± ± ± ± ± g ± ± ± 154 034 . . M ( = 0 = 0 TABLE 2 V,Gal. V,Gal. A A o 0.03 -12.31 0.02 -11.10 0.17 -12.21 0.04 -11.04 0.05 -10.11 0.03 -13.23 0.06 -10.47 0.05 -11.47 0.02 -11.94 ) 0.19 -8.62 0.05 -9.91 0.54 -7.39 ′ ± ± ± ± ± ± ± ± ± u ± ± ± nly move the cluster further from the Starburst99 model trac M ( NGC 2552 ( NGC 3738 ( from NASA/IPAC Extragalactic Database (NED) V,Gal. A V 2 A Detected Star Clusters from the Galaxies Found in the SDSS DR5 59 1.10 -12.25 93 0.00 -10.83 98 1.64 -10.96 37 0.81 -11.16 06 0.33 -10.24 68 1.51 -12.83 68 0.60 -8.22 51 0.36 -10.23 96 1.07 -10.76 65 0.17 -9.89 28 0.14 -11.62 82 0.38 -6.39 ...... ′′ ′′ ′′ ′′ ′′ ′′ ′′ ′′ ′′ ′′ ′′ ′′ 41 36 37 15 27 34 07 40 37 38 14 43 ′ ′ ′ ′ ′ ′ ′ ′ ′ ′ ′ ′ 31 31 01 00 01 00 31 00 31 00 00 31 ◦ ◦ ◦ ◦ ◦ ◦ ◦ ◦ ◦ ◦ ◦ ◦ (J2000.0) (mag) 92 54 51 54 43 50 48 50 53 50 53 50 78 54 71 50 75 54 51 50 29 50 59 54 ...... s s s s s s s s s s s s 46 46 21 24 17 20 46 18 45 18 17 45 m m m m m m m m m m m m 35 35 19 19 19 19 35 19 35 19 19 35 values of zero indicate that removing additional reddeing o h h h h h h h h h h h h V A .— Note 7 11 8 08 6 11 7 08 5 11 6 08 3 08 2 11 ID R.A. Dec. The table appears in its entirety in the electron edition. 10 11 17 08 11 08 16 08 1 18

TABLE 3 Rejection Summary

Filter Number of Clusters Irregulars Spirals

Initial Sample 4447 1738 Poor Photometry 2908 409 Rejection by color3 573 344 Old (> 20 Myr) 390 293 4.4 Low Mass (M < 10 M⊙) 136 79 7.5 High Mass (M > 10 M⊙) 1 2 Final Sample 439 611

Note. — Poor Photometry refers to objects rejected due to photometric errors > 10% or having an aperture > 1.5 times the median stellar 3 3 aperature. Rejection by color refers to objects rejected by the color routine either by emission line contaimination or lying more than 3σ and 0.50 mag from the Starburst99 models. 19 55 . 6 − 31 . 6 10 ⊙ M 31 . 6 − 07 . 6 10 ⊙ M 07 . 6 − 84 . 5 10 ⊙ M 84 . 5 − 60 . 5 10 ⊙ M 60 . 5 − 36 . 5 10 ⊙ M 36 . Spirals 5 Irregulars TABLE 4 − Cluster Recovery Rate 12 20 Myr) clusters in each galaxy. . 5 . 10 ⊙ M f young ( 12 . 5 − 88 . 4 Completeness Corrections for the Galaxy Samples 10 ⊙ M 88 . 4 − 64 . 4 10 ⊙ M 64 . 4 − 0.70 0.840.71 0.93 0.88 0.96 0.95 0.96 0.98 0.98 0.98 0.98 0.98 0.97 0.98 0.97 0.98 0.98 40 . 4 ··· ··· (Mpc) 10 . — The average for the corrections is weighted by the number o Galaxy Distance Note Average Average DDO 168 5.68 0.99 0.99 0.99 0.99 0.99 0.99 0.99 1.00 1.00 NGC 3738NGC 4190NGC 4214 6.19NGC 4449 5.53NGC 4485 2.98NGC 4656 5.63UGC 6541 9.20 0.77 8.61 0.95 3.32 0.95 0.74NGC 0.48 4571NGC 0.61 4713NGC 0.93 0.85 5457 5.80 0.99 9.54 0.97 6.71 0.86 0.71 0.82 0.94 0.94 0.89 0.51 0.97 0.72 0.97 0.93 0.87 0.93 0.98 0.98 0.94 0.74 0.99 0.89 0.98 0.95 0.94 0.96 1.00 0.97 0.95 0.89 0.99 0.96 0.98 0.97 0.94 0.96 1.00 0.97 0.97 0.97 1.00 0.97 0.99 0.98 0.97 0.97 1.00 0.97 0.98 1.00 0.99 0.97 0.99 0.99 0.96 0.98 1.00 0.99 1.00 1.00 0.99 0.98 0.98 0.98 0.95 0.97 1.00 0.99 0.93 1.00 0.99 0.98 0.96 0.99 0.96 0.97 0.94 1.00 0.99 0.97 1.00 1.00 0.98