Hierarchical Star-Formation in M33: Fundamental properties of the star-forming regions

Star-formation within galaxies appears on multiple scales, from spiral structure, to OB associations, to individual star clusters, and often sub-structure within these clusters. This multitude of scales calls for objective methods to find and classify star-forming regions, regardless of spatial size. To this end, we present an analysis of star-forming groups in the local group spiral galaxy M33, based on a new implementation of the Minimum Spanning Tree (MST) method. Unlike previous studies which limited themselves to a single spatial scale, we study star-forming structures from the effective resolution limit (~20pc) to kpc scales. We find evidence for a continuum of star-forming group sizes, from pc to kpc scales. We do not find a characteristic scale for OB associations, unlike that found in previous studies, and we suggest that the appearance of such a scale was caused by spatial resolution and selection effects. The luminosity function of the groups is found to be well represented by a power-law with an index, -2, similar to that found for clusters and GMCs. Additionally, the groups follow a similar mass-radius relation as GMCs. The size distribution of the groups is best described by a log-normal distribution and we show that within a hierarchical distribution, if a scale is selected to find structure, the resulting size distribution will have a log-normal distribution. We find an abrupt drop of the number of groups outside a galactic radius of ~4kpc, suggesting a change in the structure of the star-forming ISM, possibly reflected in the lack of GMCs beyond this radius. (abridged)Comment: 12 pages, 16 figures, accepted MNRA

Luminosity Functions of Spitzer Identified Protostars in Nine Nearby Molecular Clouds

We identify protostars in Spitzer surveys of nine star-forming molecular clouds within 1 kpc: Serpens, Perseus, Ophiuchus, Chamaeleon, Lupus, Taurus, Orion, Cep OB3, and Mon R2, which combined host over 700 protostar candidates. Our diverse cloud sample allows us to compare protostar luminosity functions in these varied environments. We combine photometry from 2MASS J, H, and Ks bands and Spitzer IRAC and MIPS 24 micron bands to create 1 - 24 micron spectral energy distributions (SEDs). Using protostars from the c2d survey with well-determined bolometric luminosities (Lbol), we derive a relationship between Lbol, L_MIR (integrated from 1 - 24 microns), and SED slope. Estimations of Lbol for protostar candidates are combined to create luminosity functions for each cloud. Contamination due to edge-on disks, reddened Class II sources, and galaxies is estimated and removed from the luminosity functions. We find that luminosity functions for high mass star forming clouds peak near 1 Lsun and show a tail extending toward luminosities above 100 Lsun. The luminosity functions of the low mass star forming clouds do not exhibit a common peak, however the combined luminosity function of these regions peaks below 1 Lsun. Finally, we examine the luminosity functions as a function of the local surface density of YSOs. In the Orion molecular cloud, we find a significant difference between the luminosity functions of protostars in regions of high and low stellar density, the former of which is biased toward more luminous sources. This may be the result of primordial mass segregation, although this interpretation is not unique. We compare our luminosity functions to those predicted by models and find that our observed luminosity functions are best matched by models which invoke competitive accretion, although we do not find strong agreement of the high mass star forming clouds with any of the models.Comment: 76 pages, 18 figures, 7 tables. Accepted for publication in the Astronomical Journa

Spitzer Imaging of the Nearby Rich Young Cluster, Cep OB3b

We map the full extent of a rich massive young cluster in the Cep OB3b association with the IRAC and MIPS instruments aboard the {\it Spitzer} Space Telescope and the ACIS instrument aboard the $\it{Chandra}$ X-Ray Observatory. At 700 pc, it is revealed to be the second nearest large ($>1000$ member), young ($< 5$ Myr) cluster known. In contrast to the nearest large cluster, the Orion Nebula Cluster, Cep OB3b is only lightly obscured and is mostly located in a large cavity carved out of the surrounding molecular cloud. Our infrared and X-ray datasets, as well as visible photometry from the literature, are used to take a census of the young stars in Cep OB3b. We find that the young stars within the cluster are concentrated in two sub-clusters; an eastern sub-cluster, near the Cep B molecular clump, and a western sub-cluster, near the Cep F molecular clump. Using our census of young stars, we examine the fraction of young stars with infrared excesses indicative of circumstellar disks. We create a map of the disk fraction throughout the cluster and find that it is spatially variable. Due to these spatial variations, the two sub-clusters exhibit substantially different average disk fractions from each other: $32$ and $50$. We discuss whether the discrepant disk fractions are due to the photodestruction of disks by the high mass members of the cluster or whether they result from differences in the ages of the sub-clusters. We conclude that the discrepant disk fractions are most likely due to differences in the ages.Comment: 48 Pages, 12 figures, 6 table

Spitzer observations of the Orion OB1 association: disk census in the low mass stars

