We present a new method to measure colors and magnitudes of the tip of the red giant branch in multiple bandpasses simultaneously by fitting an n-dimensional Gaussian to photometry of candidate tip stars. We demonstrate that this method has several advantages over traditional edge detection, particularly in regimes where the TRGB magnitude is strongly color-dependent, as is the case in the near-infrared. We apply this method to a re-reduction of a set of optical and near-IR HST data originally presented in Dalcanton et al. (2012). The re-reduction takes advantage of the increased depth and accuracy in the NIR photometry enabled by simultaneous reduction with higher resolution optical data in crowded fields (Williams et al. 2014). We compare three possible absolute calibrations of the resulting apparent TRGB measurements, one adopting the same distance moduli as in Dalcanton et al. (2012), and two based on predicted TRGB absolute magnitudes from two widely-used, modern sets of model isochrones. We find systematic offsets among the model absolute calibrations at the ~0.1 mag level, in line with previous investigations. The models also have difficulty reproducing the optical-NIR color-magnitude behavior of our measurements, making these observations a useful benchmark for future improvements.

The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a flattened shape, discrete geological units, and noticeable albedo heterogeneity. However, there is little surface color and compositional heterogeneity. No evidence for satellites, ring or dust structures, gas coma, or solar wind interactions was detected. By origin MU69 appears consistent with pebble cloud collapse followed by a low velocity merger of its two lobes.

Aims: We use new deep XMM-Newton observations of the northern disk of M 31 to trace the hot interstellar medium (ISM) in unprecedented detail and to characterise the physical properties of the X-ray emitting plasmas. Methods: We used all XMM-Newton data up to and including our new observations to produce the most detailed image yet of the hot ISM plasma in a grand design spiral galaxy such as our own. We compared the X-ray morphology to multi-wavelength studies in the literature to set it in the context of the multi-phase ISM. We performed spectral analyses on the extended emission using our new observations as they offer sufficient depth and count statistics to constrain the plasma properties. Data from the Panchromatic Hubble Andromeda Treasury were used to estimate the energy injected by massive stars and their supernovae. We compared these results to the hot gas properties. Results: The brightest emission regions were found to be correlated with populations of massive stars, notably in the 10 kpc star-forming ring. The plasma temperatures in the ring regions are ~0.2 keV up to ~0.6 keV. We suggest this emission is hot ISM heated in massive stellar clusters and superbubbles. We derived X-ray luminosities, densities, and pressures for the gas in each region. We also found large extended emission filling low density gaps in the dust morphology of the northern disk, notably between the 5 kpc and 10 kpc star-forming rings. We propose that the hot gas was heated and expelled into the gaps by the populations of massive stars in the rings. Conclusions. It is clear that the massive stellar populations are responsible for heating the ISM to X-ray emitting temperatures, filling their surroundings, and possibly driving the hot gas into the low density regions. Overall, the morphology and spectra of the hot gas in the northern disk of M 31 is similar to other galaxy disks.

We report, for the first time, the detection of the Mn-K$\alpha$ line in the Type IIb supernova (SN IIb) remnant, Cassiopeia A. Manganese ($^{55}$Mn after decay of $^{55}$Co), a neutron-rich element, together with chromium ($^{52}$Cr after decay of $^{52}$Fe), is mainly synthesized at the explosive incomplete Si burning regime. Therefore, the Mn/Cr mass ratio with its neutron excess reflects the neutronization at the relevant burning layer during the explosion. Chandra's archival X-ray data of Cassiopeia A indicate a low Mn/Cr mass ratio with values in the range 0.10--0.66, which, when compared to one-dimensional SN explosion models, requires that the electron fraction be 0.4990 $\lesssim Y_{\rm e} \lesssim$ 0.5 at the incomplete Si burning layer. An explosion model assuming a solar-metallicity progenitor with a typical explosion energy ($1 \times 10^{51}$ erg) fails to reproduce such a high electron fraction. In such models, the explosive Si-burning regime extends only to the Si/O layer established during the progenitor's hydrostatic evolution; the $Y_e$ in the Si/O layer is lower than the value required by our observational constraints. We can satisfy the observed Mn/Cr mass ratio if the explosive Si-burning regime were to extend into the O/Ne hydrostatic layer, which has a higher $Y_{\rm e}$. This would require an energetic ($> 2 \times 10^{51}$ erg) and/or asymmetric explosion of a sub-solar metallicity progenitor ($Z \lesssim 0.5Z_{\odot}$) for Cassiopeia A. The low initial metallicity can be used to rule out a single-star progenitor, leaving the possibility of a binary progenitor with a compact companion (white dwarf, neutron star or black hole). We discuss the detectability of X-rays from Bondi accretion onto such a compact companion around the explosion site. We also discuss other possible mass-loss scenarios for the progenitor system of Cassiopeia A.

A galaxy's stellar mass-to-light ratio ($M_\star/L$) is a useful tool for converting luminosity to stellar mass ($M_\star$). However, the practical utility of $M_\star/L$ inferred from stellar population synthesis (SPS) models is limited by mismatches between the real and assumed models for star formation history (SFH) and dust geometry, both of which vary within galaxies. Here, we measure spatial variations in $M_\star/L$ and their dependence on color, star formation history, and dust across the disk of M31, using a map of $M^\mathrm{CMD}_\star$ derived from color-magnitude diagrams of resolved stars in the Panchromatic Hubble Andromeda Treasury (PHAT) survey. First, we find comparable scatter in $M_\star/L$ for the optical and mid-IR, contrary to the common idea that $M_\star/L$ is less variable in the IR. Second, we confirm that $M_\star/L$ is correlated with color for both the optical and mid-IR and report color vs. $M_\star/L$ relations (CMLRs) in M31 for filters used in the Sloan Digital Sky Survey (SDSS) and Widefield Infrared Survey Explorer (WISE). Third, we show that the CMLR residuals correlate with recent SFH, such that quiescent regions are offset to higher $M_\star/L$ than star-forming regions at a fixed color. The mid-IR CMLR, however, is not linear due to the high scatter of $M_\star/L$ in star-forming regions. Finally, we find a flatter optical CMLR than any SPS-based CMLRs in the literature. We show this is an effect of dust geometry, which is typically neglected but should be accounted for when using optical data to map $M_\star/L$.

We present results from a 2-epoch HST H$\alpha$ emission line survey of the Andromeda Galaxy that overlaps the footprint of the Panchromatic Hubble Andromeda Treasury (PHAT) survey. We find 552 (542) classical Be stars and 8429 (8556) normal B-type stars in epoch # 1 (epoch # 2), yielding an overall fractional Be content of 6.15% $\pm$0.26% (5.96% $\pm$0.25%). The fractional Be content decreased with spectral sub-type from $\sim$23.6% $\pm$2.0% ($\sim$23.9% $\pm$2.0%) for B0-type stars to $\sim$3.1% $\pm$0.34% ($\sim$3.4% $\pm$0.35%) for B8-type stars in epoch # 1 (epoch # 2). We observe a clear population of cluster Be stars at early fractional main sequence lifetimes, indicating that a subset of Be stars emerge onto the ZAMS as rapid rotators. Be stars are 2.8x rarer in M31 for the earliest sub-types compared to the SMC, confirming that the fractional Be content decreases in significantly more metal rich environments (like the Milky Way and M31). However, M31 does not follow a clear trend of Be fraction decreasing with metallicity compared to the Milky Way, which may reflect that the Be phenomenon is enhanced with evolutionary age. The rate of disk-loss or disk-regeneration episodes we observe, 22% $\pm$ 2% yr$^{-1}$, is similar to that observed for seven other Galactic clusters reported in the literature, assuming these latter transient fractions scale by a linear rate. The similar number of disk-loss events (57) as disk-renewal events (43) was unexpected since disk dissipation time-scales can be $\sim$2x the typical time-scales for disk build-up phases.

{\it Chandra} X-ray observations of Kepler's supernova remnant indicate the existence of a high speed Fe-rich ejecta structure in the southwestern region. We report strong K-shell emission from Fe-peak elements (Cr, Mn, Fe, Ni), as well as Ca, in this Fe-rich structure, implying that those elements could be produced in the inner area of the exploding white dwarf. We found Ca/Fe, Cr/Fe, Mn/Fe and Ni/Fe mass ratios of 1.0--4.1\%, 1.0--4.6\%, 1--11\% and 2--30\%, respectively. In order to constrain the burning regime that could produce this structure, we compared these observed mass ratios with those in 18 one-dimensional Type Ia nucleosynthesis models (including both near-$M_{\rm Ch}$ and sub-$M_{\rm Ch}$ explosion models). The observed mass ratios agree well with those around the middle layer of incomplete Si-burning in Type Ia nucleosynthesis models with a peak temperature of $\sim$(5.0--5.3)$\times$10$^{9}$ K and a high metallicity, Z $>$ 0.0225. Based on our results, we infer the necessity for some mechanism to produce protruding Fe-rich clumps dominated by incomplete Si-burning products during the explosion. We also discuss the future perspectives of X-ray observations of Fe-rich structures in other Type Ia supernova remnants.

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.

Photodissociation regions (PDRs) are parts of the ISM consisting of predominantly neutral gas, located at the interface between H II regions and molecular clouds. The physical conditions within these regions show variations on very short spatial scales, and therefore PDRs constitute ideal laboratories for investigating the properties and evolution of dust grains. We have mapped IC 63 at high resolution from the UV to the NIR (275 nm to 1.6 $\mu$m), using the Hubble Space Telescope WFC3. Using a Bayesian SED fitting tool, we simultaneously derive a set of stellar ($T_\text{eff}$, $\log(g)$, distance) and extinction ($A_V$, $R_V$) parameters for 520 background stars. We present maps of $A_V$ and $R_V$ with a resolution of 25 arcsec based on these results. The extinction properties vary across the PDR, with values for $A_V$ between 0.5 and 1.4 mag, and a decreasing trend in $R_V$, going from 3.7 at the front of the nebula to values as low as 2.5 further in. This provides evidence for evolution of the dust optical properties. We fit two modified blackbodies to the MIR and FIR SED, obtained by combining the $A_V$ map with data from Spitzer and Herschel. We derive effective temperatures (30 K and 227 K) and the ratio of opacities at 160 $\mu$m to V band $\kappa_{160} / \kappa_V$ ($7.0 \times 10^{-4}$ and $2.9 \times 10^{-9}$) for the two dust populations. Similar fits to individual pixels show spatial variations of $\kappa_{160} / \kappa_{V}$. The analysis of our HST data, combined with these Spitzer and Herschel data, provides the first panchromatic view of dust within a PDR.

Hickson Compact Group (HCG) 16 is a prototypical compact group of galaxies, apparently in an intermediate stage of evolution, where its galaxies are losing gas to the intra-group medium (IGrM). The wide variety of ongoing process in HCG 16 make it an ideal case study for exploring which processes are likely to dominate the late stages of evolution in compact groups. In order to build a coherent picture of the evolution of this group we make use of the multi-wavelength data available, but focus particularly on HI as a tracer of interactions and evolutionary phase. We reprocess archival VLA L-band data of HCG 16 using the multi-scale CLEAN algorithm to accurately recover diffuse features. Despite the clear disruption of the HI component of HCG 16 we find that it is not globally HI deficient, even though HCG 16a and b have lost the majority of their HI and almost 50% of the group's HI is in the IGrM, including a 160 kpc long tail extending towards the SE. This study indicates that in the recent past (~1 Gyr) galaxies HCG 16a and b likely underwent major interactions that unbound gas without triggering significant star formation. This gas was then swept away by a high speed, close encounter with NGC 848. The starburst events HCG 16c and d have triggered galactic winds which, in the case of HCG 16d, appears to have disrupted its HI reservoir. The tidal features still connected to these galaxies indicate that more HI will soon be lost to the IGrM, while that which remains in the discs will likely be consumed by star formation episodes. This is expected to result in a collection of gas-poor galaxies embedded in a diffuse HI structure, which will gradually be evaporated by the UV background, resembling the final stage of the evolutionary model of compact groups. This work is accompanied by a complete workflow to reproduce the final data products, intended to improve reproducibility. (Abridged)

We present a source catalog from the first deep hard X-ray ($E>10$ keV) survey of the Small Magellanic Cloud (SMC), the NuSTAR Legacy Survey of the SMC. We observed three fields, for a total exposure time of 1 Ms, along the bar of this nearby star-forming galaxy. Fields were chosen for their young stellar and accreting binary populations. We detected 10 sources above a 3$\sigma$ significance level (4$-$25 keV) and obtained upper limits on an additional 40 sources. We reached a 3$\sigma$ limiting luminosity in the 4$-$25 keV band of $\sim$ $10^{35}$ erg s$^{-1}$, allowing us to probe fainter X-ray binary (XRB) populations than has been possible with other extragalactic NuSTAR surveys. We used hard X-ray colors and luminosities to constrain the compact-object type, exploiting the spectral differences between accreting black holes and neutron stars at $E>10$ keV. Several of our sources demonstrate variability consistent with previously observed behavior. We confirmed pulsations for seven pulsars in our 3$\sigma$ sample. We present the first detection of pulsations from a Be-XRB, SXP305 (CXO J005215.4$-$73191), with an X-ray pulse period of $305.69\pm0.16$ seconds and a likely orbital period of $\sim$1160-1180 days. Bright sources ($\gtrsim 5\times 10^{36}$ erg s$^{-1}$) in our sample have compact-object classifications consistent with their previously reported types in the literature. Lower luminosity sources ($\lesssim 5\times 10^{36}$ erg s$^{-1}$) have X-ray colors and luminosities consistent with multiple classifications. We raise questions about possible spectral differences at low luminosity between SMC pulsars and the Galactic pulsars used to create the diagnostic diagrams.

We present the star formation history (SFH) of the isolated (D~970 kpc) Local Group dwarf galaxy WLM measured from color-magnitude diagrams constructed from deep Hubble Space Telescope imaging. Our observations include a central (0.5 $r_h$) and outer field (0.7 $ r_h$) that reach below the oldest main sequence turnoff. WLM has no early dominant episode of star formation: 20% of its stellar mass formed by ~12.5 Gyr ago (z~5). It also has an SFR that rises to the present with 50% of the stellar mass within the most recent 5 Gyr (z<0.7). There is evidence of a strong age gradient: the mean age of the outer field is 5 Gyr older than the inner field despite being only 0.4 kpc apart. Some models suggest such steep gradients are associated with strong stellar feedback and dark matter core creation. The SFHs of real isolated dwarf galaxies and those from the the Feedback In Realistic Environment suite are in good agreement for $M_{\star}(z=0) \sim 10^7-10^9 M_{\odot}$, but in worse agreement at lower masses ($M_{\star}(z=0) \sim 10^5-10^7 M_{\odot}$). These differences may be explainable by systematics in the models (e.g., reionization model) and/or observations (HST field placement). We suggest that a coordinated effort to get deep CMDs between HST/JWST (crowded central fields) and WFIRST (wide-area halo coverage) is the optimal path for measuring global SFHs of isolated dwarf galaxies.

High-multiplex and deep spectroscopic follow-up of upcoming panoramic deep-imaging surveys like LSST, Euclid, and WFIRST is a widely recognized and increasingly urgent necessity. No current or planned facility at a U.S. observatory meets the sensitivity, multiplex, and rapid-response time needed to exploit these future datasets. FOBOS, the Fiber-Optic Broadband Optical Spectrograph, is a near-term fiber-based facility that addresses these spectroscopic needs by optimizing depth over area and exploiting the aperture advantage of the existing 10m Keck II Telescope. The result is an instrument with a uniquely blue-sensitive wavelength range (0.31-1.0 um) at R~3500, high-multiplex (1800 fibers), and a factor 1.7 greater survey speed and order-of-magnitude greater sampling density than Subaru's Prime Focus Spectrograph (PFS). In the era of panoramic deep imaging, FOBOS will excel at building the deep, spectroscopic reference data sets needed to interpret vast imaging data. At the same time, its flexible focal plane, including a mode with 25 deployable integral-field units (IFUs) across a 20 arcmin diameter field, enables an expansive range of scientific investigations. Its key programmatic areas include (1) nested stellar-parameter training sets that enable studies of the Milky Way and M31 halo sub-structure, as well as local group dwarf galaxies, (2) a comprehensive picture of galaxy formation thanks to detailed mapping of the baryonic environment at z~2 and statistical linking of evolving populations to the present day, and (3) dramatic enhancements in cosmological constraints via precise photometric redshifts and determined redshift distributions. In combination with Keck I instrumentation, FOBOS also provides instant access to medium-resolution spectroscopy for transient sources with full coverage from the UV to the K-band.

We present a multiwavelength examination of the supernova remnant (SNR) S26 in the nearby galaxy NGC 300 using data from Chandra X-ray Observatory, XMM-Newton X-ray Observatory, Hubble Space Telescope (HST), the Very Large Array, and the Australia Telescope Compact Array. We simultaneously fit all of the available X-ray data with a thermal plasma model and find a temperature of $0.77 \pm 0.13$ keV with a hydrogen column density of ($9.7^{+6.4}_{-4.8}$)$\times 10^{20}$ cm$^{-2}$. HST imaging allows us to measure a semimajor axis of $0.78 \pm 0.10$ arcsec ($7.5 \pm 1.0$ pc) and a semiminor axis of $0.69^{+0.14}_{-0.12}$ arcsec ($6.7^{+1.2}_{-1.4}$ pc). This precise size helps to constrain the age and velocity of the shock to be ($3.3^{+0.7}_{-0.6}$)$\times 10^{3}$ yr and $411^{+275}_{-122}$ km s$^{-1}$. We also fit photometry of the surrounding stars to infer the age and mass of the progenitor star to be $8 \pm 1$ Myr and $25^{+1}_{-5}$ M$_{\odot}$. Based on measured radio properties of the source and assuming equipartition, the estimated radio luminosity of $\sim 1.7 \times 10^{34}$ erg s$^{-1}$ over the $10^{8}-10^{11}$ Hz frequency range results in a minimum magnetic field associated with this SNR of $0.067$ mG and the minimum energy needed to power the observed synchrotron emission of $1.5 \times 10^{49}$ erg. The size and temperature of N300-S26 appear to be similar to the Galactic SNR G311.5-0.3 except that G311.5-0.3 has a significantly lower X-ray luminosity, is older, and has a slower shock velocity.

We present the detection of the first candidate colliding-wind binary (CWB) in M33, located in the giant H II region NGC 604. The source was first identified in archival {\it Chandra} imaging as a relatively soft X-ray point source, with the likely primary star determined from precise astrometric alignment between archival {\it Hubble Space Telescope} and {\it Chandra} imaging. The candidate primary star in the CWB is classified for the first time in this work as a carbon-rich Wolf-Rayet star with a likely O star companion based on spectroscopy obtained from Gemini-North. We model the X-ray spectrum using {\it Chandra} and {\it XMM-Newton} observations, and find the CWB is well-fit as a $\sim$ 1 keV thermal plasma with a median unabsorbed luminosity in the 0.5--2.0 keV band of $L_{\rm X}$ $\sim$ 3 $\times$ 10$^{35}$ erg s$^{-1}$, making this source among the brightest of CWBs observed to date. We present a long term light curve for the candidate CWB from archival {\it Chandra} and {\it XMM-Newton} observations, and discuss the constraints placed on the binary by this light curve, as well as the X-ray luminosity at maximum. Finally, we compare this candidate CWB in M33 to other well-studied, bright CWBs in the Galaxy and Magellanic Clouds, such as $\eta$ Car.

The morphology of thin stellar streams can be used to test the nature of dark matter. It is therefore crucial to extend searches for globular cluster streams to other galaxies than the Milky Way. In this paper, we investigate the current and future prospects of detecting globular cluster streams in external galaxies in resolved stars (e.g. with WFIRST) and using integrated light (e.g. with HSC, LSST and Euclid). In particular, we inject mock-streams to data from the PAndAS M31 survey, and produce simulated M31 backgrounds mimicking what WFIRST will observe in M31. Additionally, we estimate the distance limit to which globular cluster streams will be observable. Our results demonstrate that for a 1 hour (1000 sec.) exposure, using conservative estimates, WFIRST should detect globular cluster streams in resolved stars in galaxies out to distances of ~3.5 Mpc (~2 Mpc). This volume contains 199 (122) galaxies of which >90% are dwarfs. With integrated light, thin streams can be resolved out to ~100 Mpc with HSC and LSST and to ~600 Mpc with WFIRST and Euclid. The low surface brightness of the streams (typically >30 mag/arcsec$^2$), however, will make them difficult to detect, unless the streams originate from very young clusters. We emphasize that if the external galaxies do not host spiral arms or galactic bars, gaps in their stellar streams provide an ideal test case for evidence of interactions with dark matter subhalos. Furthermore, obtaining a large samples of thin stellar streams can help constrain the orbital structure and hence the potentials of external halos.

We analyze resolved stellar near-infrared photometry of 21 HST fields in M31 to constrain the impact of metallicity on the formation of carbon stars. Observations of nearby galaxies show that the carbon stars are increasingly rare at higher metallicity. Models indicate that carbon star formation efficiency drops due to the decrease in dredge-up efficiency in metal-rich thermally-pulsing Asymptotic Giant Branch (TP-AGB) stars, coupled to a higher initial abundance of oxygen. However, while models predict a metallicity ceiling above which carbon stars cannot form, previous observations have not yet pinpointed this limit. Our new observations reliably separate carbon stars from M-type TP-AGB stars across 2.6-13.7 kpc of M31's metal-rich disk using HST WFC3/IR medium-band filters. We find that the ratio of C to M stars (C/M) decreases more rapidly than extrapolations of observations in more metal-poor galaxies, resulting in a C/M that is too low by more than a factor of 10 in the innermost fields and indicating a dramatic decline in C star formation efficiency at metallicities higher than [M/H] $\approx$ -0.1 dex. The metallicity ceiling remains undetected, but must occur at metallicities higher than what is measured in M31's inner disk ([M/H] $\gtrsim$ +0.06 dex).

As galaxies evolve, they must enrich and exchange gas with the surrounding medium, but the timing of these processes and how much gas is involved remain poorly understood. In this work, we leverage metals as tracers of past gas flows to constrain the history of metal ejection and redistribution in M31. This roughly $L*$ galaxy is a unique case where spatially resolved measurements of the gas-phase and stellar metallicity, dust extinction, and neutral ISM gas content are all available, enabling a census of the present metal mass. We combine spatially resolved star formation histories from the Panchromatic Hubble Andromeda Treasury survey with a metal production model to calculate the history of metal production in M31. We find that $1.8\times10^9 \, M_\odot$ of metals, or 62\% of the metal mass formed within $r < 19 \,\mathrm{kpc}$, is missing from the disk in our fiducial model, implying that the M31 disk has experienced significant gaseous outflows over its lifetime. Under a conservative range of model assumptions, we find that between 3\% and 88\% of metals have been lost ($1.9\times10^7 - 6.4\times10^9 \, M_\odot$), so metals are missing even when all model parameters are chosen to favor metal retention. We show that the missing metal mass could be harbored in M31's CGM if the majority of the metals reside in a hot gas phase. Finally, we find that some metal mass produced in the past 1.5 Gyr in the central $\sim5 \,\mathrm{kpc}$ has likely been redistributed to larger radii within the disk.

We determine the ages of the young, resolved stellar populations at the locations of 237 optically-identified supernova remnants in M83. These age distributions put constraints on the progenitor masses of the supernovae that produced 199 of the remnants. The other 38 show no evidence for having a young progenitor and are therefore good Type Ia SNR candidates. Starting from Hubble Space Telescope broadband imaging, we measured resolved stellar photometry of seven archival WFC3/UVIS fields in F336W, F438W, and F814W. We generate color-magnitude diagrams of the stars within 50 pc of each SNR and fit them with stellar evolution models to obtain the population ages. From these ages we infer the progenitor mass that corresponds to the lifetime of the most prominent age that is $<$50 Myr. In this sample, there are 47 SNRs with best-fit progenitor masses $>$15 M$_{\odot}$, and 5 of these are $>$15 M$_{\odot}$ at 84% confidence. This is the largest collection of high-mass progenitors to date, including our highest-mass progenitor inference found so far, with a constraint of $<$8 Myr. Overall, the distribution of progenitor masses has a power-law index of $-3.0^{+0.2}_{-0.7}$, steeper than Salpeter initial mass function ($-2.35$). It remains unclear whether the reason for the low number of high-mass progenitors is due to the difficulty of finding and measuring such objects or because only a fraction of very massive stars produce supernovae.

This is a white paper submitted in response to the call from the Astro2020 Decadal Survey Committee. We outline the scientific progress that will be made in the next few decades in the study of supernova remnants in the X-ray band, using observatories like Athena, Lynx, and AXIS.

This white paper is motivated by open questions in star formation, which can be uniquely addressed by high resolution X-ray imaging and require an X-ray observatory with large collecting area along good spectral resolution. A complete census of star-forming regions in X-rays, combined with well matched infrared (IR) data, will advance our understanding of disk survival times and dissipation mechanisms. In addition, we will be able to directly observe the effects of X-ray irradiation on circumstellar grain growth to compare with grain evolution models in both high- and low-UV environments. X-rays are native to stars at all phases of star formation and affect planet-forming disks especially through flares. Moreover, X-rays trace magnetic fields which weave through the flares, providing a unique, non-gravitational feedback mechanism between disk and star. Finally, the bright X-ray emission emanating from hot plasma associated with massive stars can have large scale impacts on the topology of star-forming regions and their interface with the interstellar medium (ISM).

The Great Andromeda Galaxy (M31) is the nexus of the near-far galaxy evolution connection and a principal data point for near-field cosmology. Due to its proximity (780 kpc), M31 can be resolved into individual stars like the Milky Way (MW). Unlike the MW, we have the advantage of a global view of M31, enabling M31 to be observed with techniques that also apply to more distant galaxies. Moreover, recent evidence suggests that M31 may have survived a major merger within the last several Gyr, shaping the morphology of its stellar halo and triggering a starburst, while leaving the stellar disk largely intact. The MW and M31 thus provide complementary opportunities for in-depth studies of the disks, halos, and satellites of L* galaxies. Our understanding of the M31 system will be transformed in the 2020s if they include wide field facilities for both photometry (HST-like sensitivity and resolution) and spectroscopy (10-m class telescope, >1 sq. deg. field, highly multiplexed, R~ 3000 to 6000). We focus here on the power of these facilities to constrain the past, present, and future merger history of M31, via chemo-dynamical analyses and star formation histories of phase-mixed stars accreted at early times, as well as stars in surviving tidal debris features, M31's extended disk, and intact satellite galaxies that will eventually be tidally incorporated into the halo. This will yield an unprecedented view of the hierarchical formation of the M31 system and the subhalos that built it into the L* galaxy we observe today.

What are the most important conditions and processes governing the growth of stellar-origin compact objects? The identification of compact object type as either black hole (BH) or neutron star (NS) is fundamental to understanding their formation and evolution. To date, time-domain determination of compact object type remains a relatively untapped tool. Measurement of orbital periods, pulsations, and bursts will lead to a revolution in the study of the demographics of NS and BH populations, linking source phenomena to accretion and galaxy parameters (e.g., star formation, metallicity). To perform these measurements over sufficient parameter space, a combination of a wide-field (>5000 deg^2) transient X-ray monitor over a dynamic energy range (~1-100 keV) and an X-ray telescope for deep surveys with <5 arcsec PSF half-energy width (HEW) angular resolution are required. Synergy with multiwavelength data for characterizing the underlying stellar population will transform our understanding of the time domain properties of transient sources, helping to explain details of supernova explosions and gravitational wave event rates.

