We develop a data-driven model to map stellar parameters (effective temperature, surface gravity and metallicity) accurately and precisely to broad-band stellar photometry. This model must, and does, simultaneously constrain the passband-specific dust reddening vector in the Milky Way. The model uses a neural network to learn the (de-reddened) absolute magnitude in one band and colors across many bands, given stellar parameters from spectroscopic surveys and parallax constraints from Gaia. To demonstrate the effectiveness of this approach, we train our model on a dataset with spectroscopic labels from LAMOST, APOGEE and GALAH, Gaia parallaxes, and optical and near-infrared photometry from Gaia, Pan-STARRS 1, 2MASS and WISE. Testing the model on these datasets leads to an excellent fit and a precise -- and by construction accurate -- prediction of the color-magnitude diagrams in many bands. This flexible approach rigorously links spectroscopic and photometric surveys, and also results in an improved, temperature-dependent reddening vector. As such, it provides a simple and accurate method for predicting photometry in stellar evolutionary models. Our model will form a basis to infer stellar properties, distances and dust extinction from photometric data, which should be of great use in 3D mapping of the Milky Way. Our trained model may be obtained at https://doi.org/10.5281/zenodo.3902382.

In a cosmological setting, the disc of a galaxy is expected to continuously experience gravitational torques and perturbations from a variety of sources, which can cause the disc to wobble, flare and warp. Specifically, the study of galactic warps and their dynamical nature can potentially reveal key information on the formation history of galaxies and the mass distribution of their halos. Our Milky Way presents a unique case study for galactic warps, thanks to the detailed knowledge of its stellar distribution and kinematics. Using a simple model of how the warp's orientation is changing with time, we here measure the precession rate of the Milky Way's warp using 12 million giant stars from Gaia Data Release 2, finding that it is precessing at $10.86 \pm 0.03_{stat} \pm 3.20_{sys}$ km/s/kpc in the direction of Galactic rotation, about one third the angular velocity at the Sun's position in the Galaxy. The direction and magnitude of the warp's precession rate favour the scenario that the warp is the results of a recent or ongoing encounter with a satellite galaxy, rather than the relic of the ancient assembly history of the Galaxy.

We present a mock stellar catalog, matching in volume, depth and data model the content of the planned Gaia early data release 3 (Gaia EDR3). We have generated our catalog (GeDR3mock) using galaxia, a tool to sample stars from an underlying Milky Way (MW) model or from N-body data. We used an updated Besan\c{c}on Galactic model together with the latest PARSEC stellar evolutionary tracks, now also including white dwarfs. We added the Magellanic clouds and realistic open clusters with internal rotation. We empirically modelled uncertainties based on Gaia DR2 (GDR2) and scaled them according to the longer baseline in Gaia EDR3. The apparent magnitudes were reddened according to a new selection of 3D extinction maps. To help with the Gaia selection function we provide all-sky magnitude limit maps in G and BP for a few relevant GDR2 subsets together with the routines to produce these maps for user-defined subsets. We supplement the catalog with photometry and extinctions in non-Gaia bands. The catalog is available in the Virtual Observatory and can be queried just like the actual Gaia EDR3 will be. We highlight a few capabilities of the Astronomy Data Query Language (ADQL) with educative catalog queries. We use the data extracted from those queries to compare GeDR3mock to GDR2, which emphasises the importance of adding observational noise to the mock data. Since the underlying truth, e.g. stellar parameters, is know in GeDR3mock, it can be used to construct priors as well as mock data tests for parameter estimation. All code, models and data used to produce GeDR3mock are linked and contained in galaxia_wrap, a python package, representing a fast galactic forward model, able to project MW models and N-body data into realistic Gaia observables.

The extended main sequence turn offs (eMSTOs) of several young to intermediate age clusters are examined in the Magellanic Clouds and the Milky Way. We explore the effects of extended star formation (eSF) and a range of stellar rotation rates on the behavior of the color-magnitude diagram (CMD), paying particular attention to the MSTO. We create synthetic stellar populations based on MESA stellar models to simulate observed Hubble Space Telescope and Gaia star cluster data. We model the effect of rotation as a non-parametric distribution, allowing for maximum flexibility. In our models the slow rotators comprise the blueward, and fast rotators the redward portion of the eMSTO. We simulate data under three scenarios: non-rotating eSF, a range of rotation rates with a single age, and a combination of age and rotation effects. We find that two of the five clusters (the youngest and oldest) favor an age spread, but these also achieve the overall worst fits of all clusters. The other three clusters show comparable statistical evidence between rotation and an age spread. In all five cases, a rotation rate distribution alone is capable of qualitatively matching the observed eMSTO structure. In future work, we aim to compare our predicted Vsin(i) with observations in order to better constrain the physics related to stellar rotation.

We construct a supervised classifier based on Gaussian Mixture Models to probabilistically classify objects in Gaia data release 2 (GDR2) using only photometric and astrometric data in that release. The model is trained empirically to classify objects into three classes -- star, quasar, galaxy -- for G<=14.5 mag down to the Gaia magnitude limit of G=21.0 mag. Galaxies and quasars are identified for the training set by a cross-match to objects with spectroscopic classifications from the Sloan Digital Sky Survey. Stars are defined directly from GDR2. When allowing for the expectation that quasars are 500 times rarer than stars, and galaxies 7500 times rarer than stars (the class imbalance problem), samples classified with a threshold probability of 0.5 are predicted to have purities of 0.43 for quasars and 0.28 for galaxies, and completenesses of 0.58 and 0.72 respectively. The purities can be increased up to 0.60 by adopting a higher threshold. Not accounting for this expected low frequency of extragalactic objects (the class prior) would give both erroneously optimistic performance predictions and severely impure samples. Applying our model to all 1.20 billion objects in GDR2 with the required features, we classify 2.3 million objects as quasars and 0.37 million objects as galaxies (with individual probabilities above 0.5). The small number of galaxies is due to the strong bias of the satellite detection algorithm and on-ground data selection against extended objects. We infer the true number of quasars and galaxies -- as these classes are defined by our training set -- to be 690,000 and 110,000 respectively (+/- 50%). The aim of this work is to see how well extragalactic objects can be classified using only GDR2 data. Better classifications should be possible with the low resolution spectroscopy (BP/RP) planned for GDR3.

