(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

We present the detection of very extended stellar populations around the Large Magellanic Cloud (LMC) out to R~21 degrees, or ~18.5 kpc at the LMC distance of 50 kpc, as detected in the Survey of the MAgellanic Stellar History (SMASH) performed with the Dark Energy Camera on the NOAO Blanco 4m Telescope. The deep (g~24) SMASH color magnitude diagrams (CMDs) clearly reveal old (~9 Gyr), metal-poor ([Fe/H]=-0.8 dex) main-sequence stars at a distance of 50 kpc. The surface brightness of these detections is extremely low with our most distant detection having 34 mag per arcsec squared in g-band. The SMASH radial density profile breaks from the inner LMC exponential decline at ~13-15 degrees and a second component at larger radii has a shallower slope with power-law index of -2.2 that contributes ~0.4% of the LMC's total stellar mass. In addition, the SMASH densities exhibit large scatter around our best-fit model of ~70% indicating that the envelope of stellar material in the LMC periphery is highly disturbed. We also use data from the NOAO Source catalog to map the LMC main-sequence populations at intermediate radii and detect a steep dropoff in density on the eastern side of the LMC (at R~8 deg) as well as an extended structure to the far northeast. These combined results confirm the existence of a very extended, low-density envelope of stellar material with disturbed shape around the LMC. The exact origin of this structure remains unclear but the leading options include a classical accreted halo or tidally stripped outer disk material.

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image ~14,000 square degrees of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamic observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering an inference-based catalog which includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4um, 4.6um, 12um and 22um) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

We present a study of the three-dimensional (3D) structure of the Large Magellanic Cloud (LMC) using ~2.2 million red clump (RC) stars selected from the Survey of the MAgellanic Stellar History (SMASH). To correct for line-of-sight dust extinction, the intrinsic RC color and magnitude and their radial dependence are carefully measured by using internal nearly dust-free regions. These are then used to construct an accurate 2D reddening map (165 square degrees with ~10 arcmin resolution) of the LMC disk and the 3D spatial distribution of RC stars. An inclined disk model is fit to the 2D distance map yielding a best-fit inclination angle i = 25.86(+0.73,-1.39) degrees with random errors of +\-0.19 degrees, line-of-nodes position angle theta = 149.23(+6.43,-8.35) degrees with random errors of +/-0.49 degrees. These angles vary with galactic radius, indicating that the LMC disk is warped and twisted likely due to the repeated tidal interactions with the Small Magellanic Cloud (SMC). For the first time, our data reveal a significant warp in the southwest of the outer disk starting at rho ~ 7 degrees that departs from the defined LMC plane up to ~4 kpc towards the SMC, suggesting that it originated from a strong interaction with the SMC. In addition, the inner disk encompassing the off-centered bar appears to be tilted up to 5-15 degrees relative to the rest of the LMC disk. These findings on the outer warp and the tilted bar are consistent with the predictions from the Besla et al. (2012) simulation of a recent direct collision with the SMC.

Most massive stars end their lives as Red Supergiants (RSGs), a short-lived evolution phase when they are known to pulsate with varying amplitudes. The RSG period-luminosity (PL) relation has been measured in the Milky Way, the Magellanic Clouds and M33 for about 120 stars in total. Using over 1500 epochs of R-band monitoring from the Palomar Transient Factory (PTF) survey over a five-year period, we study the variability of 255 spectroscopically cataloged RSGs in M31. We find that all RGSs brighter than M_K~ -10 mag (log(L/L_sun)>4.8) are variable at dm_R>0.05 mag. Our period analysis finds 63 with significant pulsation periods. Using the periods found and the known values of M_K for these stars, we derive the RSG PL relation in M31 and show that it is consistent with those derived earlier in other galaxies of different metallicities. We also detect, for the first time, a sequence of likely first-overtone pulsations. Comparison to stellar evolution models from MESA confirms the first overtone hypothesis and indicates that the variable stars in this sample have 12 M_sun<M<24 M_sun. As these RSGs are the immediate progenitors to Type II-P core-collapse supernovae (SNe), we also explore the implication of their variability in the initial-mass estimates for SN progenitors based on archival images of the progenitors. We find that this effect is small compared to the present measurement errors.

