We measure the mass function for a sample of 840 young star clusters with ages between 10-300 Myr observed by the Panchromatic Hubble Andromeda Treasury (PHAT) survey in M31. The data show clear evidence of a high-mass truncation: only 15 clusters more massive than $10^4$ $M_{\odot}$ are observed, compared to $\sim$100 expected for a canonical $M^{-2}$ pure power-law mass function with the same total number of clusters above the catalog completeness limit. Adopting a Schechter function parameterization, we fit a characteristic truncation mass of $M_c = 8.5^{+2.8}_{-1.8} \times 10^3$ $M_{\odot}$. While previous studies have measured cluster mass function truncations, the characteristic truncation mass we measure is the lowest ever reported. Combining this M31 measurement with previous results, we find that the cluster mass function truncation correlates strongly with the characteristic star formation rate surface density of the host galaxy, where $M_c \propto$ $\langle \Sigma_{\mathrm{SFR}} \rangle^{\sim1.1}$. We also find evidence that suggests the observed $M_c$-$\Sigma_{\mathrm{SFR}}$ relation also applies to globular clusters, linking the two populations via a common formation pathway. If so, globular cluster mass functions could be useful tools for constraining the star formation properties of their progenitor host galaxies in the early Universe.

We use the Panchromatic Hubble Andromeda Treasury (PHAT) survey dataset to perform spatially resolved measurements of star cluster formation efficiency ($\Gamma$), the fraction of stellar mass formed in long-lived star clusters. We use robust star formation history and cluster parameter constraints, obtained through color-magnitude diagram analysis of resolved stellar populations, to study Andromeda's cluster and field populations over the last $\sim$300 Myr. We measure $\Gamma$ of 4-8% for young, 10-100 Myr old populations in M31. We find that cluster formation efficiency varies systematically across the M31 disk, consistent with variations in mid-plane pressure. These $\Gamma$ measurements expand the range of well-studied galactic environments, providing precise constraints in an HI-dominated, low intensity star formation environment. Spatially resolved results from M31 are broadly consistent with previous trends observed on galaxy-integrated scales, where $\Gamma$ increases with increasing star formation rate surface density ($\Sigma_{\mathrm{SFR}}$). However, we can explain observed scatter in the relation and attain better agreement between observations and theoretical models if we account for environmental variations in gas depletion time ($\tau_{\mathrm{dep}}$) when modeling $\Gamma$, accounting for the qualitative shift in star formation behavior when transitioning from a H$_2$-dominated to a HI-dominated interstellar medium. We also demonstrate that $\Gamma$ measurements in high $\Sigma_{\mathrm{SFR}}$ starburst systems are well-explained by $\tau_{\mathrm{dep}}$-dependent fiducial $\Gamma$ models.

We present a sample of 11 M31 Cepheids in stellar clusters, derived from the overlap of the Panchromatic Hubble Andromeda Treasury (PHAT) cluster catalog and the Pan-STARRS1 (PS1) disk Cepheid catalog. After identifying the PS1 Cepheids in the HST catalog, we calibrate the PS1 mean magnitudes using the higher resolution HST photometry, revealing up to 1 magnitude offsets due to crowding effects in the ground-based catalog. We measure ages of the clusters by performing single stellar population fits to their color-magnitude diagrams (CMDs) excluding their Cepheids. From these cluster age measurements, we derive an empirical period-age relation which agrees well with the existing literature values. By confirming this relation for M31 Cepheids, we justify its application in high-precision pointwise age estimation across M31.