The mechanisms of the formation and evolution of the Milky Way are encoded in the orbits, chemistry and ages of its stars. With the 4MOST MIlky way Disk And BuLgE Low-Resolution Survey (4MIDABLE-LR) we aim to study kinematic and chemical substructures in the Milky Way disc and bulge region with samples of unprecedented size out to larger distances and greater precision than conceivable with Gaia alone or any other ongoing or planned survey. Gaia gives us the unique opportunity for target selection based almost entirely on parallax and magnitude range, hence increasing the efficiency in sampling larger Milky Way volumes with well-defined and effective selection functions. Our main goal is to provide a detailed chrono-chemo-kinematical extended map of our Galaxy and the largest Gaia follow-up down to $G = 19$ magnitudes (Vega). The complex nature of the disc components (for example, large target densities and highly structured extinction distribution in the Milky Way bulge and disc area), prompted us to develop a survey strategy with five main sub-surveys that are tailored to answer the still open questions about the assembly and evolution of our Galaxy, while taking full advantage of the Gaia data.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

The 4MOST instrument is a multi-object spectrograph that will address Galactic and extragalactic science cases simultaneously by observing targets from a large number of different surveys within each science exposure. This parallel mode of operation and the survey nature of 4MOST require some distinct 4MOST-specific operational features within the overall operations model of ESO. The main feature is that the 4MOST Consortium will deliver, not only the instrument, but also contractual services to the user community, which is why 4MOST is also described as a facility. This white paper concentrates on information particularly useful to answering the forthcoming Call for Letters of Intent.

We use Gaia DR2 astrometric and photometric data, published radial velocities and MESA models to infer distances, orbits, surface gravities, and effective temperatures for all ultra metal-poor stars ($\FeH<-4.0$ dex) available in the literature. Assuming that these stars are old ($>11\Gyr$) and that they are expected to belong to the Milky Way halo, we find that these 42 stars (18 dwarf stars and 24 giants or sub-giants) are currently within $\sim20\kpc$ of the Sun and that they map a wide variety of orbits. A large fraction of those stars remains confined to the inner parts of the halo and was likely formed or accreted early on in the history of the Milky Way, while others have larger apocentres ($>30\kpc$), hinting at later accretion from dwarf galaxies. Of particular interest, we find evidence that a significant fraction of all known UMP stars ($\sim26$\%) are on prograde orbits confined within $3\kpc$ of the Milky Way plane ($J_z < 100 \kms \kpc$). One intriguing interpretation is that these stars belonged to the massive building block(s) of the proto-Milky Way that formed the backbone of the Milky Way disc. Alternatively, they might have formed in the early disc and have been dynamically heated, or have been brought into the Milky Way by one or more accretion events whose orbit was dragged into the plane by dynamical friction before disruption. The combination of the exquisite Gaia DR2 data and surveys of the very metal-poor sky opens an exciting era in which we can trace the very early formation of the Milky Way.

Using Gaia DR2 astrometry, we map the kinematic signature of the Galactic stellar warp out to a distance of 7 kpc from the Sun. Combining Gaia DR2 and 2MASS photometry, we identify, via a probabilistic approach, 599 494 upper main sequence stars and 12 616 068 giants without the need for individual extinction estimates. The spatial distribution of the upper main sequence stars clearly shows segments of the nearest spiral arms. The large-scale kinematics of both the upper main sequence and giant populations show a clear signature of the warp of the Milky Way, apparent as a gradient of 5-6 km/s in the vertical velocities from 8 to 14 kpc in Galactic radius. The presence of the signal in both samples, which have different typical ages, suggests that the warp is a gravitationally induced phenomenon.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

For the vast majority of stars in the second Gaia data release, reliable distances cannot be obtained by inverting the parallax. A correct inference procedure must instead be used to account for the nonlinearity of the transformation and the asymmetry of the resulting probability distribution. Here we infer distances to essentially all 1.33 billion stars with parallaxes published in the second \gaia\ data release. This is done using a weak distance prior that varies smoothly as a function of Galactic longitude and latitude according to a Galaxy model. The irreducible uncertainty in the distance estimate is characterized by the lower and upper bounds of an asymmetric confidence interval. Although more precise distances can be estimated for a subset of the stars using additional data (such as photometry), our goal is to provide purely geometric distance estimates, independent of assumptions about the physical properties of, or interstellar extinction towards, individual stars. We analyse the characteristics of the catalogue and validate it using clusters. The catalogue can be queried on the Gaia archive using ADQL at http://gea.esac.esa.int/archive/ and downloaded from http://www.mpia.de/~calj/gdr2_distances.html .

The early Universe presented a star formation environment that was almost devoid of heavy elements. The lowest metallicity stars thus provide a unique window into the earliest Galactic stages, but are exceedingly rare and difficult to find. Here we present the discovery of an ultra-metal-poor star, Pristine_221.8781+9.7844, using narrow-band Ca H&K photometry from the Pristine survey. Follow-up medium and high-resolution spectroscopy confirms the ultra-metal-poor nature of Pristine_221.8781+9.7844 ([Fe/H] = -4.66 +/- 0.13 in 1D LTE) with an enhancement of 0.3-0.4 dex in alpha-elements relative to Fe, and an unusually low carbon abundance. We derive an upper limit of A(C) = 5.6, well below typical A(C) values for such ultra metal-poor stars. This makes Pristine_221.8781+9.7844 one of the most metal-poor stars; in fact, it is very similar to the most metal-poor star known (SDSS J102915+172927). The existence of a class of ultra metal-poor stars with low(er) carbon abundances suggest that there must have been several formation channels in the early Universe through which long-lived, low-mass stars were formed.

Passing stars may play an important role in the evolution of our solar system. We search for close stellar encounters to the Sun among all 7.2 million stars in Gaia-DR2 that have six-dimensional phase space data. We characterize encounters by integrating their orbits through a Galactic potential and propagating the correlated uncertainties via a Monte Carlo resampling. After filtering to remove spurious data, we find 694 stars that have median (over uncertainties) closest encounter distances within 5 pc, all occurring within 15 Myr from now. 26 of these have at least a 50% chance of coming closer than 1 pc (and 7 within 0.5 pc), all but one of which are newly discovered here. We further confirm some and refute several other previously-identified encounters, confirming suspicions about their data. The closest encounter in the sample is Gl 710, which has a 95% probability of coming closer than 0.08 pc (17 000 AU). Taking mass estimates from Gaia astrometry and multiband photometry for essentially all encounters, we find that Gl 710 also has the largest impulse on the Oort cloud. Using a Galaxy model, we compute the completeness of the Gaia-DR2 encountering sample as a function of perihelion time and distance. Only 15% of encounters within 5 pc occurring within +/- 5 Myr of now have been identified, mostly due to the lack of radial velocities for faint and/or cool stars. Accounting for the incompleteness, we infer the present rate of encounters within 1 pc to be 19.7 +/- 2.2 per Myr, a quantity expected to scale quadratically with the encounter distance out to at least several pc. Spuriously large parallaxes in our sample from imperfect filtering would tend to inflate both the number of encounters found and this inferred rate. The magnitude of this effect is hard to quantify.