The unprecedented volume and rate of transient events that will be discovered by the Large Synoptic Survey Telescope (LSST) demands that the astronomical community update its followup paradigm. Alert-brokers -- automated software system to sift through, characterize, annotate and prioritize events for followup -- will be critical tools for managing alert streams in the LSST era. The Arizona-NOAO Temporal Analysis and Response to Events System (ANTARES) is one such broker. In this work, we develop a machine learning pipeline to characterize and classify variable and transient sources only using the available multiband optical photometry. We describe three illustrative stages of the pipeline, serving the three goals of early, intermediate and retrospective classification of alerts. The first takes the form of variable vs transient categorization, the second, a multi-class typing of the combined variable and transient dataset, and the third, a purity-driven subtyping of a transient class. While several similar algorithms have proven themselves in simulations, we validate their performance on real observations for the first time. We quantitatively evaluate our pipeline on sparse, unevenly sampled, heteroskedastic data from various existing observational campaigns, and demonstrate very competitive classification performance. We describe our progress towards adapting the pipeline developed in this work into a real-time broker working on live alert streams from time-domain surveys.

Ongoing and future surveys with repeat imaging in multiple bands are producing (or will produce) time-spaced measurements of brightness, resulting in the identification of large numbers of variable sources in the sky. A large fraction of these are periodic variables: compilations of these are of scientific interest for a variety of purposes. Unavoidably, the data-sets from many such surveys not only have sparse sampling, but also have embedded frequencies in the observing cadence that beat against the natural periodicities of any object under investigation. Such limitations can make period determination ambiguous and uncertain. For multi-band data sets with asynchronous measurements in multiple pass-bands, we want to maximally utilize the information on periodicity in a manner that is agnostic of differences in the light curve shapes across the different channels. Given large volumes of data, computational efficiency is also at a premium. This paper develops and presents a computationally economic method for determining periodicity which combines the results from two different classes of period determination algorithms. The underlying principles are illustrated through examples. The effectiveness of this approach for combining asynchronously sampled measurements in multiple observables that share an underlying fundamental frequency is also demonstrated.

The Large and Small Magellanic Clouds (LMC and SMC) are unique local laboratories for studying the formation and evolution of small galaxies in exquisite detail. The Survey of the MAgellanic Stellar History (SMASH) is an NOAO community DECam survey of the Clouds mapping 480 square degrees (distributed over ~2400 square degrees at ~20% filling factor) to ~24th mag in ugriz with the goal of identifying broadly distributed, low surface brightness stellar populations associated with the stellar halos and tidal debris of the Clouds. SMASH will also derive spatially-resolved star formation histories covering all ages out to large radii from the MCs that will further complement our understanding of their formation. Here, we present a summary of the survey, its data reduction, and a description of the first public Data Release (DR1). The SMASH DECam data have been reduced with a combination of the NOAO Community Pipeline, PHOTRED, an automated PSF photometry pipeline based mainly on the DAOPHOT suite, and custom calibration software. The attained astrometric precision is ~15 mas and the accuracy is ~2 mas with respect to the Gaia DR1 astrometric reference frame. The photometric precision is ~0.5-0.7% in griz and ~1% in u with a calibration accuracy of ~1.3% in all bands. The median 5 sigma point source depths in ugriz bands are 23.9, 24.8, 24.5, 24.2, 23.5 mag. The SMASH data already have been used to discover the Hydra II Milky Way satellite, the SMASH 1 old globular cluster likely associated with the LMC, and very extended stellar populations around the LMC out to R~18.4 kpc. SMASH DR1 contains measurements of ~100 million objects distributed in 61 fields. A prototype version of the NOAO Data Lab provides data access, including a data discovery tool, SMASH database access, an image cutout service, and a Jupyter notebook server with example notebooks for exploratory analysis.

We characterize the absolute magnitudes and colors of RR Lyrae stars in the globular cluster M5 in the ugriz filter system of the Dark Energy Camera (DECam). We provide empirical Period-Luminosity (P-L) relationships in all 5 bands based on 47 RR Lyrae stars of the type ab and 14 stars of the type c. The P-L relationships were found to be better constrained for the fundamental mode RR Lyrae stars in the riz passbands, with dispersion of 0.03, 0.02 and 0.02 magnitudes, respectively. The dispersion of the color at minimum light was found to be small, supporting the use of this parameter as a means to obtain accurate interstellar extinctions along the line of sight up to the distance of the RR Lyrae star. We found a trend of color at minimum light with pulsational period that, if taken into account, brings the dispersion in color at minimum light to < 0.016 magnitudes for the (r-i), (i-z), and (r-z) colors. These calibrations will be very useful for using RR Lyrae stars from DECam observations as both standard candles for distance determinations and color standards for reddening measurements.