Supernova remnants (SNRs) offer the means to study SN explosions, dynamics, and shocks at sub-parsec scales. X-ray observations probe the hot metals synthesized in the explosion and the TeV electrons accelerated by the shocks, and thus they are key to test recent, high-fidelity three-dimensional SN simulations. In this white paper, we discuss the major advances possible with X-ray spectro-imaging at arcsecond scales, with a few eV spectral resolution and a large effective area. These capabilities would revolutionize SN science, offering a three-dimensional view of metals synthesized in explosions and enabling population studies of SNRs in Local Group galaxies. Moreover, this future X-ray mission could detect faint, narrow synchrotron filaments and shock precursors that will constrain the diffusive shock acceleration process.

Supernova remnants (SNRs) offer the means to study SN explosions, dynamics, and shocks at sub-parsec scales. X-ray observations probe the hot metals synthesized in the explosion and the TeV electrons accelerated by the shocks, and thus they are key to test recent, high-fidelity three-dimensional SN simulations. In this white paper, we discuss the major advances possible with X-ray spectro-imaging at arcsecond scales, with a few eV spectral resolution and a large effective area. These capabilities would revolutionize SN science, offering a three-dimensional view of metals synthesized in explosions and enabling population studies of SNRs in Local Group galaxies. Moreover, this future X-ray mission could detect faint, narrow synchrotron filaments and shock precursors that will constrain the diffusive shock acceleration process.

Much of the baryonic matter in the Universe, including the most active and luminous sources, are best studied in the X-ray band. Key advances in X-ray optics and detectors have paved the way for the Advanced X-ray Imaging Satellite (AXIS), a Probe-class mission that is a major improvement over Chandra, which has generated a steady stream of important discoveries for the past 2 decades. AXIS can be launched in the late 2020s and will transform our understanding in several major areas of astrophysics, including the growth and fueling of supermassive black holes, galaxy formation and evolution, the microphysics of cosmic plasmas, the time-variable universe, and a wide variety of cutting-edge studies. Relative to Chandra, the AXIS PSF is nearly twice as sharp on-axis; its field of view for subarcsecond imaging 70 times larger by area; its effective area at 1 keV is 10 times larger. The low-Earth orbit ensures a low and stable detector background, resulting in 50 times greater sensitivity than Chandra for extended sources. AXIS has a rapid repointing response with operations similar to Swift, but is 100 times more sensitive for time-domain science. These capabilities open up a vast discovery space and complement the next generation of astronomical observatories. A high-spectral-resolution mission (Athena) operating at the same time as a high-angular-resolution mission (AXIS) greatly increases the range of scientific discovery. AXIS will use lightweight X-ray optics made of thin single-crystal silicon mirrors developed at NASA Goddard. The detector array builds on a long legacy of X-ray CCD and provides improved photon localization, much faster readout time, and broader energy band. The estimated mission costs are consistent with the $1B Probe mission cost guideline.

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution.

The stellar initial mass function (IMF) is a fundamental property of star formation, offering key insight into the physics driving the process as well as informing our understanding of stellar populations, their by-products, and their impact on the surrounding medium. While the IMF appears to be fairly uniform in the Milky Way disk, it is not yet known how the IMF might behave across a wide range of environments, such as those with extreme gas temperatures and densities, high pressures, and low metallicities. We discuss new opportunities for measuring the IMF in such environments in the coming decade with JWST, WFIRST, and thirty-meter class telescopes. For the first time, we will be able to measure the high-mass slope and peak of the IMF via direct star counts for massive star clusters across the Milky Way and Local Group, providing stringent constraints for star formation theory and laying the groundwork for understanding distant and unresolved stellar systems.

Clumpy structures are a common feature in X-ray images of young Type Ia supernova remnants (SNRs). Although the precise origin of such clumps remains unclear there are three generic possibilities: clumpiness imposed during the explosion, hydrodynamic instabilities that act during the remnant's evolution, and pre-existing structures in the ambient medium. In this article we focus on discriminating between clumping distributions that arise from the explosion and those from the remnant's evolution using existing 3D hydrodynamical simulations. We utilize the genus statistic for this discrimination, applying it to the simulations and {\it Chandra} X-ray observations of the well-known SN Ia remnant of SN 1572 (Tycho's SNR). The genus curve of Tycho's SNR strongly indicates a skewed non-Gaussian distribution of the ejecta clumps and is similar to the genus curve for the simulation with initially clumped ejecta. In contrast, the simulation of perfectly smooth ejecta where clumping arises from the action of hydrodynamic instabilities produced a genus curve that is similar to a random Gaussian field, but disagrees strongly with the genus curve of the observed image. Our results support a scenario in which the observed structure of SN Ia remnants arises from initial clumpiness in the explosion.

The Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.

Metal Evolution, Transport, and Abundance in the LMC (METAL) is a large Cycle 24 program on the Hubble Space Telescope aimed at measuring dust extinction properties and interstellar depletions in the Large Magellanic Cloud (LMC) at half-solar metallicity. The 101-orbit program is comprised of COS and STIS spectroscopy toward 33 LMC massive stars between 1150 A and 3180 A, and parallel WFC3 imaging in 7 NUV-NIR filters. The fraction of silicon in the gas-phase (depletion) obtained from the spectroscopy decreases with increasing hydrogen column density. Depletion patterns for silicon differ between the Milky Way, LMC, and SMC, with the silicon depletion level offsetting almost exactly the metallicity differences, leading to constant gas-phase abundances in those galaxies for a given hydrogen column density. The silicon depletion correlates linearly with the absolute-to-selective extinction, R$_V$, indicating a link between gas depletion and dust grain size. Extinction maps are derived from the resolved stellar photometry in the parallel imaging, which can be compared to FIR images from Herschel and Spitzer to estimate the emissivity of dust at LMC metallicity. The full METAL sample of depletions, UV extinction curves, and extinction maps will inform the abundance, size, composition, and optical properties of dust grains in the LMC, comprehensively improve our understanding of dust properties, and the accuracy with which dust-based gas masses, star formation rates and histories in nearby and high-redshift galaxies are estimated. This overview paper describes the goals, design, data reduction, and initial results of the METAL survey.

The Wide Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) is a high-contrast imager and integral field spectrograph that will enable the study of exoplanets and circumstellar disks at visible wavelengths. Ground-based high-contrast instrumentation has fundamentally limited performance at small working angles, even under optimistic assumptions for 30m-class telescopes. There is a strong scientific driver for better performance, particularly at visible wavelengths. Future flagship mission concepts aim to image Earth analogues with visible light flux ratios of more than 10^10. CGI is a critical intermediate step toward that goal, with a predicted 10^8-9 flux ratio capability in the visible. CGI achieves this through improvements over current ground and space systems in several areas: (i) Hardware: space-qualified (TRL9) deformable mirrors, detectors, and coronagraphs, (ii) Algorithms: wavefront sensing and control; post-processing of integral field spectrograph, polarimetric, and extended object data, and (iii) Validation of telescope and instrument models at high accuracy and precision. This white paper, submitted to the 2018 NAS Exoplanet Science Strategy call, describes the status of key CGI technologies and presents ways in which performance is likely to evolve as the CGI design matures.

We have compiled the most complete census of High-Mass X-ray Binaries (HMXBs) in the Small Magellanic Cloud with the aim to investigate the formation efficiency of young accreting binaries in its low metallicity environment. In total, we use 127 X-ray sources with detections in our \chandra X-ray Visionary Program (XVP), supplemented by 14 additional (likely and confirmed) HMXBs identified by \cite{2016A&A...586A..81H} that fall within the XVP area, but are not either detected in our survey (9 sources) or matched with any XVP source that has at least one OB counterpart in the OGLE-III catalog (5 sources). Specifically, we examine the number ratio of the HMXBs [N(HMXBs)] to {\it (a)} the number of OB stars, {\it (b)} the local star-formation rate (SFR), and {\it (c)} the stellar mass produced during the specific star-formation burst, all as a function of the age of their parent stellar populations. Each of these indicators serves a different role, but in all cases we find that the HMXB formation efficiency increases as a function of time (following a burst of star formation) up to $\sim$40--60\,Myr, and then gradually decreases. The peak formation efficiency N(HMXB)/SFR is (49 $\pm$ 14) $[10^{-5}~{\rm M_{\odot}/yr}]^{-1}$, in good agreement with previous estimates of the average formation efficiency in the broad $\sim$20--60\,Myr age range. The frequency of HMXBs is a factor of 8$\times$ higher than at $\sim$10\,Myr, and 4$\times$ higher than at $\sim$260\,Myr, i.e. at earlier and later epochs, respectively.

We report new Chandra observations of one of the few Galactic supernova remnants whose X-ray spectrum is dominated by nonthermal synchrotron radiation, G330.2+1.0. We find that between 2006 and 2017, some parts of the shell have expanded by about 1%, giving a free-expansion (undecelerated) age of about 1000 yr, and implying shock velocities there of 9000 km/s for a distance of 5 kpc. Somewhat slower expansion is seen elsewhere around the remnant periphery, in particular in compact knots. Because some deceleration must have taken place, we infer that G330.2+1.0 is less than about 1000 yr old. Thus, G330.2+1.0 is one of only four Galactic core-collapse remnants of the last millennium. The large size, low brightness, and young age require a very low ambient density, suggesting expansion in a stellar-wind bubble. We suggest that in the east, where some thermal emission is seen and expansion velocities are much slower, the shock has reached the edge of the cavity. The high shock velocities can easily accelerate relativistic electrons to X-ray-emitting energies. A few small regions show highly significant brightness changes by 10% to 20%, both brightening and fading, a phenomenon previously observed in only two supernova remnants, indicating strong and/or turbulent magnetic fields.

X-ray observations of supernova remnants (SNRs) in the last decade have shown that the presence of recombining plasmas is somewhat common in a certain type of object. The SNR W49B is the youngest, hottest, and most highly ionized among such objects and hence provides crucial information about how the recombination phase is reached during the early evolutionary phase of SNRs. In particular, spectral properties of radiative recombination continuum (RRC) from Fe are the key for constraining the detailed plasma conditions. Here we present imaging and spectral studies of W49B with Nuclear Spectroscopic Telescope Array (NuSTAR), utilizing the highest-ever sensitivity to the Fe RRC at > 8.8keV. We confirm that the Fe RRC is the most prominent at the western part of the SNR because of the lowest electron temperature (~ 1.2 keV) achieved there. Our spatially-resolved spectral analysis reveals a positive correlation between the electron temperature and the recombination timescale with a uniform initial temperature of ~ 4 keV, which is consistent with the rapid adiabatic cooling scenario as an origin of the overionization. This work demonstrates NuSTAR's suitability for studies of thermal emission, in addition to hard nonthermal X-rays, from young and middle-aged SNRs.

We analyze the kinematics of Andromeda's disk as a function of stellar age by using photometry from the Panchromatic Hubble Andromeda Treasury (PHAT) survey and spectroscopy from the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) survey. We use HI 21-cm and CO ($\rm J=1 \rightarrow 0$) data to examine the difference between the deprojected rotation velocity of the gas and that of the stars. We divide the stars into four stellar age bins, from shortest lived to longest lived: massive main sequence stars (0.03 Gyr), more luminous intermediate mass asymptotic giant branch (AGB) stars (0.4 Gyr), less luminous intermediate mass AGB stars (2 Gyr), and low mass red giant branch stars (4 Gyr). There is a clear correlation between the offset of the stellar and the gas rotation velocity, or the asymmetric drift: the longer lived populations lag farther behind the gas than short lived populations. We also examine possible causes of the substructure in the rotation curves and find that the most significant cause of scatter in the rotation curves comes from the tilted ring model being an imperfect way to account for the multiple warps in Andromeda's disk.

We present nearly simultaneous NuSTAR and XMM-Newton observations of the nearby (832 kpc) ultraluminous X-ray source (ULX) M33 X-8. M33 X-8 has a 0.3-10 keV luminosity of LX ~ 1.4 x 10^39 erg/s, near the boundary of the "ultraluminous" classification, making it an important source for understanding the link between typical Galactic X-ray binaries and ULXs. Past studies have shown that the 0.3-10 keV spectrum of X-8 can be characterized using an advection-dominated accretion disk model. We find that when fitting to our NuSTAR and XMM-Newton observations, an additional high-energy (>10 keV) Comptonization component is required, which allows us to rule out single advection-dominated disk and classical sub-Eddington models. With our new constraints, we analyze XMM-Newton data taken over the last 17 years to show that small (~30%) variations in the 0.3-10 keV flux of M33 X-8 result in spectral changes similar to those observed for other ULXs. The two most likely phenomenological scenarios suggested by the data are degenerate in terms of constraining the nature of the accreting compact object (i.e., black hole versus neutron star). We further present a search for pulsations using our suite of data; however, no clear pulsations are detected. Future observations designed to observe M33 X-8 at different flux levels across the full 0.3-30 keV range would significantly improve our constraints on the nature of this important source.

We report on NuSTAR observations of the mixed morphology supernova remnant (SNR) W49B, focusing on its nonthermal emission. Whereas radio observations as well as recent gamma-ray observations evidenced particle acceleration in this SNR, nonthermal X-ray emission has not been reported so far. With the unprecedented sensitivity of NuSTAR in the hard X-ray band, we detect a significant power-law-like component extending up to $\sim 20~{\rm keV}$, most probably of nonthermal origin. The newly discovered component has a photon index of $\Gamma =1.4^{+1.0}_{-1.1}$ with an energy flux between 10 and 20 keV of $(3.3 \pm 0.7) \times 10^{-13}~{\rm erg}~{\rm cm}^{-2}~{\rm s}^{-1}$. The emission mechanism is discussed based on the NuSTAR data combined with those in other wavelengths in the literature. The NuSTAR data, in terms both of the spectral slope and of the flux, are best interpreted as nonthermal electron bremsstrahlung. If this scenario is the case, then the NuSTAR emission provides a new probe to sub-relativistic particles accelerated in the SNR.

We present the results from the Hitomi Soft Gamma-ray Detector (SGD) observation of the Crab nebula. The main part of SGD is a Compton camera, which in addition to being a spectrometer, is capable of measuring polarization of gamma-ray photons. The Crab nebula is one of the brightest X-ray / gamma-ray sources on the sky, and, the only source from which polarized X-ray photons have been detected. SGD observed the Crab nebula during the initial test observation phase of Hitomi. We performed the data analysis of the SGD observation, the SGD background estimation and the SGD Monte Carlo simulations, and, successfully detected polarized gamma-ray emission from the Crab nebula with only about 5 ks exposure time. The obtained polarization fraction of the phase-integrated Crab emission (sum of pulsar and nebula emissions) is (22.1 $\pm$ 10.6)% and, the polarization angle is 110.7$^o$ + 13.2 / $-$13.0$^o$ in the energy range of 60--160 keV (The errors correspond to the 1 sigma deviation). The confidence level of the polarization detection was 99.3%. The polarization angle measured by SGD is about one sigma deviation with the projected spin axis of the pulsar, 124.0$^o$ $\pm$0.1$^o$.

We carried out new observations of two fields in the northern ring of M31 with XMM-Newton with two exposures of 100 ks each and obtained a complete list of X-ray sources down to a sensitivity limit of ~7 x 10^34 erg s^-1 (0.5 - 2.0 keV). The major objective of the observing programme was the study of the hot phase of the ISM in M31. The analysis of the diffuse emission and the study of the ISM is presented in a separate paper. We analysed the spectral properties of all detected sources using hardness ratios and spectra if the statistics were high enough. We also checked for variability. We cross-correlated the source list with the source catalogue of a new survey of the northern disc of M31 carried out with Chandra and Hubble (Panchromatic Hubble Andromeda Treasury, PHAT) as well as with other existing catalogues. We detected a total of 389 sources, including 43 foreground stars and candidates and 50 background sources. Based on the comparison to the Chandra/PHAT survey, we classify 24 hard X-ray sources as new candidates for X-ray binaries (XRBs). In total, we identified 34 XRBs and candidates and 18 supernova remnants (SNRs) and candidates. Three of the four brightest SNRs show emission mainly below 2 keV, consistent with shocked ISM. The spectra of two of them also require an additional component with a higher temperature. The SNR [SPH11] 1535 has a harder spectrum and might suggest that there is a pulsar-wind nebula inside the SNR. We find five new sources showing clear time variability. We also studied the spectral properties of the transient source SWIFT J004420.1+413702, which shows significant variation in flux over a period of seven months (June 2015 to January 2016) and associated change in absorption. Based on the likely optical counterpart detected in the Chandra/PHAT survey, the source is classified as a low-mass X-ray binary.

We present a second epoch of {\it Chandra} observations of the Type Ia LMC SNR 0509-68.7 (N103B) obtained in 2017. When combined with the earlier observations from 1999, we have a 17.4-year baseline with which we can search for evidence of the remnant's expansion. Although the lack of strong point source detections makes absolute image alignment at the necessary accuracy impossible, we can measure the change in the diameter and the area of the remnant, and find that it has expanded by an average velocity of 4170 (2860, 5450) km s$^{-1}$. This supports the picture of this being a young remnant; this expansion velocity corresponds to an undecelerated age of 850 yr, making the real age somewhat younger, consistent with results from light echo studies. Previous infrared observations have revealed high densities in the western half of the remnant, likely from circumstellar material, so it is likely that the real expansion velocity is lower on that side of the remnant and higher on the eastern side. A similar scenario is seen in Kepler's SNR. N103B joins the rare class of Magellanic Cloud SNRs with measured proper motions.

The failed supernova N6946-BH1 likely formed a black hole (BH); we age-date the surrounding population and infer an age and initial mass for the progenitor of this BH formation candidate. First, we use archival Hubble Space Telescope imaging to extract broadband photometry of the resolved stellar populations surrounding this event. Using this photometry, we fit stellar evolution models to the color-magnitude diagrams to measure the recent star formation history (SFH). Modeling the photometry requires an accurate distance; therefore, we measure the tip of the red giant branch (TRGB) and infer a distance modulus of $29.47 \pm 0.079$ to NGC~6946, or a metric distance of $7.83 \pm 0.29$ Mpc. To estimate the stellar population's age, we convert the SFH and uncertainties into a probabilistic distribution for the progenitor's age. The region in the immediate vicinity of N6946-BH1 exhibits the youngest and most vigorous star formation for several hundred pc. This suggests that the progenitor is not a runaway star. From these measurements, we infer an age for the BH progenitor of $10.6^{+14.5}_{-5.9}$ Myr. Assuming that the progenitor evolved effectively as a single star, this corresponds to an initial mass of $17.9^{+29.9}_{-7.6}$ $M_{\odot}$. Previous spectral energy distribution (SED) modeling of the progenitor suggests a mass of $\sim$27 $M_{\odot}$. Formally, the SED-derived mass falls within our narrowest 68\% confidence interval; however, $91\%$ of the probability distribtuion function we measure lies below that mass, putting some tension between the age and the direct-imaging results.

The X-ray source populations within galaxies are typically difficult to identify and classify from X-ray data alone. We are able to break through this barrier by combining deep new Chandra ACIS-I observations with extensive Hubble Space Telescope imaging from the PHAT survey of the M31 disk. We detect 373 X-ray sources down to 0.35-8.0 keV flux of 10$^{-15}$ erg cm$^{-2}$ s$^{-1}$ over 0.4 square degrees, 170 of which are reported for the first time. We identify optical counterpart candidates for 188 of the 373 sources, after using the HST data to correct the absolute astrometry of our Chandra imaging to 0.1$"$. While 58 of these 188 are associated with point sources potentially in M31, over half (107) of the counterpart candidates are extended background galaxies, 5 are star clusters, 12 are foreground stars, and 6 are supernova remnants. Sources with no clear counterpart candidate are most likely to be undetected background galaxies and low-mass X-ray binaries in M31. The 58 point sources that are not consistent with foreground stars are bright enough that they could be high mass stars in M31; however, all but 8 have optical colors inconsistent with single stars, suggesting that many could be background galaxies or binary counterparts. For point-like counterparts, we examine the star formation history of the surrounding stellar populations to look for a young component that could be associated with a high mass X-ray binary. For the 40 point-like counterpart candidates associated with young populations, we find that their age distribution has two peaks at 15-20 Myr and 40-50 Myr. If we only consider the 8 counterpart candidates with typical high-mass main sequence optical star colors, their age distribution peaks mimic those of the sample of 40. Finally, we find that intrinsic faintness, and not extinction, is the main limitation for finding further counterpart candidates.

Using the star formation histories (SFHs) near 94 supernova remnants (SNRs), we infer the progenitor mass distribution for core-collapse supernovae. We use Bayesian inference and model each SFH with multiple bursts of star formation (SF), one of which is assumed to be associated with the SNR. Assuming single-star evolution, the minimum mass of CCSNe is $7.33^{+0.02}_{-0.16}$ $\text{M}_\odot$, the slope of the progenitor mass distribution is $\alpha = -2.96^{+0.45}_{-0.25}$, and the maximum mass is greater than $\text{M}_\textrm{max} > 59$ $\text{M}_\odot$ with a 68% confidence. While these results are consistent with previous inferences, they also provide tighter constraints. The progenitor distribution is somewhat steeper than a Salpeter initial mass function ($\alpha$ = -2.35). This suggests that either SNR catalogs are biased against the youngest SF regions, or the most massive stars do not explode as easily as lower mass stars. If SNR catalogs are biased, it will most likely affect the slope but not the minimum mass. The uncertainties are dominated by three primary sources of uncertainty, the SFH resolution, the number of SF bursts, and the uncertainty on SF rate in each age bin. We address the first two of these uncertainties, with an emphasis on multiple bursts. The third will be addressed in future work.

Nearby galaxy surveys have long classified X-ray binaries (XRBs) by the mass category of their donor stars (high-mass and low-mass). The NuSTAR observatory, which provides imaging data at E $>10$ keV, has enabled the classification of extragalactic XRBs by their compact object type: neutron star (NS) or black hole (BH). We analyzed NuSTAR/Chandra/XMM-Newton observations from a NuSTAR-selected sample of 12 galaxies within 5 Mpc having stellar masses ($M_{\star}$) $10^{7-11}$ $M_{\odot}$ and star formation rates (SFR) $\approx0.01-15$ $M_{\odot}$ yr$^{-1}$. We detect 128 NuSTAR sources to a sensitivity of $\approx10^{38}$ erg s$^{-1}$. Using NuSTAR color-intensity and color-color diagrams we classify 43 of these sources as candidate NS and 47 as candidate BH. We further subdivide BH by accretion states (soft, intermediate, and hard) and NS by weak (Z/Atoll) and strong (accreting pulsar) magnetic field. Using 8 normal (Milky Way-type) galaxies in the sample, we confirm the relation between SFR and galaxy X-ray point source luminosity in the 4-25 and 12-25 keV energy bands. We also constrain galaxy X-ray point source luminosity using the relation $L_{\rm{X}}=\alpha M_{\star}+\beta\text{SFR}$, finding agreement with previous work. The XLF of all sources in the 4-25 and 12-25 keV energy bands matches with the $\alpha=1.6$ slope for high-mass XRBs. We find that NS XLFs suggest a decline beginning at the Eddington limit for a 1.4 $M_{\odot}$ NS, whereas the BH fraction shows an approximate monotonic increase in the 4-25 and 12-25keV energy bands. We calculate the overall ratio of BH to NS to be $\approx1$ for 4-25 keV and $\approx2$ for 12-25 keV.

Double detonations in double white dwarf (WD) binaries undergoing unstable mass transfer have emerged in recent years as one of the most promising Type Ia supernova (SN Ia) progenitor scenarios. One potential outcome of this "dynamically driven double-degenerate double-detonation" (D^6) scenario is that the companion WD survives the explosion and is flung away with a velocity equal to its > 1000 km/s pre-SN orbital velocity. We perform a search for these hypervelocity runaway WDs using Gaia's second data release. In this paper, we discuss seven candidates followed up with ground-based instruments. Three sources are likely to be some of the fastest known stars in the Milky Way, with total Galactocentric velocities between 1000 and 3000 km/s, and are consistent with having previously been companion WDs in pre-SN Ia systems. However, although the radial velocity of one of the stars is > 1000 km/s, the radial velocities of the other two stars are puzzlingly consistent with 0. The combined five-parameter astrometric solutions from Gaia and radial velocities from follow-up spectra yield tentative 6D confirmation of the D^6 scenario. The past position of one of these stars places it within a faint, old SN remnant, further strengthening the interpretation of these candidates as hypervelocity runaways from binary systems that underwent SNe Ia.

We have identified 55 candidate high-mass X-ray binaries (HMXBs) in M33 using available archival {\it HST} and {\it Chandra} imaging to find blue stars associated with X-ray positions. We use the {\it HST} photometric data to model the color-magnitude diagrams in the vicinity of each candidate HMXB to measure a resolved recent star formation history (SFH), and thus a formation timescale, or age for the source. Taken together, the SFHs for all candidate HMXBs in M33 yield an age distribution that suggests preferred formation timescales for HMXBs in M33 of $<$ 5 Myr and $\sim$ 40 Myr after the initial star formation episode. The population at 40 Myr is seen in other Local Group galaxies, and can be attributed to a peak in formation efficiency of HMXBs with neutron stars as compact objects and B star secondary companions. This timescale is preferred as neutron stars should form in abundance from $\sim$ 8 M$_{\odot}$ core-collapse progenitors on these timescales, and B stars are shown observationally to be most actively losing mass around this time. The young population at $<$ 5 Myr has not be observed in other Local Group HMXB population studies, but may be attributed to a population of very massive progenitors forming black holes very early on. We discuss these results in the context of massive binary evolution, and the implications for compact object binaries and gravitational wave sources.

We present 15 high mass X-ray binary (HMXB) candidates in the disk of M31 for which we are able to infer compact object type, spectral type of the donor star, and age using multiwavelength observations from NuSTAR, Chandra, and the Hubble Space Telescope (HST). The hard X-ray colors and luminosities from NuSTAR permit the tentative classification of accreting X-ray binary systems by compact object type, distinguishing black hole from neutron star systems. We find hard state black holes, pulsars, and non-magnetized neutron stars associated with optical point source counterparts with similar frequency. We also find nine non-magnetized neutron stars coincident with globular clusters and an equal number of pulsars with and without point source optical counterparts. We perform spectral energy distribution (SED) fitting for the most likely optical counterparts to the HMXB candidates, finding 7 likely high mass stars and 1 possible red Helium burning star. The remaining 7 HMXB optical counterparts have poor SED fits, so their companion stars remain unclassified. Using published star formation histories, we find that the majority of HMXB candidates --- X-ray sources with UV-bright point source optical counterpart candidates --- are found in regions with star formation bursts less than 50 Myr ago, with 3 associated with young stellar ages (<10 Myr). This is consistent with similar studies of HMXB populations in the Magellanic Clouds, M33, NGC 300, and NGC 2403.

SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$\mu$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$\sigma$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx are available at http://spherex.caltech.edu. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. Though these three themes are undoubtedly compelling, they are far from exhausting the scientific output of SPHEREx. Indeed, SPHEREx would create a unique all-sky spectral database including spectra of very large numbers of astronomical and solar system targets, including both extended and diffuse sources. These spectra would enable a wide variety of investigations, and the SPHEREx team is dedicated to making the data available to the community to enable these investigations, which we refer to as Legacy Science. To that end, we have sponsored two workshops for the general scientific community to identify the most interesting Legacy Science themes and to ensure that the SPHEREx data products are responsive to their needs. In February of 2016, some 50 scientists from all fields met in Pasadena to develop these themes and to understand their implications for the SPHEREx mission. The 2016 workshop highlighted many synergies between SPHEREx and other contemporaneous astronomical missions, facilities, and databases. Consequently, in January 2018 we convened a second workshop at the Center for Astrophysics in Cambridge to focus specifically on these synergies. This white paper reports on the results of the 2018 SPHEREx workshop.