The second Gaia data release (Gaia-DR2) contains, beyond the astrometry, three-band photometry for 1.38 billion sources. We have used these three broad bands to infer stellar effective temperatures, Teff, for all sources brighter than G=17 mag with Teff in the range 3000-10 000 K (161 million sources). Using in addition the parallaxes, we infer the line-of-sight extinction, A_G, and the reddening, E[BP-RP], for 88 million sources. Together with a bolometric correction we derive luminosity and radius for 77 million sources. These quantities as well as their estimated uncertainties are part of Gaia-DR2. Here we describe the procedures by which these quantities were obtained, including the underlying assumptions, comparison with literature estimates, and the limitations of our results. Typical accuracies are of order 324 K (Teff), 0.46 mag (A_G), 0.23 mag (E[BP-RP]), 15% (luminosity), and 10% (radius). Being based on only a small number of observable quantities and limited training data, our results are necessarily subject to some extreme assumptions that can lead to strong systematics in some cases (not included in the aforementioned accuracy estimates). One aspect is the non-negativity contraint of our estimates, in particular extinction. Yet in several regions of parameter space our results show very good performance, for example for red clump stars and solar analogues. Large uncertainties render the extinctions less useful at the individual star level, but they show good performance for ensemble estimates. We identify regimes in which our parameters should and should not be used and we define a "clean" sample. Despite the limitations, this is the largest catalogue of uniformly-inferred stellar parameters to date. More precise and detailed astrophysical parameters based on the full BP/RP spectrophotometry are planned as part of the third Gaia data release.

We present a mock catalog of Milky Way stars, matching in volume and depth the content of the Gaia data release 2 (GDR2). We generated our catalog using Galaxia, a tool to sample stars from a Besancon Galactic model, together with a realistic 3D dust extinction map. The catalog mimicks the complete GDR2 data model and contains most of the entries in the Gaia source catalog: 5-parameter astrometry, 3-band photometry, radial velocities, stellar parameters, and associated scaled nominal uncertainty estimates. In addition, we supplemented the catalog with extinctions and photometry for non-Gaia bands. This catalog can be used to prepare GDR2 queries in a realistic runtime environment, and it can serve as a Galactic model against which to compare the actual GDR2 data in the space of observables. The catalog is hosted through the virtual observatory GAVO's Heidelberg data center service and thus can be queried using ADQL as for GDR2 data.

We present a map of the three-dimensional (3D) distribution of dust in the Orion complex. Orion is the closest site of high-mass star formation, making it an excellent laboratory for studying the interstellar medium and star formation. We use data from the Gaia-TGAS catalogue combined with photometry from 2MASS and WISE to get the distances and extinctions of individual stars in the vicinity of the Orion complex. We adopt the non-parametric method of Rezaei Kh. et al. (2017) to infer the probability distribution function of the dust densities at arbitrary points throughout the region. We map the dust distribution towards different parts of the Orion complex and find the distance and depth of the cloud compatible with other recent works which show the ability of the method to be applicable on the local molecular clouds to map their 3D dust distribution. We also demonstrate the danger of only using colours of stars to derive their extinctions without considering further physical constraints like the colour-magnitude diagram (CMD).

Observational tests of stellar and Galactic chemical evolution call for the joint knowledge of a star's physical parameters, detailed element abundances, and precise age. For cool main-sequence (MS) stars the abundances of many elements can be measured from spectroscopy, but ages are very hard to determine. The situation is different if the MS star has a white dwarf (WD) companion and a known distance, as the age of such a binary system can then be determined precisely from the photometric properties of the cooling WD. As a pilot study for obtaining precise age determinations of field MS stars, we identify nearly one hundred candidates for such wide binary systems: a faint WD whose GPS1 proper motion matches that of a brighter MS star in Gaia/TGAS with a good parallax ($\sigma_\varpi/\varpi\le 0.05$). We model the WD's multi-band photometry with the BASE-9 code using this precise distance (assumed to be common for the pair) and infer ages for each binary system. The resulting age estimates are precise to $\le 10\%$ ($\le 20\%$) for $42$ ($67$) MS-WD systems. Our analysis more than doubles the number of MS-WD systems with precise distances known to date, and it boosts the number of such systems with precise age determination by an order of magnitude. With the advent of the Gaia DR2 data, this approach will be applicable to a far larger sample, providing ages for many MS stars (that can yield detailed abundances for over 20 elements), especially in the age range 2 to 8\,\Gyr, where there are only few known star clusters.

The Pristine survey is a narrow-band, photometric survey focused around the wavelength region of the Ca II H & K absorption lines, designed to efficiently search for extremely metal-poor stars. In this work, we use the first results of a medium-resolution spectroscopic follow-up to refine the selection criteria for finding extremely metal-poor stars ($\textrm{[Fe/H]} \leq -3.0$) in the Pristine survey. We consider methods by which stars can be selected from available broad-band and infrared photometry plus the additional Pristine narrow-band photometry. The spectroscopic sample presented in this paper consists of 205 stars in the magnitude range $14 < V < 18$. Applying the photometric selection criteria cuts the sample down to 149 stars, and from these we report a success rate of 70% for finding stars with $\textrm{[Fe/H]} \leq -2.5$ and 22% for finding stars with $\textrm{[Fe/H]} \leq -3.0$. These statistics compare favourably with other surveys that search for extremely metal-poor stars, namely an improvement by a factor of $\sim 4-5$ for recovering stars with $\textrm{[Fe/H]} \leq -3.0$. In addition, Pristine covers a fainter magnitude range than its predecessors, and can thus probe deeper into the Galactic halo.

We present the Pristine survey, a new narrow-band photometric survey focused on the metallicity-sensitive Ca H & K lines and conducted in the northern hemisphere with the wide-field imager MegaCam on the Canada-France-Hawaii Telescope (CFHT). This paper reviews our overall survey strategy and discusses the data processing and metallicity calibration. Additionally we review the application of these data to the main aims of the survey, which are to gather a large sample of the most metal-poor stars in the Galaxy, to further characterise the faintest Milky Way satellites, and to map the (metal-poor) substructure in the Galactic halo. The current Pristine footprint comprises over 1,000 deg2 in the Galactic halo ranging from b~30 to 78 and covers many known stellar substructures. We demonstrate that, for SDSS stellar objects, we can calibrate the photometry at the 0.02-magnitude level. The comparison with existing spectroscopic metallicities from SDSS/SEGUE and LAMOST shows that, when combined with SDSS broad-band g and i photometry, we can use the CaHK photometry to infer photometric metallicities with an accuracy of ~0.2 dex from [Fe/H]=-0.5 down to the extremely metal-poor regime ([Fe/H]<-3.0). After the removal of various contaminants, we can efficiently select metal-poor stars and build a very complete sample with high purity. The success rate of uncovering [Fe/H]SEGUE<-3.0 stars among [Fe/H]Pristine<-3.0 selected stars is 24% and 85% of the remaining candidates are still very metal poor ([Fe/H]<-2.0). We further demonstrate that Pristine is well suited to identify the very rare and pristine Galactic stars with [Fe/H]<-4.0, which can teach us valuable lessons about the early Universe.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