We have undertaken the largest systematic study of the high-mass stellar initial mass function (IMF) to date using the optical color-magnitude diagrams (CMDs) of 85 resolved, young (4 Myr < t < 25 Myr), intermediate mass star clusters (10^3-10^4 Msun), observed as part of the Panchromatic Hubble Andromeda Treasury (PHAT) program. We fit each cluster's CMD to measure its mass function (MF) slope for stars >2 Msun. For the ensemble of clusters, the distribution of stellar MF slopes is best described by $\Gamma=+1.45^{+0.03}_{-0.06}$ with a very small intrinsic scatter. The data also imply no significant dependencies of the MF slope on cluster age, mass, and size, providing direct observational evidence that the measured MF represents the IMF. This analysis implies that the high-mass IMF slope in M31 clusters is universal with a slope ($\Gamma=+1.45^{+0.03}_{-0.06}$) that is steeper than the canonical Kroupa (+1.30) and Salpeter (+1.35) values. Using our inference model on select Milky Way (MW) and LMC high-mass IMF studies from the literature, we find $\Gamma_{\rm MW} \sim+1.15\pm0.1$ and $\Gamma_{\rm LMC} \sim+1.3\pm0.1$, both with intrinsic scatter of ~0.3-0.4 dex. Thus, while the high-mass IMF in the Local Group may be universal, systematics in literature IMF studies preclude any definitive conclusions; homogenous investigations of the high-mass IMF in the local universe are needed to overcome this limitation. Consequently, the present study represents the most robust measurement of the high-mass IMF slope to date. We have grafted the M31 high-mass IMF slope onto widely used sub-solar mass Kroupa and Chabrier IMFs and show that commonly used UV- and Halpha-based star formation rates should be increased by a factor of ~1.3-1.5 and the number of stars with masses >8 Msun are ~25% fewer than expected for a Salpeter/Kroupa IMF. [abridged]

We construct a stellar cluster catalog for the Panchromatic Hubble Andromeda Treasury (PHAT) survey using image classifications collected from the Andromeda Project citizen science website. We identify 2,753 clusters and 2,270 background galaxies within ~0.5 deg$^2$ of PHAT imaging searched, or ~400 kpc$^2$ in deprojected area at the distance of the Andromeda galaxy (M31). These identifications result from 1.82 million classifications of ~20,000 individual images (totaling ~7 gigapixels) by tens of thousands of volunteers. We show that our crowd-sourced approach, which collects >80 classifications per image, provides a robust, repeatable method of cluster identification. The high spatial resolution Hubble Space Telescope images resolve individual stars in each cluster and are instrumental in the factor of ~6 increase in the number of clusters known within the survey footprint. We measure integrated photometry in six filter passbands, ranging from the near-UV to the near-IR. PHAT clusters span a range of ~8 magnitudes in F475W (g-band) luminosity, equivalent to ~4 decades in cluster mass. We perform catalog completeness analysis using >3000 synthetic cluster simulations to determine robust detection limits and demonstrate that the catalog is 50% complete down to ~500 solar masses for ages <100 Myr. We include catalogs of clusters, background galaxies, remaining unselected candidates, and synthetic cluster simulations, making all information publicly available to the community. The catalog published here serves as the definitive base data product for PHAT cluster science, providing a census of star clusters in an L$^*$ spiral galaxy with unmatched sensitivity and quality.

We investigate the clustering of early-type stars younger than 300 Myr on galactic scales in M31. Based on the stellar photometric catalogs of the Panchromatic Hubble Andromeda Treasury program that also provides stellar parameters derived from the individual energy distributions, our analysis is focused on the young stars in three star-forming regions, located at galactocentric distances of about 5, 10, and 15 kpc, corresponding to the inner spiral arms, the ring structure, and the outer arm, respectively. We apply the two-point correlation function to our selected sample to investigate the clustering behavior of these stars across different time- and length-scales. We find that young stellar structure survives across the whole extent of M31 longer than 300 Myr. Stellar distribution in all regions appears to be self-similar, with younger stars being systematically more strongly clustered than the older, which are more dispersed. The observed clustering is interpreted as being induced by turbulence, the driving source for which is probably gravitational instabilities driven by the spiral arms, which are stronger closer to the galactic centre.