The Wide Field Infrared Survey Telescope (WFIRST) was the top ranked large space mission in the 2010 New Worlds, New Horizons decadal survey, and it was formed by merging the science programs of 3 different mission concepts, including the Microlensing Planet Finder (MPF) concept (Bennett \etal\ 2010). The WFIRST science program (Spergel \etal\ 2015) consists of a general observer program, a wavefront controlled technology program, and two targeted science programs: a program to study dark energy, and a statistical census of exoplanets with a microlensing survey, which uses nearly one quarter of WFIRST's observing time in the current design reference mission. The New Worlds, New Horizons (decadal survey) midterm assessment summarizes the science case for the WFIRST exoplanet microlensing survey with this statement: "WFIRST's microlensing census of planets beyond 1 AU will perfectly complement Kepler's census of compact systems, and WFIRST will also be able to detect free-floating planets unbound from their parent stars\rlap."

We age-date the stellar populations associated with 12 historic nearby core-collapse supernovae (CCSNe) and 2 supernova impostors, and from these ages, we infer their initial masses and associated uncertainties. To do this, we have obtained new HST imaging covering these CCSNe. Using these images, we measure resolved stellar photometry for the stars surrounding the locations of the SNe. We then fit the color-magnitude distributions of this photometry with stellar evolution models to determine the ages of any young existing populations present. From these age distributions, we infer the most likely progenitor mass for all of the SNe in our sample. We find ages between 4 and 50 Myr, corresponding to masses from 7.5 to 59 solar masses. There were no SNe that lacked a young population within 50~pc. Our sample contains 4 type Ib/c SNe; their masses have a wide range of values, suggesting that the progenitors of stripped-envelope SNe are binary systems. Both impostors have masses constrained to be $\lesssim$7.5 solar masses. In cases with precursor imaging measurements, we find that age-dating and precursor imaging give consistent progenitor masses. This consistency implies that, although the uncertainties for each technique are significantly different, the results of both are reliable to the measured uncertainties. We combine these new measurements with those from our previous work and find that the distribution of 25 core-collapse SNe progenitor masses is consistent with a standard Salpeter power-law mass function, no upper mass cutoff, and an assumed minimum mass for core-collapse of 7.5~M$_{\odot}$.

Since its discovery in 2008, the Andromeda galaxy nova M31N 2008-12a has been observed in eruption every single year. This unprecedented frequency indicates an extreme object, with a massive white dwarf and a high accretion rate, which is the most promising candidate for the single-degenerate progenitor of a type-Ia supernova known to date. The previous three eruptions of M31N 2008-12a have displayed remarkably homogeneous multi-wavelength properties: (i) From a faint peak, the optical light curve declined rapidly by two magnitudes in less than two days; (ii) Early spectra showed initial high velocities that slowed down significantly within days and displayed clear He/N lines throughout; (iii) The supersoft X-ray source (SSS) phase of the nova began extremely early, six days after eruption, and only lasted for about two weeks. In contrast, the peculiar 2016 eruption was clearly different. Here we report (i) the considerable delay in the 2016 eruption date, (ii) the significantly shorter SSS phase, and (iii) the brighter optical peak magnitude (with a hitherto unobserved cusp shape). Early theoretical models suggest that these three different effects can be consistently understood as caused by a lower quiescence mass-accretion rate. The corresponding higher ignition mass caused a brighter peak in the free-free emission model. The less-massive accretion disk experienced greater disruption, consequently delaying re-establishment of effective accretion. Without the early refueling, the SSS phase was shortened. Observing the next few eruptions will determine whether the properties of the 2016 outburst make it a genuine outlier in the evolution of M31N 2008-12a.

We discuss the technical challenges we faced and the techniques we used to overcome them when reducing the PHAT photometric data set on the Amazon Elastic Compute Cloud (EC2). We first describe the architecture of our photometry pipeline, which we found particularly efficient for reducing the data in multiple ways for different purposes. We then describe the features of EC2 that make this architecture both efficient to use and challenging to implement. We describe the techniques we adopted to process our data, and suggest ways these techniques may be improved for those interested in trying such reductions in the future. Finally, we summarize the output photometry data products, which are now hosted publicly in two places in two formats. They are in simple fits tables in the high-level science products on MAST, and on a queryable database available through the NOAO Data Lab.

We present results from the Hitomi X-ray observation of a young composite-type supernova remnant (SNR) G21.5$-$0.9, whose emission is dominated by the pulsar wind nebula (PWN) contribution. The X-ray spectra in the 0.8-80 keV range obtained with the Soft X-ray Spectrometer (SXS), Soft X-ray Imager (SXI) and Hard X-ray Imager (HXI) show a significant break in the continuum as previously found with the NuSTAR observation. After taking into account all known emissions from the SNR other than the PWN itself, we find that the Hitomi spectra can be fitted with a broken power law with photon indices of $\Gamma_1=1.74\pm0.02$ and $\Gamma_2=2.14\pm0.01$ below and above the break at $7.1\pm0.3$ keV, which is significantly lower than the NuSTAR result ($\sim9.0$ keV). The spectral break cannot be reproduced by time-dependent particle injection one-zone spectral energy distribution models, which strongly indicates that a more complex emission model is needed, as suggested by recent theoretical models. We also search for narrow emission or absorption lines with the SXS, and perform a timing analysis of PSR J1833$-$1034 with the HXI and SGD. No significant pulsation is found from the pulsar. However, unexpectedly, narrow absorption line features are detected in the SXS data at 4.2345 keV and 9.296 keV with a significance of 3.65 $\sigma$. While the origin of these features is not understood, their mere detection opens up a new field of research and was only possible with the high resolution, sensitivity and ability to measure extended sources provided by an X-ray microcalorimeter.

We present moderately deep (125 ks) {\it XMM-Newton} observations of supernova remnant G330.2$+$1.0. This remnant is one of only a few known that fall into "synchrotron-dominated" category, with the emission almost entirely dominated by a nonthermal continuum. Previous X-ray observations could only characterize the spectra of a few regions. Here, we examine the spectra from fourteen regions surrounding the entire rim, finding that the spectral properties of the nonthermal emission do not vary significantly in any systematic way from one part of the forward shock to another, unlike several other remnants of this class. We confirm earlier findings that the power-law index, $\Gamma$, ranges from about 2.1-2.5, while the absorbing column density is generally between 2.0-2.6 $\times 10^{22}$ cm$^{-2}$. Fits with the {\it srcut} model find values of the roll-off frequency in the range of 10$^{17.1} - 10^{17.5}$ Hz, implying energies of accelerated electrons of $\sim 100$ TeV. These values imply a high shock velocity of $\sim 4600$ km s$^{-1}$, favoring a young age of the remnant. Diffuse emission from the interior is nonthermal in origin as well, and fits to these regions yield similar values to those along the rim, also implying a young age. Thermal emission is present in the east, and the spectrum is consistent with a $\sim 650$ km s$^{-1}$ shock wave encountering interstellar or circumstellar material with a density of $\sim 1$ cm$^{-3}$.

The present paper investigates the temperature structure of the X-ray emitting plasma in the core of the Perseus cluster using the 1.8--20.0 keV data obtained with the Soft X-ray Spectrometer (SXS) onboard the Hitomi Observatory. A series of four observations were carried out, with a total effective exposure time of 338 ks and covering a central region $\sim7'$ in diameter. The SXS was operated with an energy resolution of $\sim$5 eV (full width at half maximum) at 5.9 keV. Not only fine structures of K-shell lines in He-like ions but also transitions from higher principal quantum numbers are clearly resolved from Si through Fe. This enables us to perform temperature diagnostics using the line ratios of Si, S, Ar, Ca, and Fe, and to provide the first direct measurement of the excitation temperature and ionization temperature in the Perseus cluster. The observed spectrum is roughly reproduced by a single temperature thermal plasma model in collisional ionization equilibrium, but detailed line ratio diagnostics reveal slight deviations from this approximation. In particular, the data exhibit an apparent trend of increasing ionization temperature with increasing atomic mass, as well as small differences between the ionization and excitation temperatures for Fe, the only element for which both temperatures can be measured. The best-fit two-temperature models suggest a combination of 3 and 5 keV gas, which is consistent with the idea that the observed small deviations from a single temperature approximation are due to the effects of projection of the known radial temperature gradient in the cluster core along the line of sight. Comparison with the Chandra/ACIS and the XMM-Newton/RGS results on the other hand suggests that additional lower-temperature components are present in the ICM but not detectable by Hitomi SXS given its 1.8--20 keV energy band.

The Hitomi SXS spectrum of the Perseus cluster, with $\sim$5 eV resolution in the 2-9 keV band, offers an unprecedented benchmark of the atomic modeling and database for hot collisional plasmas. It reveals both successes and challenges of the current atomic codes. The latest versions of AtomDB/APEC (3.0.8), SPEX (3.03.00), and CHIANTI (8.0) all provide reasonable fits to the broad-band spectrum, and are in close agreement on best-fit temperature, emission measure, and abundances of a few elements such as Ni. For the Fe abundance, the APEC and SPEX measurements differ by 16%, which is 17 times higher than the statistical uncertainty. This is mostly attributed to the differences in adopted collisional excitation and dielectronic recombination rates of the strongest emission lines. We further investigate and compare the sensitivity of the derived physical parameters to the astrophysical source modeling and instrumental effects. The Hitomi results show that an accurate atomic code is as important as the astrophysical modeling and instrumental calibration aspects. Substantial updates of atomic databases and targeted laboratory measurements are needed to get the current codes ready for the data from the next Hitomi-level mission.

We present Hitomi observations of N132D, a young, X-ray bright, O-rich core-collapse supernova remnant in the Large Magellanic Cloud (LMC). Despite a very short observation of only 3.7 ks, the Soft X-ray Spectrometer (SXS) easily detects the line complexes of highly ionized S K and Fe K with 16-17 counts in each. The Fe feature is measured for the first time at high spectral resolution. Based on the plausible assumption that the Fe K emission is dominated by He-like ions, we find that the material responsible for this Fe emission is highly redshifted at ~800 km/s compared to the local LMC interstellar medium (ISM), with a 90% credible interval of 50-1500 km/s if a weakly informative prior is placed on possible line broadening. This indicates (1) that the Fe emission arises from the supernova ejecta, and (2) that these ejecta are highly asymmetric, since no blue-shifted component is found. The S K velocity is consistent with the local LMC ISM, and is likely from swept-up ISM material. These results are consistent with spatial mapping that shows the He-like Fe concentrated in the interior of the remnant and the S tracing the outer shell. The results also show that even with a very small number of counts, direct velocity measurements from Doppler-shifted lines detected in extended objects like supernova remnants are now possible. Thanks to the very low SXS background of ~1 event per spectral resolution element per 100 ks, such results are obtainable during short pointed or slew observations with similar instruments. This highlights the power of high-spectral-resolution imaging observations, and demonstrates the new window that has been opened with Hitomi and will be greatly widened with future missions such as the X-ray Astronomy Recovery Mission (XARM) and Athena.

We report a Hitomi observation of IGR J16318-4848, a high-mass X-ray binary system with an extremely strong absorption of N_H~10^{24} cm^{-2}. Previous X-ray studies revealed that its spectrum is dominated by strong fluorescence lines of Fe as well as continuum emission. For physical and geometrical insight into the nature of the reprocessing material, we utilize the high spectroscopic resolving power of the X-ray microcalorimeter (the soft X-ray spectrometer; SXS) and the wide-band sensitivity by the soft and hard X-ray imager (SXI and HXI) aboard Hitomi. Even though photon counts are limited due to unintended off-axis pointing, the SXS spectrum resolves Fe K{\alpha_1} and K{\alpha_2} lines and puts strong constraints on the line centroid and width. The line width corresponds to the velocity of 160^{+300}_{-70} km s^{-1}. This represents the most accurate, and smallest, width measurement of this line made so far from any X-ray binary, much less than the Doppler broadening and shift expected from speeds which are characteristic of similar systems. Combined with the K-shell edge energy measured by the SXI and HXI spectra, the ionization state of Fe is estimated to be in the range of Fe I--IV. Considering the estimated ionization parameter and the distance between the X-ray source and the absorber, the density and thickness of the materials are estimated. The extraordinarily strong absorption and the absence of a Compton shoulder component is confirmed. These characteristics suggest reprocessing materials which are distributed in a narrow solid angle or scattering primarily with warm free electrons or neutral hydrogen.

The origin of the narrow Fe-K{\alpha} fluorescence line at 6.4 keV from active galactic nuclei has long been under debate; some of the possible sites are the outer accretion disk, the broad line region, a molecular torus, or interstellar/intracluster media. In February-March 2016, we performed the first X-ray microcalorimeter spectroscopy with the Soft X-ray Spectrometer (SXS) onboard the Hitomi satellite of the Fanaroff-Riley type I radio galaxy NGC 1275 at the center of the Perseus cluster of galaxies. With the high energy resolution of ~5 eV at 6 keV achieved by Hitomi/SXS, we detected the Fe-K{\alpha} line with ~5.4 {\sigma} significance. The velocity width is constrained to be 500-1600 km s$^{-1}$ (FWHM for Gaussian models) at 90% confidence. The SXS also constrains the continuum level from the NGC 1275 nucleus up to ~20 keV, giving an equivalent width ~20 eV of the 6.4 keV line. Because the velocity width is narrower than that of broad H{\alpha} line of ~2750 km s$^{-1}$, we can exclude a large contribution to the line flux from the accretion disk and the broad line region. Furthermore, we performed pixel map analyses on the Hitomi/SXS data and image analyses on the Chandra archival data, and revealed that the Fe-K{\alpha} line comes from a region within ~1.6 kpc from the NGC 1275 core, where an active galactic nucleus emission dominates, rather than that from intracluster media. Therefore, we suggest that the source of the Fe-K{\alpha} line from NGC 1275 is likely a low-covering fraction molecular torus or a rotating molecular disk which probably extends from a pc to hundreds pc scale in the active galactic nucleus system.

We review the observations of dust emission in supernova rem- nants (SNRs) and supernovae (SNe). Theoretical calculations suggest that SNe, particularly core-collapse, should make significant quantities of dust, perhaps as much as a solar mass. Observations of extragalactic SNe have yet to find anywhere near this amount, but this may be the result of observa- tional limitations. SN 1987A, in the process of transitioning from a SN to an SNR, does show signs of a significant amount of dust forming in its ejecta, but whether this dust will survive the passage of the reverse shock to be injected into the ISM is unknown. IR observations of SNRs have not turned up significant quantities of dust, and the dust that is observed is generally swept-up by the forward shock, rather than created in the ejecta. Because the shock waves also destroy dust in the ISM, we explore the question of whether SNe might be net destroyers, rather than net creators of dust in the universe.

