Search

We present evidence for $\gamma$-ray emission from a stacked population of 39 high-latitude globular clusters (GCs) not detected in the Fermi Point Source Catalog, likely attributable to populations of millisecond pulsars within them. In this work, we use 13 years of data collected by the Large Area Telescope aboard the Fermi Gamma-Ray Space Telescope to search for a cumulative signal from undetected GCs and compared them to control fields (CFs), selected to match the celestial distribution of the target clusters so as to distinguish the $\gamma$-ray signal from background emission. The joint likelihood distribution of the GCs has a significant separation ($\sim4\sigma$) from that of the CFs. We also investigate correlations between detected cluster luminosities and other cluster properties such as distance, the number of millisecond pulsars associated with each cluster, and stellar encounter rate but find no significant relationships., Comment: 9 pages, 9 figures

A complete accounting of nearby objects -- from the highest-mass white dwarf progenitors down to low-mass brown dwarfs -- is now possible, thanks to an almost complete set of trigonometric parallax determinations from Gaia, ground-based surveys, and Spitzer follow-up. We create a census of objects within a Sun-centered sphere of 20-pc radius and check published literature to decompose each binary or higher-order system into its separate components. The result is a volume-limited census of $\sim$3,600 individual star formation products useful in measuring the initial mass function across the stellar ($<8 M_\odot$) and substellar ($\gtrsim 5 M_{Jup}$) regimes. Comparing our resulting initial mass function to previous measurements shows good agreement above 0.8$M_\odot$ and a divergence at lower masses. Our 20-pc space densities are best fit with a quadripartite power law, $\xi(M) = dN/dM \propto M^{-\alpha}$ with long-established values of $\alpha = 2.3$ at high masses ($0.55 < M < 8.00 M_\odot$) and $\alpha = 1.3$ at intermediate masses ($0.22 < M < 0.55 M_\odot$), but at lower masses we find $\alpha = 0.25$ for $0.05 < M <0.22 M_\odot$ and $\alpha = 0.6$ for $0.01 < M < 0.05 M_\odot$. This implies that the rate of production as a function of decreasing mass diminishes in the low-mass star/high-mass brown dwarf regime before increasing again in the low-mass brown dwarf regime. Correcting for completeness, we find a star to brown dwarf number ratio of, currently, 4:1, and an average mass per object of 0.41 $M_\odot$., Comment: Accepted for publication in The Astrophysical Journal Supplement Series. 123 pages with four ancillary files

The existence of a vast nova shell surrounding the prototypical dwarf nova Z Camelopardalis (Z Cam) proves that some old novae undergo metamorphosis to appear as dwarf novae thousands of years after a nova eruption. The expansion rates of ancient nova shells offer a way to constrain both the time between nova eruptions and the time for post-nova mass transfer rates to decrease significantly, simultaneously testing nova thermonuclear runaway models and hibernation theory. Previous limits on the expansion rate of part of the Z Cam shell constrain the inter-eruption time between Z Cam nova events to be $>$ 1300 years. Deeper narrow-band imaging of the ejecta of Z Cam with the Condor Array Telescope now reveals very low surface brightness areas of the remainder of the shell. A second, even fainter shell is also detected, concentric with and nearly three times the size of the "inner" shell. This is the first observational support of the prediction that concentric shells must surround the frequently-erupting novae of relatively massive white dwarfs. The Condor images extend our Z Cam imaging baseline to 15 years, yielding the inner shell's expansion rate as $v = 83 \pm 37$ km s$^{-1}$ at 23 degrees South of West, in excellent agreement with our 2012 prediction. This velocity corresponds to an approximate age of $t = 2672^{-817}_{+2102}$ yr. While consistent with the suggestion that the most recent nova eruption of Z Cam was the transient recorded by Chinese Imperial astrologers in the year 77 BCE, the age uncertainty is still too large to support or disprove a connection with Z Cam., Comment: 13 pages, 9 figures

M87 has been monitored with a cadence of 5 days over a 9 month-long span through the near-ultraviolet (NUV:F275W) and optical (F606W) filters of the Wide Field Camera 3 (WFC3) of the $\textit{Hubble Space Telescope}$. This unprecedented dataset yields the NUV and optical light and color curves of 94 M87 novae, characterizing the outburst and decline properties of the largest extragalactic nova dataset in the literature (after M31 and M81). We test and confirm nova modelers' prediction that recurrent novae cannot erupt more frequently that once every 45 days; show that there are zero rapidly recurring novae in the central $\sim$ 1/3 of M87 with recurrence times $ < $ 130 days; demonstrate that novae closely follow the K-band light of M87 to within a few arcsec of the galaxy nucleus; show that nova NUV light curves are as heterogeneous as their optical counterparts, and usually peak 5 to 30 days after visible light maximum; determine our observations' annual detection completeness to be 71 - 77\%; and measure the rate Rnova of nova eruptions in M87 as $352_{-37}^{+37}$/yr. The corresponding luminosity-specific classical nova rate for this galaxy is $7.91_{-1.20}^{+1.20}/yr/10^{10}L_\odot,_{K}$. These rates confirm that ground-based observations of extragalactic novae miss most faint, fast novae and those near the centers of galaxies. An annual M87 nova rate of 300 or more seems inescapable. A luminosity-specific nova rate of $\sim$ $7 - 10/yr/10^{10}L_\odot,_{K}$ in ${\it all}$ types of galaxies is indicated by the data available in 2023., Comment: Accepted/In Press in ApJS; 3 Tables, 108 Figures, 180 pages

Accurately quantifying the rates dM/dt at which massive stars lose mass is essential to any understanding of their evolution. All dM/dt estimates to date assume wind clumping factors; not allowing for clumping leads to overestimates of dM/dt and underestimates of lifetimes and masses when these stars explode as supernovae. Mid-IR spectroscopy suggested that the wind of the nearest Wolf-Rayet star, Gamma2 Vel, is resolved with a Full Width at 10 per cent intensity of 0.5 arcsec, or 171 AU at the 342 pc distance of the star. As the Zorro speckle imager on Gemini-South is capable of 0.02 arcsec resolution, we have used it to image Gamma2 Vel at two orbital phases (0.30 and 0.44) with two narrowband and two intermediate-band filters in an attempt to resolve its wind. Our observations demonstrate that the wind of Gamma2 Vel may be resolved as a 0.07 arcsec westward elongation through an 832 nm filter at orbital phase 0.3 . If confirmed, this is the smallest scale (24 AU) at which a WR star wind asymmetry has been directly imaged. Similar imaging at multiple phases is needed to determine if the asymmetry is due to stochastic wind clumping, co-rotating interaction regions or colliding-wind, cone-shaped shocks., Comment: 7 pages, 4 Figures, in press in MNRAS