We present a framework to simultaneously constrain the values and uncertainties of the strength of convective core overshooting, metallicity, extinction, distance, and age in stellar populations. We then apply the framework to archival Hubble Space Telescope observations of six stellar clusters in the Large Magellanic Cloud that have reported ages between ~1-2.5 Gyr. Assuming a canonical value of the strength of core convective overshooting, we recover the well-known age-metallicity correlation, and additional correlations between metallicity and extinction and metallicity and distance. If we allow the strength of core overshooting to vary, we find that for intermediate-aged stellar clusters, the measured values of distance and extinction are negligibly effected by uncertainties of core overshooting strength. However, cluster age and metallicity may have disconcertingly large systematic shifts when core overshooting strength is allowed to vary by more than +/- 0.05 Hp. Using the six stellar clusters, we combine their posterior distribution functions to obtain the most probable core overshooting value, 0.500 +0.016 -0.134 Hp, which is in line with canonical values.

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.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

Pulsating stars, such as Cepheids, Miras, and RR Lyrae stars, are important distance indicators and calibrators of the "cosmic distance ladder", and yet their period-luminosity-metallicity (PLZ) relations are still constrained using simple statistical methods that cannot take full advantage of available data. To enable optimal usage of data provided by the Gaia mission, we present a probabilistic approach that simultaneously constrains parameters of PLZ relations and uncertainties in Gaia parallax measurements. We demonstrate this approach by constraining PLZ relations of type $ab$ RR Lyrae stars in near-infrared W1 and W2 bands, using Tycho-Gaia Astrometric Solution (TGAS) parallax measurements for a sample of $\approx100$ type $ab$ RR Lyrae stars located within 2.5 kpc of the Sun. The fitted PLZ relations are consistent with previous studies, and in combination with other data, deliver distances precise to 6% (once various sources of uncertainty are taken into account). To a precision of 0.05 mas ($1\sigma$), we do not find a statistically significant offset in TGAS parallaxes for this sample of distant RR Lyrae stars (median parallax of 0.8 mas and distance of 1.4 kpc). With only minor modifications, our probabilistic approach can be used to constrain PLZ relations of other pulsating stars, and we intend to apply it to Cepheid and Mira stars in the near future.

We explore to which extent stars within Galactic disk open clusters resemble each other in the high-dimensional space of their photospheric element abundances, and contrast this with pairs of field stars. Our analysis is based on abundances for 20 elements, homogeneously derived from APOGEE spectra (with carefully quantified uncertainties, with a median value of $\sim 0.03$ dex). We consider 90 red giant stars in seven open clusters and find that most stars within a cluster have abundances in most elements that are indistinguishable (in a $\chi^2$-sense) from those of the other members, as expected for stellar birth siblings. An analogous analysis among pairs of $>1000$ field stars shows that highly significant abundance differences in the 20-dimensional space can be established for the vast majority of these pairs, and that the APOGEE-based abundance measurements have high discriminating power. However, pairs of field stars whose abundances are indistinguishable even at 0.03~dex precision exist: $\sim 0.3$ percent of all field star pairs, and $\sim 1.0$ percent of field star pairs at the same (solar) metallicity [Fe/H] = $0 \pm 0.02$. Most of these pairs are presumably not birth siblings from the same cluster, but rather doppelganger. Our analysis implies that 'chemical tagging' in the strict sense, identifying birth siblings for typical disk stars through their abundance similarity alone, will not work with such data. However, our approach shows that abundances have extremely valuable information for probabilistic chemo-orbital modeling and combined with velocities, we have identified new cluster members from the field.

We present synthetic far- and near-ultraviolet (FUV and NUV) maps of M31, both with and without dust reddening. These maps were constructed from spatially-resolved star formation histories (SFHs) derived from optical Hubble Space Telescope imaging of resolved stars, taken as part of the Panchromatic Hubble Andromeda Treasury program. We use stellar population synthesis modeling to generate synthetic UV maps with projected spatial resolution of $\sim$100 pc ($\sim$24 arcseconds) The predicted UV flux agrees well with the observed flux, with median ratios between the modeled and observed flux of $\log_{10}(f^{syn}/f^{obs}) = 0.03\pm0.24$ and $-0.03\pm0.16$ in the FUV and NUV, respectively. This agreement is particularly impressive given that we used only optical photometry to construct these UV maps. We use the dust-free maps to examine properties of obscured flux and star formation by comparing our reddened and dust-free FUV flux maps with the observed FUV and FUV+24{\mu}m flux to examine the fraction of obscured flux. The synthetic flux maps require that $\sim$90% of the FUV flux in M31 is obscured by dust, while the GALEX-based methods suggest that $\sim$70% of the flux is obscured. This increase in the obscured flux estimate is driven by significant differences between the dust-free synthetic FUV flux and that derived when correcting the observed FUV for dust with 24{\mu}m observations. The difference is further illustrated when we compare the SFRs derived from the FUV+24{\mu}m flux with the 100 Myr average SFR from the SFHs. The 24{\mu}m-corrected FUV flux underestimates the SFR by a factor of $\sim$2.3 - 2.5. [abridged]

We present a non-parametric model for inferring the three-dimensional (3D) distribution of dust density in the Milky Way. Our approach uses the extinction measured towards stars at different locations in the Galaxy at approximately known distances. Each extinction measurement is proportional to the integrated dust density along its line-of-sight. Making simple assumptions about the spatial correlation of the dust density, we can infer the most probable 3D distribution of dust across the entire observed region, including along sight lines which were not observed. This is possible because our model employs a Gaussian Process to connect all lines-of-sight. We demonstrate the capability of our model to capture detailed dust density variations using mock data as well as simulated data from the Gaia Universe Model Snapshot. We then apply our method to a sample of giant stars observed by APOGEE and Kepler to construct a 3D dust map over a small region of the Galaxy. Due to our smoothness constraint and its isotropy, we provide one of the first maps which does not show the "fingers of god" effect.