Extending the earlier measurements reported in Hitomi collaboration (2016, Nature, 535, 117), we examine the atmospheric gas motions within the central 100~kpc of the Perseus cluster using observations obtained with the Hitomi satellite. After correcting for the point spread function of the telescope and using optically thin emission lines, we find that the line-of-sight velocity dispersion of the hot gas is remarkably low and mostly uniform. The velocity dispersion reaches maxima of approximately 200~km~s$^{-1}$ toward the central active galactic nucleus (AGN) and toward the AGN inflated north-western `ghost' bubble. Elsewhere within the observed region, the velocity dispersion appears constant around 100~km~s$^{-1}$. We also detect a velocity gradient with a 100~km~s$^{-1}$ amplitude across the cluster core, consistent with large-scale sloshing of the core gas. If the observed gas motions are isotropic, the kinetic pressure support is less than 10\% of the thermal pressure support in the cluster core. The well-resolved optically thin emission lines have Gaussian shapes, indicating that the turbulent driving scale is likely below 100~kpc, which is consistent with the size of the AGN jet inflated bubbles. We also report the first measurement of the ion temperature in the intracluster medium, which we find to be consistent with the electron temperature. In addition, we present a new measurement of the redshift to the brightest cluster galaxy NGC~1275.

Thanks to its high spectral resolution (~5 eV at 6 keV), the Soft X-ray Spectrometer (SXS) on board Hitomi enables us to measure the detailed structure of spatially resolved emission lines from highly ionized ions in galaxy clusters for the first time. In this series of papers, using the SXS we have measured the velocities of gas motions, metallicities and the multi-temperature structure of the gas in the core of the Perseus cluster. Here, we show that when inferring physical properties from line emissivities in systems like Perseus, the resonant scattering (RS) effect should be taken into account. In the Hitomi waveband, RS mostly affects the FeXXV He$\alpha$ line ($w$) - the strongest line in the spectrum. The flux measured by Hitomi in this line is suppressed by a factor ~1.3 in the inner ~30 kpc, compared to predictions for an optically thin plasma; the suppression decreases with the distance from the center. The $w$ line also appears slightly broader than other lines from the same ion. The observed distortions of the $w$ line flux, shape and distance dependence are all consistent with the expected effect of the resonant scattering in the Perseus core. By measuring the ratio of fluxes in optically thick ($w$) and thin (FeXXV forbidden, He$\beta$, Ly$\alpha$) lines, and comparing these ratios with predictions from Monte Carlo radiative transfer simulations, the velocities of gas motions have been obtained. The results are consistent with the direct measurements of gas velocities from line broadening described elsewhere in this series, although the systematic and statistical uncertainties remain significant. Further improvements in the predictions of line emissivities in plasma models, and deeper observations with future X-ray missions will enable RS measurements to provide powerful constraints on the amplitude and anisotropy of clusters gas motions.

The 2008 discovery of an eruption of M31N 2008-12a began a journey on which the true nature of this remarkable recurrent nova continues to be revealed. M31N 2008-12a contains a white dwarf close to the Chandrasekhar limit, accreting at a high rate from its companion, and undergoes thermonuclear eruptions which are observed yearly and may even be twice as frequent. In this paper, we report on Hubble Space Telescope Space Telescope Imaging Spectrograph ultraviolet spectroscopy taken within days of the predicted 2015 eruption, coupled with Keck spectroscopy of the 2013 eruption. Together, this spectroscopy permits the reddening to be constrained to E(B-V) = 0.10 +/- 0.03. The UV spectroscopy reveals evidence for highly ionized, structured, and high velocity ejecta at early times. No evidence for neon is seen in these spectra however, but it may be that little insight can be gained regarding the composition of the white dwarf (CO vs ONe).

The recurrent nova M31N 2008-12a experiences annual eruptions, contains a near-Chandrasekhar mass white dwarf, and has the largest mass accretion rate in any nova system. In this paper, we present Hubble Space Telescope (HST) WFC3/UVIS photometry of the late decline of the 2015 eruption. We couple these new data with archival HST observations of the quiescent system and Keck spectroscopy of the 2014 eruption. The late-time photometry reveals a rapid decline to a minimum luminosity state, before a possible recovery / re-brightening in the run-up to the next eruption. Comparison with accretion disk models supports the survival of the accretion disk during the eruptions, and uncovers a quiescent disk mass accretion rate of the order of $10^{-6}\,M_\odot\,\mathrm{yr}^{-1}$, which may rise beyond $10^{-5}\,M_\odot\,\mathrm{yr}^{-1}$ during the super-soft source phase - both of which could be problematic for a number of well-established nova eruption models. Such large accretion rates, close to the Eddington limit, might be expected to be accompanied by additional mass loss from the disk through a wind and even collimated outflows. The archival HST observations, combined with the disk modeling, provide the first constraints on the mass donor; $L_\mathrm{donor}=103^{+12}_{-11}\,L_\odot$, $R_\mathrm{donor}=14.14^{+0.46}_{-0.47}\,R_\odot$, and $T_\mathrm{eff, donor}=4890\pm110$ K, which may be consistent with an irradiated M31 red-clump star. Such a donor would require a system orbital period $\gtrsim5$ days. Our updated analysis predicts that the M31N 2008-12a WD could reach the Chandrasekhar mass in < 20 kyr.

We use Hubble Space Telescope (HST) observations of red clump stars taken as part of the Small Magellanic Cloud Investigation of Dust and Gas Evolution (SMIDGE) program to measure the average dust extinction curve in a ~ 200 pc x 100 pc region in the southwest bar of the Small Magellanic Cloud (SMC). The rich information provided by our 8-band ultra-violet through near-infrared photometry allows us to model the color-magnitude diagram of the red clump accounting for the extinction curve shape, a log-normal distribution of $A_{V}$, and the depth of the stellar distribution along the line of sight. We measure an extinction curve with $R_{475} = A_{475}/(A_{475}-A_{814})$ = 2.65 $\pm$ 0.11. This measurement is significantly larger than the equivalent values of published Milky Way $R_{V}$ = 3.1 ($R_{475} = 1.83$) and SMC Bar $R_{V}$ = 2.74 ($R_{475} = 1.86$) extinction curves. Similar extinction curve offsets in the Large Magellanic Cloud (LMC) have been interpreted as the effect of large dust grains. We demonstrate that the line-of-sight depth of the SMC (and LMC) introduces an apparent "gray" contribution to the extinction curve inferred from the morphology of the red clump. We show that no gray dust component is needed to explain extinction curve measurements when a full-width half-max depth of 10 $\pm$ 2 kpc in the stellar distribution of the SMC (5 $\pm$ 1 kpc for the LMC) is considered, which agrees with recent studies of Magellanic Cloud stellar structure. The results of our work demonstrate the power of broad-band HST imaging for simultaneously constraining dust and galactic structure outside the Milky Way.

We map the star formation history across M31 by fitting stellar evolution models to color-magnitude diagrams of each 83$''{\times}$83$''$ (0.3$\times$1.4 kpc, deprojected) region of the PHAT survey outside of the innermost 6$'{\times}$12$'$ portion. We find that most of the star formation occurred prior to $\sim$8 Gyr ago, followed by a relatively quiescent period until $\sim$4 Gyr ago, a subsequent star formation episode about 2 Gyr ago and a return to relative quiescence. There appears to be little, if any, structure visible for populations with ages older than 2 Gyr, suggesting significant mixing since that epoch. Finally, assuming a Kroupa IMF from 0.1$-$100 M$_{\odot}$, we find that the total amount of star formation over the past 14 Gyr in the area over which we have fit models is 5${\times}$10$^{10}$ M$_{\odot}$. Fitting the radial distribution of this star formation and assuming azimuthal symmetry, (1.5$\pm$0.2)${\times}$10$^{11}$ M$_{\odot}$ of stars have formed in the M31 disk as a whole, (9$\pm$2)${\times}$10$^{10}$ M$_{\odot}$ of which has likely survived to the present after accounting for evolutionary effects. This mass is about one fifth of the total dynamical mass of M31.

To search for giant X-ray pulses correlated with the giant radio pulses (GRPs) from the Crab pulsar, we performed a simultaneous observation of the Crab pulsar with the X-ray satellite Hitomi in the 2 -- 300 keV band and the Kashima NICT radio observatory in the 1.4 -- 1.7 GHz band with a net exposure of about 2 ks on 25 March 2016, just before the loss of the Hitomi mission.The timing performance of the Hitomi instruments was confirmed to meet the timing requirement and about 1,000 and 100 GRPs were simultaneously observed at the main and inter-pulse phases, respectively, and we found no apparent correlation between the giant radio pulses and the X-ray emission in either the main or inter-pulse phases.All variations are within the 2 sigma fluctuations of the X-ray fluxes at the pulse peaks, and the 3 sigma upper limits of variations of main- or inter- pulse GRPs are 22\% or 80\% of the peak flux in a 0.20 phase width, respectively, in the 2 -- 300 keV band.The values become 25\% or 110\% for main or inter-pulse GRPs, respectively, when the phase width is restricted into the 0.03 phase.Among the upper limits from the Hitomi satellite, those in the 4.5-10 keV and the 70-300 keV are obtained for the first time, and those in other bands are consistent with previous reports.Numerically, the upper limits of main- and inter-pulse GRPs in the 0.20 phase width are about (2.4 and 9.3) $\times 10^{-11}$ erg cm$^{-2}$, respectively. No significant variability in pulse profiles implies that the GRPs originated from a local place within the magnetosphere and the number of photon-emitting particles temporally increases.However, the results do not statistically rule out variations correlated with the GRPs, because the possible X-ray enhancement may appear due to a $>0.02$\% brightening of the pulse-peak flux under such conditions.

We have carried out a study of the X-ray properties of the supernova remnant (SNR) population in M33 with XMM-Newton, comprising deep observations of 8 fields in M33 covering all of the area within the D$_{25}$ contours, and with a typical luminosity of 7.1$\times$10$^{34}$ erg s$^{-1}$ (0.2-2.0 keV) . Here we report our work to characterize the X-ray properties of the previously identified SNRs in M33, as well as our search for new X-ray detected SNRs. With our deep observations and large field of view we have detected 105 SNRs at the 3$\sigma$ level, of which 54 SNRs are newly detected in X-rays, and three are newly discovered SNRs. Combining XMM-Newton data with deep Chandra survey data allows detailed spectral fitting of 15 SNRs, for which we have measured temperatures, ionization timescales, and individual abundances. This large sample of SNRs allows us to construct an X-ray luminosity function, and compare its shape to luminosity functions from host galaxies of differing metallicities and star formation rates to look for environmental effects on SNR properties. We conclude that while metallicity may play a role in SNR population characteristics, differing star formation histories on short timescales, and small-scale environmental effects appear to cause more significant differences between X-ray luminosity distributions. In addition, we analyze the X-ray detectability of SNRs, and find that in M33 SNRs with higher [SII]/H$\alpha$ ratios, as well as those with smaller galactocentric distances, are more detectable in X-rays.

We present results of integral field spectroscopy of Balmer-dominated shocks in the LMC supernova remnant (SNR) N103B, carried out using the Wide Field Integral Spectrograph (WiFeS) on the 2.3 m telescope at Siding Spring Observatory in Australia. Existing X-ray studies of N103B have indicated a Type Ia supernova origin. Radiative shock emission from clumpy material surrounding the SNR may result from interaction of the forward shock with relic stellar wind material, possibly implicating a thermonuclear explosion in a single-degenerate binary system. The recently discovered Balmer-dominated shocks mark the impact of the forward shock with low density, partially neutral CSM gas, and form a partial shell encircling clumps of material exhibiting radiative shocks. The WiFeS spectra of N103B cover the range 3500 - 7000 Angstroms and reveal broad Halpha profiles up to 2360 km/s in width along the northern rim, and both Halpha and Hbeta broad profiles approximately 1300 km/s in width along the southern rim. Our fits to the Halpha line profiles indicate that in addition to from the usual broad and narrow emission components, a third component of intermediate width exists in the Balmer-dominated shocks, ranging from around 125 km/s up to 225 km/s in width. Recent Balmer-dominated shock models have predicted such an intermediate-width component will be generated in a fast neutral precursor. We derive a Sedov age of approximately 685 years for N103B from the Balmer-dominated spectra, consistent with the young age of 380 - 860 years estimated from light echo studies.

We used four-year baseline HST/WFC3 IR observations of the Galactic Centre in the F153M band (1.53 micron) to identify variable stars in the central ~2.3'x2.3' field. We classified 3845 long-term (periods from months to years) and 76 short-term (periods of a few days or less) variables among a total sample of 33070 stars. For 36 of the latter ones, we also derived their periods (<3 days). Our catalog not only confirms bright long period variables and massive eclipsing binaries identified in previous works, but also contains many newly recognized dim variable stars. For example, we found \delta Scuti and RR Lyrae stars towards the Galactic Centre for the first time, as well as one BL Her star (period < 1.3 d). We cross-correlated our catalog with previous spectroscopic studies and found that 319 variables have well-defined stellar types, such as Wolf-Rayet, OB main sequence, supergiants and asymptotic giant branch stars. We used colours and magnitudes to infer the probable variable types for those stars without accurately measured periods or spectroscopic information. We conclude that the majority of unclassified variables could potentially be eclipsing/ellipsoidal binaries and Type II Cepheids. Our source catalog will be valuable for future studies aimed at constraining the distance, star formation history and massive binary fraction of the Milky Way nuclear star cluster.

Because of strong and spatially highly variable interstellar extinction and extreme source crowding, the faint (K>15) stellar population in the Milky Way's nuclear cluster is still poorly studied. RR Lyrae stars provide us with a tool to estimate the mass of the oldest, relative dim stellar population. Recently, we analyzed HST/WFC3/IR observations of the central 2.3'x2.3' of the Milky Way and found 21 variable stars with periods between 0.2 and 1d. Here, we present a further comprehensive analysis of these stars. The period-luminosity relationship of RR Lyrae is used to derive their extinctions and distances. Using multiple approaches, we classify our sample as four RRc, four RRab and three candidates, ten binaries. Especially, the four RRabs show sawtooth light curves and fall exactly onto the Oosterhoff I division in the Bailey diagram. Compared to the RRabs reported by Minniti et al, 2016, our new RRabs have higher extinction (A_K>1.8) and should be closer to the Galactic Centre. The extinction and distance of one RRab match those for the nuclear star cluster given in previous works. We perform simulations and find that after correcting for incompleteness, there could be no more than 40 RRabs within the nuclear star cluster and in our field-of-view. Through comparing with the known globular clusters of the Milky Way, we estimate that if there exists an old, metal-poor (-1.5<[Fe/H]<-1) stellar population in the Milky Way nuclear star cluster on a scale of 5x5pc, then it contributes at most 4.7x10^5 solar mass, i.e., ~18% of the stellar mass.

The Crab nebula originated from a core-collapse supernova (SN) explosion observed in 1054 A.D. When viewed as a supernova remnant (SNR), it has an anomalously low observed ejecta mass and kinetic energy for an Fe-core collapse SN. Intensive searches were made for a massive shell that solves this discrepancy, but none has been detected. An alternative idea is that the SN1054 is an electron-capture (EC) explosion with a lower explosion energy by an order of magnitude than Fe-core collapse SNe. In the X-rays, imaging searches were performed for the plasma emission from the shell in the Crab outskirts to set a stringent upper limit to the X-ray emitting mass. However, the extreme brightness of the source hampers access to its vicinity. We thus employed spectroscopic technique using the X-ray micro-calorimeter onboard the Hitomi satellite. By exploiting its superb energy resolution, we set an upper limit for emission or absorption features from yet undetected thermal plasma in the 2-12 keV range. We also re-evaluated the existing Chandra and XMM-Newton data. By assembling these results, a new upper limit was obtained for the X-ray plasma mass of <~ 1Mo for a wide range of assumed shell radius, size, and plasma temperature both in and out of the collisional equilibrium. To compare with the observation, we further performed hydrodynamic simulations of the Crab SNR for two SN models (Fe-core versus EC) under two SN environments (uniform ISM versus progenitor wind). We found that the observed mass limit can be compatible with both SN models if the SN environment has a low density of <~ 0.03 cm-3 (Fe core) or <~ 0.1 cm-3 (EC) for the uniform density, or a progenitor wind density somewhat less than that provided by a mass loss rate of 10-5 Mo yr-1 at 20 km s-1 for the wind environment.

We present the first three-dimensional measurements of the velocity of various ejecta knots in Tycho's supernova remnant, known to result from a Type Ia explosion. Chandra X-ray observations over a 12-year baseline from 2003 to 2015 allow us to measure the proper motion of nearly 60 "tufts" of Si-rich ejecta, giving us the velocity in the plane of the sky. For the line of sight velocity, we use two different methods: a non-equilibrium ionization model fit to the strong Si and S lines in the 1.2-2.8 keV regime, and a fit consisting of a series of Gaussian lines. These methods give consistent results, allowing us to determine the red or blue shift of each of the knots. Assuming a distance of 3.5 kpc, we find total velocities that range from 2400 to 6600 km s$^{-1}$, with a mean of 4430 km s$^{-1}$. We find several regions where the ejecta knots have overtaken the forward shock. These regions have proper motions in excess of 6000 km s$^{-1}$. Some Type Ia supernova explosion models predict a velocity asymmetry in the ejecta. We find no such velocity asymmetries in Tycho, and discuss our findings in light of various explosion models, favoring those delayed detonation models with relatively vigorous and symmetrical deflagrations. Finally, we compare measurements with models of the remnant's evolution that include both smooth and clumpy ejecta profiles, finding that both ejecta profiles can be accommodated by the observations.

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.

The requirements-driven OSIRIS-REx Camera Suite (OCAMS) acquires images essential to collecting a sample from the surface of Bennu. During proximity operations, these images document the presence of satellites and plumes, record spin state, enable an accurate digital terrain model of the shape of the asteroid and identify any surface hazards. They confirm the presence of sampleable regolith on the surface, observe the sampling event itself, and image the sample head in order to verify its readiness to be stowed. They document the history of Bennu as an example of early solar system material, as a microgravity body with a planetesimal size-scale, and as a carbonaceous object. OCAMS is fitted with three cameras. The MapCam records point-source color images on approach to the asteroid in order to connect ground-based point-source observations of Bennu to later higher-resolution surface spectral imaging. The SamCam documents the sample site before, during, and after it is disturbed by the sample mechanism. The PolyCam, using its focus mechanism, observes the sample site at sub-centimeter resolutions, revealing surface texture and morphology. While their imaging requirements divide naturally between the three cameras, they preserve a strong degree of functional overlap. OCAMS and the other spacecraft instruments allow the OSIRIS-REx mission to collect a sample from a microgravity body on the same visit during which it was first optically acquired from long range, a useful capability as humanity explores near-Earth, Main-Belt and Jupiter Trojan asteroids.

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

We report the identification of a bright hard X-ray source dominating the M31 bulge above 25 keV from a simultaneous NuSTAR-Swift observation. We find that this source is the counterpart to Swift J0042.6+4112, which was previously detected in the Swift BAT All-sky Hard X-ray Survey. This Swift BAT source had been suggested to be the combined emission from a number of point sources; our new observations have identified a single X-ray source from 0.5 to 50 keV as the counterpart for the first time. In the 0.5-10 keV band, the source had been classified as an X-ray binary candidate in various Chandra and XMM studies; however, since it was not clearly associated with Swift J0042.6+4112, the previous E < 10 keV observations did not generate much attention. This source has a spectrum with a soft X-ray excess (kT~ 0.2 keV) plus a hard spectrum with a power law of Gamma ~ 1 and a cutoff around 15-20 keV, typical of the spectral characteristics of accreting pulsars. Unfortunately, any potential pulsation was undetected in the NuSTAR data, possibly due to insufficient photon statistics. The existing deep HST images exclude high-mass (>3 Msun) donors at the location of this source. The best interpretation for the nature of this source is an X-ray pulsar with an intermediate-mass (<3 Msun) companion or a symbiotic X-ray binary. We discuss other possibilities in more detail.

We present the spectral analysis of a 200~ks observation of the broad-line radio galaxy 3C~120 performed with the high energy transmission grating (HETG) spectrometer on board the \emph{Chandra} X-ray Observatory. We find (i) a neutral absorption component intrinsic to the source with column density of $\text{log}N_H = 20.67\pm0.05$~cm$^{-2}$, (ii) no evidence for a warm absorber with an upper limit on the column density of just $\text{log}N_H < 19.7$~cm$^{-2}$ assuming the typical ionization parameter log$\xi$$\simeq$2.5~erg~s$^{-1}$~cm, the warm absorber may instead be replaced by (iii) a hot emitting gas with temperature $kT \simeq 0.7$~keV observed as soft X-ray emission from ionized Fe L-shell lines which may originate from a kpc scale shocked bubble inflated by the AGN wind or jet with a shock velocity of about 1,000~km~s$^{-1}$ determined by the emission line width, (iv) a neutral Fe K$\alpha$ line and accompanying emission lines indicative of a Compton-thick cold reflector with low reflection fraction $R\simeq0.2$, suggesting a large opening angle of the torus, (v) a highly ionized Fe~XXV emission feature indicative of photoionized gas with ionization parameter log$\xi$$=$$3.75^{+0.27}_{-0.38}$~erg~s$^{-1}$~cm and a column density of $\text{log}N_H > 22$~cm$^{-2}$ localized within $\sim$2~pc from the X-ray source, and (vi) possible signatures for a highly ionized disk wind. Together with previous evidence for intense molecular line emission, these results indicate that 3C~120 is likely a late state merger undergoing strong AGN feedback.

High-resolution X-ray spectroscopy with Hitomi was expected to resolve the origin of the faint unidentified E=3.5 keV emission line reported in several low-resolution studies of various massive systems, such as galaxies and clusters, including the Perseus cluster. We have analyzed the Hitomi first-light observation of the Perseus cluster. The emission line expected for Perseus based on the XMM-Newton signal from the large cluster sample under the dark matter decay scenario is too faint to be detectable in the Hitomi data. However, the previously reported 3.5 keV flux from Perseus was anomalously high compared to the sample-based prediction. We find no unidentified line at the reported high flux level. Taking into account the XMM measurement uncertainties for this region, the inconsistency with Hitomi is at a 99% significance for a broad dark-matter line and at 99.7% for a narrow line from the gas. We do not find anomalously high fluxes of the nearby faint K line or the Ar satellite line that were proposed as explanations for the earlier 3.5 keV detections. We do find a hint of a broad excess near the energies of high-n transitions of Sxvi (E=3.44 keV rest-frame) -- a possible signature of charge exchange in the molecular nebula and another proposed explanation for the unidentified line. While its energy is consistent with XMM pn detections, it is unlikely to explain the MOS signal. A confirmation of this interesting feature has to wait for a more sensitive observation with a future calorimeter experiment.

We confirm that the object DDO216-A1 is a substantial globular cluster at the center of Local Group galaxy DDO216 (the Pegasus dwarf irregular), using Hubble Space Telescope ACS imaging. By fitting isochrones, we find the cluster metallicity to be [M/H] = -1.6 +/-0.2, for reddening E(B-V) = 0.16 +/-0.02; the best-fit age is 12.3 +/-0.8 Gyr. There are ~30 RR Lyrae variables in the cluster; the magnitude of the fundamental mode pulsators gives a distance modulus of 24.77 +/-0.08 - identical to the host galaxy. The ratio of overtone to fundamental mode variables and their mean periods make DDO216-A1 an Oosterhoff Type I cluster. We find an I-band central surface brightness 20.85 +/-0.17 F814W mag per square arcsecond, a half-light radius of 3.1 arcsec (13.4 pc), and an absolute magnitude M814 = -7.90 +/-0.16 (approximately 10^5 solar masses). King models fit to the cluster give the core radius and concentration index, r_c = 2.1" +/-0.9" and c = 1.24 +/-0.39. The cluster is an "extended" cluster somewhat typical of some dwarf galaxies and the outer halo of the Milky Way. The cluster is projected <30 pc south of the center of DDO216, unusually central compared to most dwarf galaxy globular clusters. Analytical models of dynamical friction and tidal destruction suggest that it probably formed at a larger distance, up to ~1 kpc, and migrated inward. DDO216 has an unexceptional cluster specific frequency, S_N = 10. DDO216 is the lowest-luminosity Local Group galaxy to host a 10^5 solar mass globular cluster, and the only transition-type (dSph/dIrr) in the Local Group with a globular.

In young starburst galaxies, the X-ray population is expected to be dominated by the relics of the most massive and short-lived stars, black-hole and neutron-star high mass X-ray binaries (XRBs). In the closest such galaxy, IC 10, we have made a multi-wavelength census of these objects. Employing a novel statistical correlation technique, we have matched our list of 110 X-ray point sources, derived from a decade of Chandra observations, against published photometric data. We report an 8 sigma correlation between the celestial coordinates of the two catalogs, with 42 X-ray sources having an optical counterpart. Applying an optical color-magnitude selection to isolate blue supergiant (SG) stars in IC 10, we find 16 matches. Both cases show a statistically significant overabundance versus the expectation value for chance alignments. The blue objects also exhibit systematically higher fx/fv ratios than other stars in the same magnitude range. Blue SG-XRBs include a major class of progenitors of double-degenerate binaries, hence their numbers are an important factor in modeling the rate of gravitational wave sources. We suggest that the anomalous features of the IC 10 stellar population are explained if the age of the IC 10 starburst is close to the time of the peak of interaction for massive binaries.

We have monitored the Cassiopeia dwarf galaxy (IC 10) in a series of 10 Chandra ACIS-S observations to capture its variable and transient X-ray source population, which is expected to be dominated by High Mass X-ray Binaries (HMXBs). We present a sample of 21 X-ray sources that are variable between observations at the 3 sigma level, from a catalog of 110 unique point sources. We find 4 transients (flux variability ratio greater than 10) and a further 8 objects with ratio > 5. The observations span years 2003 - 2010 and reach a limiting luminosity of >10$^{35}$ erg/s, providing sensitivity to X-ray binaries in IC 10 as well as flare stars in the foreground Milky Way. The nature of the variable sources is investigated from light-curves, X-ray spectra, energy quantiles, and optical counterparts. The purpose of this study is to discover the composition of the X-ray binary population in a young starburst environment. IC 10 provides a sharp contrast in stellar population age (<10 My) when compared to the Magellanic Clouds (40-200 My) where most of the known HMXBs reside. We find 10 strong HMXB candidates, 2 probable background Active Galactic Nuclei, 4 foreground flare-stars or active binaries, and 5 not yet classifiable sources. Complete classification of the sample requires optical spectroscopy for radial velocity analysis and deeper X-ray observations to obtain higher S/N spectra and search for pulsations. A catalog and supporting dataset are provided.

We present a sample of 224 stars that emit H$\alpha$ (H$\alpha$ stars) in the Andromeda galaxy (M31). The stars were selected from $\sim$ 5000 spectra, collected as part of the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo survey using Keck II/DEIMOS. We used six-filter Hubble Space Telescope photometry from the Panchromatic Hubble Andromeda Treasury survey to classify and investigate the properties of the H$\alpha$ stars. We identified five distinct categories of H$\alpha$ star: B-type main sequence (MS) stars, `transitioning'-MS (T-MS) stars, red core He burning (RHeB) stars, non-C-rich asymptotic giant branch (AGB) stars, and C-rich AGB stars. We found $\sim$ 12 per cent of B-type stars exhibit H$\alpha$ emission (Be stars). The frequency of Be to all B stars is known to vary with the metallicity of their environment. Comparing this proportion of Be stars with other environments around the Local Group, the result could indicate that M31 is more metal rich than the Milky Way. We predict that the 17 T-MS H$\alpha$ stars are Be stars evolving off the MS with fading H$\alpha$ emission. We separated RHeB from AGB H$\alpha$ stars. We conclude that the 61 RHeB and AGB stars are likely to be Long Period Variables. We found that $\sim$ 14 per cent of C-rich AGB stars (C stars) emit H$\alpha$, which is an upper limit for the ratio of C-rich Miras to C stars. This catalogue of H$\alpha$ stars will be useful to constrain stellar evolutionary models, calibrate distance indicators for intermediate age populations, and investigate the properties of M31.

We have obtained three epochs of Chandra ACIS-I observations (totaling $\sim$184 ks) of the nearby spiral galaxy NGC~300 to study the logN-logS distributions of its X-ray point source population down to $\sim$2$\times$10$^{-15}$ erg s$^{-1}$ cm$^{-2}$ in the 0.35-8 keV band (equivalent to $\sim$10$^{36}$ erg s$^{-1}$). The individual epoch logN-logS distributions are best described as the sum of a background AGN component, a simple power law, and a broken power law, with the shape of the logN-logS distributions sometimes varying between observations. The simple power law and AGN components produce a good fit for "persistent" sources (i.e., with fluxes that remain constant within a factor of $\sim$2). The differential power law index of $\sim$1.2 and high fluxes suggest that the persistent sources intrinsic to NGC~300 are dominated by Roche lobe overflowing low mass X-ray binaries. The variable X-ray sources are described by a broken power law, with a faint-end power law index of $\sim$1.7, a bright-end index of $\sim$2.8-4.9, and a break flux of $\sim$8$\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$ ($\sim$4$\times10^{36}$ erg s$^{-1}$), suggesting they are mostly outbursting, wind-fed high mass X-ray binaries, although the logN-logS distribution of variable sources likely also contains low-mass X-ray binaries. We generate model logN-logS distributions for synthetic X-ray binaries and constrain the distribution of maximum X-ray fluxes attained during outburst. Our observations suggest that the majority of outbursting X-ray binaries occur at sub-Eddington luminosities, where mass transfer likely occurs through direct wind accretion at $\sim$1-3% of the Eddington rate.

We explore the nucleus of the nearby 10$^9M_{\odot}$~early-type galaxy (ETGs), NGC~404, using Hubble Space Telescope (HST)/STIS spectroscopy and WFC3 imaging. We first present evidence for nuclear variability in UV, optical, and infrared filters over a time period of 15~years. This variability adds to the already substantial evidence for an accreting black hole at the center of NGC~404. We then redetermine the dynamical black hole mass in NGC~404 including modeling of the nuclear stellar populations. We combine HST/STIS spectroscopy with WFC3 images to create a local color-$M/L$~relation derived from stellar population modeling of the STIS data. We then use this to create a mass model for the nuclear region. We use Jeans modeling to fit this mass model to adaptive optics (AO) stellar kinematic observations from Gemini/NIFS. From our stellar dynamical modeling, we find a 3$\sigma$ upper limit on the black hole mass of $1.5\times10^5M_{\odot}$. Given the accretion evidence for a black hole, this upper limit makes NGC~404 the lowest mass central black hole with dynamical mass constraints. We find that the kinematics of H$_2$ emission line gas show evidence for non-gravitational motions preventing the use of gas dynamical modeling to constrain the black hole mass. Our stellar population modeling also reveals that the central, counter-rotating region of the nuclear cluster is dominated by $\sim$1~Gyr old populations.

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]

The recent star formation history (SFH) in the outer disk of NGC 300 is presented through the analysis of color magnitude diagrams (CMDs). We analyze resolved stellar photometry by creating CMDs from four Hubble Space Telescope fields containing a combination of images from the Advanced Camera for Surveys and the UVIS imager aboard the Wide Field Camera 3. From the best models of these CMDs, we derive the SFH in order to extract the young stellar component for the past 200 Myrs. We find that the young stellar disk of NGC 300 is unbroken out to at least ${\sim}$8 scale lengths (including an upper limit out to ${\sim}$10 scale lengths) with $r{_s}$ = 1.4 ${\pm}$ 0.1 kpc, which is similar to the total stellar surface brightness profile. This unbroken profile suggests that NGC 300 is undisturbed, similar to the isolated disk galaxy NGC 2403. We compare the environments of NGC 300, NGC 2403, and M33 along with the properties of the gas and stellar disks. We find that the disturbed HI outer disk morphology is not accompanied by a break in the young stellar disk. This may indicate that processes which affect the outer HI morphology may not leave an imprint on the young stellar disk.

The Andromeda Galaxy recurrent nova M31N 2008-12a had been observed in eruption ten times, including yearly eruptions from 2008-2014. With a measured recurrence period of $P_\mathrm{rec}=351\pm13$ days (we believe the true value to be half of this) and a white dwarf very close to the Chandrasekhar limit, M31N 2008-12a has become the leading pre-explosion supernova type Ia progenitor candidate. Following multi-wavelength follow-up observations of the 2013 and 2014 eruptions, we initiated a campaign to ensure early detection of the predicted 2015 eruption, which triggered ambitious ground and space-based follow-up programs. In this paper we present the 2015 detection; visible to near-infrared photometry and visible spectroscopy; and ultraviolet and X-ray observations from the Swift observatory. The LCOGT 2m (Hawaii) discovered the 2015 eruption, estimated to have commenced at Aug. $28.28\pm0.12$ UT. The 2013-2015 eruptions are remarkably similar at all wavelengths. New early spectroscopic observations reveal short-lived emission from material with velocities $\sim13000$ km s$^{-1}$, possibly collimated outflows. Photometric and spectroscopic observations of the eruption provide strong evidence supporting a red giant donor. An apparently stochastic variability during the early super-soft X-ray phase was comparable in amplitude and duration to past eruptions, but the 2013 and 2015 eruptions show evidence of a brief flux dip during this phase. The multi-eruption Swift/XRT spectra show tentative evidence of high-ionization emission lines above a high-temperature continuum. Following Henze et al. (2015a), the updated recurrence period based on all known eruptions is $P_\mathrm{rec}=174\pm10$ d, and we expect the next eruption of M31N 2008-12a to occur around mid-Sep. 2016.

We present far-infrared images of the Galactic oxygen-rich supernova remnant (SNR) G292.0+1.8, acquired with the PACS and SPIRE instruments of the Herschel Space Observatory. We find that the SNR shell is detected in the PACS blue (100 micron) band, but not in the red (160 micron) band, broadly consistent with results from AKARI observations. There is no discernible emission from G292.0+1.8 in SPIRE imagery at 250, 350 and 500 micron. Comparing the 100 micron emission to that observed with Spitzer, at 24 and 70 micron, we find a very similar appearance for G292.0+1.8 at all three wavelengths. The IR emission is dominated by dust from non-radiative circumstellar shocks. In addition, the radiatively shocked O-rich clump known as the 'Spur' on the eastern side of G292.0+1.8, is clearly detected in the PACS blueimages, with marginal detection in the red. Fitting the existing 14-40 micron IRS spectra of the Spur together with photometric measurements from 70 micron MIPS and 100 micron PACS photometry, we place an upper limit of <~ 0.04 solar masses, of ejecta dust mass in the Spur, under the most conservative assumption that the ejecta dust has a temperature of 15 K. Modeling the dust continuum in the IRS spectra at four positions around the rim, we estimate postshock densities ranging from n_p = 3.5 cm^-3 to 11 cm^-3. The integrated spectrum of the entire SNR, dominated by swept up circumstellar dust, can be fit with a two-component dust model with a silicate component at 62 K and graphite component at 40 K for a total dust mass of 0.023 solar masses.

Clusters of galaxies are the most massive gravitationally-bound objects in the Universe and are still forming. They are thus important probes of cosmological parameters and a host of astrophysical processes. Knowledge of the dynamics of the pervasive hot gas, which dominates in mass over stars in a cluster, is a crucial missing ingredient. It can enable new insights into mechanical energy injection by the central supermassive black hole and the use of hydrostatic equilibrium for the determination of cluster masses. X-rays from the core of the Perseus cluster are emitted by the 50 million K diffuse hot plasma filling its gravitational potential well. The Active Galactic Nucleus of the central galaxy NGC1275 is pumping jetted energy into the surrounding intracluster medium, creating buoyant bubbles filled with relativistic plasma. These likely induce motions in the intracluster medium and heat the inner gas preventing runaway radiative cooling; a process known as Active Galactic Nucleus Feedback. Here we report on Hitomi X-ray observations of the Perseus cluster core, which reveal a remarkably quiescent atmosphere where the gas has a line-of-sight velocity dispersion of 164+/-10 km/s in a region 30-60 kpc from the central nucleus. A gradient in the line-of-sight velocity of 150+/-70 km/s is found across the 60 kpc image of the cluster core. Turbulent pressure support in the gas is 4% or less of the thermodynamic pressure, with large scale shear at most doubling that estimate. We infer that total cluster masses determined from hydrostatic equilibrium in the central regions need little correction for turbulent pressure.

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.

Located in the Small Magellanic Cloud (SMC), SXP214 is an X-ray pulsar in a high mass X-ray binary system with a Be-star companion. A recent survey of the SMC under a Chandra X-ray Visionary program found the source in a transition when the X-ray flux was on a steady rise. The Lomb-Scargle periodogram revealed a pulse period of 211.49 +/- 0.42 s, which is significantly (>5sigma) shorter than the previous measurements with XMM-Newton and RXTE. This implies that the system has gone through sudden spin-up episodes recently. The pulse profile shows a sharp eclipse-like feature with a modulation amplitude of >95%. The linear rise of the observed X-ray luminosity from <~2x to 7x10^35 erg s^-1 is correlated with steady softening of the X-ray spectrum, which can be described by the changes in the local absorption from N_H ~ 10^24 to <~10^20 cm^-2 for an absorbed power-law model. The soft X-ray emission below 2 keV was absent in the early part of the observation when only the pulsating hard X-ray component was observed, whereas at later times both soft and hard X-ray components were observed pulsating. A likely explanation is that the neutron star was initially hidden in the circumstellar disk of the companion, and later came out of the disk with the accreted material that continued fueling the observed pulsation.

We present newly obtained X-ray and radio observations of Tycho's supernova remnant using {\it Chandra} and the Karl G. Jansky Very Large Array in 2015 and 2013/14, respectively. When combined with earlier epoch observations by these instruments, we now have time baselines for expansion measurements of the remnant of 12-15 year in the X-rays and 30 year in the radio. The remnant's large angular size allows for proper motion measurements at many locations around the periphery of the blast wave. We find, consistent with earlier measurements, a clear gradient in the expansion velocity of the remnant, despite its round shape. The proper motions on the western and southwestern sides of the remnant are about a factor of two higher than those in the east and northeast. We showed in an earlier work that this is related to an offset of the explosion site from the geometric center of the remnant due to a density gradient in the ISM, and using our refined measurements reported here, we find that this offset is $\sim 23"$ towards the northeast. An explosion center offset in such a circular remnant has implications for searches for progenitor companions in other remnants.

Thermally-Pulsing Asymptotic Giant Branch (TP-AGB) stars are relatively short lived (less than a few Myr), yet their cool effective temperatures, high luminosities, efficient mass-loss and dust production can dramatically effect the chemical enrichment histories and the spectral energy distributions of their host galaxies. The ability to accurately model TP-AGB stars is critical to the interpretation of the integrated light of distant galaxies, especially in redder wavelengths. We continue previous efforts to constrain the evolution and lifetimes of TP-AGB stars by modeling their underlying stellar populations. Using Hubble Space Telescope (HST) optical and near-infrared photometry taken of 12 fields of 10 nearby galaxies imaged via the ACS Nearby Galaxy Survey Treasury and the near-infrared HST/SNAP follow-up campaign, we compare the model and observed TP-AGB luminosity functions as well as the number ratio of TP-AGB to red giant branch stars. We confirm the best-fitting mass-loss prescription, introduced by Rosenfield et al. 2014, in which two different wind regimes are active during the TP-AGB, significantly improves models of many galaxies that show evidence of recent star formation. This study extends previous efforts to constrain TP-AGB lifetimes to metallicities ranging -1.59 < [Fe/H] < -0.56 and initial TP-AGB masses up to ~ 4 Msun, which include TP-AGB stars that undergo hot-bottom burning.