Accurate determination of the rates of nova eruptions in different kinds of galaxies give us strong constraints on those galaxies' underlying white dwarf and binary populations, and those stars' spatial distributions. Until 2016, limitations inherent in ground-based surveys of external galaxies - and dust extinction in the Milky Way - significantly hampered the determination of those rates and how much they differ between different types of galaxies. Infrared Galactic surveys and dense cadence Hubble Space Telescope (HST)-based surveys are overcoming these limitations, leading to sharply increased nova-in-galaxy rates relative to those previously claimed. Here we present 14 nova candidates that were serendipitously observed during a year-long HST survey of the massive spiral galaxy M51 (the "Whirlpool Galaxy"). We use simulations based on observed nova light curves to model the incompleteness of the HST survey in unprecedented detail, determining a nova detection efficiency $\epsilon = 20.3$ percent. The survey's M51 area coverage, combined with $\epsilon$, indicates a conservative M51 nova rate of $172^{+46}_{-37}$ novae yr$^{-1}$, corresponding to a luminosity-specific nova rate (LSNR) of $\sim10.4^{+2.8}_{-2.2}$ novae yr$^{-1}$/$10^{10} L_{\odot,K}$. Both these rates are approximately an order of magnitude higher than those estimated by ground-based studies, contradicting claims of universal low nova rates in all types of galaxies determined by low cadence, ground-based surveys. They demonstrate that, contrary to theoretical models, the HST-determined LSNR in a giant elliptical galaxy (M87) and a giant spiral galaxy (M51) likely do not differ by an order of magnitude or more, and may in fact be quite similar.

We present a far-ultraviolet (FUV) study of the globular cluster M30 (NGC 7099). The images were obtained using the Advanced Camera for Surveys (ACS/SBC, F150LP, FUV) and the Wide Field Planetary Camera 2 (WFPC2, F300W, UV) on board the Hubble Space Telescope (HST). We compare the catalogue of FUV objects to ten known X-ray sources and find six confident matches of two cataclysmic variables (CVs), one RS CVn, one red giant with strong FUV emission and two sources only detected in the FUV. We also searched for variable sources in our dataset and found a total of seven blue stragglers (BSs), four horizontal branch (HB) stars, five red giant branch stars, 28 main sequence stars and four gap objects that demonstrated variability. One BS star is a known W-UMa contact binary, one of the gap objects is a known CV identified in this work to be a dwarf nova, and the three other gap sources are weak variables. The periods and positions of two of the variable HB stars match them to two previously known RR Lyrae variables of types RRab and RRc., Comment: 13 pages, 12 figures, accepted for publication to MNRAS

We present a far-ultraviolet (FUV) study of the globular cluster M30 (NGC 7099). The images were obtained using the Advanced Camera for Surveys (ACS/SBC, F150LP, FUV) and the Wide Field Planetary Camera 2 (WFPC2, F300W, UV) which were both on board the Hubble Space Telescope (HST). The FUV-UV colour-magnitude diagram (CMD) shows a main sequence (MS) turnoff at FUV $\approx$ 22 mag and FUV-UV $\approx$ 3 mag. The MS extends 4 mag below the turnoff, and a prominent horizontal branch (HB) and blue straggler (BS) sequence can be seen. A total of 1218 MS stars, 185 red giant branch stars, 47 BS stars and 41 HB stars are identified, along with 78 sources blueward of the MS which consist of white dwarfs (WDs) and objects in the gap between the WDs and the MS that include potential cataclysmic variable (CV) candidates. The radial distribution of the BS population is concentrated towards the cluster centre, indicating that mass segregation has occurred. The blue and red sub-populations of the double BS sequence appear mixed in the ultraviolet CMD, and no significant central concentration of CV candidates is seen in this cluster., Comment: 11 pages, 7 figures. Accepted for publication in MNRAS

Gravitationally bound companions to stars enable determinations of their masses, and offer clues to their formation, evolution and dynamical histories. So motivated, we have carried out a speckle imaging survey of eight of the nearest and brightest Wolf-Rayet (WR) stars to directly measure the frequency of their resolvable companions, and to search for much fainter companions than hitherto possible. We found one new, close companion to each of WR 113, WR 115 and WR 120 in the separation range 0.2" - 1.2". Our results provide more evidence that similar-brightness, close companions to WR stars are common. More remarkably, they also demonstrate that the predicted, but much fainter and thus elusive companions to WR stars are now within reach of modern speckle cameras on 8m class telescopes by finding the first example. The new companion to WR 113 is just 1.16" distant from it, and is 8 magnitudes fainter than the WR star. The empirical probability of a chance line-of-sight of the faint companion at the position of WR 113 is < 0.5%, though we cannot yet prove or disprove if the two stars are gravitationally bound. If these three new detections are physical companions we suggest, based on their narrowband magnitudes, colors, reddenings and GAIA distances that the companions to WR113, WR 115 and WR 120 are an F-type dwarf, an early B-type dwarf, and a WNE-type WR star, respectively., Comment: 14 pages, 4 figures, accepted for publication in MNRAS

The orbital periods of most eclipsing cataclysmic binaries are not undergoing linear secular decreases of order a few parts per billion as expected from simple theory. Instead, they show several parts per million increases and decreases on timescales of years to decades, ascribed to magnetic effects in their donors, triple companions, or both. To directly test the triple companion hypothesis, we carried out a speckle imaging survey of six of the nearest and brightest cataclysmic variables. We found no main sequence companions earlier than spectral types M4V in the separation range 0.02" - 1.2", corresponding to projected linear separations of 2 - 100 AU, and periods of 3 - 1000 years. We conclude that main sequence triple companions to CVs are not very common, but cannot rule out the presence of the faintest M dwarfs or close brown dwarf companions., Comment: 8 pages, 1 figure, accepted in MNRAS

Some massive, merging black holes (BH) may be descendants of binary O stars. The evolution and mass transfer between these O stars determines the spins of their progeny BH. These will be measurable with future gravitational wave detectors, incentivizing the measurement of the spins of O stars in binaries. We previously measured the spins of O stars in Galactic Wolf-Rayet (WR) + O binaries. Here we measure the vsini of four LMC and two SMC O stars in WR + O binaries to determine whether lower metallicity might affect the spin rates. We find that the O stars in Galactic and Magellanic WR + O binaries display average vsini = 258 +/- 18 km/s and 270 +/- 15 km/s, respectively. Two LMC O stars measured on successive nights show significant line width variability, possibly due to differing orbital phases exhibiting different parts of the O stars illuminated differently by their WR companions. Despite this variability, the vsini are highly super-synchronous but distinctly subcritical for the O stars in all these binaries; thus we conclude that an efficient mechanism for shedding angular momentum from O stars in WR + O binaries must exist. This mechanism, probably related to Roche lobe overflow-created dynamo magnetic fields, prevents nearly 100% breakup spin rates, as expected when RLOF operates, as it must, in these stars. A Spruit-Tayler dynamo and O star wind might be that mechanism., Comment: 8 pages, 13 figures, accepted in MNRAS