In the Milky Way, the thick disk can be defined using individual stellar abundances, kinematics, or age; or geometrically, as stars high above the mid-plane. In nearby galaxies, where only a geometric definition can be used, thick disks appear to have large radial scale-lengths, and their red colors suggest that they are uniformly old. The Milky Way's geometrically thick disk is also radially extended, but it is far from chemically uniform: alpha-enhanced stars are confined within the inner Galaxy. In simulated galaxies, where old stars are centrally concentrated, geometrically thick disks are radially extended, too. Younger stellar populations flare in the simulated disks' outer regions, bringing those stars high above the mid-plane. The resulting geometrically thick disks therefore show a radial age gradient, from old in their central regions to younger in their outskirts. Based on our age estimates for a large sample of giant stars in the APOGEE survey, we can now test this scenario for the Milky Way. We find that the geometrically-defined thick disk in the Milky Way has indeed a strong radial age gradient: the median age for red clump stars goes from ~9 Gyr in the inner disk to 5 Gyr in the outer disk. We propose that at least some nearby galaxies could also have thick disks that are not uniformly old, and that geometrically thick disks might be complex structures resulting from different formation mechanisms in their inner and outer parts.

We present the Bayesian Extinction And Stellar Tool (BEAST), a probabilistic approach to modeling the dust extinguished photometric spectral energy distribution of an individual star while accounting for observational uncertainties common to large resolved star surveys. Given a set of photometric measurements and an observational uncertainty model, the BEAST infers the physical properties of the stellar source using stellar evolution and atmosphere models and constrains the line of sight extinction using a newly developed mixture model that encompasses the full range of dust extinction curves seen in the Local Group. The BEAST is specifically formulated for use with large multi-band surveys of resolved stellar populations. Our approach accounts for measurement uncertainties and any covariance between them due to stellar crowding (both systematic biases and uncertainties in the bias) and absolute flux calibration, thereby incorporating the full information content of the measurement. We illustrate the accuracy and precision possible with the BEAST using data from the Panchromatic Hubble Andromeda Treasury. While the BEAST has been developed for this survey, it can be easily applied to similar existing and planned resolved star surveys.

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 show that the masses of red giant stars can be well predicted from their photospheric carbon and nitrogen abundances, in conjunction with their spectroscopic stellar labels log g, Teff, and [Fe/H]. This is qualitatively expected from mass-dependent post main sequence evolution. We here establish an empirical relation between these quantities by drawing on 1,475 red giants with asteroseismic mass estimates from Kepler that also have spectroscopic labels from APOGEE DR12. We assess the accuracy of our model, and find that it predicts stellar masses with fractional r.m.s. errors of about 14% (typically 0.2 Msun). From these masses, we derive ages with r.m.s errors of 40%. This empirical model allows us for the first time to make age determinations (in the range 1-13 Gyr) for vast numbers of giant stars across the Galaxy. We apply our model to 52,000 stars in APOGEE DR12, for which no direct mass and age information was previously available. We find that these estimates highlight the vertical age structure of the Milky Way disk, and that the relation of age with [alpha/M] and metallicity is broadly consistent with established expectations based on detailed studies of the solar neighbourhood.

We present a study of spatial variations in the metallicity of old red giant branch stars in the Andromeda galaxy. Photometric metallicity estimates are derived by interpolating isochrones for over seven million stars in the Panchromatic Hubble Andromeda Treasury (PHAT) survey. This is the first systematic study of stellar metallicities over the inner 20 kpc of Andromeda's galactic disk. We see a clear metallicity gradient of $-0.020\pm0.004$ dex/kpc from $\sim4-20$ kpc assuming a constant RGB age. This metallicity gradient is derived after correcting for the effects of photometric bias and completeness and dust extinction and is quite insensitive to these effects. The unknown age gradient in M31's disk creates the dominant systematic uncertainty in our derived metallicity gradient. However, spectroscopic analyses of galaxies similar to M31 show that they typically have small age gradients that make this systematic error comparable to the 1$\sigma$ error on our metallicity gradient measurement. In addition to the metallicity gradient, we observe an asymmetric local enhancement in metallicity at radii of 3-6 kpc that appears to be associated with Andromeda's elongated bar. This same region also appears to have an enhanced stellar density and velocity dispersion.

We map the distribution of dust in M31 at 25pc resolution, using stellar photometry from the Panchromatic Hubble Andromeda Treasury. We develop a new mapping technique that models the NIR color-magnitude diagram (CMD) of red giant branch (RGB) stars. The model CMDs combine an unreddened foreground of RGB stars with a reddened background population viewed through a log-normal column density distribution of dust. Fits to the model constrain the median extinction, the width of the extinction distribution, and the fraction of reddened stars. The resulting extinction map has >4 times better resolution than maps of dust emission, while providing a more direct measurement of the dust column. There is superb morphological agreement between the new map and maps of the extinction inferred from dust emission by Draine et al. 2014. However, the widely-used Draine & Li (2007) dust models overpredict the observed extinction by a factor of ~2.5, suggesting that M31's true dust mass is lower and that dust grains are significantly more emissive than assumed in Draine et al. (2014). The discrepancy we identify is consistent with similar findings in the Milky Way by the Planck Collaboration (2015), but has a more complex dependence on parameters from the Draine & Li (2007) dust models. We also show that the discrepancy with the Draine et al. (2014) map is lowest where the interstellar radiation field has a harder spectrum than average. We discuss possible improvements to the CMD dust mapping technique, and explore further applications.

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.

The Ophiuchus stream is a recently discovered stellar tidal stream in the Milky Way. We present high-quality spectroscopic data for 14 stream member stars obtained using the Keck and MMT telescopes. We confirm the stream as a fast moving ($v_{los}\sim290$ km s$^{-1}$), kinematically cold group ($\sigma_{v_{los}}\lesssim1$ km s$^{-1}$) of $\alpha$-enhanced and metal-poor stars (${\rm [\alpha/Fe]\sim0.4}$ dex, ${\rm [Fe/H]\sim-2.0}$ dex). Using a probabilistic technique, we model the stream simultaneously in line-of-sight velocity, color-magnitude, coordinate, and proper motion space, and so determine its distribution in 6D phase-space. We find that that the stream extends in distance from 7.5 to 9 kpc from the Sun; it is 50 times longer than wide, merely appearing highly foreshortened in projection. The analysis of the stellar population contained in the stream suggests that it is $\sim12$ Gyr old, and that its initial stellar mass was $\sim2\times10^4$ $M_{\odot}$ (or at least $\gtrsim7\times10^3$ $M_{\odot}$). Assuming a fiducial Milky Way potential, we fit an orbit to the stream which matches the observed phase-space distribution, except for some tension in the proper motions: the stream has an orbital period of $\sim350$ Myr, and is on a fairly eccentric orbit ($e\sim0.66$) with a pericenter of $\sim3.5$ kpc and an apocenter of $\sim17$ kpc. The phase-space structure and stellar population of the stream show that its progenitor must have been a globular cluster that was disrupted only $\sim240$ Myr ago. We do not detect any significant overdensity of stars along the stream that would indicate the presence of a progenitor, and conclude that the stream is all that is left of the progenitor.