The Arizona-NOAO Temporal Analysis and Response to Events System (ANTARES) is a joint effort of NOAO and the Department of Computer Science at the University of Arizona to build prototype software to process alerts from time-domain surveys, especially LSST, to identify those alerts that must be followed up immediately. Value is added by annotating incoming alerts with existing information from previous surveys and compilations across the electromagnetic spectrum and from the history of past alerts. Comparison against a knowledge repository of properties and features of known or predicted kinds of variable phenomena is used for categorization. The architecture and algorithms being employed are described.

We present initial results from a program aimed at establishing a network of hot DA white dwarfs to serve as spectrophotometric standards for present and future wide-field surveys. These stars span the equatorial zone and are faint enough to be conveniently observed throughout the year with large-aperture telescopes. Spectra of these white dwarfs are analyzed to generate a non-local-thermodynamic-equilibrium (NLTE) model atmosphere normalized to HST colors, including adjustments for wavelength-dependent interstellar extinction. Once established, this standard star network will serve ground-based observatories in both hemispheres as well as space-based instrumentation from the UV to the near IR. We demonstrate the effectiveness of this concept and show how two different approaches to the problem using somewhat different assumptions produce equivalent results. We discuss lessons learned and the resulting corrective actions applied to our program.

We report the discovery of one RR Lyrae star in the ultra--faint satellite galaxy Hydra II based on time series photometry in the g, r and i bands obtained with the Dark Energy Camera at Cerro Tololo Interamerican Observatory, Chile. The RR Lyrae star has a mean magnitude of $i = 21.30\pm 0.04$ which translates to a heliocentric distance of $151\pm 8$ kpc for Hydra II; this value is $\sim 13\%$ larger than the estimate from the discovery paper based on the average magnitude of several blue horizontal branch star candidates. The new distance implies a slightly larger half-light radius of $76^{+12}_{-10}$ pc and a brighter absolute magnitude of $M_V = -5.1 \pm 0.3$, which keeps this object within the realm of the dwarf galaxies. The pulsational properties of the RR Lyrae star ($P=0.645$ d, $\Delta g = 0.68$ mag) suggest Hydra II may be a member of the intermediate Oosterhoff or Oosterhoff II group. A comparison with other RR Lyrae stars in ultra--faint systems indicates similar pulsational properties among them, which are different to those found among halo field stars and those in the largest of the Milky Way satellites. We also report the discovery of 31 additional short period variables in the field of view (RR Lyrae, SX Phe, eclipsing binaries, and a likely anomalous cepheid). However, given their magnitudes and large angular separation from Hydra II, they must be field stars not related to Hydra II.

We explore the assembly history of the M31 bulge within a projected major-axis radius of 180" (~680 pc) by studying its stellar populations in Hubble Space Telescope WFC3 and ACS observations. Colors formed by comparing near-ultraviolet vs. optical bands are found to become bluer with increasing major-axis radius, which is opposite to that predicted if the sole sources of near-ultraviolet light were old extreme horizontal branch stars with a negative radial gradient in metallicity. Spectral energy distribution fits require a metal-rich intermediate-age stellar population (300 Myr to 1 Gyr old, ~solar metallicity) in addition to the dominant old population. The radial gradients in age and metallicity of the old stellar population are consistent with those in previous works. For the intermediate-age population, we find an increase in age with radius and a mass fraction that increases up to 2% at 680 pc away from the center. We exclude contamination from the M31 disk and/or halo as the main origin for this population. Our results thus suggest that intermediate-age stars exist beyond the central 5" (19 pc) of M31 and contribute ~1% of the total stellar mass in the bulge. These stars could be related to the secular growth of the M31 bulge.