Herbig Ae/Be stars are intermediate-mass pre-main sequence stars surrounded by circumstellar dust disks. Some are observed to produce jets, whose appearance as a sequence of shock fronts (knots) suggests a past episodic outflow variability. This "jet fossil record" can be used to reconstruct the outflow history. We present the first optical to near-infrared (NIR) VLT/X-shooter spectra of the jet from the Herbig Ae star HD 163296. We determine physical conditions in the knots, as well as their kinematic "launch epochs". Knots are formed simultaneously on either side of the disk, with a regular interval of ~16 yr. The velocity dispersion versus jet velocity and the energy input are comparable in both lobes. However, the mass loss rate, velocity, and shock conditions are asymmetric. We find Mjet/Macc ~ 0.01-0.1, consistent with magneto-centrifugal jet launching models. No evidence for dust is found in the high-velocity jet, suggesting it is launched within the sublimation radius (<0.5 au). The jet inclination measured from proper motions and radial velocities confirms it is perpendicular to the disk. A tentative relation is found between the structure of the jet and the photometric variability of the source. Episodes of NIR brightening were previously detected and attributed to a dusty disk wind. We report for the first time significant optical fadings lasting from a few days up to a year, coinciding with the NIR brightenings. These are likely caused by dust lifted high above the disk plane; this supports the disk wind scenario. The disk wind is launched at a larger radius than the high-velocity atomic jet, although their outflow variability may have a common origin. No significant relation between outflow and accretion variability could be established. Our findings confirm that this source undergoes periodic ejection events, which may be coupled with dust ejections above the disk plane.

We present a probabilistic approach for inferring the parameters of the present day power-law stellar mass function (MF) of a resolved young star cluster. This technique (a) fully exploits the information content of a given dataset; (b) accounts for observational uncertainties in a straightforward way; (c) assigns meaningful uncertainties to the inferred parameters; (d) avoids the pitfalls associated with binning data; and (e) is applicable to virtually any resolved young cluster, laying the groundwork for a systematic study of the high mass stellar MF (M > 1 Msun). Using simulated clusters and Markov chain Monte Carlo sampling of the probability distribution functions, we show that estimates of the MF slope, {\alpha}, are unbiased and that the uncertainty, {\Delta}{\alpha}, depends primarily on the number of observed stars and stellar mass range they span, assuming that the uncertainties on individual masses and the completeness are well-characterized. Using idealized mock data, we compute the lower limit precision on {\alpha} and provide an analytic approximation for {\Delta}{\alpha} as a function of the observed number of stars and mass range. We find that ~ 3/4 of quoted literature uncertainties are smaller than the theoretical lower limit. By correcting these uncertainties to the theoretical lower limits, we find the literature studies yield <{\alpha}>=2.46 with a 1-{\sigma} dispersion of 0.35 dex. We verify that it is impossible for a power-law MF to obtain meaningful constraints on the upper mass limit of the IMF. We show that avoiding substantial biases in the MF slope requires: (1) including the MF as a prior when deriving individual stellar mass estimates; (2) modeling the uncertainties in the individual stellar masses; and (3) fully characterizing and then explicitly modeling the completeness for stars of a given mass. (abridged)

The apparent age and mass of a stellar cluster can be strongly affected by stochastic sampling of the stellar initial mass function, when inferred from the integrated color of low mass clusters (less than ~10^4 solar masses). We use simulated star clusters to show that these effects are minimized when the brightest, rapidly evolving stars in a cluster can be resolved, and the light of the fainter, more numerous unresolved stars can be analyzed separately. When comparing the light from the less luminous cluster members to models of unresolved light, more accurate age estimates can be obtained than when analyzing the integrated light from the entire cluster under the assumption that the initial mass function is fully populated. We show the success of this technique first using simulated clusters, and then with a stellar cluster in M31. This method represents one way of accounting for the discrete, stochastic sampling of the stellar initial mass function in less massive clusters and can be leveraged in studies of clusters throughout the Local Group and other nearby galaxies.