Stationary density waves rotating at a constant pattern speed $\Omega_{\rm P}$ would produce age gradients across spiral arms. We test whether such age gradients are present in M81 by deriving the recent star formation histories (SFHs) of 20 regions around one of M81's grand-design spiral arms. For each region, we use resolved stellar populations to determine the SFH by modeling the observed color-magnitude diagram (CMD) constructed from archival Hubble Space Telescope (HST) F435W and F606W imaging. Although we should be able to detect systematic time delays in our spatially-resolved SFHs, we find no evidence of star formation propagation across the spiral arm. Our data therefore provide no convincing evidence for a stationary density wave with a single pattern speed in M81, and instead favor the scenario of kinematic spiral patterns that are likely driven by tidal interactions with the companion galaxies M82 and NGC 3077.

We measure the recent star formation history (SFH) across M31 using optical images taken with the \texit{Hubble Space Telescope} as part of the Panchromatic Hubble Andromeda Treasury (PHAT). We fit the color-magnitude diagrams in ~9000 regions that are ~100 pc $\times$ 100 pc in projected size, covering a 0.5 square degree area (~380 kpc$^2$, deprojected) in the NE quadrant of M31. We show that the SFHs vary significantly on these small spatial scales but that there are also coherent galaxy-wide fluctuations in the SFH back to ~500 Myr, most notably in M31's 10-kpc star-forming ring. We find that the 10-kpc ring is at least 400 Myr old, showing ongoing star formation over the past ~500 Myr. This indicates the presence of molecular gas in the ring over at least 2 dynamical times at this radius. We also find that the ring's position is constant throughout this time, and is stationary at the level of 1 km/s, although there is evidence for broadening of the ring due to diffusion of stars into the disk. Based on existing models of M31's ring features, the lack of evolution in the ring's position makes a purely collisional ring origin highly unlikely. We find that the global SFR has been fairly constant over the last ~500 Myr, though it does show a small increase at 50 Myr that is 1.3 times the average SFR over the past 100 Myr. During the last ~500 Myr, ~60% of all SF occurs in the 10-kpc ring. Finally, we find that in the past 100 Myr, the average SFR over the PHAT survey area is $0.28\pm0.03$ M$_\odot$ yr$^{-1}$ with an average deprojected intensity of $7.3 \times 10^{-4}$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$, which yields a total SFR of ~0.7 M$_\odot$ yr$^{-1}$ when extrapolated to the entire area of M31's disk. This SFR is consistent with measurements from broadband estimates. [abridged]

We have identified a major global enhancement of star formation in the inner M31 disk that occurred between 2-4 Gyr ago, producing $\sim$60% of the stellar mass formed in the past 5 Gyr. The presence of this episode in the inner disk was discovered by modeling the optical resolved star color-magnitude diagrams of low extinction regions in the main disk of M31 (3$<$R$<$20 kpc) as part of the Panchromatic Hubble Andromeda Treasury. This measurement confirms and extends recent measurements of a widespread star formation enhancement of similar age in the outer disk, suggesting that this burst was both massive and global. Following the galaxy-wide burst, the star formation rate of M31 has significantly declined. We briefly discuss possible causes for these features of the M31 evolutionary history, including interactions with M32, M33 and/or a merger.

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.

The [CII] 158 micron line is one of the strongest emission lines observed in star-forming galaxies, and has been empirically measured to correlate with the star formation rate (SFR) globally and on ~kpc scales. However, due to the multi-phase origins of [CII], one might expect this relation to break down at small scales. We investigate the origins of [CII] emission by examining high spatial resolution observations of [CII] in M31, with the Survey of Lines in M31 (SLIM). We present five ~700x700 pc (3"x3") Fields mapping the [CII] emission, Halpha emission, combined with ancillary infrared (IR) data. We spatially separate star-forming regions from diffuse gas and dust emission on ~50 pc scales. We find that the [CII] - SFR correlation holds even at these scales, although the relation typically has a flatter slope than found at larger (~kpc) scales. While the Halpha emission in M31 is concentrated in the SFR regions, we find that a significant amount (~20-90%) of the [CII] emission comes from outside star-forming regions, and that the total IR (TIR) emission has the highest diffuse fraction of all SFR tracers. We find a weak correlation of the [CII]/TIR to dust color in each Field, and find a large scale trend of increasing [CII]/TIR with galactocentric radius. The differences in the relative diffuse fractions of [CII], Halpha and IR tracers are likely caused by a combination of energetic photon leakage from HII regions and heating by the diffuse radiation field arising from older (B-star) stellar populations. However, we find that by averaging our measurements over ~kpc scales, these effects are minimized, and the relation between [CII] and SFR found in other nearby galaxy studies is retrieved.

Using Hubble Space Telescope (HST) photometry to measure star formation histories, we age-date the stellar populations surrounding supernova remnants (SNRs) in M31 and M33. We then apply stellar evolution models to the ages to infer the corresponding masses for their supernova progenitor stars. We analyze 33 M33 SNR progenitors and 29 M31 SNR progenitors in this work. We then combine these measurements with 53 previously published M31 SNR progenitor measurements to bring our total number of progenitor mass estimates to 115. To quantify the mass distributions, we fit power laws of the form $dN/dM \propto M^{-\alpha}$. Our new, larger sample of M31 progenitors follows a distribution with $\alpha = 4.4\pm 0.4$, and the M33 sample follows a distribution with $\alpha = 3.8^{+0.4}_{-0.5}$. Thus both samples are consistent within the uncertainties, and the full sample across both galaxies gives $\alpha = 4.2\pm 0.3$. Both the individual and full distributions display a paucity of massive stars when compared to a Salpeter initial mass function (IMF), which we would expect to observe if all massive stars exploded as SN that leave behind observable SNR. If we instead fix $\alpha = 2.35$ and treat the maximum mass as a free parameter, we find $M_{max} \sim 35-45M_{sun}$, indicative of a potential maximum cutoff mass for SN production. Our results suggest that either SNR surveys are biased against finding objects in the youngest (<10 Myr old) regions, or the highest mass stars do not produce SNe.