We present stellar metallicities in Leo I, Leo II, IC 1613, and Phoenix dwarf galaxies derived from medium (F390M) and broad (F555W, F814W) band photometry using the Wide Field Camera 3 (WFC3) instrument aboard the Hubble Space Telescope. We measured metallicity distribution functions (MDFs) in two ways, 1) matching stars to isochrones in color-color diagrams, and 2) solving for the best linear combination of synthetic populations to match the observed color-color diagram. The synthetic technique reduces the effect of photometric scatter, and produces MDFs 30-50 % narrower than the MDFs produced from individually matched stars. We fit the synthetic and individual MDFs to analytical chemical evolution models (CEM) to quantify the enrichment and the effect of gas flows within the galaxies. Additionally, we measure stellar metallicity gradients in Leo I and II. For IC 1613 and Phoenix our data do not have the radial extent to confirm a metallicity gradient for either galaxy. We find the MDF of Leo I (dwarf spheroidal) to be very peaked with a steep metal rich cutoff and an extended metal poor tail, while Leo II (dwarf spheroidal), Phoenix (dwarf transition) and IC 1613 (dwarf irregular) have wider, less peaked MDFs than Leo I. A simple CEM is not the best fit for any of our galaxies, therefore we also fit the `Best Accretion Model' of Lynden-Bell 1975. For Leo II, IC 1613 and Phoenix we find similar accretion parameters for the CEM, even though they all have different effective yields, masses, star formation histories and morphologies. We suggest that the dynamical history of a galaxy is reflected in the MDF, where broad MDFs are seen in galaxies that have chemically evolved in relative isolation and narrowly peaked MDFs are seen in galaxies that have experienced more complicated dynamical interactions concurrent with their chemical evolution.

We present the discovery of a new dwarf galaxy, Hydra II, found serendipitously within the data from the ongoing Survey of the MAgellanic Stellar History (SMASH) conducted with the Dark Energy Camera on the Blanco 4m Telescope. The new satellite is compact (r_h = 68 +/- 11 pc) and faint (M_V = -4.8 +/- 0.3), but well within the realm of dwarf galaxies. The stellar distribution of HydraII in the color-magnitude diagram is well-described by a metal-poor ([Fe/H] = -2.2) and old (13 Gyr) isochrone and shows a distinct blue horizontal branch, some possible red clump stars, and faint stars that are suggestive of blue stragglers. At a heliocentric distance of 134 +/- 10 kpc, Hydra II is located in a region of the Galactic halo that models have suggested may host material from the leading arm of the Magellanic Stream. A comparison with N-body simulations hints that the new dwarf galaxy could be or could have been a satellite of the Magellanic Clouds.

The Arizona-NOAO Temporal Analysis and Response to Events System (ANTARES) is a joint project of the National Optical Astronomy Observatory and the Department of Computer Science at the University of Arizona. The goal is to build the software infrastructure necessary to process and filter alerts produced by time-domain surveys, with the ultimate source of such alerts being the Large Synoptic Survey Telescope (LSST). The ANTARES broker will add value to alerts by annotating them with information from external sources such as previous surveys from across the electromagnetic spectrum. In addition, the temporal history of annotated alerts will provide further annotation for analysis. These alerts will go through a cascade of filters to select interesting candidates. For the prototype, `interesting' is defined as the rarest or most unusual alert, but future systems will accommodate multiple filtering goals. The system is designed to be flexible, allowing users to access the stream at multiple points throughout the process, and to insert custom filters where necessary. We describe the basic architecture of ANTARES and the principles that will guide development and implementation.

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.

The science goals for ground-based large-area surveys, such as the Dark Energy Survey, Pan-STARRS, and the Large Synoptic Survey Telescope, require calibration of broadband photometry that is stable in time and uniform over the sky to precisions of a per cent or better. This performance will need to be achieved with data taken over the course of many years, and often in less than ideal conditions. This paper describes a strategy to achieve precise internal calibration of imaging survey data taken in less than photometric conditions, and reports results of an observational study of the techniques needed to implement this strategy. We find that images of celestial fields used in this case study with stellar densities of order one per arcmin-squared and taken through cloudless skies can be calibrated with relative precision of 0.5 per cent (reproducibility). We report measurements of spatial structure functions of cloud absorption observed over a range of atmospheric conditions, and find it possible to achieve photometric measurements that are reproducible to 1 per cent in images that were taken through cloud layers that transmit as little as 25 per cent of the incident optical flux (1.5 magnitudes of extinction). We find, however, that photometric precision below 1 per cent is impeded by the thinnest detectable cloud layers. We comment on implications of these results for the observing strategies of future surveys.