We present the results of a joint Swift-NuSTAR spectroscopy campaign on M31. We focus on the five brightest globular cluster X-ray sources in our fields. Two of these had previously been argued to be black hole candidates on the basis of apparent hard-state spectra at luminosities above those for which neutron stars are in hard states. We show that these two sources are likely to be Z-sources (i.e. low magnetic field neutron stars accreting near their Eddington limits), or perhaps bright atoll sources (low magnetic field neutron stars which are just a bit fainter than this level) on the basis of simultaneous Swift and NuSTAR spectra which cover a broader range of energies. These new observations reveal spectral curvature above 6-8 keV that would be hard to detect without the broader energy coverage the NuSTAR data provide relative to Chandra and XMM-Newton. We show that the other three sources are also likely to be bright neutron star X-ray binaries, rather than black hole X-ray binaries. We discuss why it should already have been realized that it was unlikely that these objects were black holes on the basis of their being persistent sources, and we re-examine past work which suggested that tidal capture products would be persistently bright X-ray emitters. We discuss how this problem is likely due to neglecting disk winds in older work that predict which systems will be persistent and which will be transient.

A precise measurement of shock velocities is crucial for constraining the mechanism and efficiency of cosmic-ray (CR) acceleration at supernova remnant (SNR) shock fronts. The northeastern rim of the SNR RCW 86 is thought to be a particularly efficient CR acceleration site, owing to the recent result in which an extremely high shock velocity of ~6000 km/s was claimed (Helder et al. 2009). Here we revisit the same SNR rim with the Chandra X-ray Observatory, 11 years after the first observation. This longer baseline than previously available allows us to determine a more accurate proper motion of the nonthermal X-ray filament, revealing a much lower velocity of 3000 \pm 340 km/s (and even slower at a brighter region). Although the value has dropped to a half of that from the previous X-ray measurement, it is still higher than the mean velocity of the H-alpha filaments in this region (~1200 km/s). This discrepancy implies that the filaments bright in nonthermal X-rays and H-alpha emission trace different velocity components, and thus a CR pressure constrained by combining the X-ray kinematics and the H-alpha spectroscopy can easily be overestimated. We also measure the proper motion of the thermal X-ray filament immediately to the south of the nonthermal one. The inferred velocity (720 \pm 360 km/s) is significantly lower than that of the nonthermal filament, suggesting the presence of denser ambient material, possibly a wall formed by a wind from the progenitor, which has drastically slowed down the shock.

We present new observations of the "supernova impostor" SN 2010da using the Chandra X-ray Observatory and the Hubble Space Telescope. During the initial 2010 outburst, the 0.3-10 keV luminosity was observed by Swift to be $\sim5\times10^{38}$ erg s$^{-1}$ and faded by a factor of $\sim$25 in a four month period. Our two new Chandra observations show a factor of $\sim$10 increase in the 0.35-8 keV X-ray flux, from $\sim$4$\times10^{36}$ erg s$^{-1}$ to $4\times10^{37}$ erg s$^{-1}$ in $\sim$6 months, and the X-ray spectrum is consistent in both observations with a power law photon index of $\Gamma\sim0$. We find evidence of X-ray spectral state changes: when SN 2010da is in a high-luminosity state, the X-ray spectrum is harder ($\Gamma\sim0$) compared to the low-luminosity state ($\Gamma\sim1.2\pm0.8$). Using our Hubble observations, we fit the color magnitude diagram of the coeval stellar population to estimate a time since formation of the SN 2010da progenitor system of $\lesssim$5 Myr. Our observations are consistent with SN 2010da being a high-mass X-ray binary (HMXB) composed of a neutron star and a luminous blue variable-like companion, although we cannot rule out the possibility that SN 2010da is an unusually X-ray bright massive star. The $\lesssim$5 Myr age is consistent with the theoretically predicted delay time between the formation of a massive binary and the onset of the HMXB phase. It is possible that the initial 2010 outburst marked the beginning of X-ray production in the system, making SN 2010da possibly the first massive progenitor binary ever observed to evolve into an HMXB.

We report on the proper motions of Balmer-dominated filaments in Kepler's supernova remnant using high resolution images obtained with the Hubble Space Telescope at two epochs separated by about 10 years. We use the improved proper motion measurements and revised values of shock velocities to derive a distance to Kepler of 5.1 [+0.8, -0.7] kpc. The main shock around the northern rim of the remnant has a typical speed of 1690 km/s and is encountering material with densities of about 8 cm^-3. We find evidence for the variation of shock properties over small spatial scales, including differences in the driving pressures as the shock wraps around a curved cloud surface. We find that the Balmer filaments ahead of the ejecta knot on the northwest boundary of the remnant are becoming fainter and more diffuse. We also find that the Balmer filaments associated with circumstellar material in the interior regions of the remnant are due to shocks with significantly lower velocities and that the brightness variations among these filaments trace the density distribution of the material, which may have a disk-like geometry.

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.

The Pluto system was recently explored by NASA's New Horizons spacecraft, making closest approach on 14 July 2015. Pluto's surface displays diverse landforms, terrain ages, albedos, colors, and composition gradients. Evidence is found for a water-ice crust, geologically young surface units, surface ice convection, wind streaks, volatile transport, and glacial flow. Pluto's atmosphere is highly extended, with trace hydrocarbons, a global haze layer, and a surface pressure near 10 microbars. Pluto's diverse surface geology and long-term activity raise fundamental questions about how small planets remain active many billions of years after formation. Pluto's large moon Charon displays tectonics and evidence for a heterogeneous crustal composition, its north pole displays puzzling dark terrain. Small satellites Hydra and Nix have higher albedos than expected.

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.

Several young supernova remnants exhibit thin X-ray bright rims of synchrotron radiation at their forward shocks. Thin rims require strong magnetic field amplification beyond simple shock compression if rim widths are only limited by electron energy losses. But, magnetic field damping behind the shock could produce similarly thin rims with less extreme field amplification. Variation of rim width with energy may thus discriminate between competing influences on rim widths. We measured rim widths around Tycho's supernova remnant in 5 energy bands using an archival 750 ks Chandra observation. Rims narrow with increasing energy and are well described by either loss-limited or damped scenarios, so X-ray rim width-energy dependence does not uniquely specify a model. But, radio counterparts to thin rims are not loss-limited and better reflect magnetic field structure. Joint radio and X-ray modeling favors magnetic damping in Tycho's SNR with damping lengths ~1--5% of remnant radius and magnetic field strengths ~50--400 $\mu$G assuming Bohm diffusion. X-ray rim widths are ~1% of remnant radius, somewhat smaller than inferred damping lengths. Electron energy losses are important in all models of X-ray rims, suggesting that the distinction between loss-limited and damped models is blurred in soft X-rays. All loss-limited and damping models require magnetic fields $\gtrsim$ 20 $\mu$G, affirming the necessity of magnetic field amplification beyond simple compression.

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 explore the ratio (C/M) of carbon-rich to oxygen-rich thermally pulsing asymptotic giant branch(TP-AGB) stars in the disk of M31 using a combination of moderate-resolution optical spectroscopy from the Spectroscopic Landscape of Andromeda's Stellar Halo (SPLASH) survey and six-filter Hubble Space Telescope photometry from the Panchromatic Hubble Andromeda Treasury (PHAT) survey.Carbon stars were identified spectroscopically. Oxygen-rich M-stars were identifed using three different photometric definitions designed to mimic, and thus evaluate, selection techniques common in the literature. We calculate the C/M ratio as a function of galactocentric radius, present-day gas-phase oxygen abundance, stellar metallicity, age (via proxy defined as the ratio of TP-AGB stars to red giant branch, RGB, stars), and mean star formation rate over the last 400 Myr. We find statistically significant correlations between log(C/M) and all parameters. These trends are consistent across different M-star selection methods, though the fiducial values change. Of particular note is our observed relationship between log(C/M) and stellar metallicity, which is fully consistent with the trend seen across Local Group satellite galaxies. The fact that this trend persists in stellar populations with very different star formation histories indicates that the C/M ratio is governed by stellar properties alone.

We present comprehensive X-ray point source catalogs of NGC~55, NGC~2403, and NGC~4214 as part of the Chandra Local Volume Survey. The combined archival observations have effective exposure times of 56.5 ks, 190 ks, and 79 ks for NGC~55, NGC~2403, and NGC~4214, respectively. When combined with our published catalogs for NGC 300 and NGC 404, our survey contains 629 X-ray sources total down to a limiting unabsorbed luminosity of $\sim5\times10^{35}$ erg s$^{-1}$ in the 0.35-8 keV band in each of the five galaxies. We present X-ray hardness ratios, spectral analysis, radial source distributions, and an analysis of the temporal variability for the X-ray sources detected at high significance. To constrain the nature of each X-ray source, we carried out cross-correlations with multi-wavelength data sets. We searched overlapping Hubble Space Telescope observations for optical counterparts to our X-ray detections to provide preliminary classifications for each X-ray source as a likely X-ray binary, background AGN, supernova remnant, or foreground star.

Using Hubble Space Telescope Advanced Camera for Surveys (HST/ACS) and Wide Field Camera 3 (WFC3) observations from the Panchromatic Hubble Andromeda Treasury (PHAT), we present new period-luminosity (P-L) relations for Cepheid variables in M31. Cepheids from several ground-based studies are identified in the PHAT photometry to derive new P-L and Wesenheit P-L relations in the near infrared and visual filters. We derive a distance modulus to M31 of 24.51 +/- 0.08 in the IR bands and 24.32 +/- 0.09 in the visual bands, including the first P-L relations in the F475W and F814W filters for M31. Our derived visual and IR distance moduli disagree at slightly more than a 1{\sigma} level. Differences in the P-L relations between ground-based and HST observations are investigated for a subset of Cepheids. We find a significant discrepancy between ground-based and HST P-L relations with the same Cepheids, suggesting adverse effects from photometric contamination in ground-based Cepheid observations. Additionally, a statistically significant radial trend in the P-L relation is found which does not appear to be explained by metallicity.

We present photometric observations from the {\it Stratospheric Observatory for Infrared Astronomy (SOFIA)} at 11.1 $\mu$m of the Type IIn supernova (SN IIn) 2010jl. The SN is undetected by {\it SOFIA}, but the upper limits obtained, combined with new and archival detections from {\it Spitzer} at 3.6 \& 4.5 $\mu$m allow us to characterize the composition of the dust present. Dust in other Type IIn SNe has been shown in previous works to reside in a circumstellar shell of material ejected by the progenitor system in the few millenia prior to explosion. Our model fits show that the dust in the system shows no evidence for the strong, ubiquitous 9.7 $\mu$m feature from silicate dust, suggesting the presence of carbonaceous grains. The observations are best fit with 0.01-0.05 $\msun$ of carbonaceous dust radiating at a temperature of $\sim 550-620$ K. The dust composition may reveal clues concerning the nature of the progenitor system, which remains ambiguous for this subclass. Most of the single star progenitor systems proposed for SNe IIn, such as luminous blue variables, red supergiants, yellow hypergiants, and B[e] stars, all clearly show silicate dust in their pre-SN outflows. However, this post-SN result is consistent with the small sample of SNe IIn with mid-IR observations, none of which show signs of emission from silicate dust in their IR spectra.

NGC~300 X-1 is a Wolf Rayet + black hole binary that exhibits periodic decreases in X-ray flux. We present two new observations of NGC~300 X-1 from the Chandra X-ray Observatory (totaling $\sim$130 ks) along with ACS imaging data from the Hubble Space Telescope. We observe significant short-term variability in the X-ray emission that is inconsistent with an occultation by the donor star, but is consistent with structure in the outer accretion disk or the wind of the donor star. We simultaneously fit a partially-covered disk blackbody and Comptonized corona model to the eclipse egress and non-eclipsing portions of the X-ray spectrum. We find that the only model parameters that varied between the eclipse egress and non-eclipsing portions of the spectra were the partial covering fraction ($\sim$86% during eclipse egress and $\sim$44% during non-eclipse) and absorbing column ($\sim$12.3$\times10^{22}$ cm$^{-2}$ during eclipse egress, compared to $\sim$1.4$\times10^{22}$ cm$^{-2}$ during non-eclipse). The X-ray spectra are consistent with the movement of the X-ray source through the dense stellar winds of the companion star. From our new HST imaging, we find the WR star within the X-ray error circle, along with additional optical sources including an AGB star and an early-type main sequence star. Finally, we use our egress measurement to rephase previous radial velocity measurements reported in the literature, and find evidence that the velocities are strongly affected by the ionization of the wind by the compact object. Thus, we argue the inferred mass of the black hole may not be reliable.

We present results from two {\it Chandra} observations of the 16.6 d X-ray/$\gamma$-ray high-mass binary 1FGL J1018.6$-$5856 located at the center of the supernova remnant G284.3$-$1.8. The binary spectra, separated by 0.25 in binary phase, are fit with an absorbed power-law model with $\Gamma\approx1.7-1.8$ for both observations (the flux during the second observation is a factor of 1.7 smaller). In the high-resolution ACIS-I image we found a hint of extended emission $\approx2''$--$3''$ southeast of the binary, significant at the $3\sigma$ level. Binary evolution codes reproduce the system's observed properties with two massive stars with an initial 18-day period, undergoing mass transfer and leaving behind a heavy $\approx2$ $\msun$ neutron star. The initial mass of the progenitor star in this scenario is 27 $\pm 4$ $\msun$. \emph{Chandra} and {\it XMM-Newton} images of the remnant show it has a relatively low X-ray surface brightness. The two brightest regions of extended X-ray emission, with luminosities $\sim10^{33}$ erg s$^{-1}$ for $d=5$ kpc, lie in the northern and western portions and show significantly different spectra. The northern patch is consistent with shocked ISM, with a low temperature and long ionization timescale. However, the western patch is dominated by ejecta, and shows significantly enhanced Mg content relative to other ejecta products. The abundance ratios inferred resemble those from the Large Magellanic Cloud remnant N49B. To our knowledge, this is only the second case of such Mg-rich ejecta found in an SNR. Nucleosynthesis models for core-collapse SNe predict Mg-rich ejecta from very massive progenitors of $>25$ $\msun$.

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 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]

Recent studies have shown that starburst dwarf galaxies have steeply rising rotation curves in their inner parts, pointing to a close link between the intense star formation and a centrally concentrated mass distribution (baryons and dark matter). More quiescent dwarf irregulars typically have slowly rising rotation curves, although some "compact" irregulars with steep, inner rotation curves exist. We analyze archival Hubble Space Telescope images of two nearby "compact" irregular galaxies (NGC 4190 and NGC 5204), which were selected solely on the basis of their dynamical properties and their proximity. We derive their recent star-formation histories by fitting color-magnitude diagrams of resolved stellar populations, and find that the star-formation properties of both galaxies are consistent with those of known starburst dwarfs. Despite the small sample, this strongly reinforces the notion that the starburst activity is closely related to the inner shape of the potential well.

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 obtained a deep 8-field XMM-Newton mosaic of M33 covering the galaxy out to the D$_{25}$ isophote and beyond to a limiting 0.2--4.5 keV unabsorbed flux of 5$\times$10$^{-16}$ erg cm$^{-2}$ s$^{-1}$ (L${>}$4$\times$10$^{34}$ erg s$^{-1}$ at the distance of M33). These data allow complete coverage of the galaxy with high sensitivity to soft sources such as diffuse hot gas and supernova remnants. Here we describe the methods we used to identify and characterize 1296 point sources in the 8 fields. We compare our resulting source catalog to the literature, note variable sources, construct hardness ratios, classify soft sources, analyze the source density profile, and measure the X-ray luminosity function. As a result of the large effective area of XMM-Newton below 1 keV, the survey contains many new soft X-ray sources. The radial source density profile and X-ray luminosity function for the sources suggests that only $\sim$15% of the 391 bright sources with L${>}$3.6$\times$10$^{35}$ erg s$^{-1}$ are likely to be associated with M33, and more than a third of these are known supernova remnants. The log(N)--log(S) distribution, when corrected for background contamination, is a relatively flat power-law with a differential index of 1.5, which suggests many of the other M33 sources may be high-mass X-ray binaries. Finally, we note the discovery of an interesting new transient X-ray source, which we are unable to classify.

We explore the quenching of low-mass galaxies (10^4 < Mstar < 10^8 Msun) as a function of lookback time using the star formation histories (SFHs) of 38 Local Group dwarf galaxies. The SFHs were derived from analyzing color-magnitude diagrams of resolved stellar populations in archival Hubble Space Telescope/Wide Field Planetary Camera 2 imaging. We find: (1) Lower mass galaxies quench earlier than higher mass galaxies; (2) Inside of virial radius there is no correlation between a satellite's current proximity to a massive host and its quenching epoch; (3) There are hints of systematic differences in quenching times of M31 and Milky Way (MW) satellites, although the sample sample size and uncertainties in the SFHs of M31 dwarfs prohibit definitive conclusions. Combined with literature results, we qualitatively consider the redshift evolution (z=0-1) of the quenched galaxy fraction over ~7 dex in stellar mass (10^4 < Mstar < 10^11.5 Msun). The quenched fraction of all galaxies generally increases toward the present, with both the lowest and highest mass systems exhibiting the largest quenched fractions at all redshifts. In contrast, galaxies between Mstar ~ 10^8-10^10 Msun have the lowest quenched fractions. We suggest that such intermediate-mass galaxies are the least efficient at quenching. Finally, we compare our quenching times with predictions for infall times of low-mass galaxies associated with the MW. We find that some of the lowest-mass satellites (e.g., CVn II, Leo IV) may have been quenched before infall while higher mass satellites (e.g., Leo I, Fornax) typically quench ~1-4 Gyr after infall.

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]

Despite decades of intense efforts, many fundamental aspects of Type Ia supernova (SNe Ia) remain elusive. One of the major open questions is whether the mass of the exploding white dwarf (WD) is close to the Chandrasekhar limit. Here we report the detection of strong K-shell emission from stable Fe-peak elements in the Suzaku X-ray spectrum of the Type Ia supernova remnant (SNR) 3C 397. The high Ni/Fe and Mn/Fe mass ratios (0.11-0.24 and 0.018-0.033, respectively) in the hot plasma component that dominates the K-shell emission lines indicate a degree of neutronization in the SN ejecta which can only be achieved by electron captures in the dense cores of exploding WDs with a near-Chandrasekhar mass. This suggests a single-degenerate origin for 3C 397, since Chandrasekhar mass progenitors are expected naturally if the WD accretes mass slowly from a companion. Together with other results supporting the double-degenerate scenario, our work adds to the mounting evidence that both progenitor channels make a significant contribution to the SN Ia rate in star-forming galaxies.

The stellar kinematics of galactic disks are key to constraining disk formation and evolution processes. In this paper, for the first time, we measure the stellar age-velocity dispersion correlation in the inner 20 kpc (3.5 disk scale lengths) of M31 and show that it is dramatically different from that in the Milky Way. We use optical Hubble Space Telescope/Advanced Camera for Surveys photometry of 5800 individual stars from the Panchromatic Hubble Andromeda Treasury (PHAT) survey and Keck/DEIMOS radial velocity measurements of the same stars from the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) survey. We show that the average line-of-sight velocity dispersion is a steadily increasing function of stellar age exterior to R=10 kpc, increasing from 30 km/s for the young upper main sequence stars to 90 km/s for the old red giant branch stars. This monotonic increase implies that a continuous or recurring process contributed to the evolution of the disk. Both the slope and normalization of the dispersion vs. age relation are significantly larger than in the Milky Way, allowing for the possibility that the disk of M31 has had a more violent history than the disk of the Milky Way, more in line with cosmological predictions. We also find evidence for an inhomogeneous distribution of stars from a second kinematical component in addition to the dominant disk component. One of the largest and hottest high-dispersion patches is present in all age bins, and may be the signature of the end of the long bar.

We have analyzed new HST/ACS and HST/WFC3 imaging in F475W and F814W of two previously-unobserved fields along the M31 minor axis to confirm our previous constraints on the shape of M31's inner stellar halo. Both of these new datasets reach a depth of at least F814W$<$27 and clearly detect the blue horizontal branch (BHB) of the field as a distinct feature of the color-magnitude diagram. We measure the density of BHB stars and the ratio of BHB to red giant branch stars in each field using identical techniques to our previous work. We find excellent agreement with our previous measurement of a power-law for the 2-D projected surface density with an index of 2.6$^{+0.3}_{-0.2}$ outside of 3 kpc, which flattens to $\alpha <$1.2 inside of 3 kpc. Our findings confirm our previous suggestion that the field BHB stars in M31 are part of the halo population. However, the total halo profile is now known to differ from this BHB profile, which suggests that we have isolated the metal-poor component. This component appears to have an unbroken power-law profile from 3-150 kpc but accounts for only about half of the total halo stellar mass. Discrepancies between the BHB density profile and other measurements of the inner halo are therefore likely due to the different profile of the metal-rich halo component, which is not only steeper than the profile of the met al-poor component, but also has a larger core radius. These profile differences also help to explain the large ratio of BHB/RGB stars in our observations.

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 evolution and lifetimes of thermally pulsating asymptotic giant branch (TP-AGB) stars suffer from significant uncertainties. We present a detailed framework for constraining model luminosity functions of TP-AGB stars using resolved stellar populations. We show an example of this method that compares various TP-AGB mass-loss prescriptions that differ in their treatments of mass loss before the onset of dust-driven winds (pre-dust). We find that models with more efficient pre-dust driven mass loss produce results consistent with observations, as opposed to more canonical mass-loss models. Efficient pre-dust driven mass-loss predicts for [Fe/H] < -1.2, lower mass TP-AGB stars (M < 1 Msun) must have lifetimes less than about 1.2 Myr.

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 have measured stellar photometry with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) and Advanced Camera for Surveys (ACS) in near ultraviolet (F275W, F336W), optical (F475W, F814W), and near infrared (F110W, F160W) bands for 117 million resolved stars in M31. As part of the Panchromatic Hubble Andromeda Treasury (PHAT) survey, we measured photometry with simultaneous point spread function fitting across all bands and at all source positions after precise astrometric image alignment (<5-10 milliarcsecond accuracy). In the outer disk, the photometry reaches a completeness-limited depth of F475W~28, while in the crowded, high surface brightness bulge, the photometry reaches F475W~25. We find that simultaneous photometry and optimized measurement parameters significantly increase the detection limit of the lowest resolution filters (WFC3/IR) providing color-magnitude diagrams that are up to 2.5 magnitudes deeper when compared with color-magnitude diagrams from WFC3/IR photometry alone. We present extensive analysis of the data quality including comparisons of luminosity functions and repeat measurements, and we use artificial star tests to quantify photometric completeness, uncertainties and biases. We find that largest sources of systematic error in the photometry are due to spatial variations in the point spread function models and charge transfer efficiency corrections. This stellar catalog is the largest ever produced for equidistant sources, and is publicly available for download by the community.

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 present new optical spectroscopy of 20 candidate counterparts of 17 X-ray sources in the direction of the M31 disc. By comparing the X-ray catalogue from the XMM-Newton survey of M31 with star catalogues from the Local Group Galaxy Survey, we chose counterpart candidates based on optical colour and X-ray hardness. We have discovered 17 counterpart candidates with spectra containing stellar features. Eight of these are early-type stars of O or B type in M31, with hard X-ray spectra, making them good HMXB candidates. Three of these eight exhibit emission lines, which we consider to be the strongest HMXB candidates. In addition, our spectra reveal two likely Galactic cataclysmic variables, one foreground M star, two probable LMXBs related to M31 globular clusters, one emission line region with an embedded Wolf-Rayet star, and one newly-discovered supernova remnant. Finally, two of the sources have stellar spectra with no features indicative of association with an X-ray source.

Several young supernova remnants (SNRs), including SN1006, emit synchrotron X-rays in narrow filaments, hereafter thin rims, along their periphery. The widths of these rims imply 50 to 100 $\mu$G fields in the region immediately behind the shock, far larger than expected for the interstellar medium compressed by unmodified shocks, assuming electron radiative losses limit rim widths. However, magnetic-field damping could also produce thin rims. Here we review the literature on rim width calculations, summarizing the case for magnetic-field amplification. We extend these calculations to include an arbitrary power-law dependence of the diffusion coefficient on energy, $D \propto E^{\mu}$. Loss-limited rim widths should shrink with increasing photon energy, while magnetic-damping models predict widths almost independent of photon energy. We use these results to analyze Chandra observations of SN 1006, in particular the southwest limb. We parameterize the full widths at half maximum (FWHM) in terms of energy as FWHM $\propto E^{m_E}_{\gamma}$. Filament widths in SN1006 decrease with energy; $m_E \sim -0.3$ to $-0.8$, implying magnetic field amplification by factors of 10 to 50, above the factor of 4 expected in strong unmodified shocks. For SN 1006, the rapid shrinkage rules out magnetic damping models. It also favors short mean free paths (small diffusion coefficients) and strong dependence of $D$ on energy ($\mu \ge 1$).

We report results from Spitzer observations of SNR 0509-68.7, also known as N103B, a young Type Ia supernova remnant in the Large Magellanic Cloud that shows interaction with a dense medium in its western hemisphere. Our images show that N103B has strong IR emission from warm dust in the post-shock environment. The post-shock gas density we derive, 45 cm$^{-3}$, is much higher than in other Type Ia remnants in the LMC, though a lack of spatial resolution may bias measurements towards regions of higher than average density. This density is similar to that in Kepler's SNR, a Type Ia interacting with a circumstellar medium. Optical images show H$\alpha$ emission along the entire periphery of the western portion of the shock, with [O III] and [S II] lines emitted from a few dense clumps of material where the shock has become radiative. The dust is silicate in nature, though standard silicate dust models fail to reproduce the "18 $\mu$m" silicate feature that peaks instead at 17.3 $\mu$m. We propose that the dense material is circumstellar material lost from the progenitor system, as with Kepler. If the CSM interpretation is correct, this remnant would become the second member, along with Kepler, of a class of Type Ia remnants characterized by interaction with a dense CSM hundreds of years post-explosion. A lack of N enhancement eliminates symbiotic AGB progenitors. The white dwarf companion must have been relatively unevolved at the time of the explosion.

The evolution and lifetimes of thermally pulsating asymptotic giant branch (TP-AGB) stars suffer from significant uncertainties. In this work, we analyze the numbers and luminosity functions of TP-AGB stars in six quiescent, low metallicity ([Fe/H] $\lesssim -0.86$) galaxies taken from the ANGST sample, using HST photometry in both optical and near-infrared filters. The galaxies contain over 1000 TP-AGB stars (at least 60 per field). We compare the observed TP-AGB luminosity functions and relative numbers of TP-AGB and RGB stars, to models generated from different suites of TP-AGB evolutionary tracks after adopting star formation histories (SFH) derived from the HST deep optical observations. We test various mass-loss prescriptions that differ in their treatments of mass-loss before the onset of dust-driven winds (pre-dust). These comparisons confirm that pre-dust mass-loss is important, since models that neglect pre-dust mass-loss fail to explain the observed TP-AGB/RGB ratio or the luminosity functions. In contrast, models with more efficient pre-dust mass-loss produce results consistent with observations. We find that for [Fe/H]$\lesssim-0.86$, lower mass TP-AGB stars ($M\lesssim1M_\odot$) must have lifetimes of $\sim0.5$ Myr and higher masses ($M\lesssim 3M_\odot$) must have lifetimes $\lesssim 1.2$ Myr. In addition, assuming our best-fitting mass-loss prescription, we show that the third dredge up has no significant effect on TP-AGB lifetimes in this mass and metallicity range.