We present an analysis of the second epoch HST WFC3 F110W near-Infrared (NIR) imaging data of the globular cluster M4. The new dataset suggests that one of the previously suggested four brown dwarf candidates in this cluster is indeed a high-probability cluster member. The position of this object in the NIR colour magnitude diagrams (CMDs) is in the white dwarf/brown dwarf area. The source is too faint to be a low-mass main sequence star, but, according to theoretical considerations, also most likely somewhat too bright to be a bona-fide brown dwarf. Since we know that the source is a cluster member, we determined a new optical magnitude estimate at the position the source should have in the optical image. This new estimate places the source closer to the white dwarf sequence in the optical-NIR CMD and suggests that it might be a very cool (T_eff < 4500 K) white dwarf at the bottom of the white dwarf cooling sequence in M4, or a white dwarf/brown dwarf binary. We cannot entirely exclude the possibility that the source is a very massive, bright brown dwarf, or a very low-mass main sequence star, however, we conclude that we still have not convincingly detected a brown dwarf in a globular cluster, but we expect to be very close to the start of the brown dwarf cooling sequence in this cluster. We also note that the main sequence ends at F110W approx.22.5 mag in the proper-motion cleaned CMDs, where completeness is still high., Comment: 13 pages, including 7 figures. Accepted for publication in MNRAS

The giant elliptical galaxy M87 has been imaged over 30 consecutive days in 2001, 60 consecutive days in 2005-2006, and every 5 days over a 265 day span in 2016-2017 with the Hubble Space Telescope, leading to the detection of 137 classical novae throughout M87. We have identified 2134 globular clusters (GC) in M87 in these images, and carried out searches of the clusters for classical novae erupting in or near them. One GC CN was detected in the 2001 data, while zero novae were found during the 2005-2006 observations. Four candidate GC novae were (barely) detected in visible light during the 2016-2017 observations, but none of the four were seen in near-ultraviolet light, leading us to reject them. Combining these results with our detection of one M87 GC nova out of a total of 137 detected CN, we conclude that such novae may be overabundant relative to the field, but small number statistics dominate this (and all other) searches. A definitive determination of GC CN overabundance (or not) will require much larger samples which LSST should provide in the coming decade., Comment: 11 pages, 2 figures

Late-time observations of Type Ia supernovae (SNe Ia), >900 days after explosion, have shown that this type of SN does not suffer an "IR catastrophe" at 500 days as previously predicted. Instead, several groups have observed a slow-down in the optical light curves of these SNe. A few reasons have been suggested for this slow-down, from a changing fraction of positrons reprocessed by the expanding ejecta, through a boost of energy from slow radioactive decay chains such as 57Co --> 57Fe, to atomic "freeze-out." Discovering which of these (or some other) heating mechanisms is behind the slow-down will directly impact studies of SN Ia progenitors, explosion models, and nebular-stage physics. Recently, Graur et al. (2018) suggested a possible correlation between the shape of the late-time light curves of four SNe Ia and their stretch values, which are proxies for their intrinsic luminosities. Here, we present Hubble Space Telescope observations of the SN Ia ASASSN-14lp at ~850-960 days past maximum light. With a stretch of s = 1.15 +/- 0.05, it is the most luminous normal SN Ia observed so far at these late times. We rule out contamination by light echoes and show that the late-time, optical light curve of ASASSN-14lp is flatter than that of previous SNe Ia observed at late times. This result is in line with-and strengthens-the Graur et al. (2018) correlation, but additional SNe are needed to verify it., Comment: Accepted for publication in ApJ; 9 pages, 4 figures, 1 table

The late-time light curves of Type Ia supernovae (SNe Ia), observed $>900$ days after explosion, present the possibility of a new diagnostic for SN Ia progenitor and explosion models. First, however, we must discover what physical process (or combination of processes) leads to the slow-down of the late-time light curve relative to a pure $^{56}$Co decay, as observed in SNe 2011fe, 2012cg, and 2014J. We present Hubble Space Telescope observations of SN 2015F, taken $\approx 600-1040$ days past maximum light. Unlike those of the three other SNe Ia, the light curve of SN 2015F remains consistent with being powered solely by the radioactive decay of $^{56}$Co. We fit the light curves of these four SNe Ia in a consistent manner and measure possible correlations between the light curve stretch - a proxy for the intrinsic luminosity of the SN - and the parameters of the physical model used in the fit (e.g., the mass ratio of $^{56}$Co and $^{57}$Co produced in the explosion, or the time at which freeze-out sets in). We propose a new, late-time Phillips-like correlation between the stretch of the SNe and the shape of their late-time light curves, which we parametrize as the difference between their pseudo-bolometric luminosities at 600 and 900 days: $\Delta L_{900} = {\rm log}(L_{600}/L_{900})$. This model-independent correlation provides a new way to test which physical process lies behind the slow-down of SN Ia light curves $>900$ days after explosion, and, ultimately, fresh constraints on the various SN Ia progenitor and explosion models., Comment: Accepted for publication in ApJ (14 pages, 6 figures, 3 tables)

We have used HST/COS to obtain a series of 4 FUV (915-2148A) spectroscopic observations of the prototypical dwarf nova U Geminorum during its cooling following a two-week outburst. Our FUV spectral analysis of the data indicates that the white dwarf (WD) cools from a temperature of 41,500 K, 15 days after the peak of the outburst, to 36,250 K, 56 days after the peak of the outburst, assuming a massive WD (log(g)=8.8) and a distance of 100.4 pc. These results are self-consistent with a 1.1 solar mass WD with a 5,000 km radius. The spectra show many absorption lines of but no emission features. We find supra-solar abundances of nitrogen confirming the anomalous high N/C ratio. The FUV lightcurve reveals a 5% modulation with the orbital phase, showing dips near phase 0.25 and 0.75, where the spectra exhibit an increase in the depth of some absorption lines and in particular strong absorption lines from Si, Al, and Ar. The phase dependence we observe is consistent with material overflowing the disk rim at the hot spot, reaching a maximum elevation near phase 0.75, falling back at smaller radii near phase 0.5 where it bounces off the disk surface and again rising above the disk near phase 0.25. There is a large scatter in the absorption lines' velocities, especially for the silicon lines, while the carbon lines seem to match more closely the orbital velocity of the WD. This indicates that many absorption lines are affected by- or form in- the overflowing stream material veiling the WD, making the analysis of the WD spectra more difficult., Comment: Accepted for publication in The Astrophysical Journal, This is the arXiv/astroph copy which might be slightly different than the published ApJ version due to formatting