We quantified and calibrated the metallicity and temperature sensitivities of colors derived from nine Wide Field Camera 3 (WFC3) filters aboard the Hubble Space Telescope using Dartmouth isochrones and Kurucz atmospheres models. The theoretical isochrone colors were tested and calibrated against observations of five well studied galactic clusters: M92, NGC 6752, NGC 104, NGC 5927, and NGC 6791, all of which have spectroscopically determined metallicities spanning -2.30 < [Fe/H] < +0.4. We found empirical corrections to the Dartmouth isochrone grid for each of the following color magnitude diagrams (CMD) (F555W--F814W, F814W), (F336W-F555W, F814W), (F390M-F555W, F814W) and (F390W-F555W, F814W). Using the empirical corrections we tested the accuracy and spread of the photometric metallicities assigned from CMDs and color-color diagrams (which are necessary to break the age-metallicity degeneracy). Testing three color-color diagrams [(F336W-F555W),(F390M-F555W),(F390W-F555W), vs (F555W-F814W)], we found the colors (F390M-F555W) and (F390W-F555W), to be the best suited to measure photometric metallicities. The color (F390W-F555W) requires much less integration time, but generally produces wider metallicity distributions, and, at very-low metallicity, the MDF from (F390W-F555W) is ~60% wider than that from (F390M-F555W). Using the calibrated isochrones we recovered the overall cluster metallicity to within ~0.1 dex in [Fe/H] when using CMDs (i.e. when the distance, reddening and ages are approximately known). The measured metallicity distribution function (MDF) from color-color diagrams show this method measures metallicities of stellar clusters of unknown age and metallicity with an accuracy of ~0.2 - 0.5 dex using F336W--F555W, ~0.15 - 0.25 dex using F390M-F555W, and ~0.2 - 0.4 dex with F390W-F555W, with the larger uncertainty pertaining to the lowest metallicity range.

This document summarizes the results of a community-based discussion of the potential science impact of the Mayall+BigBOSS highly multiplexed multi-object spectroscopic capability. The KPNO Mayall 4m telescope equipped with the DOE- and internationally-funded BigBOSS spectrograph offers one of the most cost-efficient ways of accomplishing many of the pressing scientific goals identified for this decade by the "New Worlds, New Horizons" report. The BigBOSS Key Project will place unprecedented constraints on cosmological parameters related to the expansion history of the universe. With the addition of an open (publicly funded) community access component, the scientific impact of BigBOSS can be extended to many important astrophysical questions related to the origin and evolution of galaxies, stars, and the IGM. Massive spectroscopy is the critical missing ingredient in numerous ongoing and planned ground- and space-based surveys, and BigBOSS is unique in its ability to provide this to the US community. BigBOSS data from community-led projects will play a vital role in the education and training of students and in maintaining US leadership in these fields of astrophysics. We urge the NSF-AST division to support community science with the BigBOSS multi-object spectrograph through the period of the BigBOSS survey in order to ensure public access to the extraordinary spectroscopic capability.

We present a multi-wavelength photometric study of ~15,000 resolved stars in the nearby spiral galaxy M83 (NGC5236, D=4.61Mpc) based on Hubble Space Telescope Wide Field Camera 3 observations using four filters: F336W, F438W, F555W, and F814W. We select 50 regions (an average size of 260 pc by 280 pc) in the spiral arm and inter-arm areas of M83, and determine the age distribution of the luminous stellar populations in each region. This is accomplished by correcting for extinction towards each individual star by comparing its colors with predictions from stellar isochrones. We compare the resulting luminosity weighted mean ages of the luminous stars in the 50 regions with those determined from several independent methods, including the number ratio of red-to-blue supergiants, morphological appearance of the regions, surface brightness fluctuations, and the ages of clusters in the regions. We find reasonably good agreement between these methods. We also find that young stars are much more likely to be found in concentrated aggregates along spiral arms, while older stars are more dispersed. These results are consistent with the scenario that star formation is associated with the spiral arms, and stars form primarily in star clusters and then disperse on short timescales to form the field population. The locations of Wolf-Rayet stars are found to correlate with the positions of many of the youngest regions, providing additional support for our ability to accurately estimate ages. We address the effects of spatial resolution on the measured colors, magnitudes, and age estimates. While individual stars can occasionally show measurable differences in the colors and magnitudes, the age estimates for entire regions are only slightly affected.

Using observations obtained with the Wide Field Camera 3 (WFC3) on board the Hubble Space Telescope (HST), we have studied the properties of the stellar populations in the central regions of 30 Dor, in the Large Magellanic Cloud. The observations clearly reveal the presence of considerable differential extinction across the field. We characterise and quantify this effect using young massive main sequence stars to derive a statistical reddening correction for most objects in the field. We then search for pre-main sequence (PMS) stars by looking for objects with a strong (> 4 sigma) Halpha excess emission and find about 1150 of them over the entire field. Comparison of their location in the Hertzsprung-Russell diagram with theoretical PMS evolutionary tracks for the appropriate metallicity reveals that about one third of these objects are younger than ~4Myr, compatible with the age of the massive stars in the central ionising cluster R136, whereas the rest have ages up to ~30Myr, with a median age of ~12Myr. This indicates that star formation has proceeded over an extended period of time, although we cannot discriminate between an extended episode and a series of short and frequent bursts that are not resolved in time. While the younger PMS population preferentially occupies the central regions of the cluster, older PMS objects are more uniformly distributed across the field and are remarkably few at the very centre of the cluster. We attribute this latter effect to photoevaporation of the older circumstellar discs caused by the massive ionising members of R136.