Using resolved stellar photometry measured from archival HST imaging, we generate color-magnitude diagrams of the stars within 50 pc of the locations of historic core-collapse supernovae that took place in galaxies within 8 Mpc. We fit these color-magnitude distributions with stellar evolution models to determine the best-fit age distribution of the young population. We then translate these age distributions into probability distributions for the progenitor mass of each SNe. The measurements are anchored by the main-sequence stars surrounding the event, making them less sensitive to assumptions about binarity, post-main-sequence evolution, or circumstellar dust. We demonstrate that, in cases where the literature contains masses that have been measured from direct imaging, our measurements are consistent with (but less precise than) these measurements. Using this technique, we constrain the progenitor masses of 17 historic SNe, 11 of which have no previous estimates from direct imaging. Our measurements still allow the possibility that all SNe progenitor masses are <20 M_sun. However, the large uncertainties for the highest-mass progenitors also allow the possibility of no upper-mass cutoff.

We search for signatures of reionization in the star formation histories (SFHs) of 38 Local Group dwarf galaxies (10$^4$ $<$ M$_{\star}$ $<$ 10$^9$ M$_{\odot}$). The SFHs are derived from color-magnitude diagrams using archival Hubble Space Telescope/Wide Field Planetary Camera 2 imaging. Only five quenched galaxies (And V, And VI, And XIII, Leo IV, Hercules) are consistent with forming the bulk of their stars before reionization, when full uncertainties are considered. Observations of 13 of the predicted `true fossils' identified by Bovill & Ricotti show that only two (Hercules and Leo IV) indicate star formation quenched by reionization. However, both are within the virial radius of the Milky Way and evidence of tidal disturbance complicates this interpretation. We argue that the late-time gas capture scenario posited by Ricotti for the low mass, gas-rich, and star-forming fossil candidate Leo T is observationally indistinguishable from simple gas retention. Given the ambiguity between environmental effects and reionization, the best reionization fossil candidates are quenched low mass field galaxies (e.g., KKR 25).

We present uniformly measured star formation histories (SFHs) of 40 Local Group dwarf galaxies based on color-magnitude diagram (CMD) analysis from archival Hubble Space Telescope imaging. We demonstrate that accurate SFHs can be recovered from CMDs that do not reach the oldest main sequence turn-off (MSTO), but emphasize that the oldest MSTO is critical for precisely constraining the earliest epochs of star formation. We find that: (1) the average lifetime SFHs of dwarf spheroidals (dSphs) can be approximated by an exponentially declining SFH with $\tau$ $\sim$ 5 Gyr; (2) lower luminosity dSphs are less likely to have extended SFHs than more luminous dSphs; (3) the average SFHs of dwarf irregulars (dIrrs), transition dwarfs (dTrans), and dwarf ellipticals (dEs) can be approximated by the combination of an exponentially declining SFH ($\tau$ $\sim$ 3-4 Gyr) for lookback ages $>$ 10-12 Gyr ago and a constant SFH thereafter; (4) the observed fraction of stellar mass formed prior to z=2 ranges considerably (80\% for galaxies with M $<$ 10$^5$ M$_{\odot}$ to 30\% for galaxies with M$>$10$^7$ M$_{\odot}$) and is largely explained by environment; (5) the distinction between "ultra-faint" and "classical" dSphs is arbitrary; (6) LG dIrrs formed a significantly higher fraction of stellar mass prior to z=2 than the SDSS galaxies from Leiter 2012 and the SFHs from the abundance matching models of Behroozi et al. 2013. This may indicate higher than expected star-formation efficiencies at early times in low mass galaxies. Finally, we provide all the SFHs in tabulated electronic format for use by the community.

We have used optical observations of resolved stars from the Panchromatic Hubble Andromeda Treasury (PHAT) to measure the recent (< 500 Myr) star formation histories (SFHs) of 33 FUV-bright regions in M31. The region areas ranged from ~$10^4$ to $10^6$ pc$^2$, which allowed us to test the reliability of FUV flux as a tracer of recent star formation on sub-kpc scales. The star formation rates (SFRs) derived from the extinction-corrected observed FUV fluxes were, on average, consistent with the 100-Myr mean SFRs of the SFHs to within the 1$\sigma$ scatter. Overall, the scatter was larger than the uncertainties in the SFRs and particularly evident among the smallest regions. The scatter was consistent with an even combination of discrete sampling of the initial mass function and high variability in the SFHs. This result demonstrates the importance of satisfying both the full-IMF and the constant-SFR assumptions for obtaining precise SFR estimates from FUV flux. Assuming a robust FUV extinction correction, we estimate that a factor of 2.5 uncertainty can be expected in FUV-based SFRs for regions smaller than $10^5$ pc$^2$, or a few hundred pc. We also examined ages and masses derived from UV flux under the common assumption that the regions are simple stellar populations (SSPs). The SFHs showed that most of the regions are not SSPs, and the age and mass estimates were correspondingly discrepant from the SFHs. For those regions with SSP-like SFHs, we found mean discrepancies of 10 Myr in age and a factor of 3 to 4 in mass. It was not possible to distinguish the SSP-like regions from the others based on integrated FUV flux.

We report on spectra of two positions in the XA region of the Cygnus Loop supernova remnant obtained with the InfraRed Spectrograph on the Spitzer Space Telescope. The spectra span the 10-35 micron wavelength range, which contains a number of collisionally excited forbidden lines. These data are supplemented by optical spectra obtained at the Whipple Observatory and an archival UV spectrum from the International Ultraviolet Explorer. Coverage from the UV through the IR provides tests of shock wave models and tight constraints on model parameters. Only lines from high ionization species are detected in the spectrum of a filament on the edge of the remnant. The filament traces a 180 km/s shock that has just begun to cool, and the oxygen to neon abundance ratio lies in the normal range found for Galactic H II regions. Lines from both high and low ionization species are detected in the spectrum of the cusp of a shock-cloud interaction, which lies within the remnant boundary. The spectrum of the cusp region is matched by a shock of about 150 km/s that has cooled and begun to recombine. The post-shock region has a swept-up column density of about 1.3E18 cm^-2, and the gas has reached a temperature of 7000 to 8000 K. The spectrum of the Cusp indicates that roughly half of the refractory silicon and iron atoms have been liberated from the grains. Dust emission is not detected at either position.

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.

We introduce a deep (670 ks) X-ray survey of the SN1006 remnant from Chandra, plus a deep H-alpha image from the 4m telescope at CTIO. Comparison with Chandra images from 2003 gives the first measurement of X-ray proper motions around the entire rim. We find that the expansion velocity varies significantly with azimuth: the highest velocity of ~7400 km/s (almost 2.5 times that in the NW) is found along the SE rim, where both kinematics and spectra indicate that most of the X-rays stem from undecelerated ejecta. Asymmetries in the distribution of ejecta are seen on a variety of spatial scales. Si-rich ejecta are especially prominent in the SE quadrant, while O and Mg are more uniformly distributed, indicating large-scale asymmetries arising from the explosion itself. Ne emission is strongest in a sharp filament just behind the primary shock along the NW rim, where the preshock density is highest. Here the Ne is likely interstellar, while Ne within the shell may include a contribution from ejecta. Within the interior of the projected shell we find a few isolated "bullets" of what appear to be supernova ejecta that are immediately preceded by bowshocks seen in H-alpha--features we interpret as ejecta knots that have reached relatively dense regions of the surrounding ISM. Recent 3-dimensional hydrodynamic models for SN Ia display small-scale features that resemble the ones seen in X-rays in SN1006; an origin in the explosion itself or from subsequent instabilities both remain viable options. We have expanded the search for precursor X-ray emission ahead of a synchrotron-dominated shock front. Profiles along both NE and SW rims require that a precursor be thinner than about 3 arcsec and fainter than about 5% of the post-shock peak. These limits suggest that the magnetic field is amplified by a factor of 7 or more in a narrow precursor region, promoting diffusive particle acceleration.

The B-Machine Telescope is the culmination of several years of development, construction, characterization and observation. The telescope is a departure from standard polarization chopping of correlation receivers to a half wave plate technique. Typical polarimeters use a correlation receiver to chop the polarization signal to overcome the $1/f$ noise inherent in HEMT amplifiers. B-Machine uses a room temperature half wave plate technology to chop between polarization states and measure the polarization signature of the CMB. The telescope has a demodulated $1/f$ knee of 5 mHz and an average sensitivity of 1.6 $\mathrm{mK}\sqrt{\mathrm{s}}$. This document examines the construction, characterization, observation of astronomical sources, and data set analysis of B-Machine. Preliminary power spectra and sky maps with large sky coverage for the first year data set are included.

Destriping is a well-established technique for removing low-frequency correlated noise from Cosmic Microwave Background (CMB) survey data. In this paper we present a destriping algorithm tailored to data from a polarimeter, i.e. an instrument where each channel independently measures the polarization of the input signal. We also describe a fully parallel implementation in Python released as Free Software and analyze its results and performance on simulated datasets, both the design case of signal and correlated noise, and with additional systematic effects. Finally we apply the algorithm to 30 days of 37.5 GHz polarized microwave data gathered from the B-Machine experiment, developed at UCSB. The B-Machine data and destriped maps are made publicly available. The purpose is the development of a scalable software tool to be applied to the upcoming 12 months of temperature and polarization data from LATTE (Low frequency All sky TemperaTure Experiment) at 8 GHz and to even larger datasets.

We characterize the bulge, disk, and halo subcomponents in the Andromeda galaxy (M31) over the radial range 4 < R_proj < 225 kpc. The cospatial nature of these subcomponents renders them difficult to disentangle using surface brightness (SB) information alone, especially interior to ~20 kpc. Our new decomposition technique combines information from the luminosity function (LF) of over 1.5 million bright (20 < m_814W < 22) stars from the Panchromatic Hubble Andromeda Treasury (PHAT) survey, radial velocities of over 5000 red giant branch stars in the same magnitude range from the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) survey, and integrated I-band SB profiles from various sources. We use an affine-invariant Markov chain Monte Carlo algorithm to fit an appropriate toy model to these three data sets. The bulge, disk, and halo SB profiles are modeled as a Sersic, exponential, and cored power-law, respectively, and the LFs are modeled as broken power-laws. We present probability distributions for each of 32 parameters describing the SB profiles and LFs of the three subcomponents. We find that the number of stars with a disk-like LF is ~5% larger than the the number with disk-like (dynamically cold) kinematics, suggesting that some stars born in the disk have been dynamically heated to the point that they are kinematically indistinguishable from halo members. This is the first kinematical evidence for a "kicked-up disk" halo population in M31. The fraction of kicked-up disk stars is consistent with that found in simulations. We also find evidence for a radially varying disk LF, consistent with a negative metallicity gradient in the stellar disk.

Dust grains are sputtered away in the hot gas behind shock fronts in supernova remnants, gradually enriching the gas phase with refractory elements. We have measured emission in C IV $\lambda$1550 from C atoms sputtered from dust in the gas behind a non-radiative shock wave in the northern Cygnus Loop. Overall, the intensity observed behind the shock agrees approximately with predictions from model calculations that match the Spitzer 24 micron and the X-ray intensity profiles. Thus these observations confirm the overall picture of dust destruction in SNR shocks and the sputtering rates used in models. However, there is a discrepancy in that the CIV intensity 10" behind the shock is too high compared to the intensities at the shock and 25" behind it. Variations in the density, hydrogen neutral fraction and the dust properties over parsec scales in the pre-shock medium limit our ability to test dust destruction models in detail.

We use medium-band near-infrared (NIR) Hubble Space Telescope WFC3 photometry with model NIR spectra of Asymptotic Giant Branch (AGB) stars to develop a new tool for efficiently distinguishing carbon-rich (C-type) AGB stars from oxygen-rich (M-type) AGB stars in galaxies at the edge of and outside the Local Group. We present the results of a test of this method on a region of the inner disk of M31, where we find a surprising lack of C stars, contrary to the findings of previous C star searches in other regions of M31. We find only 1 candidate C star (plus up to 6 additional, less certain C stars candidates), resulting in an extremely low ratio of C to M stars (C/M = (3.3(+20,-0.1))x10^-4) that is 1-2 orders of magnitude lower than other C/M estimates in M31. The low C/M ratio is likely due to the high metallicity in this region which impedes stars from achieving C/O > 1 in their atmospheres. These observations provide stringent constraints evolutionary models of metal-rich AGB stars and suggest that there is a metallicity threshold above which M stars are unable to make the transition to C stars, dramatically affecting AGB mass loss and dust production and, consequently, the observed global properties of metal-rich galaxies.

Turbulent neutral hydrogen (HI) line widths are often thought to be driven primarily by star formation (SF), but the timescale for converting SF energy to HI kinetic energy is unclear. As a complication, studies on the connection between Hi line widths and SF in external galaxies often use broadband tracers for the SF rate, which must implicitly assume that SF histories (SFHs) have been constant over the timescale of the tracer. In this paper, we compare measures of HI energy to time-resolved SFHs in a number of nearby dwarf galaxies. We find that HI energy surface density is strongly correlated only with SF that occurred 30-40 Myr ago. This timescale corresponds to the approximate lifetime of the lowest mass supernova progenitors (~ 8 Msun). This analysis suggests that the coupling between SF and the neutral ISM is strongest on this timescale, due either to an intrinsic delay between the release of the peak energy from SF or to the coherent effects of many SNe during this interval. At \Sigma_SFR > 10^-3 Msun yr^-1 kpc^-2, we find a mean coupling efficiency between SF energy and HI energy of \epsilon = 0.11 +/- 0.04 using the 30-40 Myr timescale. However, unphysical efficiencies are required in lower \Sigma_SFR systems, implying that SF is not the primary driver of HI kinematics at \Sigma_SFR < 10^-3 Msun yr^-1 kpc^-2.

We use empirical star formation histories (SFHs), measured from HST-based resolved star color-magnitude diagrams, as input into population synthesis codes to model the broadband spectral energy distributions (SEDs) of ~50 nearby dwarf galaxies (6.5 < log M/M_* < 8.5, with metallicities ~10% solar). In the presence of realistic SFHs, we compare the modeled and observed SEDs from the ultraviolet (UV) through near-infrared (NIR) and assess the reliability of widely used UV-based star formation rate (SFR) indicators. In the FUV through i bands, we find that the observed and modeled SEDs are in excellent agreement. In the Spitzer 3.6micron and 4.5micron bands, we find that modeled SEDs systematically over-predict observed luminosities by up to ~0.2 dex, depending on treatment of the TP-AGB stars in the synthesis models. We assess the reliability of UV luminosity as a SFR indicator, in light of independently constrained SFHs. We find that fluctuations in the SFHs alone can cause factor of ~2 variations in the UV luminosities relative to the assumption of a constant SFH over the past 100 Myr. These variations are not strongly correlated with UV-optical colors, implying that correcting UV-based SFRs for the effects of realistic SFHs is difficult using only the broadband SED. Additionally, for this diverse sample of galaxies, we find that stars older than 100 Myr can contribute from <5% to100% of the present day UV luminosity, highlighting the challenges in defining a characteristic star formation timescale associated with UV emission. We do find a relationship between UV emission timescale and broadband UV-optical color, though it is different than predictions based on exponentially declining SFH models. Our findings have significant implications for the comparison of UV-based SFRs across low-metallicity populations with diverse SFHs.

We have examined resolved stellar photometry from HST imaging surrounding 18 high-mass X-ray binary (HMXB) candidates in NGC 300 and NGC 2403 as determined from combined Chandra/HST analysis. We have fit the color-magnitude distribution of the surrounding stars with stellar evolution models. All but one region in NGC 300 and two in NGC 2403 contain a population with an age between 20 and 70 Myr. One of the candidates is the ultraluminous X-ray source (ULX) in NGC 2403, which we associate with a 60 Myr old population. These age distributions provide additional evidence that 16 of these 18 candidates are HMXBs. Furthermore, our results suggest that the most common HMXB age in these galaxies is 40-55 Myr. This preferred age is similar to observations of HMXBs in the Small Magellanic Cloud, providing new evidence of this formation timescale, but in higher metallicity populations. We suggest that this preferred HMXB age is the result of the fortuitous combination of two physical effects. First, this is the age of a population when the greatest rate of core-collapse events should be occurring, maximizing neutron star production. Second, this is the age when B stars are most likely to be actively losing mass. We also discuss our results in the context of HMXB feedback in galaxies, confirming HMXBs as a potentially important source of energy for the interstellar medium in low-mass galaxies.

Many of the most exciting questions in astrophysics and cosmology, including the majority of observational probes of dark energy, rely on an understanding of the nonlinear regime of structure formation. In order to fully exploit the information available from this regime and to extract cosmological constraints, accurate theoretical predictions are needed. Currently such predictions can only be obtained from costly, precision numerical simulations. This paper is the third in a series aimed at constructing an accurate calibration of the nonlinear mass power spectrum on Mpc scales for a wide range of currently viable cosmological models, including dark energy. The first two papers addressed the numerical challenges, and the scheme by which an interpolator was built from a carefully chosen set of cosmological models. In this paper we introduce the "Coyote Univers"' simulation suite which comprises nearly 1,000 N-body simulations at different force and mass resolutions, spanning 38 wCDM cosmologies. This large simulation suite enables us to construct a prediction scheme, or emulator, for the nonlinear matter power spectrum accurate at the percent level out to k~1 h/Mpc. We describe the construction of the emulator, explain the tests performed to ensure its accuracy, and discuss how the central ideas may be extended to a wider range of cosmological models and applications. A power spectrum emulator code is released publicly as part of this paper.

{\it Spitzer} images of Tycho's supernova remnant in the mid-infrared reveal limb-brightened emission from the entire periphery of the shell and faint filamentary structures in the interior. As with other young remnants, this emission is produced by dust grains, warmed to $\sim 100$ K in the post-shock environment by collisions with energetic electrons and ions. The ratio of the 70 to 24 $\mu$m fluxes is a diagnostic of the dust temperature, which in turn is a sensitive function of the plasma density. We find significant variations in the 70/24 flux ratio around the periphery of Tycho's forward shock, implying order-of-magnitude variations in density. While some of these are likely localized interactions with dense clumps of the interstellar medium, we find an overall gradient in the ambient density surrounding Tycho, with densities 3-10 times higher in the NE than in the SW. This large density gradient is qualitatively consistent with the variations in the proper motion of the shock observed in radio and X-ray studies. Overall, the mean ISM density around Tycho is quite low ($\sim 0.1-0.2$ cm$^{-3}$), consistent with the lack of thermal X-ray emission observed at the forward shock. We perform two-dimensional hydrodynamic simulations of a Type Ia SN expanding into a density gradient in the ISM, and find that the overall round shape of the remnant is still easily acheivable, even for explosions into significant gradients. However, this leads to an offset of the center of the explosion from the geometric center of the remnant of up to 20%, although lower values of 10% are preferred. The best match with hydrodynamical simulations is achieved if Tycho is located at a large (3-4 kpc) distance in a medium with a mean preshock density of $\sim 0.2$ cm$^{-3}$. Such preshock densities are obtained for highly ($\ga 50$%) porous ISM grains.

The Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) observed 10 transits of XO-2b over a period of three years. We present measurements which confirm previous estimates of system parameters like the normalized semi-major axis (a/R_{*}), stellar density (\rho_{*}), impact parameter (b) and orbital inclination (i_{orb}). Our errors on system parameters like a/R_{*} and \rho_{*} have improved by ~40% compared to previous best ground-based measurements. Our study of the transit times show no evidence for transit timing variations and we are able to rule out co-planar companions with masses \ge 0.20 \mearth\ in low order mean motion resonance with XO-2b. We also explored the stability of the XO-2 system given various orbital configurations of a hypothetical planet near the 2:1 mean motion resonance. We find that a wide range of orbits (including Earth-mass perturbers) are both dynamically stable and produce observable TTVs. We find that up to 51% of our stable simulations show TTVs that are smaller than the typical transit timing errors (~20 sec) measured for XO-2b, and hence remain undetectable.

We present detailed analysis of color-magnitude diagrams of NGC2403, obtained from a deep (m<28) Hubble Space Telescope (HST) Wide Field Planetary Camera 2 observation of the outer disk of NGC2403, supplemented by several shallow (m<26) HST Advanced Camera for Surveys fields. We derive the spatially resolved star formation history of NGC2403 out to 11 disk scale lengths. In the inner portions of the galaxy, we compare the recent star formation rates (SFRs) we derive from the resolved stars with those measured using GALEX FUV + Spitzer 24-micron fluxes, finding excellent agreement between the methods. Our measurements also show that the radial gradient in recent SFR mirror s the disk exponential profile to 11 scale lengths with no break, extending to SFR densities a factor of 100 lower than those that can be measured with GALEX and Spitzer (2x10^{-6} M_{\sun} yr^{-1} kpc^{-2}). Furthermore, we find that the cumulative stellar mass of the disk was formed at similar times at all radii. We compare these characteristics of NGC2403 to those of its "morphological twins," NGC300 and M33, showing that the structure and age distributions of the NGC2403 disk are more similar to those of the relatively isolated system NGC300 than to those of the Local Group analog M33. We also discuss the environments and HI morphologies of these three nearby galaxies, comparing them to integrated light studies of larger samples of more distant galaxy disks. Taken together, the physical properties and evolutionary history of NGC2403 suggest that the galaxy has had no close encounters with other M81 group members and may be falling into the group for the first time.

We present transit observations of the WASP-2 exoplanet system by the Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) program. Model fitting to these data allows us to improve measurements of the hot-Jupiter exoplanet WASP-2b and its orbital parameters by a factor of ~2 over prior studies; we do not find evidence for transit depth variations. We do find reduced chi^2 values greater than 1.0 in the observed minus computed transit times. A sinusoidal fit to the residuals yields a timing semi-amplitude of 32 seconds and a period of 389 days. However, random rearrangements of the data provide similar quality fits, and we cannot with certainty ascribe the timing variations to mutual exoplanet interactions. This inconclusive result is consistent with the lack of incontrovertible transit timing variations (TTVs) observed in other hot-Jupiter systems. This outcome emphasizes that unique recognition of TTVs requires dense sampling of the libration cycle (e.g. continuous observations from space-based platforms). However, even in systems observed with the Kepler spacecraft, there is a noted lack of transiting companions and TTVs in hot-Jupiter systems. This result is more meaningful, and indicates that hot-Jupiter systems, while they are easily observable from the ground, do not appear to be currently configured in a manner favorable to the detection of TTVs. The future of ground-based TTV studies may reside in resolving secular trends, and/or implementation at extreme quality observing sites to minimize atmospheric red noise.

The 2013 Planetary Science Decadal Survey identified a detailed investigation of the Trojan asteroids occupying Jupiter's L4 and L5 Lagrange points as a priority for future NASA missions. Observing these asteroids and measuring their physical characteristics and composition would aid in identification of their source and provide answers about their likely impact history and evolution, thus yielding information about the makeup and dynamics of the early Solar System. We present a conceptual design for a mission to the Jovian Trojan asteroids: the Trojan ASteroid Tour, Exploration, and Rendezvous (TASTER) mission, that is consistent with the NASA New Frontiers candidate mission recommended by the Decadal Survey and the final result of the 2011 NASA-JPL Planetary Science Summer School. Our proposed mission includes visits to two Trojans in the L4 population: a 500 km altitude fly-by of 1999 XS143, followed by a rendezvous with and detailed observations of 911 Agamemnon at orbital altitudes of 1000 - 100 km over a 12 month nominal science data capture period. Our proposed instrument payload - wide- and narrow-angle cameras, a visual and infrared mapping spectrometer, and a neutron/gamma ray spectrometer - would provide unprecedented high-resolution, regional-to-global datasets for the target bodies, yielding fundamental information about the early history and evolution of the Solar System. Although our mission design was completed as part of an academic exercise, this study serves as a useful starting point for future Trojan mission design studies. In particular, we identify and discuss key issues that can make large differences in the complex trade-offs required when designing a mission to the Trojan asteroids.

We report results of infrared imaging and spectroscopic observations of the SN 1006 remnant, carried out with the Spitzer Space Telescope. The 24 micron image from MIPS clearly shows faint filamentary emission along the northwest rim of the remnant shell, nearly coincident with the Balmer filaments that delineate the present position of the expanding shock. The 24 micron emission traces the Balmer filaments almost perfectly, but lies a few arcsec within, indicating an origin in interstellar dust heated by the shock. Subsequent decline in the IR behind the shock is presumably due largely to grain destruction through sputtering. The emission drops far more rapidly than current models predict, however, even for a higher proportion of small grains than would be found closer to the Galactic plane. The rapid drop may result in part from a grain density that has always been lower -- a relic effect from an earlier epoch when the shock was encountering a lower density -- but higher grain destruction rates still seem to be required. Spectra from three positions along the NW filament from the IRS instrument all show only a featureless continuum, consistent with thermal emission from warm dust. The dust-to-gas mass ratio in the pre-shock interstellar medium is lower than that expected for the Galactic ISM -- as has also been observed in the analysis of IR emission from other SNRs but whose cause remains unclear. As with other SN Ia remnants, SN 1006 shows no evidence for dust grain formation in the supernova ejecta.

We present a comprehensive X-ray point source catalog of NGC 404 obtained as part of the Chandra Local Volume Survey. A new, 97 ks Chandra ACIS-S observation of NGC 404 was combined with archival observations for a total exposure of ~123 ks. Our survey yields 74 highly significant X-ray point sources and is sensitive to a limiting unabsorbed luminosity of ~6x10^35 erg s^-1 in the 0.35-8 keV band. To constrain the nature of each X-ray source, cross-correlations with multi-wavelength data were generated. We searched overlapping HST observations for optical counterparts to our X-ray detections, but find only two X-ray sources with candidate optical counterparts. We find 21 likely low mass X-ray binaries (LMXBs), although this number is a lower limit due to the difficulties in separating LMXBs from background AGN. The X-ray luminosity functions (XLFs) in both the soft and hard energy bands are presented. The XLFs in the soft band (0.5-2 keV) and the hard band (2-8 keV) have a limiting luminosity at the 90% completeness limit of 10^35 erg s^-1 and 10^36 erg s^-1, respectively, significantly lower than previous X-ray studies of NGC 404. We find the XLFs to be consistent with those of other X-ray populations dominated by LMXBs. However, the number of luminous (>10^37 erg s^-1) X-ray sources per unit stellar mass in NGC 404 is lower than is observed for other galaxies. The relative lack of luminous XRBs may be due to a population of LMXBs with main sequence companions formed during an epoch of elevated star formation ~0.5 Gyr ago.

The Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) observed eleven transits of TrES-3b over two years in order to constrain system parameters and look for transit timing and depth variations. We describe an updated analysis protocol for APOSTLE data, including the reduction pipeline, transit model and Markov Chain Monte Carlo analyzer. Our estimates of the system parameters for TrESb are consistent with previous estimates to within the 2\sigma\ confidence level. We improved the errors (by 10--30%) on system parameters like the orbital inclination ($i_{\text{orb}}$), impact parameter (b) and stellar density (\rho$_{\star}$) compared to previous measurements. The near-grazing nature of the system, and incomplete sampling of some transits, limited our ability to place reliable uncertainties on individual transit depths and hence we do not report strong evidence for variability. Our analysis of the transit timing data show no evidence for transit timing variations and our timing measurements are able to rule out Super-Earth and Gas Giant companions in low order mean motion resonance with TrES-3b.

We report the results of an X-ray proper motion measurement for the NW rim of SN1006, carried out by comparing Chandra observations from 2001 and 2012. The NW limb has predominantly thermal X-ray emission, and it is the only location in SN1006 with significant optical emission: a thin, Balmer-dominated filament. For most of the NW rim, the proper motion is about 0.30 arcsec/yr, essentially the same as has been measured from the H-alpha filament. Isolated regions of the NW limb are dominated by nonthermal emission, and here the proper motion is much higher, 0.49 arcsec/yr, close to the value measured in X-rays along the much brighter NE limb, where the X-rays are overwhelmingly nonthermal. At the 2.2 kpc distance to SN1006, the proper motions imply shock velocities of about 3000 km/s and 5000 km/s in the thermal and nonthermal regions, respectively. A lower velocity behind the H-alpha filament is consistent with the picture that SN1006 is encountering denser gas in the NW, as is also suggested by its overall morphology. In the thermally-dominated portion of the X-ray shell, we also see an offset in the radial profiles at different energies; the 0.5-0.6 keV peak dominated by O VII is closer to the shock front than that of the 0.8-3 keV emission--due to the longer times for heavier elements to reach ionization states where they produce strong X-ray emission.

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.

We map the spatial distribution of recent star formation over a few x 100 Myr timescales in fifteen starburst dwarf galaxies using the location of young blue helium burning stars identified from optically resolved stellar populations in archival Hubble Space Telescope observations. By comparing the star formation histories from both the high surface brightness central regions and the diffuse outer regions, we measure the degree to which the star formation has been centrally concentrated during the galaxies' starbursts, using three different metrics for the spatial concentration. We find that the galaxies span a full range in spatial concentration, from highly centralized to broadly distributed star formation. Since most starbursts have historically been identified by relatively short timescale star formation tracers (e.g., Halpha emission), there could be a strong bias towards classifying only those galaxies with recent, centralized star formation as starbursts, while missing starbursts that are spatially distributed.

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.

We attempt to constrain the shape of M31's inner stellar halo by tracing the surface density of blue horizontal branch (BHB) stars at galactocentric distances ranging from 2 kpc to 35 kpc. Our measurements make use of resolved stellar photometry from a section of the Panchromatic Hubble Andromeda Treasury (PHAT) survey, supplemented by several archival Hubble Space Telescope observations. We find that the ratio of BHB to red giant stars is relatively constant outside of 10 kpc, suggesting that the BHB is as reliable a tracer of the halo population as the red giant branch. In the inner halo, we do not expect BHB stars to be produced by the high metallicity bulge and disk, making BHB stars a good candidate to be a reliable tracer of the stellar halo to much smaller galactocentric distances. If we assume a power-law profile r^(-\alpha) for the 2-D projected surface density BHB distribution, we obtain a high-quality fit with a 2-D power-law index of \alpha=2.6^{+0.3}_{-0.2} outside of 3 kpc, which flattens to \alpha<1.2 inside of 3 kpc. This slope is consistent with previous measurements but is anchored to a radial baseline that extends much farther inward. Finally, assuming azimuthal symmetry and a constant mass-to-light ratio, the best-fitting profile yields a total halo stellar mass of 2.1^{+1.7}_{-0.4} x 10^9 M_sun. These properties are comparable with both simulations of stellar halo formation formed by satellite disruption alone, and with simulations that include some in situ formation of halo stars.

We present the source catalog of a new Chandra ACIS-I observation of NGC 300 obtained as part of the Chandra Local Volume Survey. Our 63 ks exposure covers ~88% of the D25 isophote (R~6.3 kpc) and yields a catalog of 95 X-ray point sources detected at high significance to a limiting unabsorbed 0.35-8 keV luminosity of ~10^36 erg s^-1. Sources were cross-correlated with a previous XMM-Newton catalog, and we find 75 "X-ray transient candidate" sources that were detected by one observatory, but not the other. We derive an X-ray scale length of 1.7+/-0.2 kpc and a recent star formation rate of 0.12 Msun yr^-1, in excellent agreement with optical observations. Deep, multi-color imaging from the Hubble Space Telescope, covering ~32% of our Chandra field, was used to search for optical counterparts to the X-ray sources, and we have developed a new source classification scheme to determine which sources are likely X-ray binaries, supernova remnants, and background AGN candidates. Finally, we present the X-ray luminosity functions (XLFs) at different X-ray energies, and we find the total NGC 300 X-ray point source population to be consistent with other late-type galaxies hosting young stellar populations (<50 Myr). We find XLF of sources associated with older stellar populations has a steeper slope than the XLF of X-ray sources coinciding with young stellar populations, consistent with theoretical predictions.

Variable NaI absorption lines have been reported in a number of type Ia supernovae (SNeIa). The presence of this circumstellar material suggests that cataclysmic variables (CVs) with a giant donor star may be the progenitors of these SNeIa (Patat et al. 2007). We present echelle spectra of the CV QU Carinae which strengthen the connection between CVs of the V Sge class, the Accretion Wind Evolution scenario, variable wind features, variable NaI absorption, and SNIa. This thread not only provides insight into the spectral peculiarities of QU Car, but also links SNeIa as a class with their parent systems.

We present the analysis of supernova remnants (SNRs) and candidates in M31 identified in the XMM-Newton large programme survey of M31. SNRs are among the bright X-ray sources in a galaxy. They are good indicators of recent star formation activities of a galaxy and of the interstellar environment in which they evolve. By combining the X-ray data of sources in M31 with optical data as well as with optical and radio catalogues, we aim to compile a complete, revised list of SNRs emitting X-rays in M31 detected with XMM-Newton, study their luminosity and spatial distribution, and understand the X-ray spectrum of the brightest SNRs. We analysed the X-ray spectra of the twelve brightest SNRs and candidates using XMM-Newton data. The four brightest sources allowed us to perform a more detailed spectral analysis and the comparison of different models to describe their spectrum. For all M31 large programme sources we searched for optical counterparts on the Ha, [Sii], and [Oiii] images of the Local Group Galaxy Survey. We confirm 21 X-ray sources as counterparts of known SNRs. In addition, we identify five new X-ray sources as X-ray and optically emitting SNRs. Seventeen sources are no longer considered as SNR candidates. We have thus created a list of 26 X-ray SNRs and 20 candidates in M31 based on their X-ray, optical, and radio emission, which is the most recent complete list of X-ray SNRs in M31. The brightest SNRs have X-ray luminosities of up to 8 x 10^36 erg/s in the 0.35 - 2.0 keV band.

As part of the Panchromatic Hubble Andromeda Treasury (PHAT) multi-cycle program, we observed a 12' \times 6.5' area of the bulge of M31 with the WFC3/UVIS filters F275W and F336W. From these data we have assembled a sample of \sim4000 UV-bright, old stars, vastly larger than previously available. We use updated Padova stellar evolutionary tracks to classify these hot stars into three classes: Post-AGB stars (P-AGB), Post-Early AGB (PE-AGB) stars and AGB-manqu\'e stars. P-AGB stars are the end result of the asymptotic giant branch (AGB) phase and are expected in a wide range of stellar populations, whereas PE-AGB and AGB-manqu\'e (together referred to as the hot post-horizontal branch; HP-HB) stars are the result of insufficient envelope masses to allow a full AGB phase, and are expected to be particularly prominent at high helium or {\alpha} abundances when the mass loss on the RGB is high. Our data support previous claims that most UV-bright sources in the bulge are likely hot (extreme) horizontal branch stars (EHB) and their progeny. We construct the first radial profiles of these stellar populations, and show that they are highly centrally concentrated, even more so than the integrated UV or optical light. However, we find that this UV-bright population does not dominate the total UV luminosity at any radius, as we are detecting only the progeny of the EHB stars that are the likely source of the UVX. We calculate that only a few percent of MS stars in the central bulge can have gone through the HP-HB phase and that this percentage decreases strongly with distance from the center. We also find that the surface density of hot UV-bright stars has the same radial variation as that of low-mass X-ray binaries. We discuss age, metallicity, and abundance variations as possible explanations for the observed radial variation in the UV-bright population.

Characterization of the relatively poorly-understood progenitor systems of Type Ia supernovae is of great importance in astrophysics, particularly given the important cosmological role that these supernovae play. Kepler's Supernova Remnant, the result of a Type Ia supernova, shows evidence for an interaction with a dense circumstellar medium (CSM), suggesting a single-degenerate progenitor system. We present 7.5-38 $\mu$m infrared (IR) spectra of the remnant, obtained with the {\it Spitzer Space Telescope}, dominated by emission from warm dust. Broad spectral features at 10 and 18 $\mu$m, consistent with various silicate particles, are seen throughout. These silicates were likely formed in the stellar outflow from the progenitor system during the AGB stage of evolution, and imply an oxygen-rich chemistry. In addition to silicate dust, a second component, possibly carbonaceous dust, is necessary to account for the short-wavelength IRS and IRAC data. This could imply a mixed chemistry in the atmosphere of the progenitor system. However, non-spherical metallic iron inclusions within silicate grains provide an alternative solution. Models of collisionally-heated dust emission from fast shocks ($>$ 1000 km s$^{-1}$) propagating into the CSM can reproduce the majority of the emission associated with non-radiative filaments, where dust temperatures are $\sim 80-100$ K, but fail to account for the highest temperatures detected, in excess of 150 K. We find that slower shocks (a few hundred km s$^{-1}$) into moderate density material ($n_{0} \sim 50-250$ cm$^{-3}$) are the only viable source of heating for this hottest dust. We confirm the finding of an overall density gradient, with densities in the north being an order of magnitude greater than those in the south.

Using star formation histories derived from optically resolved stellar populations in nineteen nearby starburst dwarf galaxies observed with the Hubble Space Telescope, we measure the stellar mass surface densities of stars newly formed in the bursts. By assuming a star formation efficiency (SFE), we then calculate the inferred gas surface densities present at the onset of the starbursts. Assuming a SFE of 1%, as is often assumed in normal star-forming galaxies, and assuming that the gas was purely atomic, translates to very high HI surface densities (~10^2-10^3 Msun pc^-2), which are much higher than have been observed in dwarf galaxies. This implies either higher values of SFE in these dwarf starburst galaxies or the presence of significant amounts of H_2 in dwarfs (or both). Raising the assumed SFEs to 10% or greater (in line with observations of more massive starbursts associated with merging galaxies), still results in HI surface densities higher than observed in 10 galaxies. Thus, these observations appear to require that a significant fraction of the gas in these dwarf starbursts galaxies was in the molecular form at the onset of the bursts. Our results imply molecular gas column densities in the range 10^19-10^21 cm^-2 for the sample. In those galaxies where CO observations have been made, these densities correspond to values of the CO-H_2 conversion factor (X_CO) in the range >3-80x10^20 cm^-2 (K km s^-1)^-1, or up to 40x greater than Galactic X_CO values.

The combination of large size, high stellar density, high metallicity, and Sersic surface brightness profile of the spheroidal component of the Andromeda galaxy (M31) within R_proj ~ 20 kpc suggest that it is unlike any subcomponent of the Milky Way. In this work we capitalize on our proximity to and external view of M31 to probe the kinematical properties of this "inner spheroid." We employ a Markov chain Monte Carlo (MCMC) analysis of resolved stellar kinematics from Keck/DEIMOS spectra of 5651 red giant branch stars to disentangle M31's inner spheroid from its stellar disk. We measure the mean velocity and dispersion of the spheroid in each of five spatial bins after accounting for a locally cold stellar disk as well as the Giant Southern Stream and associated tidal debris. For the first time, we detect significant spheroid rotation (v_rot ~ 50 km/s) beyond R_proj ~ 5 kpc. The velocity dispersion decreases from about 140 km/s at R_proj = 7 kpc to 120 km/s at R_proj = 14 kpc, consistent to 2 sigma with existing measurements and models. We calculate the probability that a given star is a member of the spheroid and find that the spheroid has a significant presence throughout the spatial extent of our sample. Lastly, we show that the flattening of the spheroid is due to velocity anisotropy in addition to rotation. Though this suggests that the inner spheroid of M31 more closely resembles an elliptical galaxy than a typical spiral galaxy bulge, it should be cautioned that our measurements are much farther out (2 - 14 r_eff) than for the comparison samples.

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.

The Panchromatic Hubble Andromeda Treasury (PHAT) is an on-going HST Multicycle Treasury program to image ~1/3 of M31's star forming disk in 6 filters, from the UV to the NIR. The full survey will resolve the galaxy into more than 100 million stars with projected radii from 0-20 kpc over a contiguous 0.5 square degree area in 828 orbits, producing imaging in the F275W and F336W filters with WFC3/UVIS, F475W and F814W with ACS/WFC, and F110W and F160W with WFC3/IR. The resulting wavelength coverage gives excellent constraints on stellar temperature, bolometric luminosity, and extinction for most spectral types. The photometry reaches SNR=4 at F275W=25.1, F336W=24.9, F475W=27.9, F814W=27.1, F110W=25.5, and F160W=24.6 for single pointings in the uncrowded outer disk; however, the optical and NIR data are crowding limited, and the deepest reliable magnitudes are up to 5 magnitudes brighter in the inner bulge. All pointings are dithered and produce Nyquist-sampled images in F475W, F814W, and F160W. We describe the observing strategy, photometry, astrometry, and data products, along with extensive tests of photometric stability, crowding errors, spatially-dependent photometric biases, and telescope pointing control. We report on initial fits to the structure of M31's disk, derived from the density of RGB stars, in a way that is independent of the assumed M/L and is robust to variations in dust extinction. These fits also show that the 10 kpc ring is not just a region of enhanced recent star formation, but is instead a dynamical structure containing a significant overdensity of stars with ages >1 Gyr. (Abridged)

Using high spatial resolution HST WFC3 and ACS imaging of resolved stellar populations, we constrain the contribution of thermally-pulsing asymptotic giant branch (TP-AGB) stars and red helium burning (RHeB) stars to the 1.6 um near-infrared (NIR) luminosities of 23 nearby galaxies. The TP-AGB phase contributes as much as 17% of the integrated F160W flux, even when the red giant branch is well populated. The RHeB population contribution can match or even exceed the TP-AGB contribution, providing as much as 21% of the integrated F160W light. The NIR mass-to-light (M/L) ratio should therefore be expected to vary significantly due to fluctuations in the star formation rate over timescales from 25 Myr to several Gyr. We compare our observational results to predictions based on optically derived star formation histories and stellar population synthesis (SPS) models, including models based on the Padova isochrones (used in popular SPS programs). The SPS models generally reproduce the expected numbers of TP-AGB stars in the sample. The same SPS models, however, give a larger discrepancy in the F160W flux contribution from the TP-AGB stars, over-predicting the flux by a weighted mean factor of 2.3 +/-0.8. This larger offset is driven by the prediction of modest numbers of high luminosity TP-AGB stars at young (<300 Myrs) ages. The best-fit SPS models simultaneously tend to under-predict the numbers and fluxes of stars on the RHeB sequence, typically by a factor of 2.0+/-0.6 for galaxies with significant numbers of RHeBs. Coincidentally, over-prediction of the TP-AGB and under-prediction of the RHeBs result in a NIR M/L ratio largely unchanged for a rapid star formation rate. However, the NIR-to-optical flux ratio of galaxies could be significantly smaller than AGB-rich models would predict, an outcome that has been observed in some intermediate redshift post-starburst galaxies. (Abridged)

We present the star formation history (SFH) of the faintest known star-forming galaxy, Leo T, based on imaging taken with the Hubble Space Telescope (HST) Wide Field Planetary Camera 2 (WFPC2). The HST/WFPC2 color-magnitude diagram (CMD) of Leo T is exquisitely deep, extending ~ 2 magnitudes below the oldest main sequence turnoff, permitting excellent constraints on star formation at all ages. We use a maximum likelihood CMD fitting technique to measure the SFH of Leo T assuming three different sets of stellar evolution models: Padova (solar-scaled metallicity) and BaSTI (both solar-scaled and alpha-enhanced metallicities). The resulting SFHs are remarkably consistent at all ages, indicating that our derived SFH is robust to the choice of stellar evolution model. From the lifetime SFH of Leo T, we find that 50% of the total stellar mass formed prior to z ~ 1 (7.6 Gyr ago). Subsequent to this epoch, the SFH of Leo T is roughly constant until the most recent ~ 25 Myr, where the SFH shows an abrupt drop. This decrease could be due to a cessation of star formation or stellar initial mass function sampling effects, but we are unable to distinguish between the two scenarios. Overall, our measured SFH is consistent with previously derived SFHs of Leo T. However, the HST-based solution provides improved age resolution and reduced uncertainties at all epochs. The SFH, baryonic gas fraction, and location of Leo T are unlike any of the other recently discovered faint dwarf galaxies in the Local Group, and instead bear strong resemblance to gas-rich dwarf galaxies (irregular or transition), suggesting that gas-rich dwarf galaxies may share common modes of star formation over a large range of stellar mass (~ 10^5-10^9 Msun).

Using Hubble Space Telescope (HST) photometry, we characterize the age of the stellar association in the vicinity of supernova (SN) 2011dh and use it to infer the zero-age main sequence mass (M_{ZAMS}) of the progenitor star. We find two distinct and significant star formation events with ages of <6 and 17^{+3}_{-4} Myrs, and the corresponding M_{ZAMS} are >29 and 13^{+2}_{-1} M_{Sun}, respectively. These two bursts represent 18^{+4}_{-9}% (young) and 64^{+10}_{-14}% (old) of the total star formation in the last 50 Myrs. Adopting these fractions as probabilities suggests that the most probable M_{ZAMS} is 13^{+2}_{-1} M_{Sun}. These results are most sensitive to the luminosity function along the well-understood main sequence and are less sensitive to uncertain late-stage stellar evolution. Therefore, they stand even if the progenitor suffered disruptive post-main-sequence evolution (e.g. eruptive mass loss or binary Roche-lobe overflow). Progenitor identification will help to further constrain the appropriate population. Even though pre-explosion images show a yellow supergiant (YSG) at the site of the SN, panchromatic SN light curves suggest a more compact star as the progenitor. In spite of this, our results suggest an association between the YSG and the SN. Not only was the star located at the SN site, but reinforcing an association, the star's bolometric luminosity is consistent with the final evolutionary stage of the 17 Myr old star burst. If the YSG disappears, then M_{ZAMS}=13^{+2}_{-1} M_{Sun}, but if it persists, then our results allow the possibility that the progenitor was an unseen star of >29 M_{Sun}.

We present near-infrared (NIR) color-magnitude diagrams (CMDs) for the resolved stellar populations within 26 fields of 23 nearby galaxies (<4 Mpc), based on F110W and F160W images from Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST). The CMDs sample both old dormant and young star-forming populations. We match key NIR CMD features with their counterparts in optical CMDs, and identify the red core Helium burning (RHeB) sequence as a significant contributor to the NIR flux in stellar populations younger than a few 100 Myrs old, suggesting that star formation can drive surprisingly rapid variations in the NIR mass-to-light ratio. The NIR luminosity of star forming galaxies is therefore not necessarily proportional to the stellar mass. We note that these individual bright RHeB stars may be misidentified as old stellar clusters in low resolution imaging. We also discuss the CMD location of asymptotic giant branch (AGB) stars, and the separation of AGB sub-populations using a combination of optical and NIR colors. We empirically calibrate the NIR magnitude of the tip of the red giant branch (TRGB) as a function of color, allowing this widely adopted filter to be used for distance measurements. We find a clear trend between NIR RGB color and metallicity. However, it appears unlikely that the slope of the NIR RGB can be used as a metallicity indicator in extragalactic systems with comparable data. Finally, we discuss scattered light in the WFC3, which becomes significant for exposures taken close to a bright earth limb.

We consider the effects of non-constant star formation histories (SFHs) on Halpha and GALEX far ultra-violet (FUV) star formation rate (SFR) indicators. Under the assumption of a fully populated Chabrier IMF, we compare the distribution of Halpha-to-FUV flux ratios from ~ 1500 simple, periodic model SFHs with observations of 185 galaxies from the Spitzer Local Volume Legacy survey. We find a set of SFH models that are well matched to the data, such that more massive galaxies are best characterized by nearly constant SFHs, while low mass systems experience bursts amplitudes of ~ 30 (i.e., an increase in the SFR by a factor of 30 over the SFR during the inter-burst period), burst durations of tens of Myr, and periods of ~ 250 Myr; these SFHs are broadly consistent with the increased stochastic star formation expected in systems with lower SFRs. We analyze the predicted temporal evolution of galaxy stellar mass, R-band surface brightness, Halpha-derived SFR, and blue luminosity, and find that they provide a reasonable match to observed flux distributions. We find that our model SFHs are generally able to reproduce both the observed systematic decline and increased scatter in Halpha-to-FUV ratios toward low mass systems, without invoking other physical mechanisms. We also compare our predictions with those from the Integrated Galactic IMF theory with a constant SFR. We find that while both predict a systematic decline in the observed ratios, only the time variable SFH models are capable of producing the observed population of low mass galaxies ($M_{*}$ < 10$^{7}$ Msun) with normal Halpha-to-FUV ratios. These results demonstrate that a variable IMF alone has difficulty explaining the observed scatter in the Halpha-to-FUV ratios. We conclude by considering the limitations of the model SFHs, and discuss the use of additional empirical constraints to improve future SFH modeling efforts.

We compare the cumulative star formation histories (SFHs) of Local Group (LG) dwarf galaxies with those in the volume-limited ACS Nearby Galaxy Survey Treasury (ANGST) sample (D < 4 Mpc), in order to understand how typical the LG dwarf galaxies are relative to those in the nearby universe. The SFHs were derived in a uniform manner from high quality optical color-magnitude diagrams constructed from Hubble Space Telescope imaging. We find that the {\it mean} cumulative SFHs of the LG dwarfs are comparable to the mean cumulative SFHs of the ANGST sample for the three different morphological types (dwarf spheroidals/ellipticals: dSph/dE; dwarf irregulars: dI; transition dwarfs: dTrans). We also discuss effects such as population gradients and systematic uncertainties in the stellar models that may influence the derived SFHs. Both the ANGST and Local Group dwarf galaxies show a consistent and strong morphology-density relationship, emphasizing the importance of environment in the evolution of dwarf galaxies. Specifically, we confirm that dIs are found at lower densities and higher luminosities than dSphs, within this large sample. We also find that dTrans are located in similar environments to those occupied by dwarf irregular galaxies, but have systematically lower luminosities that are more comparable to those of dwarf spheroidals. The similarity of the SFHs and morphology-density relationships of the LG and ANGST dwarf galaxies suggests that the LG dwarfs are a good representation of dwarf galaxies in the local universe.

We present the results of a 100 ks {\it Chandra} observation of the NGC 404 nuclear region. The long exposure and excellent spatial resolution of {\it Chandra} has enabled us to critically examine the nuclear environment of NGC 404, which is known to host a nuclear star cluster and potentially an intermediate-mass black hole (on the order of a few times $10^5$ \Msun). We find two distinct X-ray sources: a hard, central point source coincident with the optical and radio centers of the galaxy, and a soft extended region that is coincident with areas of high H$\alpha$ emission and likely recent star formation. When we fit the 0.3-8 keV spectra of each region separately, we find the hard nuclear point source to be dominated by a power law (\PL = 1.88), while the soft off-nuclear region is best fit by a thermal plasma model ($kT$ = 0.67 keV). We therefore find evidence for both a power law component and hot gas in the nuclear region of NGC 404. We estimate the 2-10 keV luminosity to be 1.3$^{+0.8}_{-0.5}\times10^{37}$ erg s$^{-1}$. A low-level of diffuse X-ray emission is detected out to $\sim$15\as ($\sim$0.2 kpc) from the nucleus. We compare our results to the observed relationships between power law photon index and Eddington ratio for both X-ray binaries and low luminosity active galaxies and find NGC 404 to be consistent with other low luminosity active galaxies. We therefore favor the conclusion that NGC 404 harbors an intermediate-mass black hole accreting at a very low level.

We derive the optical luminosity, colors, and ratios of the blue and red helium burning (HeB) stellar populations from archival Hubble Space Telescope observations of nineteen starburst dwarf galaxies and compare them with theoretical isochrones from Padova stellar evolution models across metallicities from Z=0.001 to 0.009. We find that the observational data and the theoretical isochrones for both blue and red HeB populations overlap in optical luminosities and colors and the observed and predicted blue to red HeB ratios agree for stars older than 50 Myr over the time bins studied. These findings confirm the usefulness of applying isochrones to interpret observations of HeB populations. However, there are significant differences, especially for the red HeB population. Specifically we find: (1) offsets in color between the observations and theoretical isochrones of order 0.15 mag (0.5 mag) for the blue (red) HeB populations brighter than M_V ~ -4 mag, which cannot be solely due to differential extinction; (2) blue HeB stars fainter than M_V ~ -3 mag are bluer than predicted; (3) the slope of the red HeB sequence is shallower than predicted by a factor of ~3; and (4) the models overpredict the ratio of the most luminous blue to red HeB stars corresponding to ages <50 Myr. Additionally, we find that for the more metal-rich galaxies in our sample (Z> 0.5 Zsolar) the red HeB stars overlap with the red giant branch stars in the color magnitude diagrams, thus reducing their usefulness as indicators of star formation for ages >100 Myr.

We report results from a multi-wavelength analysis of the Galactic SNR RCW 86, the proposed remnant of the supernova of 185 A.D. We report new infrared observations from {\it Spitzer} and {\it WISE}, where the entire shell is detected at 24 and 22 $\mu$m. We fit the infrared flux ratios with models of collisionally heated ambient dust, finding post-shock gas densities in the non-radiative shocks of 2.4 and 2.0 cm$^{-3}$ in the SW and NW portions of the remnant, respectively. The Balmer-dominated shocks around the periphery of the shell, large amount of iron in the X-ray emitting ejecta, and lack of a compact remnant support a Type Ia origin for this remnant. From hydrodynamic simulations, the observed characteristics of RCW 86 are successfully reproduced by an off-center explosion in a low-density cavity carved by the progenitor system. This would make RCW 86 the first known case of a Type Ia supernova in a wind-blown bubble. The fast shocks ($> 3000$ km s$^{-1}$) observed in the NE are propagating in the low-density bubble, where the shock is just beginning to encounter the shell, while the slower shocks elsewhere have already encountered the bubble wall. The diffuse nature of the synchrotron emission in the SW and NW is due to electrons that were accelerated early in the lifetime of the remnant, when the shock was still in the bubble. Electrons in a bubble could produce gamma-rays by inverse-Compton scattering. The wind-blown bubble scenario requires a single-degenerate progenitor, which should leave behind a companion star.

Aims. M 31 has been completely imaged for the first time obtaining a luminosity lower limit ~10^35 erg s-1 in the 0.2-4.5 keV band. Our XMM-Newton EPIC survey combines archival observations with observations taken between June 2006 and February 2008. The main goal of the paper is a study of the X-ray source population of M 31. Methods. An X-ray catalogue of 1897 sources was created; 914 were detected for the first time. Source classification and identification were based on X-ray hardness ratios, spatial extent of the sources, and by cross correlating with catalogues in the X-ray, optical, infrared and radio wavelengths. We also analysed the long-term variability of the X-ray sources. Results. Fourteen of the 30 supersoft source candidates represent supersoft emission of optical novae. Many of the 25 supernova remnants (SNRs) and 31 SNR candidates lie within the 10 kpc dust ring and other star forming regions in M 31; this connection between SNRs and star forming regions implies that most of the remnants originate in type II supernovae. The brightest sources in X-rays in M 31 belong to the class of X-ray binaries (XRBs). Ten low mass XRBs (LMXBs) and 26 LMXB candidates were identified based on their temporal variability. In addition 36 LMXBs and 17 LMXB candidates were identified due to correlations with globular clusters and globular cluster candidates. From optical and X-ray colour-colour diagrams, possible high mass XRB (HMXB) candidates were selected. Two of these candidates have an X-ray spectrum as is expected for an HMXB containing a neutron star primary. Conclusions. While our survey has greatly improved our understanding of the X-ray source populations in M 31, at this point 65% of the sources can still only be classified as "hard" sources. Deeper observations in X-ray and at other wavelengths would help classify these sources. ABRIDGED

We present a measurement of the age distribution of stars residing in spiral disks and dwarf galaxies. We derive a complete star formation history of the ~140 Mpc^3 covered by the volume-limited sample of galaxies in the Advanced Camera for Surveys (ACS) Nearby Galaxy Survey Treasury (ANGST). The total star formation rate density history is dominated by the large spirals in the volume, although the sample consists mainly of dwarf galaxies. Our measurement shows a factor of ~3 drop at z~2, in approximate agreement with results from other measurement techniques. While our results show that the overall star formation rate density has decreased since z~1, the measured rates during this epoch are higher than those obtained from other measurement techniques. This enhanced recent star formation rate appears to be largely due to an increase in the fraction of star formation contained in low-mass disks at recent times. Finally, our results indicate that despite the differences at recent times, the epoch of formation of ~50% of the stellar mass in dwarf galaxies was similar to that of ~50% of the stellar mass in large spiral galaxies (z>~2), despite the observed galaxy-to-galaxy diversity among the dwarfs.

We present the results of a 63 ks {\it Chandra} observation of the "supernova impostor" SN 2010da four months after it was first observed on 25 May 2010. We detect an X-ray source at $\sim7\sigma$ confidence coincident with the optical location of the SN 2010da outburst. Our imaging analysis has revealed a hard central point source, surrounded by soft diffuse emission extending as far as $\sim$8\asn north of the central source. The diffuse emission has a hardness ratio, 0.35-2 keV luminosity ($\sim6\times10^{35}$ erg s$^{-1}$), and size ($\sim20$ pc) consistent with that of a supernova remnant, although the low number of counts prohibits detailed spectral modeling. Our best-fit spectral model for the hard central source is a black body ($kT = 1.79^{+0.66}_{-0.43}$ keV) with no evidence for intrinsic absorption beyond the Galactic column. We estimate the 0.3-10 keV luminosity to be 1.7$^{+0.2}_{-0.5}\times10^{37}$ erg s$^{-1}$, a factor of $\sim$25 decrease since the initial outburst four months previously. The high X-ray luminosity and slow fading rate is not consistent with this object being a single massive star undergoing an outburst; instead, we favor the scenario where the massive star powering the SN 2010da optical transient is part of a wind-fed supergiant X-ray binary system with a compact companion powering the observed X-ray emission.

We present deep Hubble Space Telescope Wide Field Planetary Camera 2 (WFPC2) optical observations obtained as part of the ACS Nearby Galaxy Survey Treasury (ANGST) as well as early release Wide Field Camera 3 (WFC3) ultra-violet and infrared observations of the nearby dwarf starbursting galaxy NGC 4214. Our data provide a detailed example of how covering such a broad range in wavelength provides a powerful tool for constraining the physical properties of stellar populations. The deepest data reach the ancient red clump at M_F814W -0.2. All of the optical data reach the main sequence turnoff for stars younger than ~300 Myr, and the blue He burning sequence for stars younger than 500 Myr. The full CMD-fitting analysis shows that all three fields in our data set are consistent with ~75% of the stellar mass being older than 8 Gyr, in spite of showing a wide range in star formation rates at the present day. Thus, our results suggest that the scale length of NGC 4214 has remained relatively constant for many Gyr. As previously noted by others, we also find the galaxy has recently ramped up production, consistent with its bright UV luminosity and its population of UV-bright massive stars. In the central field we find UV point sources with F336W magnitudes as bright as -9.9. These are as bright as stars with masses of at least 52-56 M_sun and ages near 4 Myr in stellar evolution models. Assuming a standard IMF, our CMD is well-fitted by an increase in star formation rate beginning 100 Myr ago. The stellar populations of this late-type dwarf are compared with those of NGC 404, an early-type dwarf that is also the most massive galaxy in its local environment. The late-type dwarf appears to have a similar high fraction of ancient stars, suggesting that these dominant galaxies may form at early epochs even if they have low total mass and very different present-day morphologies.

This study presents the final source catalog of the Chandra ACIS Survey of M33 (ChASeM33). With a total exposure time of 1.4 Ms, ChASeM33 covers ~70% of the D25 isophote (R\approx4kpc) of M33 and provides the deepest, most complete, and detailed look at a spiral galaxy in X-rays. The source catalog includes 662 sources, reaches a limiting unabsorbed luminosity of ~2.4x10^(34) erg/s in the 0.35-8.0keV energy band, and contains source positions, source net counts, fluxes and significances in several energy bands, and information on source variability. The analysis challenges posed by ChASeM33 and the techniques adopted to address these challenges are discussed. To constrain the nature of the detected X-ray source, hardness ratios were constructed and spectra were fit for 254 sources, followup MMT spectra of 116 sources were acquired, and cross-correlations with previous X-ray catalogs and other multi-wavelength data were generated. Based on this effort, 183 of the 662 ChASeM33 sources could be identified. Finally, the luminosity function for the detected point sources as well as the one for the X-ray binaries in M33 is presented. The luminosity functions in the soft band (0.5-2.0 keV) and the hard band (2.0-8.0 keV) have a limiting luminosity at the 90% completeness limit of 4.0x10^(34) erg/s and 1.6x10^(35) erg/s (for D=817kpc), respectively, which is significantly lower than what was reported by previous X-ray binary population studies in galaxies more distant than M33. The resulting distribution is consistent with a dominant population of high mass X-ray binaries as would be expected for M33.

We report the serendipitous detection of a giant X-ray flare from the source 2XMM J043527.2-144301 during an XMM-Newton observation of the high latitude molecular cloud MBM20. The source has not been previously studied at any wavelength. The X-ray flux increases by a factor of more than 52 from quiescent state to peak of flare. A 2MASS counterpart has been identified (2MASS J04352724-1443017), and near-infrared colors reveal a spectral type of M8-M8.5 and a distance of (67\pm 13) pc, placing the source in front of MBM20. Spectral analysis and source luminosity are also consistent with this conclusion. The measured distance makes this object the most distant source (by about a factor of 4) at this spectral type detected in X-rays. The X-ray flare was characterized by peak X-ray luminosity of ~8.2E28 erg s-1 and integrated X-ray energy of ~2.3E32 erg. The flare emission has been characterized with a 2-temperature model with temperatures of ~10 and 46 MK (0.82 and 3.97 keV), and is dominated by the higher temperature component.

We present results from mid-IR spectroscopic observations of two young supernova remnants (SNRs) in the Large Magellanic Cloud (LMC) done with the {\it Spitzer Space Telescope}. We imaged SNRs B0509-67.5 and B0519-69.0 with {\it Spitzer} in 2005, and follow-up spectroscopy presented here confirms the presence of warm, shock heated dust, with no lines present in the spectrum. We use model fits to {\it Spitzer} IRS data to estimate the density of the postshock gas. Both remnants show asymmetries in the infrared images, and we interpret bright spots as places where the forward shock is running into material that is several times denser than elsewhere. The densities we infer for these objects depend on the grain omposition assumed, and we explore the effects of differing grain porosity on the model fits. We also analyze archival {\it XMM-Newton} RGS spectroscopic data, where both SNRs show strong lines of both Fe and Si, coming from ejecta, as well as strong O lines, which may come from ejecta or shocked ambient medium. We use model fits to IRS spectra to predict X-ray O line strengths for various grain models and values of the shock compression ratio. For 0509-67.5, we find that compact (solid) grain models require nearly all O lines in X-ray spectra to originate in reverse-shocked ejecta. Porous dust grains would lower the strength of ejecta lines relative to those arising in the shocked ambient medium. In 0519-69.0, we find significant evidence for a higher than standard compression ratio of 12, implying efficient cosmic-ray acceleration by the blast wave. A compact grain model is favored over porous grain models. We find that the dust-to-gas mass ratio of the ambient medium is significantly lower than what is expected in the ISM.

We present uniformly measured star formation histories (SFHs) of 60 nearby (D~4Mpc) dwarf galaxies based on CMDs of resolved stellar populations from images taken with HST as part of the ACS Nearby Galaxy Survey Treasury program (ANGST). This volume-limited sample contains 12 dSph/dE, 5 dwarf spiral, 28 dI, 12 transition, and 3 tidal dwarf galaxies. From the best fit SFHs we find three significant results: (1) the average dwarf galaxy formed >50% of its stars by z~2 and 60% of its stars by z~1, regardless of current morphological type; (2) the mean SFHs of dIs, dTrans, and dSphs are similar over most of cosmic time, with the clearest differences appearing during the most recent 1 Gyr; and (3) the mean values are inconsistent with simple SFH models, e.g., exponentially declining SFRs. The mean SFHs are in general agreement with the cosmic SFH, although we observe offsets near z~1 that could be evidence that low mass systems experienced delayed star formation relative to more massive galaxies. The sample shows a strong density-morphology relationship, such that the dSphs are less isolated than dIs. We find that the transition from a gas-rich to gas-poor galaxy cannot be solely due to internal mechanisms such as stellar feedback, and instead is likely the result of external mechanisms, e.g., ram pressure and tidal forces. The average transition dwarf is slightly less isolated and less gas-rich than the typical dI. Further, the transition dwarfs can be divided into two groups: interacting and gas-poor or isolated and gas-rich, suggesting two possible evolutionary pathways.

Characterization of the frequency response of coherent radiometric receivers is a key element in estimating the flux of astrophysical emissions, since the measured signal depends on the convolution of the source spectral emission with the instrument band shape. Laboratory Radio Frequency (RF) measurements of the instrument bandpass often require complex test setups and are subject to a number of systematic effects driven by thermal issues and impedance matching, particularly if cryogenic operation is involved. In this paper we present an approach to modeling radiometers bandpasses by integrating simulations and RF measurements of individual components. This method is based on QUCS (Quasi Universal Circuit Simulator), an open-source circuit simulator, which gives the flexibility of choosing among the available devices, implementing new analytical software models or using measured S-parameters. Therefore an independent estimate of the instrument bandpass is achieved using standard individual component measurements and validated analytical simulations. In order to automate the process of preparing input data, running simulations and exporting results we developed the Python package python-qucs and released it under GNU Public License. We discuss, as working cases, bandpass response modeling of the COFE and Planck Low Frequency Instrument (LFI) radiometers and compare results obtained with QUCS and with a commercial circuit simulator software. The main purpose of bandpass modeling in COFE is to optimize component matching, while in LFI they represent the best estimation of frequency response, since end-to-end measurements were strongly affected by systematic effects.

We present three transits of GJ 1214b, observed as part of the Apache Point Observatory Survey of Transit Lightcurves of Exoplanets (APOSTLE). We used APOSTLE r-band lightcurves in conjunction with previously gathered data of GJ 1214b to re-derive system parameters. By using parameters such as transit duration and ingress/egress length we are able to reduce the degeneracies between parameters in the fitted transit model, which is a preferred condition for Markov Chain Monte Carlo techniques typically used to quantify uncertainties in measured parameters. The joint analysis of this multi-wavelength dataset confirms earlier estimates of system parameters including planetary orbital period, the planet-to-star radius ratio and stellar density. We fit the photometric spectralenergy distribution of GJ 1214 to derive stellar luminosity, which we then use to derive its absolute mass and radius. From these derived stellar properties and previously published radial velocity data we were able to refine estimates of the absolute parameters for the planet GJ 1214b. Transit times derived from our study show no evidence for strong transit timing variations. Some lightcurves we present show features that we believe are due to stellar activity. During the first night we observed a rise in the out-of-eclipse flux of GJ 1214 with a characteristic fast-rise exponential decay shape commonly associated with stellar flares. On the second night we observed a minor brightening during transit, which we believe might have been caused by the planet obscuring a star-spot on the stellar disk.

In an attempt to constrain evolutionary models of the asymptotic giant branch (AGB) phase at the limit of low masses and low metallicities, we have examined the luminosity functions and number ratio between AGB and red giant branch (RGB) stars from a sample of resolved galaxies from the ACS Nearby Galaxy Survey Treasury (ANGST). This database provides HST optical photometry together with maps of completeness, photometric errors, and star formation histories for dozens of galaxies within 4 Mpc. We select 12 galaxies characterized by predominantly metal-poor populations as indicated by a very steep and blue RGB, and which do not present any indication of recent star formation in their color--magnitude diagrams. Thousands of AGB stars brighter than the tip of the RGB (TRGB) are present in the sample (between 60 and 400 per galaxy), hence the Poisson noise has little impact in our measurements of the AGB/RGB ratio. We model the photometric data with a few sets of thermally pulsing AGB (TP-AGB) evolutionary models with different prescriptions for the mass loss. This technique allows us to set stringent constraints to the TP-AGB models of low-mass metal-poor stars (with M<1.5 Msun, [Fe/H]<~-1.0). Indeed, those which satisfactorily reproduce the observed AGB/RGB ratios have TP-AGB lifetimes between 1.2 and 1.8 Myr, and finish their nuclear burning lives with masses between 0.51 and 0.55 Msun. This is also in good agreement with recent observations of white dwarf masses in the M4 old globular cluster. These constraints can be added to those already derived from Magellanic Cloud star clusters as important mileposts in the arduous process of calibrating AGB evolutionary models.

The starburst phenomenon can shape the evolution of the host galaxy and the surrounding intergalactic medium. The extent of the evolutionary impact is partly determined by the duration of the starburst, which has a direct correlation with both the amount of stellar feedback and the development of galactic winds, particularly for smaller mass dwarf systems. We measure the duration of starbursts in twenty nearby, ongoing, and "fossil" starbursts in dwarf galaxies based on the recent star formation histories derived from resolved stellar population data obtained with the Hubble Space Telescope. Contrary to the shorter times of 3-10 Myr often cited, the starburst durations we measure range from 450 - 650 Myr in fifteen of the dwarf galaxies and up to 1.3 Gyr in four galaxies; these longer durations are comparable to or longer than the dynamical timescales for each system. The same feedback from massive stars that may quench the flickering SF does not disrupt the overall burst event in our sample of galaxies. While five galaxies present fossil bursts, fifteen galaxies show ongoing bursts and thus the final durations may be longer than we report here for these systems. One galaxy shows a burst that has been ongoing for only 20 Myr; we are likely seeing the beginning of a burst event in this system. Using the duration of the starbursts, we calculate that the bursts deposited 10^(53.9)-10^(57.2) ergs of energy into the interstellar medium through stellar winds and supernovae and produced 3%-26% of the host galaxy's mass.

The power spectrum of density fluctuations is a foundational source of cosmological information. Precision cosmological probes targeted primarily at investigations of dark energy require accurate theoretical determinations of the power spectrum in the nonlinear regime. To exploit the observational power of future cosmological surveys, accuracy demands on the theory are at the one percent level or better. Numerical simulations are currently the only way to produce sufficiently error-controlled predictions for the power spectrum. The very high computational cost of (precision) N-body simulations is a major obstacle to obtaining predictions in the nonlinear regime, while scanning over cosmological parameters. Near-future observations, however, are likely to provide a meaningful constraint only on constant dark energy equation of state 'wCDM' cosmologies. In this paper we demonstrate that a limited set of only 37 cosmological models -- the "Coyote Universe" suite -- can be used to predict the nonlinear matter power spectrum at the required accuracy over a prior parameter range set by cosmic microwave background observations. This paper is the second in a series of three, with the final aim to provide a high-accuracy prediction scheme for the nonlinear matter power spectrum for wCDM cosmologies.

We use archival HST observations of resolved stellar populations to derive the star formation histories (SFHs) of 18 nearby starburst dwarf galaxies. In this first paper we present the observations, color-magnitude diagrams, and the SFHs of the 18 starburst galaxies, based on a homogeneous approach to the data reduction, differential extinction, and treatment of photometric completeness. We adopt a star formation rate (SFR) threshold normalized to the average SFR of the individual system as a metric for classifying starbursts in SFHs derived from resolved stellar populations. This choice facilitates finding not only currently bursting galaxies but also "fossil" bursts increasing the sample size of starburst galaxies in the nearby (D<8 Mpc) universe. Thirteen of the eighteen galaxies are experiencing ongoing bursts and five galaxies show fossil bursts. From our reconstructed SFHs, it is evident that the elevated SFRs of a burst are sustained for hundreds of Myr with variations on small timescales. A long >100 Myr temporal baseline is thus fundamental to any starburst definition or identification method. The longer lived bursts rule out rapid "self-quenching" of starbursts on global scales. The bursting galaxies' gas consumption timescales are shorter than the Hubble time for all but one galaxy confirming the short-lived nature of starbursts based on fuel limitations. Additionally, we find the strength of the H{\alpha} emission usually correlates with the CMD based SFR during the last 4-10 Myr. However, in four cases, the H{\alpha} emission is significantly less than what is expected for models of starbursts; the discrepancy is due to the SFR changing on timescales of a few Myr. The inherently short timescale of the H{\alpha} emission limits identifying galaxies as starbursts based on the current characteristics which may or may not be representative of the recent SFH of a galaxy.

The luminosity, light curve, post--maximum spectrum, and lack of a progenitor on deep pre-outburst images suggest that SN 2010U was a luminous, fast nova. Its outburst magnitude is consistent with that for a fast nova using the Maximum Magnitude-Rate of Decline relationship for classical novae.

Interstellar dust grains play a crucial role in the evolution of the galactic interstellar medium (ISM). Despite its importance, however, dust remains poorly understood in terms of its origin, composition, and abundance throughout the universe. Supernova remnants (SNRs) provide a laboratory for studying the evolution of dust grains, as they are one of the only environments in the universe where it is possible to observe grains being both created and destroyed. SNRs exhibit collisionally heated dust, allowing dust to serve as a diagnostic both for grain physics and for the plasma conditions in the SNR. I present theoretical models of collisionally heated dust which calculate grain emission as well as destruction rates. In these models, I incorporate physics such as nonthermal sputtering caused by grain motions through the gas, a more realistic approach to sputtering for small grains, and arbitrary grain compositions porous and composite grains. I apply these models to infrared and X-ray observations of Kepler's supernova and the Cygnus Loop in the galaxy, and SNRs 0509-67.5, 0519-69.0, and 0540-69.3 in the LMC. X-ray observations characterize the hot plasma while IR observations constrain grain properties and destruction rates. Such a multi-wavelength approach is crucial for a complete understanding of gas and dust interaction and evolution. Modeling of both X-ray and IR spectra allows disentangling of parameters such as pre and postshock gas density, as well as swept-up masses of gas and dust, and can provide constraints on the shock compression ratio. Observations also show that the dust-to-gas mass ratio in the ISM is lower by a factor of several than what is inferred by extinction studies of starlight. Future observatories, such as the James Webb Space Telescope and the International X-ray Observatory, will allow testing of models far beyond what is possible now.

We present HST/WFPC2 observations across the disk of the nearby isolated dwarf S0 galaxy NGC 404, which hosts an extended gas disk. Our deepest field reaches the red clump and main-sequence stars with ages <500 Myr. Although we detect trace amounts of star formation at times more recent than 10 Gyr for all fields, the proportion of red giant stars to asymptotic giants and main sequence stars suggests that the disk is dominated by an ancient (>10 Gyr) population. Detailed modeling of the color-magnitude diagram suggests that ~70% of the stellar mass in the NGC 404 disk formed by z~2 (10 Gyr ago) and at least ~90% formed prior to z~1 (8 Gyr ago). These results indicate that the stellar populations of the NGC 404 disk are on average significantly older than those of other nearby disk galaxies, suggesting that early and late type disks may have different long-term evolutionary histories, not simply differences in their recent star formation rates. Comparisons of the spatial distribution of the young stellar mass and FUV emission in GALEX images show that the brightest FUV regions contain the youngest stars, but that some young stars (<160 Myr) lie outside of these regions. FUV luminosity appears to be strongly affected by both age and stellar mass within individual regions. Finally, we use our measurements to infer the relationship between the star formation rate and the gas density of the disk at previous epochs. We find that most of the history of the NGC 404 disk is consistent with star formation that has decreased with the gas density according to the Schmidt law. However, 0.5-1 Gyr ago, the star formation rate was unusually low for the inferred gas density, consistent with the possibility that there was a gas accretion event that reignited star formation ~0.5 Gyr ago. Such an event could explain why this S0 galaxy hosts an extended gas disk.

(Abridged) Spitzer data at 24, 70, and 160 micron and ground-based H-alpha images are analyzed for a sample of 189 nearby star-forming and starburst galaxies to investigate whether reliable star formation rate (SFR) indicators can be defined using the monochromatic infrared dust emission centered at 70 and 160 micron. We compare recently published recipes for SFR measures using combinations of the 24 micron and observed H-alpha luminosities with those using 24 micron luminosity alone. From these comparisons, we derive a reference SFR indicator for use in our analysis. Linear correlations between SFR and the 70 and 160 micron luminosity are found for L(70)>=1.4x10^{42} erg/s and L(160)>=2x10^{42} erg/s, corresponding to SFR>=0.1-0.3 M_sun/yr. Below those two luminosity limits, the relation between SFR and 70 micron (160 micron) luminosity is non-linear and SFR calibrations become problematic. The dispersion of the data around the mean trend increases for increasing wavelength, becoming about 25% (factor ~2) larger at 70 (160) micron than at 24 micron. The increasing dispersion is likely an effect of the increasing contribution to the infrared emission of dust heated by stellar populations not associated with the current star formation. The non-linear relation between SFR and the 70 and 160 micron emission at faint galaxy luminosities suggests that the increasing transparency of the interstellar medium, decreasing effective dust temperature, and decreasing filling factor of star forming regions across the galaxy become important factors for decreasing luminosity. The SFR calibrations are provided for galaxies with oxygen abundance 12+Log(O/H)>8.1. At lower metallicity the infrared luminosity no longer reliably traces the SFR because galaxies are less dusty and more transparent.

M33 contains a large number of emission nebulae identified as supernova remnants (SNRs) based on the high [S II]:Ha ratios characteristic of shocked gas. Using Chandra data from the ChASeM33 survey with a 0.35-2 keV sensitivity of about 2 x 10**34 ergs/s, we have detected 82 of 137 SNR candidates, yielding confirmation of (or at least strongly support for) their SNR identifications. This provides the largest sample of remnants detected at optical and X-ray wavelengths in any galaxy, including the Milky Way. A spectral analysis of the seven X-ray brightest SNRs reveals that two, G98-31 and G98-35, have spectra that appear to indicate enrichment by ejecta from core-collapse supernova explosions. In general, the X-ray detected SNRs have soft X-ray spectra compared to the vast majority of sources detected along the line of sight to M33. It is unlikely that there are any other undiscovered thermally dominated X-ray SNRs with luminosities in excess of about 4 x 10**35 ergs/s in the portions of M33 covered by the ChASeM33 survey. We have used a combination of new and archival optical and radio observations to attempt to better understand why some objects are detected as X-ray sources and others are not. We have also developed a morphological classification scheme for the optically-identified SNRs, and discuss the efficacy of this scheme as a predictor of X-ray detectability. Finally, we have compared the SNRs found in M33 to those that have been observed in the Galaxy and the Magellanic Clouds. There are no close analogs of Cas A, Kepler's SNR, Tycho's SNR or the Crab Nebula in the regions of M33 surveyed, but we have found an X-ray source with a power law spectrum coincident with a small-diameter radio source that may be the first pulsar-wind nebula recognized in M33.

We present new Hubble Space Telescope (HST) observations of NGC 300 taken as part of the ACS Nearby Galaxy Survey Treasury (ANGST). Individual stars are resolved in these images down to an absolute magnitude of M_F814W = 1.0 (below the red clump). We determine the star formation history of the galaxy in 6 radial bins by comparing our observed color-magnitude diagrams (CMDs) with synthetic CMDs based on theoretical isochrones. We find that the stellar disk out to 5.4 kpc is primarily old, in contrast with the outwardly similar galaxy M33. We determine the scale length as a function of age and find evidence for inside-out growth of the stellar disk: the scale length has increased from 1.1 +/- 0.1 kpc 10 Gyr ago to 1.3 +/- 0.1 kpc at present, indicating a buildup in the fraction of young stars at larger radii. As the scale length of M33 has recently been shown to have increased much more dramatically with time, our results demonstrate that two galaxies with similar sizes and morphologies can have very different histories. With an N-body simulation of a galaxy designed to be similar to NGC 300, we determine that the effects of radial migration should be minimal. We trace the metallicity gradient as a function of time and find a present day metallicity gradient consistent with that seen in previous studies. Consistent results are obtained from archival images covering the same radial extent but differing in placement and filter combination.

We investigate the utility of the asymptotic giant branch (AGB) and the red giant branch (RGB) as probes of the star formation history (SFH) of the nearby (D=2.5 Mpc) dwarf irregular galaxy, KKH 98. Near-infrared (IR) Keck Laser Guide Star Adaptive Optics (AO) images resolve 592 IR bright stars reaching over 1 magnitude below the Tip of the Red Giant Branch. Significantly deeper optical (F475W and F814W) Hubble Space Telescope images of the same field contain over 2500 stars, reaching to the Red Clump and the Main Sequence turn-off for 0.5 Gyr old populations. Compared to the optical color magnitude diagram (CMD), the near-IR CMD shows significantly tighter AGB sequences, providing a good probe of the intermediate age (0.5 - 5 Gyr) populations. We match observed CMDs with stellar evolution models to recover the SFH of KKH 98. On average, the galaxy has experienced relatively constant low-level star formation (5 x 10^-4 Mo yr^-1) for much of cosmic time. Except for the youngest main sequence populations (age < 0.1 Gyr), which are typically fainter than the AO data flux limit, the SFH estimated from the the 592 IR bright stars is a reasonable match to that derived from the much larger optical data set. Differences between the optical and IR derived SFHs for 0.1 - 1 Gyr populations suggest that current stellar evolution models may be over-producing the AGB by as much as a factor of three in this galaxy. At the depth of the AO data, the IR luminous stars are not crowded. Therefore these techniques can potentially be used to determine the stellar populations of galaxies at significantly further distances.

We announce the discovery of a quasar behind the disk of M31, which was previously classified as a remarkable nova in our neighbour galaxy. The paper is primarily aimed at the outburst of J004457+4123 (Sharov 21), with the first part focussed on the optical spectroscopy and the improvement in the photometric database. Both the optical spectrum and the broad band spectral energy distribution of Sharov 21 are shown to be very similar to that of normal, radio-quiet type 1 quasars. We present photometric data covering more than a century and resulting in a long-term light curve that is densely sampled over the past five decades. The variability of the quasar is characterized by a ground state with typical fluctuation amplitudes of ~0.2 mag around B~20.5, superimposed by a singular flare of ~2 yr duration (observer frame) with the maximum at 1992.81 where the UV flux has increased by a factor of ~20. The total energy in the flare is at least three orders of magnitudes higher than the radiated energy of the most luminous supernovae, provided that it comes from an intrinsic process and the energy is radiated isotropically. The profile of the flare light curve appears to be in agreement with the standard predictions for a stellar tidal disruption event where a ~10 M_sun giant star was shredded in the tidal field of a ~2...5 10^8 M_sun black hole. The short fallback time derived from the light curve requires an ultra-close encounter where the pericentre of the stellar orbit is deep within the tidal disruption radius. Gravitational microlensing provides an alternative explanation, though the probability of such a high amplification event is very low.