We present an observational far-UV (FUV) and near-UV (NUV) study of the core region of the globular cluster NGC 6397. The observations were obtained with the Space Telescope Imaging Spectrograph (STIS, FUV), and the Wide Field Camera 3 (WFC3, NUV) on board the Hubble Space Telescope. Here, we focus on the UV bright stellar populations such as blue stragglers (BSs), white dwarfs (WDs) and cataclysmic variables (CVs). We present the first FUV-NUV color-magnitude diagram (CMD) for this cluster. To support our classification of the stellar populations, we compare our FUV-NUV CMD with optical data from the ACS Survey of Galactic Globular Clusters. The FUV-NUV CMD indicates 16 sources located in the WD area, and ten BSs within the 25"x 25" of the STIS FUV data. Eighteen Chandra X-ray sources are located within the FUV field of view. Thirteen of those have a NUV counterpart, of which nine sources also have a FUV counterpart. Out of those, five sources are previously suggested CVs, and indeed all five are located in the WD/CV region in our FUV-NUV CMD. Another CV only has a FUV but no NUV counterpart. We also detect a NUV (but no FUV) counterpart to the MSP located in the core of this cluster. The NUV lightcurves of the CVs and MSP show flickering behaviour typical of CVs. We found that the BSs and CVs are the most centrally concentrated population. This might be an effect of mass segregation or indicate the preferred birth place of BSs and CVs via dynamical interactions in the dense core region of GCs. HB stars are the least centrally concentrated population and absent in the innermost area of the core., Comment: 13 pages, MNRAS accepted

Here we report the recovery of the binary underlying the classical nova of 11 March 1437 A.D. whose age is independently confirmed by proper motion-dating, and show that in the 20th century it exhibits dwarf nova eruptions. The four oldest recovered classical novae are now all dwarf novae. Taken together they strongly suggest that mass transfer rates decrease by an order of magnitude or more in the centuries after a classical nova event, greatly slowing the evolution, and lengthening the lifetimes of these explosive binary stars., Comment: 30 pages, 7 Figures

Sub-subgiants are stars observed to be redder than normal main-sequence stars and fainter than normal subgiant (and giant) stars in an optical color-magnitude diagram. The red straggler stars, which lie redward of the red giant branch, may be related and are often grouped together with the sub-subgiants in the literature. These stars defy our standard theory of single-star evolution, and are important tests for binary evolution and stellar collision models. In total, we identify 65 sub-subgiants and red stragglers in 16 open and globular star clusters from the literature; 50 of these, including 43 sub-subgiants, pass our strict membership selection criteria (though the remaining sources may also be cluster members). In addition to their unique location on the color-magnitude diagram, we find that at least 58% (25/43) of sub-subgiants in this sample are X-ray sources with typical 0.5-2.5 keV luminosities of order 10^30 - 10^31 erg/s. Their X-ray luminosities and optical-to-X-ray flux ratios are similar to those of RS CVn active binaries. At least 65% (28/43) of the sub-subgiants in our sample are variables, 21 of which are known to be radial-velocity binaries. Typical variability periods are <15 days. At least 33% (14/43) of the sub-subgiants are H-alpha emitters. These observational demographics provide strong evidence that binarity is important for sub-subgiant formation. Finally, we find that the number of sub-subgiants per unit mass increases toward lower-mass clusters, such that the open clusters in our sample have the highest specific frequencies of sub-subgiants., Comment: 21 pages, 7 figures, 4 tables, accepted for publication in ApJ

Ten weeks of daily imaging of the giant elliptical galaxy M87 with the Hubble Space Telescope (HST) has yielded 41 nova light curves of unprecedented quality for extragalactic cataclysmic variables. We have recently used these light curves to demonstrate that the observational scatter in the so-called Maximum-Magnitude Rate of Decline (MMRD) relation for classical novae is so large as to render the nova-MMRD useless as a standard candle. Here we demonstrate that a modified Buscombe - de Vaucouleurs hypothesis, namely that novae with decline times t2 > 10 days converge to nearly the same absolute magnitude about two weeks after maximum light in a giant elliptical galaxy, is supported by our M87 nova data. For 13 novae with daily-sampled light curves, well determined times of maximum light in both the F606W and F814W filters, and decline times $t2 > 10 days we find that M87 novae display M(606W,15) = -6.37 +/- 0.46 and M(814W,15) = -6.11 +/- 0.43. If very fast novae with decline times t2 < 10 days are excluded, the distances to novae in elliptical galaxies with stellar binary populations similar to those of M87 should be determinable with 1 sigma accuracies of +/-20% with the above calibrations., Comment: 8 pages, 4 figures, accepted in Monthly Notices of the Royal Astronomical Society

The extensive grid of numerical simulations of nova eruptions of Yaron et al.(2005) first predicted that some classical novae might deviate significantly from the Maximum Magnitude - Rate of Decline (MMRD) relation, which purports to characterise novae as standard candles. Kasliwal et al. (2011) have announced the observational detection of an apparently new class of faint, fast classical novae in the Andromeda galaxy. These objects deviate strongly from the MMRD relationship, exactly as predicted by Yaron et al. (2005). Shara et al. (2016) recently reported the first detections of faint, fast novae in M87. These previously overlooked objects are as common in the giant elliptical galaxy M87 as they are in the giant spiral M31; they comprise about 40% of all classical nova eruptions and greatly increase the observational scatter in the MMRD relation. We use the extensive grid of nova simulations of Yaron et al. (2005) to identify the underlying causes of the existence of faint, fast novae. These are systems which have accreted, and can thus eject, only very low mass envelopes, of order 10^-7 - 10^-8 Msun, on massive white dwarfs. Such binaries include, but are not limited to, the recurrent novae. These same models predict the existence of ultrafast novae which display decline times t2 as short as five hours. We outline a strategy for their future detection., Comment: 22 pages, 6 figures; submitted to the Astrophysical Journal

We present and discuss initial selection criteria and first results in M33 from a systematic search for extragalactic symbiotic stars. We show that the presence of diffuse interstellar gas emission can significantly contaminate the spectra of symbiotic star candidates. This important effect forces upon us a more stringent working definition of an extragalactic symbiotic star. We report the first detections and spectroscopic characterisation of 12 symbiotic binaries in M33. We found that four of our systems contain carbon-rich giants. In another two of them the giant seems to be a Zr-enhanced MS star, while the remaining six objects host M-type giants. The high number ratio of C to M giants in these binaries is consistent with the low metallicity of M33. The spatial and radial velocity distributions of these new symbiotic binaries are consistent with a wide range of progenitor star ages., Comment: 13 pages, 15 figures, submitted to MNRAS