One of the key predictions of modeling from the IR excess of Herbig Ae stars is that for protoplanetary disks, where significant grain growth and settling has occurred, the dust disk has flattened to the point that it can be partially or largely shadowed by the innermost material at or near the dust sublimation radius. When the self-shadowing has already started, the outer disk is expected to be detected in scattered light only in the exceptional cases that the scale height of the dust disk at the sublimation radius is smaller than usual. High-contrast imaging combined with the IR spectral energy distribution allow us to measure the degree of flattening of the disk, as well as to determine the properties of the outer disk. We present polarimetric differential imaging in $H$ band obtained with Subaru/HiCIAO of one such system, MWC 480. The HiCIAO data were obtained at a historic minimum of the NIR excess. The disk is detected in scattered light from 0\farcs2-1\farcs0 (27.4-137AU). Together with the marginal detection of the disk from 1998 February 24 by HST/NICMOS, our data constrain the opening half angle for the disk to lie between 1.3$\leq\theta\leq 2.2^\circ$. When compared with similar measures in CO for the gas disk from the literature, the dust disk subtends only $\sim$30% of the gas disk scale height (H/R$\sim$0.03). Such a dust disk is a factor of 5-7 flatter than transitional disks, which have structural signatures that giant planets have formed.

The Panchromatic Hubble Andromeda Treasury (PHAT) survey is an on-going Hubble Space Telescope (HST) multi-cycle program to obtain high spatial resolution imaging of one-third of the M31 disk at ultraviolet through near-infrared wavelengths. In this paper, we present the first installment of the PHAT stellar cluster catalog. When completed, the PHAT cluster catalog will be among the largest and most comprehensive surveys of resolved star clusters in any galaxy. The exquisite spatial resolution achieved with HST has allowed us to identify hundreds of new clusters that were previously inaccessible with existing ground-based surveys. We identify 601 clusters in the Year 1 sample, representing more than a factor of four increase over previous catalogs within the current survey area (390 arcmin^2). This work presents results derived from the first \sim25% of the survey data; we estimate that the final sample will include \sim2500 clusters. For the Year 1 objects, we present a catalog with positions, radii, and six-band integrated photometry. Along with a general characterization of the cluster luminosities and colors, we discuss the cluster luminosity function, the cluster size distributions, and highlight a number of individually interesting clusters found in the Year 1 search.

We present thirteen epochs of near-infrared (0.8-5 micron) spectroscopic observations of the pre-transitional, "gapped" disk system in SAO 206462 (=HD 135344B). In all, six gas emission lines (including Br gamma, Pa beta, and the 0.8446 micron line of O I) along with continuum measurements made near the standard J, H, K, and L photometric bands were measured. A mass accretion rate of approximately 2 x 10^-8 solar masses per year was derived from the Br gamma and Pa beta lines. However, the fluxes of these lines varied by a factor of over two during the course of a few months. The continuum also varied, but by only ~30%, and even decreased at a time when the gas emission was increasing. The H I line at 1.083 microns was also found to vary in a manner inconsistent with that of either the hydrogen lines or the dust. Both the gas and dust variabilities indicate significant changes in the region of the inner gas and the inner dust belt that may be common to many young disk systems. If planets are responsible for defining the inner edge of the gap, they could interact with the material on time scales commensurate with what is observed for the variations in the dust, while other disk instabilities (thermal, magnetorotational) would operate there on longer time scales than we observe for the inner dust belt. For SAO 206462, the orbital period would likely be 1-3 years. If the changes are being induced in the disk material closer to the star than the gap, a variety of mechanisms (disk instabilities, interactions via planets) might be responsible for the changes seen. The He I feature is most likely due to a wind whose orientation changes with respect to the observer on time scales of a day or less. To further constrain the origin of the gas and dust emission will require multiple spectroscopic and interferometric observations on both shorter and longer time scales that have been sampled so far.