We search Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) and Wide Field Camera 3 (WFC3) broadband imaging data from the Panchromatic Hubble Andromeda Treasury (PHAT) survey to identify detections of cataloged planetary nebulae (PNe). Of the 711 PNe currently in the literature within the PHAT footprint, we find 467 detected in the broadband. For these 467 we are able to refine their astrometric accuracy from ~0."3 to 0."05. Using the resolution of HST, we are able to show that 152 objects currently in the catalogs are definitively not PNe, and we show that 32 objects thought to be extended in ground-based images are actually point-like and therefore good PN candidates. We also find one PN candidate that is marginally resolved. If this is a PN, it is up to 0.7 pc in diameter. With our new photometric data, we develop a method of measuring the level of excitation in individual PNe by comparing broadband and narrowband imaging and describe the effects of excitation on a PN's photometric signature. Using the photometric properties of the known PNe in the PHAT catalogs, we search for more PN, but do not find any new candidates, suggesting that ground-based emission-line surveys are complete in the PHAT footprint to F475W $\simeq$ 24.

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.

We measure the extinction curve in the central 200 pc of M31 at mid-ultraviolet to near-infrared wavelengths (from 1928A to 1.5{\mu}m), using Swift/UVOT and HST WFC3/ACS observations in thirteen bands. Taking advantage of the high angular resolution of the HST WFC3 and ACS detectors, we develop a method to simultaneously determine the relative extinction and the fraction of obscured starlight for five dusty complexes located in the circumnuclear region. The extinction curves of these clumps (RV =2.4-2.5) are steeper than the average Galactic one (RV =3.1), but are similar to optical and near-infrared curves recently measured toward the Galactic Bulge (RV~2.5). This similarity suggests that steep extinction curves may be common in the inner bulge of galaxies. In the ultraviolet, the extinction curves of these clumps are also unusual. We find that one dusty clump (size <2 pc) exhibits a strong UV bump (extinction at 2175A), more than three standard deviation higher than that predicted by common models. Although the high stellar metallicity of the M31 bulge indicates that there are sufficient carbon and silicon to produce large dust grains, the grains may have been destroyed by supernova explosions or past activity of the central super-massive black hole, resulting in the observed steepened extinction curve.

We present a study on the gender balance, in speakers and attendees, at the recent major astronomical conference, the American Astronomical Society meeting 223, in Washington, DC. We conducted an informal survey, yielding over 300 responses by volunteers at the meeting. Each response included gender data about a single talk given at the meeting, recording the gender of the speaker and all question-askers. In total, 225 individual AAS talks were sampled. We analyze basic statistical properties of this sample. We find that the gender ratio of the speakers closely matched the gender ratio of the conference attendees. The audience asked an average of 2.8 questions per talk. Talks given by women had a slightly higher number of questions asked (3.2$\pm$0.2) than talks given by men (2.6$\pm$0.1). The most significant result from this study is that while the gender ratio of speakers very closely mirrors that of conference attendees, women are under-represented in the question-asker category. We interpret this to be an age-effect, as senior scientists may be more likely to ask questions, and are more commonly men. A strong dependence on the gender of session chairs is found, whereby women ask disproportionately fewer questions in sessions chaired by men. While our results point to laudable progress in gender-balanced speaker selection, we believe future surveys of this kind would help ensure that collaboration at such meetings is as inclusive as possible.

We present ages and masses for 601 star clusters in M31 from the analysis of the six filter integrated light measurements from near ultraviolet to near infrared wavelengths, made as part of the Panchromatic Hubble Andromeda Treasury (PHAT). We derive the ages and masses using a probabilistic technique, which accounts for the effects of stochastic sampling of the stellar initial mass function. Tests on synthetic data show that this method, in conjunction with the exquisite sensitivity of the PHAT observations and their broad wavelength baseline, provides robust age and mass recovery for clusters ranging from $\sim 10^2 - 2 \times 10^6 M_\odot$. We find that the cluster age distribution is consistent with being uniform over the past $100$ Myr, which suggests a weak effect of cluster disruption within M31. The age distribution of older ($>100$ Myr) clusters fall towards old ages, consistent with a power-law decline of index $-1$, likely from a combination of fading and disruption of the clusters. We find that the mass distribution of the whole sample can be well-described by a single power-law with a spectral index of $-1.9 \pm 0.1$ over the range of $10^3-3 \times 10^5 M_\odot$. However, if we subdivide the sample by galactocentric radius, we find that the age distributions remain unchanged. However, the mass spectral index varies significantly, showing best fit values between $-2.2$ and $-1.8$, with the shallower slope in the highest star formation intensity regions. We explore the robustness of our study to potential systematics and conclude that the cluster mass function may vary with respect to environment.

In order to empirically determine the timescale and environmental dependence of stellar cluster disruption, we have undertaken an analysis of the unprecedented multi-pointing (seven), multi-wavelength (U, B, V, H$\alpha$, and I) Hubble Space Telescope imaging survey of the nearby, face-on spiral galaxy M83. The images are used to locate stellar clusters and stellar associations throughout the galaxy. Estimation of cluster properties (age, mass, and extinction) was done through a comparison of their spectral energy distributions with simple stellar population models. We constructed the largest catalog of stellar clusters and associations in this galaxy to-date, with ~1800 sources with masses above ~5000 M$_{\odot}$ and ages younger than ~300 Myr. In the present letter, we focus on the age distribution of the resulting clusters and associations. In particular, we explicitly test whether the age distributions are related with the ambient environment. Our results are in excellent agreement with previous studies of age distributions in the centre of the galaxy, which gives us confidence to expand out to search for similarities or differences in the other fields which sample different environments. We find that the age distribution of the clusters inside M83 varies strongly as a function of position within the galaxy, indicating a strong correlation with the galactic environment. If the age distributions are approximated as a power-law of the form $dN/dt \propto t^{\zeta}$, we find $\zeta$ values between 0 and -0.62 ( $\zeta$~-0.40 for the whole galaxy), in good agreement with previous results and theoretical predictions.

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)

Using HST photometry, we age-date 59 supernova remnants (SNRs) in the spiral galaxy M31 and use these ages to estimate zero-age main sequence masses (MZAMS) for their progenitors. To accomplish this, we create color-magnitude diagrams (CMDs) and use CMD fitting to measure the recent star formation history (SFH) of the regions surrounding cataloged SNR sites. We identify any young coeval population that likely produced the progenitor star and assign an age and uncertainty to that population. Application of stellar evolution models allows us to infer the MZAMS from this age. Because our technique is not contingent on precise location of the progenitor star, it can be applied to the location of any known SNR. We identify significant young SF around 53 of the 59 SNRs and assign progenitor masses to these, representing a factor of 2 increase over currently measured progenitor masses. We consider the remaining 6 SNRs as either probable Type Ia candidates or the result of core-collapse progenitors that have escaped their birth sites. The distribution of recovered progenitor masses is bottom heavy, showing a paucity of the most massive stars. If we assume a single power law distribution, dN/dM proportional to M^alpha, we find a distribution that is steeper than a Salpeter IMF (alpha=-2.35). In particular, we find values of alpha outside the range -2.7 to -4.4 inconsistent with our measured distribution at 95% confidence. If instead we assume a distribution that follows a Salpeter IMF up to some maximum mass, we find that values of M_max greater than 26 Msun are inconsistent with the measured distribution at 95% confidence. In either scenario, the data suggest that some fraction of massive stars may not explode. The result is preliminary and requires more SNRs and further analysis. In addition, we use our distribution to estimate a minimum mass for core collapse between 7.0 and 7.8 Msun.

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.

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 study the relationship between the field star formation and cluster formation properties in a large sample of nearby dwarf galaxies. We use optical data from the Hubble Space Telescope and from ground-based telescopes to derive the ages and masses of the young (t_age < 100Myr) cluster sample. Our data provides the first constraints on two proposed relationships between the star formation rate of galaxies and the properties of their cluster systems in the low star formation rate regime. The data show broad agreement with these relationships, but significant galaxy-to-galaxy scatter exists. In part, this scatter can be accounted for by simulating the small number of clusters detected from stochastically sampling the cluster mass function. However, this stochasticity does not fully account for the observed scatter in our data suggesting there may be true variations in the fraction of stars formed in clusters in dwarf galaxies. Comparison of the cluster formation and the brightest cluster in our sample galaxies also provide constraints on cluster destruction models.