Asynchronous polars (APs) are accreting white dwarfs (WDs) that have different WD and orbital angular velocities, unlike the rest of the known polars, which rotate synchronously (i.e., their WD and orbital angular velocities are the same). Past nova eruptions are the predicted cause of the asynchronicity, in part due to the fact that one of the APs, V1500 Cyg, was observed to undergo a nova eruption in 1975. We used the Southern African Large Telescope 10m class telescope and the MDM 2.4m Hiltner telescope to search for nova shells around three of the remaining four APs (V1432 Aql, BY Cam, and CD Ind) as well as one Intermediate Polar with a high asynchronicity (EX Hya). We found no evidence of nova shells in any of our images. We therefore cannot say that any of the systems besides V1500 Cyg had nova eruptions, but because not all post-nova systems have detectable shells, we also cannot exclude the possibility of a nova eruption occurring in any of these systems and knocking the rotation out of sync., Comment: 7 pages, 4 figures, to be published in MNRAS

The Hubble Space Telescope has imaged the central part of M87 over a 10 week span, leading to the discovery of 32 classical novae and nine fainter, likely very slow and/or symbiotic novae. In this first in a series of papers we present the M87 nova finder charts, and the light and color curves of the novae. We demonstrate that the rise and decline times, and the colors of M87 novae are uncorrelated with each other and with position in the galaxy. The spatial distribution of the M87 novae follows the light of the galaxy, suggesting that novae accreted by M87 during cannibalistic episodes are well-mixed. Conservatively using only the 32 brightest classical novae we derive a nova rate for M87: $363_{-45}^{+33}$ novae/yr. We also derive the luminosity-specific classical nova rate for this galaxy, which is $7.88_{-2.6}^{+2.3} /yr/ 10^{10}L_\odot,_{K}$. Both rates are 3-4 times higher higher than those reported for M87 in the past, and similarly higher than those reported for all other galaxies. We suggest that most previous ground-based surveys for novae in external galaxies, including M87, miss most faint, fast novae, and almost all slow novae near the centers of galaxies., Comment: 116 pages, 17 figures, submitted to the Astrophysical Journal Supplement Series Minor updates, conclusions unchanged

We present the results of initial spectrographic followup with the Very Large Telescope (UT3, Melipal) for $K_s \ge 14$ Galactic plane CIV emission-line candidates in the near-infrared (NIR). These 7 faint stars all display prominent HeI and CIV emission lines characteristic of a carbon-rich Wolf-Rayet star. They have NIR colours which are much too blue to be those of distant, classical WR stars. The magnitudes and colours are compatible with those expected for central stars of planetary nebulae, and are likely to come from massive progenitor populations, and themselves be more massive than any sample of planetary nebulae known. Our survey has identified thousands of such candidates., Comment: 9 pages, 10 figures

The black holes (BH) in merging BH-BH binaries are likely progeny of binary O stars. Their properties, including their spins, will be strongly influenced by the evolution of their progenitor O stars. The remarkable observation that many single O stars spin very rapidly can be explained if they accreted angular momentum from a mass-transferring, O-type or Wolf-Rayet companion before that star blew up as a supernova. To test this prediction, we have measured the spin rates of eight O stars in Wolf-Rayet (WR) + O binaries, increasing the total sample size of such O stars' measured spins from two to ten. Polarimetric and other determinations of these systems' sin i allow us to determine an average equatorial rotation velocity from HeI (HeII) lines of $v_e$ = 348 (173) km/s for these O stars, with individual star's $v_e$ from HeI (HeII) lines ranging from 482 (237) to 290 (91) km/s. We argue that the $\sim$ 100\% difference between HeI and HeII speeds is due to gravity darkening. Super-synchronous spins, now observed in all 10 O stars in WR + O binaries where it has been measured, are strong observational evidence that Roche lobe overflow mass transfer from a WR progenitor companion has played a critical role in the evolution of WR+OB binaries. While theory predicts that this mass transfer rapidly spins-up the O-type mass gainer to a nearly break-up rotational velocity $v_e \sim 530 $ km/s, the observed average $v_e$ of the O-type stars in our sample is 65\% that speed. This demonstrates that, even over the relatively short WR-phase timescale, tidal and/or other effects causing rotational spin-down must be efficient. A challenge to tidal synchronization theory is that the two longest-period binaries in our sample (with periods of 29.7 and 78.5 days) unexpectedly display super-synchronous rotation., Comment: 10 pages, 10 figures, submitted to MNRAS

Three decades of searches have revealed 154 Wolf-Rayet (WR) stars in M31, with 62 of WC type, 92 of WN type and zero of transition type WN/C or WC/N. In apparent contrast, about two percent of the WR stars in the Galaxy, the LMC and M33 simultaneously display strong lines of carbon and nitrogen, i.e. they are transition type WN/C or WC/N stars. We report here the serendipitous discovery of M31 WR 84-1, the first transition star in M31, located at RA = 00:43:43.61 DEC = +41:45:27.95 (J2000). We present its spectrum, classify it as WN5/WC6, and compare it with other known transition stars. The star is unresolved in Hubble Space Telescope narrowband and broadband images, while its spectrum displays strong, narrow emission lines of hydrogen, [NII], [SII] and [OIII]; this indicates a nebula surrounding the star. The radial velocity of the nebular lines is consistent with that of gas at the same position in the disc of M31. The metallicity at the 11.8 kpc galactocentric distance of M31 84-1 is approximately solar, consistent with other known transition stars. We suggest that modest numbers of reddened WR stars remain to be found in M31., Comment: 5 pages, 3 figures, submitted to Monthly Notices of the Royal Astronomical Society

Seitenzahl et al. (2009) have predicted that roughly three years after its explosion, the light we receive from a Type Ia supernova (SN Ia) will come mostly from reprocessing of electrons and X-rays emitted by the radioactive decay chain $^{57}{\rm Co}~\to~^{57}{\rm Fe}$, instead of positrons from the decay chain $^{56}{\rm Co}~\to~^{56}{\rm Fe}$ that dominates the SN light at earlier times. Using the {\it Hubble Space Telescope}, we followed the light curve of the SN Ia SN 2012cg out to $1055$ days after maximum light. Our measurements are consistent with the light curves predicted by the contribution of energy from the reprocessing of electrons and X-rays emitted by the decay of $^{57}$Co, offering evidence that $^{57}$Co is produced in SN Ia explosions. However, the data are also consistent with a light echo $\sim14$ mag fainter than SN 2012cg at peak. Assuming no light-echo contamination, the mass ratio of $^{57}$Ni and $^{56}$Ni produced by the explosion, a strong constraint on any SN Ia explosion model, is $0.043^{+0.012}_{-0.011}$, roughly twice Solar. In the context of current explosion models, this value favors a progenitor white dwarf with a mass near the Chandrasekhar limit., Comment: Updated to reflect the final version published by ApJ. For a video about the paper, see https://youtu.be/t3pUbZe8wqk