We use new WFC3 observations of the nearby grand design spiral galaxy M83 to develop two independent methods for estimating the ages of young star clusters. The first method uses the physical extent and morphology of Halpha emission to estimate the ages of clusters younger than tau ~10 Myr. It is based on the simple premise that the gas in very young (tau < few Myr) clusters is largely coincident with the cluster stars, is in a small, ring-like structure surrounding the stars in slightly older clusters (e.g., tau ~5 Myr), and is in a larger ring-like bubble for still older clusters (i.e., ~5-10 Myr). The second method is based on an observed relation between pixel-to-pixel flux variations within clusters and their ages. This method relies on the fact that the brightest individual stars in a cluster are most prominent at ages around 10 Myr, and fall below the detection limit (i.e., M_V < -3.5) for ages older than about 100 Myr. These two methods are the basis for a new morphological classification system which can be used to estimate the ages of star clusters based on their appearance. We compare previous age estimates of clusters in M83 determined from fitting UBVI Halpha measurements using predictions from stellar evolutionary models with our new morphological categories and find good agreement at the ~95% level. The scatter within categories is ~0.1 dex in log tau for young clusters (<10 Myr) and ~0.5 dex for older (>10 Myr) clusters. A by-product of this study is the identification of 22 "single-star" HII regions in M83, with central stars having ages ~4 Myr.

We present new observations of two Galactic globular clusters, PAL4 and PAL14, using the Wide-Field Channel of the Advanced Camera for Surveys (ACS) on board the Hubble Space Telescope (HST), and reanalyze archival data from a third, NGC2419. We matched our photometry of hundreds of stars in these fields from the ACS images to existing, ground-based photometry of faint sequences which were calibrated on the standard BVRI system of Landolt. These stars are significantly fainter than those generally used for HST calibration purposes, and therefore are much better matched to supporting precision photometry of ACS science targets. We were able to derive more accurate photometric transformation coefficients for the commonly used ACS broad-band filters compared to those published by Sirianni, et al. (2005), owing to the use of a factor of several more calibration stars which span a greater range of color. We find that the inferred transformations from each cluster individually do not vary significantly from the average, except for a small offset of the photometric zeropoint in the F850LP filter. Our results suggest that the published prescriptions for the time-dependent correction of CCD charge-transfer efficiency appear to work very well over the ~3.5 yr interval that spans our observations of PAL4 and PAL14 and the archived images of NGC2419.

We investigate the ionization structure of the nebular gas in M83 using the line diagnostic diagram, [O III](5007 \degA)/H{\beta} vs. [S II](6716 \deg A+6731 \deg A)/H{\alpha} with the newly available narrowband images from the Wide Field Camera 3 (WFC3) of the Hubble Space Telescope (HST). We produce the diagnostic diagram on a pixel-by-pixel (0.2" x 0.2") basis and compare it with several photo- and shock-ionization models. For the photo-ionized gas, we observe a gradual increase of the log([O III]/H{\beta}) ratios from the center to the spiral arm, consistent with the metallicity gradient, as the H II regions go from super solar abundance to roughly solar abundance from the center out. Using the diagnostic diagram, we separate the photo-ionized from the shock-ionized component of the gas. We find that the shock-ionized H{\alpha} emission ranges from ~2% to about 15-33% of the total, depending on the separation criteria used. An interesting feature in the diagnostic diagram is an horizontal distribution around log([O III]/H{\beta}) ~ 0. This feature is well fit by a shock-ionization model with 2.0 Z\odot metallicity and shock velocities in the range of 250 km/s to 350 km/s. A low velocity shock component, < 200 km/s, is also detected, and is spatially located at the boundary between the outer ring and the spiral arm. The low velocity shock component can be due to : 1) supernova remnants located nearby, 2) dynamical interaction between the outer ring and the spiral arm, 3) abnormal line ratios from extreme local dust extinction. The current data do not enable us to distinguish among those three possible interpretations. Our main conclusion is that, even at the HST resolution, the shocked gas represents a small fraction of the total ionized gas emission at less than 33% of the total. However, it accounts for virtually all of the mechanical energy produced by the central starburst in M83.