We present new Spitzer Space Telescope observations of two fields in the Orion OB1 association. We report here IRAC/MIPS observations for 115 confirmed members and 41 photometric candidates of the ~10 Myr 25 Orionis aggregate in the OB1a subassociation, and 106 confirmed members and 65 photometric candidates of the 5 Myr region located in the OB1b subassociation. The 25 Orionis aggregate shows a disk frequency of 6% while the field in the OB1b subassociation shows a disk frequency of 13%. Combining IRAC, MIPS and 2MASS photometry we place stars bearing disks in several classes: stars with optically thick disks (class II systems), stars with an inner transitional disks (transitional disk candidates) and stars with "evolved disks"; the last exhibit smaller IRAC/MIPS excesses than class II systems. In all, we identify 1 transitional disk candidate in the 25 Orionis aggregate and 3 in the OB1b field; this represents ~10% of the disk bearing stars, indicating that the transitional disk phase can be relatively fast. We find that the frequency of disks is a function of the stellar mass, suggesting a maximum around stars with spectral type M0. Comparing the infrared excess in the IRAC bands among several stellar groups we find that inner disk emission decays with stellar age, showing a correlation with the respective disk frequencies. The disk emission at the IRAC and MIPS bands in several stellar groups indicates that disk dissipation takes place faster in the inner region of the disks. Comparison with models of irradiated accretion disks, computed with several degrees of settling, suggests that the decrease in the overall accretion rate observed in young stellar groups is not sufficient to explain the weak disk emission observed in the IRAC bands for disk bearing stars with ages 5 Myr or older.Comment: Accepted in the Astrophysical Journa

The Spatial Distribution of Star Formation in the Solar Neighbourhood: Do all stars form in clusters?

We present a global study of low mass, young stellar object (YSO) surface densities in nearby (< 500 pc) star forming regions based on a comprehensive collection of Spitzer Space Telescope surveys. We show that the distribution of YSO surface densities in the solar neighbourhood is a smooth distribution, being adequately described by a lognormal function from a few to 10^3 YSOs per pc^2, with a peak at 22 stars/pc^2 and a dispersion of 0.85. We do not find evidence for multiple discrete modes of star-formation (e.g. clustered and distributed). Comparing the observed surface density distribution to previously reported surface density threshold definitions of clusters, we find that the fraction of stars in clusters is crucially dependent on the adopted definitions, ranging from 40 to 90%. However, we find that only a low fraction (< 26%) of stars are formed in dense environments where their formation/evolution (along with their circumstellar disks and/or planets) may be affected by the close proximity of their low-mass neighbours.Comment: 6 pages, 2 figures, MNRAS letters, accepte

A Structural Analysis of Star-Forming Region AFGL 490

We present Spitzer IRAC and MIPS observations of the star-forming region containing intermediate-mass young stellar object (YSO) AFGL 490. We supplement these data with near-IR 2MASS photometry and with deep SQIID observations off the central high extinction region. We have more than doubled the known membership of this region to 57 Class I and 303 Class II YSOs via the combined 1-24 um photometric catalog derived from these data. We construct and analyze the minimum spanning tree of their projected positions, isolating one locally over-dense cluster core containing 219 YSOs (60.8% of the region's members). We find this cluster core to be larger yet less dense than similarly analyzed clusters. Although the structure of this cluster core appears irregular, we demonstrate that the parsec-scale surface densities of both YSOs and gas are correlated with a power law slope of 2.8, as found for other similarly analyzed nearby molecular clouds. We also explore the mass segregation implications of AFGL 490's offset from the center of its core, finding that it has no apparent preferential central position relative to the low-mass members.Comment: 44 pages, 13 figures, accepted to Ap

Spitzer Observations of NGC 1333: A Study of Structure and Evolution in a Nearby Embedded Cluster

We present a comprehensive analysis of structure in the young, embedded cluster, NGC 1333 using members identified with Spitzer and 2MASS photometry based on their IR-excess emission. In total, 137 members are identified in this way, composed of 39 protostars and 98 more evolved pre-main sequence stars with disks. Of the latter class, four are transition/debris disk candidates. The fraction of exposed pre-main sequence stars with disks is 83% +/- 11%, showing that there is a measurable diskless pre-main sequence population. The sources in each of the Class I and Class II evolutionary states are shown to have very different spatial distributions relative to the distribution of the dense gas in their natal cloud. However, the distribution of nearest neighbor spacings among these two groups of sources are found to be quite similar, with a strong peak at spacings of 0.045 pc. Radial and azimuthal density profiles and surface density maps computed from the identified YSOs show that NGC 1333 is elongated and not strongly centrally concentrated, confirming previous claims in the literature. We interpret these new results as signs of a low velocity dispersion, extremely young cluster that is not in virial equilibrium.Comment: 59 pages, 20 figures, accepted to ApJ, verion with full resolution figures available at http://www.cfa.harvard.edu/~rgutermuth/preprints/gutermuth_ngc1333.pdf . Updated to fix astro-ph figure garblin