T Pyxidis is the only recurrent nova surrounded by knots of material ejected in previous outbursts. Following the eruption that began on 2011 April 14.29, we obtained seven epochs (from 4 to 383 days after eruption) of Hubble Space Telescope narrowband Ha images of T Pyx . The flash of radiation from the nova event had no effect on the ejecta until at least 55 days after the eruption began. Photoionization of hydrogen located north and south of the central star was seen 132 days after the beginning of the eruption. That hydrogen recombined in the following 51 days, allowing us to determine a hydrogen atom density of at least 7e5 cm^-3 - at least an order of magnitude denser than the previously detected, unresolved [NII] knots surrounding T Pyx. Material to the northwest and southeast was photoionized between 132 and 183 days after the eruption began. 99 days later that hydrogen had recombined. Both then (282 days after outburst) and 101 days later, we detected almost no trace of hydrogen emission around T Pyx. There is a large reservoir of previously unseen, cold diffuse hydrogen overlapping the previously detected, [NII] - emitting knots of T Pyx ejecta. The mass of this newly detected hydrogen is probably an order of magnitude larger than that of the [NII] knots. We also determine that there is no significant reservoir of undetected ejecta from the outer boundaries of the previously detected ejecta out to about twice that distance, near the plane of the sky. The lack of distant ejecta is consistent with the Schaefer et al (2010) scenario for T Pyx, in which the star underwent its first eruption within five years of 1866 after many millennia of quiescence, followed by the six observed recurrent nova eruptions since 1890. This lack of distant ejecta is not consistent with scenarios in which T Pyx has been erupting continuously as a recurrent nova for many centuries or millennia., Comment: 27 pages, 10 figures, submitted to the Astrophysical Journal

A new method of image subtraction is applied to images from a J, K, and narrow-band imaging survey of 300 square degrees of the plane of the Galaxy, searching for new Wolf-Rayet stars. Our survey spans 150 degrees in Galactic longitude and reaches 1 degree above and below the Galactic plane. The survey has a useful limiting magnitude of K =15 over most of the observed Galactic plane, and K =14 (due to severe crowding) within a few degrees of the Galactic center. The new image subtraction method described here (better than aperture or even point-spread-function photometry in very crowded fields) detected several thousand emission-line candidates. In June and July 2011 and 2012, we spectroscopically followed up on 333 candidates with MDM-TIFKAM and IRTF-SPEX, discovering 89 emission-line sources. These include 49 Wolf-Rayet stars, 43 of them previously unidentified, including the most distant known Galactic WR stars, more than doubling the number on the far side of the Milky Way. We also demonstrate our survey's ability to detect very faint PNe and other NIR emission objects., Comment: 22 pages, 19 figures; in press at MNRAS

We present refined color-color selection criteria for identifying Wolf-Rayet (WR) stars using available mid infrared (MIR) photometry from WISE in combination with near infrared (NIR) photometry from 2MASS. Using a sample of spectrally classified objects, we find that WR stars are well distinguished from the field stellar population in the (W1-W2) vs. (J-Ks) color-color diagram, and further distinguished from other emission line objects such as Planetary Nebulae, Be, and Cataclysmic variable stars using a combination of NIR and MIR color constraints. As proof of concept we applied the color constraints to a photometric sample in the Galactic plane, located WR star candidates, and present five new spectrally confirmed and classified WC (1) and WN (4) stars. Analysis of the 0.8-5.0 micron spectral data for a subset of known, bright WC and WN stars shows that emission lines (primarily He I) extend into the 3.0 - 5.0 micron spectral region, although their strength is greatly diminished compared to the 0.8-2.5 micron region. The WR population stands out relative to background field stars at NIR and MIR colors due to an excess continuum contribution, likely caused by free-free scattering in dense winds. Mean photometric properties of known WRs are presented and imply that reddened late-type WN and WC sources are easier to detect than earlier-type sources at larger Galactic radii. WISE W3 and W4 images of 10 Wolf Rayet stars show evidence of circumstellar shells linked to mass ejections from strong stellar winds., Comment: Accepted for publication in AJ. 29 pages, 18 figures, 8 tables

Hubble Space Telescope} images of the ejecta surrounding the nova T Pyxidis resolve the emission into more than two thousand bright knots. We simulate the dynamical evolution of the ejecta from T Pyxidis during its multiple eruptions over the last 150 years using the adaptive mesh refinement capability of the gas dynamics code Ramses. We demonstrate that the observed knots are the result of Richtmeyer-Meshkov gas dynamical instabilities (the equivalent of Rayleigh-Taylor instabilities in an accelerated medium). These instabilities are caused by the overrunning of the ejecta from the classical nova of 1866 by fast moving ejecta from the subsequent six recurrent nova outbursts. The model correctly predicts the observed expansion and dimming of the T Pyx ejecta as well as the knotty morphology. The model also predicts that deeper, high resolution imagery will show filamentary structure connecting the knots. We show reprocessed Hubble Space Telescope imagery that shows the first hints of such structure., Comment: 11 pages, 11 figures Changed typo in Author's name in abstract

M101 is an ideal target in which to test predictions of massive star birth and evolution. The large abundance gradient across M101 (a factor of 20) suggests that many more WR stars must be found in the inner parts of this galaxy than in the outer regions. Many HII regions and massive star-forming complexes have been identified in M101; they should be rich in WR stars, and surrounded by RSG stars. Finally, the Wolf-Rayet stars in M101 may be abundant enough for one to explode as a Type Ib or Ic supernova and/or GRB within a generation. The clear identification of the progenitor of a Type Ib or Ic supernova as a WR star would be a major confirmation of current stellar evolution theory. Motivated by these considerations, we have used the Hubble Space Telescope to carry out a deep, HeII optical narrowband imaging survey of the massive star populations in the ScI spiral galaxy M101. Combined with archival broadband images, we were able to image almost the entire galaxy with the unprecedented depth and resolution that only HST affords. We describe the extent of the survey and our images, as well as our data reduction procedures. A detailed study of a field east of the center of M101, containing the giant star-forming region NGC 5462, demonstrates how we find candidates, their properties and spatial distribution, and how we rule out most contaminants. The spatial distributions of the WR and RSG stars near a giant star-forming complex are strikingly different. WR stars dominate the complex core, while RSG dominate the complex halo. Future papers in this series will describe and catalog more than a thousand WR and RSG candidates that are detectable in our images, as well as spectra of many of those candidates., Comment: 33 pages, 10 Figures, submitted to the Astronomical Journal