This contribution addresses the question of whether the initial cluster mass function (ICMF) has a fundamental limit (or truncation) at high masses. The shape of the ICMF at high masses can be studied using the most massive young (<10 Myr) clusters, however this has proven difficult due to low-number statistics. In this contribution we use an alternative method based on the luminosities of the brightest clusters, combined with their ages. If a truncation is present, a generic prediction (nearly independent of the cluster disruption law adopted) is that the median age of bright clusters should be younger than that of fainter clusters. In the case of an non-truncated ICMF, the median age should be independent of cluster luminosity. Here, we present optical spectroscopy of twelve young stellar clusters in the face-on spiral galaxy NGC 2997. The spectra are used to estimate the age of each cluster, and the brightness of the clusters is taken from the literature. The observations are compared with the model expectations of Larsen (2009) for various ICMF forms and both mass dependent and mass independent cluster disruption. While there exists some degeneracy between the truncation mass and the amount of mass independent disruption, the observations favour a truncated ICMF. For low or modest amounts of mass independent disruption, a truncation mass of 5-6*10^5 Msun is estimated, consistent with previous determinations. Additionally, we investigate possible truncations in the ICMF in the spiral galaxy M83, the interacting Antennae galaxies, and the collection of spiral and dwarf galaxies present in Larsen (2009) based on photometric catalogues taken from the literature, and find that all catalogues are consistent with having a (environmentally dependent) truncation in the cluster mass functions.

The majority of clusters in the Universe have masses well below 10^5 Msun. Hence their integrated fluxes and colors can be affected by the random presence of a few bright stars introduced by stochastic sampling of the stellar mass function. Specific methods are being developed to extend the analysis of cluster SEDs into the low-mass regime. In this paper, we apply such a method to observations of star clusters, in the nearby spiral galaxy M83. We reassess ages and masses of a sample of 1242 objects for which UBVIHalpha fluxes were obtained with the HST/WFC3 images. Synthetic clusters with known properties are used to characterize the limitations of the method. The ensemble of color predictions of the discrete cluster models are in good agreement with the distribution of observed colors. We emphasize the important role of the Halpha data in the assessment of the fraction of young objects, particularly in breaking the age-extinction degeneracy that hampers an analysis based on UBVI only. We find the mass distribution of the cluster sample to follow a power-law of index -2.1 +/-0.2, and the distribution of ages a power-law of index -1.0 +/-0.2 for M > 10^3.5 Msun and ages between 10^7 and 10^9 yr. An extension of our main method, that makes full use of the probability distributions of age and mass of the individual clusters, is explored. It produces similar power-law slopes and will deserve further investigation. Although the properties derived for individual clusters significantly differ from those obtained with traditional, non-stochastic models in ~30% of the objects, the first order aspect of the age and mass distributions are similar to those obtained previously for this M83 sample in the range of overlap of the studies. We extend the power-law description to lower masses with better mass and age resolution and without most of the artifacts produced by the classical method.

Star clusters are studied widely both as benchmarks for stellar evolution models and in their own right. Cluster age and mass distributions within galaxies are probes of star formation histories, and of cluster formation and disruption processes. The vast majority of clusters in the Universe is small, and it is well known that the integrated fluxes and colors have broad probability distributions, due to small numbers of bright stars. This paper goes beyond the description of predicted probability distributions, and presents results of the analysis of cluster energy distributions in an explicitly stochastic context. The method developed is Bayesian. It provides posterior probability distributions in the age-mass-extinction space, using multi-wavelength photometric observations and a large collection of Monte-Carlo simulations of clusters of finite stellar masses. Both UBVI and UBVIK datasets are considered, and the study conducted in this paper is restricted to the solar metallicity. We first reassess and explain errors arising from the use of standard analysis methods, which are based on continuous population synthesis models: systematic errors on ages and random errors on masses are large, while systematic errors on masses tend to be smaller. The age-mass distributions obtained after analysis of a synthetic sample are very similar to those found for real galaxies in the literature. The Bayesian approach on the other hand, is very successful in recovering the input ages and masses. Taking stochastic effects into account is important, more important for instance than the choice of adding or removing near-IR data in many cases. We found no immediately obvious reason to reject priors inspired by previous (standard) analyses of cluster populations in galaxies, i.e. cluster distributions that scale with mass as M^-2 and are uniform on a logarithmic age scale.

The integrated spectro-photometric properties of star clusters are subject to large cluster-to-cluster variations. They are distributed in non trivial ways around the average properties predicted by standard population synthesis models. This results from the stochastic mass distribution of the finite (small) number of luminous stars in each cluster, stars which may be either particularly blue or particularly red. The color distributions are broad and usually far from Gaussian, especially for young and intermediate age clusters, as found in interacting galaxies. When photometric measurements of clusters are used to estimate ages and masses in conjunction with standard models, biases are to be expected. We present a Bayesian approach that explicitly accounts for stochasticity when estimating ages and masses of star clusters that cannot be resolved into stars. Based on Monte-Carlo simulations, we are starting to explore the probability distributions of star cluster properties obtained given a set of multi-wavelength photometric data.

The integrated properties of young star clusters are subject to large cluster-to-cluster variations because they depend directly on small numbers of bright stars. At ages at which red supergiants are expected to exist, these luminous but rare stars can be blamed for most of the variations. If unresolved clusters are to be used to constrain red supergiant models, methods must be developed that take these stochastic fluctuations into account. We discuss prospects and open issues in this field, based on recent work on high mass clusters in M82, and on first experiments towards a Bayesian study of cluster populations.