We are systematically surveying all known and suspected Z Cam-type dwarf novae for classical nova shells. This survey is motivated by the discovery of the largest known classical nova shell, which surrounds the archetypal dwarf nova Z Camelopardalis. The Z Cam shell demonstrates that at least some dwarf novae must have undergone classical nova eruptions in the past, and that at least some classical novae become dwarf novae long after their nova thermonuclear outbursts, in accord with the hibernation scenario of cataclysmic binaries. Here we report the detection of a fragmented "shell", 3 arcmin in diameter, surrounding the dwarf nova AT Cancri. This second discovery demonstrates that nova shells surrounding Z Cam-type dwarf novae cannot be very rare. The shell geometry is suggestive of bipolar, conical ejection seen nearly pole-on. A spectrum of the brightest AT Cnc shell knot is similar to that of the ejecta of the classical nova GK Per, and of Z Cam, dominated by [NII] emission. Galex FUV imagery reveals a similar-sized, FUV-emitting shell. We determine a distance of 460 pc to AT Cnc, and an upper limit to its ejecta mass of ~ 5 x 10^{-5} Msun, typical of classical novae., Comment: 17 pages, 6 figures, submitted to AJ. arXiv admin note: text overlap with arXiv:1205.3531

The dwarf nova Z Camelopardalis is surrounded by the largest known classical nova shell. This shell demonstrates that at least some dwarf novae must have undergone classical nova eruptions in the past, and that at least some classical novae become dwarf novae long after their nova thermonuclear outbursts. The current size of the shell, its known distance, and the largest observed nova ejection velocity set a lower limit to the time since Z Cam's last outburst of 220 years. The brightest part of the Z Cam shell's radius is currently p ~ 1690 pixels. No expansion of the radius of the brightest part of the ejecta was detected, with an upper limit of pdot < 1 pixel/3 years. This suggests that the last Z Cam eruption occurred p/pdot > 5,000 years ago. However, including the important effect of deceleration as the ejecta sweeps up interstellar matter in its snowplow phase reduces the lower limit to 1300 years. This is the first strong test of the prediction of nova thermonuclear runaway theory that the inter-outburst times of classical novae are longer than 1000 yr. The intriguing suggestion that Z Cam was a bright nova, recorded by Chinese imperial astrologers in October - November 77 BCE, is consistent with our measurements. If Z Cam was indeed the nova of 77 BCE we predict that its ejecta are currently expanding at 85 km/s, or 0.11 arcsec/yr. Detection and measurement of this rate of expansion are doable in just a few years., Comment: 6 figures, submitted to the Astrophysical Journal

We have imaged the ejecta of GK Persei (Nova Persei 1901 A.D.) with the Hubble Space Telescope (HST), revealing hundreds of cometary-like structures. One or both ends of the structures often show a brightness enhancement relative to the structures' middle sections, but there is no simple regularity to their morphologies (in contrast with the Helix nebula). Some of the structures' morphologies suggest the presence of slow-moving or stationary material with which the ejecta is colliding, while others suggest shaping from a wind emanating from GK Per itself. A detailed expansion map of the nova's ejecta was created by comparing HST images taken in successive years. WFPC2 narrowband images and STIS spectra demonstrate that the physical conditions in the ejecta vary strongly on spatial scales much smaller than those of the ejecta. Directly measuring accurate densities and compositions, and hence masses of this and other nova shells, will demand data at least as resolved spatially as those presented here. The filling factor the ejecta is < 1%, and the nova ejecta mass must be less than $10^{-4} \Msun$. A few of the nebulosities vary in brightness by up to a factor of two on timescales of one year. Finally, we present the deepest images yet obtained of a jet-like feature outside the main body of GK Per nebulosity, and the first spectrum of that feature. Dominated by strong, narrow emission lines of [NII], [OII], [OIII], and [SII], this feature is probably a shock due to ejected material running into stationary ISM, slowly moving ejecta from a previous nova episode, or circum-binary matter present before 1901. An upper limit to the mass of the jet is of order a few times $10^{-6} \Msun$. The jet might be an important, or even dominant mass sink from the binary system. The jet's faintness suggests that similar features could easily have been missed in other cataclysmic binaries., Comment: 43 pages, 17 figures

We explore the possible contribution of classical and recurrent novae from isolated white dwarfs accreting from the intracluster medium to the abundances of "second generation" globular cluster stellar populations. We show that under reasonable assumptions the helium abundances of clusters can be enhanced substantially by these novae and argue that novae should be considered as an important, and perhaps even dominant channel in the evolution of the intracluster medium. We also discuss a possible test for whether helium enhancement really is the cause of the multiple main sequences in globular clusters that is independent of the positions of stars in the color-magnitude diagram., Comment: 6 pages, accepted to MNRAS

We are continuing a J, K and narrow-band imaging survey of 300 square degrees of the plane of the Galaxy, searching for new Wolf-Rayet stars. Our survey spans 150 degrees in Galactic longitude and reaches 1 degree above and below the Galactic plane. The survey has a useful limiting magnitude of K = 15 over most of the observed Galactic plane, and K = 14 (due to severe crowding) within a few degrees of the Galactic center. Thousands of emission line candidates have been detected. In spectrographic follow-ups of 146 relatively bright WR star candidates we have re-examined 11 previously known WC and WN stars and discovered 71 new WR stars, 17 of type WN and 54 of type WC. Our latest image analysis pipeline now picks out WR stars with a 57% success rate. Star subtype assignments have been confirmed with K band spectra, and distances approximated using the method of spectroscopic parallax. Some of the new WR stars are amongst the most distant known in our Galaxy. The distribution of these new WR stars is beginning to trace the locations of massive stars along the distant spiral arms of the Milky Way., Comment: 48 pages, 25 figures, 8 Tables, Accepted by AJ 4-18-2012

Extremely luminous, red eruptive variables like RV in M31 are being suggested as exemplars of a new class of astrophysical object. Our greatly extended series of nova simulations shows that classical nova models can produce very red, luminous eruptions. In a poorly studied corner of 3-D nova parameter space (very cold, low-mass white dwarfs, accreting at very low rates) we find bona fide classical novae that are very luminous and red because they eject very slowly moving, massive envelopes. A crucial prediction of these nova models - in contrast to the predictions of merging star ("mergeburst") models - is that a hot remnant, the underlying white dwarf, will emerge after the massive ejected envelope has expanded enough to become optically thin. This blue remnant must fade on a timescale of decades - much faster than a "mergeburst", which must fade on timescales of millennia or longer. Furthermore, the cooling nova white dwarf and its expanding ejecta must become redder in the years after eruption, while a contracting mergeburst must become hotter and bluer. We predict that red novae will always brighten to L~1000 Lsun for about 1 year before rising to maximum luminosity at L~10^6 - 10^7 Lsun. The maximum luminosity attainable by a nova is likely to be L~10^7 Lsun, corresponding to M-12. In an accompanying paper we describe a fading, luminous blue candidate for the remnant of M31-RV; it is observed with HST to be compatible only with the nova model., Comment: Accepted for publication in ApJ. Comments welcome

M31-RV was an extraordinarily luminous (~10^6 Lsun) eruptive variable, displaying very cool temperatures (roughly 1000 Kelvins) as it faded. The photometric behavior of M31-RV (and several other very red novae, i.e. luminous eruptive red variables) has led to several models of this apparently new class of astrophysical object. One of the most detailed models is that of "mergebursts": hypothetical mergers of close binary stars. These are predicted to rival or exceed the brightest classical novae in luminosity, but to be much cooler and redder than classical novae, and to become slowly hotter and bluer as they age. This prediction suggests two stringent and definitive tests of the mergeburst hypothesis. First, there should always be a cool red remnant, and NOT a hot blue remnant at the site of such an outburst. Second, the inflated envelope of a mergeburst event should be slowly contracting, hence it must display a slowly rising effective temperature. We have located a luminous, UV-bright object within 0.4 arcsec (1.5 sigma of the astrometric position) of M31-RV in archival WFPC2 images taken 10 years after the outburst: it resembles an old nova. Twenty years after the outburst, the object remains much too hot to be a mergeburst. Its behavior remains consistent with that of theoretical nova models which erupt on a low mass white dwarf. Future Hubble UV and visible images could determine if the M31-RV analogs (in M85 and in M99) are also behaving like old novae., Comment: Accepted for publication in ApJ, comments welcome

Using archival RXTE data, we show that the ultracompact X-ray binary in NGC 1851 exhibits large amplitude X-ray flux varations of more than a factor of 10 on timescales of days to weeks and undergoes sustained periods of months where the time-averaged luminosty varies by factors of two. Variations of this magnitude and timescale have not been reported previously in other ultracompact X-ray binaries. Mass transfer in ultracompact binaries is thought to be driven by gravitational radiation and the predicted transfer rates are so high that the disks of ultracompact binaries with orbits as short as that of this object should not be susceptible to ionization instabilities. Therefore the variability characteristics we observe were unexpected, and need to be understood. We briefly discuss a few alternatives for producing the observed variations in light of the fact that the viscous timescale of the disk is of order a week, comparable to the shorter time scale variation that is observed but much less than the longer term variation. We also discuss the implications for interpretation of observations of extragalactic binaries if the type of variability seen in the source in NGC 1851 is typical., Comment: 7 pages, 5 figures, accepted for publication in MNRAS

We have searched for variable sources in the core region of M80, using far ultra-violet data taken with the Advanced Camera for Surveys on board the Hubble Space Telescope. We found three sources that exhibit strong signs of variability in our data. Among these is source TDK1, which we believe to be an RR Lyrae star that reached maximum brightness during our observations. The light curve shows a >3 mag FUV brightening over the course of ~5 hours, with an estimated peak brightness of ~16.7 mag, followed by a decrease to ~20 mag. Archival optical data obtained with WFPC2 confirm that TDK1 is variable in all wavebands. TDK1's SED is reasonably fit by a star with temperature T(eff)=6700K and radius R=4.2R(sun), consistent with the suggestion that it is an RR Lyrae. Based on the photometric and variability characteristics of the other two variables, we suggest that TDK2 is likely to be an SX Phoenicis star with ~55 minutes period, and TDK3 is likely another RR Lyrae. Finally, we briefly discuss the FUV counterparts to two previously known variables in M80, the classical nova T Sco and the dwarf nova DN1., Comment: 12 pages, 9 figures and 3 tables. Accepted for publication in MNRAS.

We present an updated catalogue of M31 globular clusters (GCs) based on images from the Wide Field CAMera (WFCAM) on the UK Infrared Telescope and from the Sloan Digital Sky Survey (SDSS). Our catalogue includes new, self-consistent ugriz and K-band photometry of these clusters. We discuss the difficulty of obtaining accurate photometry of clusters projected against M31 due to small scale background structure in the galaxy. We consider the effect of this on the accuracy of our photometry and provide realistic photometric error estimates. We investigate possible contamination in the current M31 GC catalogues using the excellent spatial resolution of these WFCAM images combined with the SDSS multicolour photometry. We identify a large population of clusters with very blue colours. Most of these have recently been proposed by other work as young clusters. We distinguish between these, and old clusters, in the final classifications. Our final catalogue includes 416 old clusters, 156 young clusters and 373 candidate clusters. We also investigate the structure of M31's old GCs using previously published King model fits to these WFCAM images. We demonstrate that the structure and colours of M31's old GC system are similar to those of the Milky Way. One GC (B383) is found to be significantly brighter in previous observations than observed here. We investigate all of the previous photometry of this GC and suggest that this variability appears to be genuine and short lived. We propose that the large increase in its luminosity my have been due to a classical nova in the GC at the time of the previous observations in 1989., Comment: Accepted for publication in MNRAS: 17 pages, 11 figures and 3 tables. Table 1 is available in full from http://www.astro.soton.ac.uk/~m.b.peacock/m31gc.html

The discovery of new Wolf-Rayet (WR) stars in our Galaxy via large-scale narrowband optical surveys has been severely limited by dust extinction. Recent improvements in infrared technology have made narrowband-broadband imaging surveys viable again. We report a new J, K and narrow-band imaging survey of 300 square degrees of the plane of the Galaxy, spanning 150 degrees in Galactic longitude and reaching 1 degree above and below the Galactic plane. The survey has a useful limiting magnitude of K = 15 over most of the observed Galactic plane, and K = 14 within a few degrees of the Galactic center. Thousands of emission line candidates have been detected. In spectrographic follow-ups of 173 WR star candidates we have discovered 41 new WR stars, 15 of type WN and 26 of type WC. Star subtype assignments have been confirmed with K band spectra, and distances approximated using the method of spectroscopic parallax. A few of the new WR stars are amongst the most distant known in our Galaxy. The distribution of these new WR stars is seen to follow that of previously known WR stars along the spiral arms of the Galaxy. Tentative radial velocities were also measured for most of the new WR stars., Comment: 55 pages, 23 figures, 7 tables, accepted to Astronomical Journal

We investigate the relationship between Low Mass X-ray Binaries (LMXBs) and globular clusters (GCs) using UKIRT observations of M31 and existing Chandra, XMM-Newton, and ROSAT catalogues. By fitting King models to these data we have estimated the structural parameters and stellar collision rates of 239 of its GCs. We show a highly significant trend between the presence of a LMXB and the stellar collision rate of a cluster. The stellar collision rate is found to be a stronger predictor of which clusters will host LMXBs than the host cluster mass. We argue that our results show that the stellar collision rate of the clusters is the fundamental parameter related to the production LMXBs. This is consistent with the formation of LMXBs through dynamical interactions with little direct dependence on the neutron star retention fraction or cluster mass., Comment: 5 pages, 4 figures, 1 table (available in full from http://www.astro.soton.ac.uk/~m.b.peacock/m31gc.html), accepted for publication in MNRAS Letters