AT2016dah and AT2017fyp are fairly typical Andromeda Galaxy (M31) classical novae. AT2016dah is an almost text book example of a 'very fast' declining, yet uncommon, Fe II'b' (broad-lined) nova, discovered during the rise to peak optical luminosity, and decaying with a smooth broken power-law light curve. AT2017fyp is classed as a 'fast' nova, unusually for M31, its early decline spectrum simultaneously shows properties of both Fe II and He/N spectral types - a 'hybrid'. Similarly, the light curve of AT2017fyp has a broken power-law decline but exhibits an extended flat-topped maximum. Both novae were followed in the UV and X-ray by the Neil Gehrels Swift Observatory, but no X-ray source was detected for either nova. The pair were followed photometrically and spectroscopically into their nebular phases. The progenitor systems were not visible in archival optical data, implying that the mass donors are main sequence stars. What makes AT2016dah and AT2017fyp particularly interesting is their position with respect to M31. The pair are close on the sky but are located far from the centre of M31, lying almost along the semi-minor axis of their host. Radial velocity measurements and simulations of the M31 nova population leads to the conclusion that both novae are members of the Andromeda Giant Stellar Stream (GSS). We find the probability of at least two M31 novae appearing coincident with the GSS by chance is ~1%. Therefore, we claim that these novae arose from the GSS progenitor, not M31 - the first confirmed novae discovered in a tidal steam.

The hydrogen-rich envelopes accreted by white dwarf stars from their red dwarf companions lead to thermonuclear runaways observed as classical nova eruptions peaking at up to 1 Million solar luminosities. Virtually all nova progenitors are novalike binaries exhibiting high rates of mass transfer to their white dwarfs before and after an eruption. It is a puzzle that binaries indistinguishable from novalikes, but with much lower mass transfer rates, and resulting dwarf nova outbursts, co-exist at the same orbital periods. Nova shells surrounding several dwarf novae demonstrate that at least some novae become dwarf novae between successive nova eruptions, though the mechanisms and timescales governing mass transfer rate variations are poorly understood. Here we report simulations of the multiGyr evolution of novae which self-consistently model every eruption's thermonuclear runaway, mass and angular momentum losses, feedback due to irradiation and variable mass transfer, and orbital size and period changes. The simulations reproduce the observed wide range of mass transfer rates at a given orbital period, with large and cyclic changes in white dwarf-red dwarf binaries emerging on kyr to Myr timescales. They also demonstrate that deep hibernation, (complete stoppage of mass transfer for long periods), occurs only in short-period binaries; that initially very different binaries converge to become nearly identical systems; that while almost all prenovae should be novalike binaries, dwarf novae should also occasionally be observed to give rise to novae; and that the masses of white dwarfs decrease only slightly while their red dwarf companions are consumed.

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.

Large numbers of low-to-medium resolution spectra of M-type dwarf stars from both the local Galactic disk and halo are available from various surveys. In order to fully exploit these data, we develop a template-fit method using a set of empirically assembled M dwarf/subdwarf classification templates, based on the measurements of the TiO and CaH molecular bands, to classify M dwarfs/subdwarfs by spectral type and metallicity class. We further present a pipeline to automatically determine the effective temperature, metallicity ([M/H]), alpha-element to iron abundance ratio ([alpha/Fe]), and surface gravity of M dwarfs/subdwarfs using the latest version of BT-Settl model atmospheres. We apply these methods to a set of low-to-medium resolution optical spectra of 1,544 high proper-motion (> 0.4"/yr) M dwarfs/subdwarfs. Our metallicity estimates appear to be consistent with the expected color-magnitude variation of stars relative to atmospheric composition, as our sample shows a clear stratification with respect to metallicity in the Gaia H-R diagram. Furthermore, the measured chemical parameters of the two components in 48 binary systems are in a good agreement with each other, which suggest a precision of 0.22 dex in [M/H], 0.08 dex in [alpha/Fe], and 0.16 dex in the combined index [alpha/Fe]+[M/H]. There is also a good consistency between metallicity class, obtained from the empirical classification templates, and the index [alpha/Fe]+[M/H] (~ [alpha/H]), obtained from model-fitting, which means that the more easily measured metallicity class can be used as a relatively reliable indicator of absolute alpha-element abundance, [alpha/H], in low-mass stars. Finally, we examine the distribution of our stars in the [alpha/Fe] vs. [M/H] diagram, which shows evidence of clustering in chemical abundance makeup, suggestive of discrete populations among the local disk and halo stars.

We investigate the old open cluster M67 using ultraviolet photometric data of Ultra-Violet Imaging Telescope in multi-filter far-UV bands. M67, well known for the presence of several blue straggler stars (BSS), has been put to detailed tests to understand their formation pathways. Currently, there are three accepted formation channels: mass transfer due to Roche-lobe overflow in binary systems, stellar mergers either due to dynamical collisions or through coalescence of close binaries. So far, there had not been any confirmed detection of a white dwarf (WD) companion to any of the BSSs in this cluster. Here, we present the detection of WD companions to 5 bright BSSs in M67. The multiwavelength spectral energy distributions covering 0.12 -11.5 $\mu$m range, were found to require binary spectral fits for 5 BSSs, consisting of a cool (BSS) and a hot companion. The parameters (Luminosity, Temperature, Radius and Mass) of the hot companions suggest them to be WDs with mass in the range 0.2 - 0.35 M$_{\odot}$ with T$_{eff}$ $\sim$ 11000 - 24000 K.

The old open cluster M67, populated with blue straggler stars (BSSs), is a well known test bed to study the BSS formation pathways. Here, we report the first direct detection of a white dwarf (WD) companion to a BSS in M67, using far-UV images from the \textit{Ultra Violet Imaging telescope} ({\it UVIT}) on {\it ASTROSAT}. Near-simultaneous observations in three far-UV bands combined with {\it GALEX}, {\it IUE}, ground and space based photometric data covering 0.14 -11.5 $\mu$m range for WOCS1007 were found to require a binary fit to its spectral energy distribution (SED), consisting of a BSS and a hot companion. On the other hand, a single spectral fit was found to be satisfactory for the SEDs of two other BSSs, WOCS1006 and WOCS2011, with the latter showing a deficient far-UV flux. The hot companion of WOCS1007 is found to have a T$_{eff}$ $\sim$ 13250-13750K and a radius of 0.09$\pm$0.01 R$_\odot$. A comparison with WD models suggests it to be a low mass WD ($\sim$ 0.18$M_\odot$), in agreement with the kinematic mass from the literature. As a low mass WD ($<$ 0.4$M_\odot$) necessitates formation through mass transfer (MT) in close binaries, WOCS1007 with a known period of 4.2 days along with its fast rotation, is likely to be formed by a case A or case B binary evolution.

The supersoft X-ray and optical transient ASASSN-16oh has been interpreted by Maccarone et al. (2019) as having being induced by an accretion event on a massive white dwarf, resembling a dwarf nova super-outburst. These authors argued that the supersoft X-ray spectrum had a different origin than in an atmosphere heated by shell nuclear burning, because no mass was ejected. We find instead that the event's timescale and other characteristics are typical of non-mass ejecting thermonuclear runaways, as already predicted by Shara et al. (1977) and the extensive grid of nova models by Yaron et al. (2005). We suggest that the low X-ray and bolometric luminosity in comparison to the predictions of the models of nuclear burning are due to an optically thick accretion disk, hiding most of the white dwarf surface. If this is the case, we calculated that the optical transient can be explained as a non-ejective thermonuclear event on a WD of $\simeq$1.1M$_\odot$ accreting at the rate of $\simeq3.5{-}5{\times}10^{-7}$M$_\odot$yr$^{-1}$. We make predictions that should prove whether the nature of the transient event was due to thermonuclear burning or to accretion; observational proof should be obtained in the next few years, because a new outburst should occur within $\simeq$10-15 years of the event.

We report the discovery of a hot white dwarf (WD) companion to a blue straggler star (BSS) in the globular cluster (GC) NGC 5466, based on observations from the Ultra-Violet Imaging Telescope (UVIT) on board AstroSat. The Spectral Energy Distribution (SED) of the Far-UV detected BSS NH 84 was constructed by combining the flux measurements from 4 filters of UVIT, with GALEX, GAIA and other ground-based observations. The SED of NH 84 reveals the presence of a hot companion to the BSS. The temperature and radius of the BSS (T$_{\mathrm{eff}} = 8000^{+1000}_{-250}$ K, R/R$_\odot = 1.44 \pm 0.05$) derived from Gemini spectra and SED fitting using Kurucz atmospheric models are consistent with each other. The temperature and radius of the hotter companion of NH 84 (T$_{\mathrm{eff}} = 32,000 \pm 2000$ K, R/R$_\odot = 0.021 \pm 0.007$) derived by fitting Koester WD models to the SED suggest that it is likely to be a hot WD. The radial velocity derived from the spectra along with the proper motion from GAIA DR2 confirms NH 84 to be a kinematic member of the cluster. This is the second detection of a BSS-WD candidate in a GC, and the first in the outskirts of a low density GC. The location of this BSS in NGC 5466 along with its dynamical age supports the mass-transfer pathway for BSS formation in low density environments.

We have initiated a survey aimed at locating a nearly complete sample of classical symbiotic stars (SySt) in the Magellanic Clouds. Such a sample is nearly impossible to obtain in the Milky Way, and is essential to constrain the formation, evolution and demise of these strongly interacting, evolved binary stars. We have imaged both Clouds in Halpha and He II 4686 narrow-band filters deeply enough to detect all known symbiotic stars. While He II 4686 is not present in all SySt, our method should yield a high success rate because the mimics of SySt are not as likely as true symbiotics to show this emission line. We demonstrate the viability of our method through the discovery and characterization of three new SySt in the Small Magellanic Cloud: 2MASS J00411657-7233253, 2MASS J01104404-7208464 and 2MASS J01113745-7159023. Enigmatic variability was observed in 2MASS J01113745-7159023, where changes in the amplitude of its quasi-periodic variability may suggest an enhanced mass transfer rate during a periastron passage on an elliptical orbit. 2MASS J01104404-7208464 is an ellipsoidal variable with an orbital period of 403d.

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.

For a sample of 38 Galactic globular clusters (GCs), we confront the observed distributions of blue straggler (BS) proper motions and masses (derived from isochrone fitting) from the BS catalog of Simunovic & Puzia with theoretical predictions for each of the two main competing BS formation mechanisms. These are mass transfer from an evolved donor on to a main-sequence (MS) star in a close binary system, and direct collisions involving MS stars during binary encounters. We use the \texttt{FEWBODY} code to perform simulations of single-binary and binary-binary interactions. This provides collisional velocity and mass distributions for comparison to the observed distributions. Most clusters are consistent with BSs derived from a dynamically relaxed population, supportive of the binary mass-transfer scenario. In a few clusters, including all the post-core collapse clusters in our sample, the collisional velocities provide the best fit.

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.

The MeerKAT radio telescope array, the Large Synoptic Survey Telescope (LSST), and eventually the Square Kilometer Array (SKA) will usher in a remarkable new era in astronomy, with thousands of transients being discovered and transmitted to the astronomical community in near-real-time each night. Immediate spectroscopic follow-up will be critical to understanding their early-time physics - a task to which the Southern African Large Telescope (SALT) is uniquely suited, given its southerly latitude and the 14-degree-diameter uncorrected field (patrol area) of its 10-m spherical primary mirror. A new telescope configuration is envisioned, incorporating multiple mini-trackers that range around a much larger patrol area of 35 degrees in diameter. Each mini-tracker is equipped with a small spherical aberration corrector feeding an efficient, low resolution spectrograph to perform contemporaneous follow-up observations.

In this paper, the fourth in the series, we continue our study of combinatorics in chaotic Newtonian dynamics. We focus once again on the chaotic four-body problem in Newtonian gravity assuming finite-sized particles, and interactions that produce direct collisions between any two particles. Our long-term goal is to predict the probability of a given collision event occurring over the course of an interaction, as a function of the numbers and properties of the particles. In previous papers, we varied the number of interacting particles, as well as the distributions of particle radii and masses. Here, we refine the methods developed in these preceding studies, and arrive at a final and robust methodology that can be used to study collisional dynamics in a variety of astrophysical contexts, ranging from stars in star clusters, galaxies in galaxy groups and clusters and even the collisional growth of planetesimals in protoplanetary disks. We further present and refine the concept of a Collision Rate Diagram (CRD), the primary tool we use to quantify the relative rates for different collision scenarios to occur. The agreement between our final theoretical model and the results of numerical scattering simulations is excellent.

Models have long predicted that the frequency-averaged masses of white dwarfs in Galactic classical novae are twice as large as those of field white dwarfs. Only a handful of dynamically well-determined nova white dwarf masses have been published, leaving the theoretical predictions poorly tested. The recurrence time distributions and mass accretion rate distributions of novae are even more poorly known. To address these deficiencies, we have combined our extensive simulations of nova eruptions with the Strope et al (2010) and Schaefer et al (2010) databases of outburst characteristics of Galactic classical and recurrent novae to determine the masses of 92 white dwarfs in novae. We find that the mean mass (frequency averaged mean mass) of 82 Galactic classical novae is 1.06 (1.13) Msun, while the mean mass of 10 recurrent novae is 1.31 Msun. These masses, and the observed nova outburst amplitude and decline time distributions allow us to determine the long-term mass accretion rate distribution of classical novae. Remarkably, that value is just 1.3 x 10^{-10} Msun/yr, which is an order of magnitude smaller than that of cataclysmic binaries in the decades before and after classical nova eruptions. This predicts that old novae become low mass transfer rate systems, and hence dwarf novae, for most of the time between nova eruptions. We determine the mass accretion rates of each of the 10 known Galactic RN, finding them to be in the range 10^{-7} - 10^{-8} $ Msun/yr. We are able to predict the recurrence time distribution of novae and compare it with the predictions of population synthesis models.

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.

Despite 50 years of extensive binary research we have to conclude that the Roche lobe overflow/mass transfer process that governs close binary evolution is still poorly understood. It is the scope of the present paper to lift a tip of the veil by studying the spin-up and spin-down processes of the O-type components of WR+O binaries. We critically analyze the available observational data of rotation speeds of the O-type components in WR+O binaries. By combining a binary evolutionary code and a formalism that describes the effects of tides in massive stars with an envelope in radiative equilibrium, we compute the corresponding rotational velocities during the Roche lobe overflow of the progenitor binaries. In all the studied WR+O binaries, we find that the O-type stars were affected by accretion of matter during the RLOF of the progenitor. This means that common envelope evolution which excludes any accretion onto the secondary O-star, has not played an important role to explain the WR+O binaries. Moreover, although it is very likely that the O-type star progenitors were spun-up by the mass transfer, many ended the RLOF/mass transfer phase with a rotational velocity that is significantly smaller than the critical rotation speed. This may indicate that during the mass transfer phase there is a spin-down process which is of the same order as, although significantly less than that of the spin-up process. We propose a Spruit-Tayler type dynamo spin-down suggested in the past to explain the rotation speeds of the mass gainers in long-period Algols.

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.

Here we report that the most rapidly recurring nova, M31N 2008-12a, which erupts annually, is surrounded by a "nova super-remnant" which demonstrates that M31N 2008-12a has erupted with high frequency for millions of years.

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.

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.

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

We present the spectroscopic evolution of AT 2017gfo, the optical counterpart of the first binary neutron star (BNS) merger detected by LIGO and Virgo, GW170817. While models have long predicted that a BNS merger could produce a kilonova (KN), we have not been able to definitively test these models until now. From one day to four days after the merger, we took five spectra of AT 2017gfo before it faded away, which was possible because it was at a distance of only 39.5 Mpc in the galaxy NGC 4993. The spectra evolve from blue ($\sim6400$K) to red ($\sim3500$K) over the three days we observed. The spectra are relatively featureless --- some weak features exist in our latest spectrum, but they are likely due to the host galaxy. However, a simple blackbody is not sufficient to explain our data: another source of luminosity or opacity is necessary. Predictions from simulations of KNe qualitatively match the observed spectroscopic evolution after two days past the merger, but underpredict the blue flux in our earliest spectrum. From our best-fit models, we infer that AT 2017gfo had an ejecta mass of $0.03M_\odot$, high ejecta velocities of $0.3c$, and a low mass fraction $\sim10^{-4}$ of high-opacity lanthanides and actinides. One possible explanation for the early excess of blue flux is that the outer ejecta is lanthanide-poor, while the inner ejecta has a higher abundance of high-opacity material. With the discovery and follow-up of this unique transient, combining gravitational-wave and electromagnetic astronomy, we have arrived in the multi-messenger era.

We report on SALT low resolution optical spectroscopy and optical/IR photometry undertaken with other SAAO telescopes (MASTER-SAAO and IRSF) of the kilonova AT 2017gfo (aka SSS17a) in the galaxy NGC4993 during the first 10 days of discovery. This event has been identified as the first ever electromagnetic counterpart of a gravitational wave event, namely GW170817, which was detected by the LIGO and Virgo gravitational wave observatories. The event is likely due to a merger of two neutron stars, resulting in a kilonova explosion. SALT was the third telescope to obtain spectroscopy of AT 2017gfo and the first spectrum, 1.2 d after the merger, is quite blue and shows some broad features, but no identifiable spectral lines and becomes redder over time. We compare the spectral and photometric evolution with recent kilonova simulations and conclude that they are in qualitative agreement for post-merger wind models with proton: nucleon ratios of $Y_e$ = 0.25$-$0.30. The blue colour of the first spectrum is consistent with the lower opacity of the Lathanide-free r-process elements in the ejecta. Differences between the models and observations are likely due to the choice of system parameters combined with the absence of atomic data for more elements in the ejecta models.

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.

Horizontal branch (HB) morphology is a complex multiple-parameter problem. Besides the metallicity, two other leading parameters are the mass loss rate (MLR) and the initial He abundance of the HB progenitors. Using the STAREV stellar evolution code, we produce a wide array of Extreme Horizontal Branch (EHB) stars and also examine their post-HB evolution. EHB stars are produced in our calculations by the so called `delayed (late) hot core flash' scenario. The MLR is increased on the red giant branch (RGB) to the extent that, prior to reaching core flash conditions, only a very thin H-rich envelope remains and helium ignition takes place at hotter positions on the HRD. We perform an extensive, self-consistent parameter study, covering populations I and II (Z=0.0001-0.03), for both normal initial helium abundances and He-enriched models (up to Y=0.40). For completeness of the study and in order to point to complete trends, we chose NOT to cut out several combinations (or results of) that may extend beyond realistic limits. We present results and complete evolutionary tracks for the covered parameter space, showing in particular that: a) Increased He abundance ON ITS OWN -- without having a significant-enough MLR on the RGB -- DOES NOT lead to the production of EHB stars; however, b) The bluest (hottest) HB positions do result from the COMBINED effect of He-enhancement and increased MLR; c) The general trend is that the effective temperature on the HB increases with decreasing metallicity, but there is an indication for a halt, or even a reversal of this trend, as Z further decreases below $10^{-3}$; d) EHB stars can serve as major contributors to the UV flux emanating from their host system. Thus, the present comprehensive study both complements and lends support to previous, more restricted studies of the HB phase, and adds results for unexplored regions of the parameter space.

Deep narrow-band HST imaging of the iconic spiral galaxy M101 has revealed over a thousand new Wolf Rayet (WR) candidates. We report spectrographic confirmation of 10 HeII emission line sources hosting 15 WR stars. We find WR stars present at both sub- and super-solar metalicities with WC stars favouring more metal-rich regions compared to WN stars. We investigate the association of WR stars with HII regions using archival HST imaging and conclude that the majority of WR stars are in or associated with HII regions. Of the 10 emission lines sources, only one appears to be unassociated with a star-forming region. Our spectroscopic survey provides confidence that our narrow-band photometric candidates are in fact bonafide WR stars, which will allow us to characterise the progenitors of any core-collapse supernovae that erupt in the future in M101.

If accreting white dwarfs (WD) in binary systems are to produce type Ia supernovae (SNIa), they must grow to nearly the Chandrasekhar mass and ignite carbon burning. Proving conclusively that a WD has grown substantially since its birth is a challenging task. Slow accretion of hydrogen inevitably leads to the erosion, rather than the growth of WDs. Rapid hydrogen accretion does lead to growth of a helium layer, due to both decreased degeneracy and the inhibition of mixing of the accreted hydrogen with the underlying WD. However, until recently, simulations of helium-accreting WDs all claimed to show the explosive ejection of a helium envelope once it exceeded $\sim 10^{-1}\, \rm M_{\odot}$. Because CO WDs cannot be born with masses in excess of $\sim 1.1\, \rm M_{\odot}$, any such object, in excess of $\sim 1.2\, \rm M_{\odot}$, must have grown substantially. We demonstrate that the WD in the symbiotic nova RS Oph is in the mass range 1.2-1.4\,M$_{\odot}$. We compare UV spectra of RS Oph with those of novae with ONe WDs, and with novae erupting on CO WDs. The RS Oph WD is clearly made of CO, demonstrating that it has grown substantially since birth. It is a prime candidate to eventually produce an SNIa.

In this paper, the third in the series, we continue our study of combinatorics in chaotic Newtonian dynamics. We study the chaotic four-body problem in Newtonian gravity assuming finite-sized particles, and we focus on interactions that produce direct collisions between any two stars. Our long-term goal is to construct an equation that gives the probability of a given collision event occurring over the course of the interaction, as a function of the total encounter energy and angular momentum as well as the numbers and properties of the particles. In previous papers, we varied the number of interacting particles and the distribution of particle radii, for all equal mass particles. Here, we focus on the effects of different combinations of particle masses. We develop an analytic formalism for calculating the time-scales for different collision scenarios to occur. Our analytic time-scales reproduce the simulated time-scales when gravitational focusing is included. We present a method for calculating the relative rates for different types of collisions to occur, assuming two different limits for the particle orbits; radial and tangential. These limits yield relative collision probabilities that bracket the probabilities we obtain directly from numerical scattering experiments, and are designed to reveal important information about the (time-averaged) trajectories of the particles as a function of the interaction parameters. Finally, we present a Collision Rate Diagram (CRD), which directly compares the predictions of our analytic rates to the simulations and quantifies the quality of the agreement. The CRD will facilitate refining our analytic collision rates in future work, as we expand in to the remaining parameter space.

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.

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.

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.

The Z Cam-type dwarf nova AT Cnc displays a classical nova (CN) shell, demonstrating that mass transfer in cataclysmic binaries decreases substantially after a CN eruption. The hibernation scenario of cataclysmic binaries predicts such a decrease, on a timescale of a few centuries. In order to measure the time since AT Cnc's last CN eruption, we have measured the radial velocities of a hundred clumps in its ejecta with SITELLE, CFHT's recently commissioned imaging Fourier transform spectrometer. These range from -455 to +490 km/s. Coupled with the known distance to AT Cnc of 460 pc (Shara 2012), the size of AT Cnc's shell, and a simple model of nova ejecta deceleration, we determine that the last CN eruption of this system occurred $330_{-90}^{+135}$ years ago. This is the most rapid transition from a high mass transfer rate, novalike variable to a low mass transfer rate, dwarf nova yet measured, and in accord with the hibernation scenario of cataclysmic binaries. We conclude by noting the similarity in deduced outburst date (within a century of 1686 CE) of AT Cnc with a "guest star" reported in the constellation Cancer by Korean observers in 1645 CE.

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.

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.

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.

In this paper, we study the chaotic four-body problem in Newtonian gravity. Assuming point particles and total encounter energies $\le$ 0, the problem has three possible outcomes. We describe each outcome as a series of discrete transformations in energy space, using the diagrams first presented in Leigh \& Geller (2012; see the Appendix). Furthermore, we develop a formalism for calculating probabilities for these outcomes to occur, expressed using the density of escape configurations per unit energy, and based on the Monaghan description originally developed for the three-body problem. We compare this analytic formalism to results from a series of binary-binary encounters with identical point particles, simulated using the \texttt{FEWBODY} code. Each of our three encounter outcomes produces a unique velocity distribution for the escaping star(s). Thus, these distributions can potentially be used to constrain the origins of dynamically-formed populations, via a direct comparison between the predicted and observed velocity distributions. Finally, we show that, for encounters that form stable triples, the simulated single star escape velocity distributions are the same as for the three-body problem. This is also the case for the other two encounter outcomes, but only at low virial ratios. This suggests that single and binary stars processed via single-binary and binary-binary encounters in dense star clusters should have a unique velocity distribution relative to the underlying Maxwellian distribution (provided the relaxation time is sufficiently long), which can be calculated analytically.

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.

We consider the origins of enigmatic stellar populations in four Local Group galactic nuclei, specifically the Milky Way, M31, M32 and M33. These are centrally concentrated blue stars, found in three out of the four nuclear star clusters (NSCs) considered here. Their origins are unknown, but could include blue straggler (BS) stars, extended horizontal branch stars and young recently formed stars. Here, we calculate order-of-magnitude estimates for various collision rates, as a function of the host NSC environment and distance from the cluster centre. These rates are sufficiently high that BSs, formed via collisions between main sequence (MS) stars, could contribute non-negligibly ($\sim$ 1-10% in mass) to every surface brightness profile, with the exception of the Milky Way. Stellar evolution models show that the envelopes of red giant branch (RGB) stars must be nearly completely stripped to significantly affect their photometric appearance, which requires multiple collisions. Hence, the collision rates for individual RGB stars are only sufficiently high in the inner $\lesssim$ 0.1 pc of M31 and M32 for RGB destruction to occur. Collisions between white dwarfs and MS stars, which should ablate the stars, could offer a steady and significant supply of gas in every NSC in our sample. The gas could either fragment to form new stars, or accrete onto old MS stars already present. Thus, collisional processes could contribute significantly to the observed blue excesses in M31 and M33; future studies should be aimed at better constraining theoretical predictions to compliment existing and future observational data.

We present the detection of a variable star with an $18.05$ minute period in far-ultraviolet (FUV) images of the globular cluster NGC 1851 taken with the Hubble Space Telescope (HST). A candidate optical counterpart lies on the red horizontal branch or the asymptotic giant branch star of the cluster, but it is statistically possible that this is a chance superposition. This interpretation is supported by optical spectroscopt obtained with HST/STIS: the spectrum contains none of the strong emission lines that would be expected if the object was a symbiotic star (i.e. a compact accretor fed by a giant donor). We therefore consider two other possibilities for the nature of FUV variable: (i) an intermediate polar (i.e. a compact binary containing an accreting magnetic white dwarf), or (ii) an AM CVn star (i.e. an interacting double-degenerate system). In the intermediate polar scenario, the object is expected to be an X-ray source. However, no X-rays are detected at its location in $\simeq 65$~ksec of {\em Chandra} imaging, which limits the X-ray luminosity to $L_X \leqslant 10^{32}$~erg~s$^{-1}$. We therefore favour the AM CVn interpretation, but a FUV spectrum is needed to distinguish conclusively between the two possibilities. If the object is an AM CVn binary, it would be the first such system known in any globular cluster.

What is the difference between a long-lived unstable (or quasi-stable) multiple star system and a bona fide star cluster? In this paper, we present a possible framework to address this question, by studying the distributions of disruption times for chaotic gravitational encounters as a function of the number of interacting particles. To this end, we perform a series of numerical scattering experiments with the \texttt{FEWBODY} code, to calculate the distributions of disruption times as a function of both the particle number N and the virial coefficient k. The subsequent distributions are fit with a physically-motivated function, consisting of an initial exponential decay followed by a very slowly decreasing tail at long encounter times due to long-lived quasi-stable encounters. We find three primary features characteristic of the calculated distributions of disruption times. These are: (1) the system half-life increases with increasing particle number, (2) the fraction of long-lived quasi-stable encounters increases with increasing particle number and (3) both the system half-life and the fraction of quasi-stable encounters increase with decreasing virial coefficient. We discuss the significance of our results for collisional dynamics, and consider the extrapolation of our results to larger-N systems. We suggest that this could potentially offer a clear and unambiguous distinction between star clusters and (unstable or quasi-stable) multiple star systems. Although we are limited by very small-number statistics, our results tentatively suggest that (for our assumptions) this transition occurs at a critical particle number of order 100.

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.

Nova LMC 2009a is confirmed as a Recurrent Nova (RN) from positional coincidence with nova LMC 1971b. The observational data set is one of the most comprehensive for any Galactic or extragalactic RN: optical and near-IR photometry from outburst until over 6 years later; optical spectra for the first 6 months, and Swift satellite Ultraviolet and X-ray observations from 9 days to almost 1 year post-outburst. We find $M_V = -8.4\pm0.8_{\mathrm{r}}\pm0.7_{\mathrm{s}}$ and expansion velocities between 1000 and 4000 km s$^{-1}$. Coronal line emission before day 9 indicates shocks in the ejecta. Strengthening of He II $\lambda$4686 preceded the emergence of the Super-Soft Source (SSS) in X-rays at $\sim63-70$ days, which was initially very variable. Periodic modulations, $P=1.2$ days, most probably orbital in nature, were evident in the UV and optical from day 43. Subsequently, the SSS shows an oscillation with the same period but with a delay of 0.28P. The progenitor system has been identified; the secondary is most likely a sub-giant feeding a luminous accretion disk. Properties of the SSS infer a white dwarf (WD) mass $1.1 \mathrm{M}_\odot \lesssim M_{\rm WD} \lesssim 1.3 \mathrm{M}_\odot$. If the accretion occurs at constant rate, $\dot{\it{M}}_{\rm acc} \simeq 3.6^{+4.7}_{-2.5} \times 10^{-7} \mathrm{M}_\odot$ yr$^{-1}$ is needed, consistent with nova models for an inter-eruption interval of 38 years, low outburst amplitude, progenitor position in the color-magnitude diagram, and spectral energy distribution at quiescence. We note striking similarities between LMC 2009a and the Galactic nova KT Eri, suggesting that KT Eri is a candidate RN.

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.

Explaining the origin and evolution of exoplanetary "hot Jupiters" remains a significant challenge. One possible mechanism for their production is planet-planet interactions, which produces hot Jupiters from planets born far from their host stars but near their dynamical stability limits. In the much more likely case of planets born far from their dynamical stability limits, can hot Jupiters can be formed in star clusters? Our N-body simulations of planetary systems inside star clusters answer this question in the affirmative, and show that hot Jupiter formation is not a rare event. We detail three case studies of the dynamics-induced births of hot Jupiters on highly eccentric orbits that can only occur inside star clusters. The hot Jupiters' orbits bear remarkable similarities to those of some of the most extreme exoplanets known: HAT-P-32 b, HAT-P-2 b, HD 80606 b and GJ 876 d. If stellar perturbations formed these hot Jupiters then our simulations predict that these very hot, inner planets are often accompanied by much more distant gas giants in highly eccentric orbits.

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.

Preliminary results of the ongoing search for symbiotic binary stars in the Local Group of Galaxies are presented and discussed.

Can a white dwarf, accreting hydrogen-rich matter from a non-degenerate companion star, ever exceed the Chandrasekhar mass and explode as a type Ia supernova? We explore the range of accretion rates that allow a white dwarf (WD) to secularly grow in mass, and derive limits on the accretion rate and on the initial mass that will allow it to reach $1.4M_\odot$ --- the Chandrasekhar mass. We follow the evolution through a long series of hydrogen flashes, during which a thick helium shell accumulates. This determines the effective helium mass accretion rate for long-term, self-consistent evolutionary runs with helium flashes. We find that net mass accumulation always occurs despite helium flashes. Although the amount of mass lost during the first few helium shell flashes is a significant fraction of that accumulated prior to the flash, that fraction decreases with repeated helium shell flashes. Eventually no mass is ejected at all during subsequent flashes. This unexpected result occurs because of continual heating of the WD interior by the helium shell flashes near its surface. The effect of heating is to lower the electron degeneracy throughout the WD, and especially in the outer layers. This key result yields helium burning that is quasi-steady state, instead of explosive. We thus find a remarkably large parameter space within which long-term, self-consistent simulations show that a WD can grow in mass and reach the Chandrasekhar limit, despite its helium flashes.

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.

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.

We present our wide field Halpha+N[II] observations of 15 cataclysmic variables to search for remnant nova shells. Such shells have been found around other cataclysmic variables that were hitherto not known as novae. Our candidates were selected as objects in the period regime of high-mass transfer systems that experience - at least occasionally - low mass transfer rates. The fact that we find no indication of a nova shell in any of these systems allows us to set a lower limit of 13000 years to the recurrence time of these objects.

We present the first spectrum of the most luminous infrared star in M33, and use it to demonstrate that the object is almost certainly a binary composed of a massive O star and a dust-enshrouded Red Hypergiant. This is the most luminous symbiotic binary ever discovered. Its radial velocity is an excellent match to that of the hydrogen gas in the disk of M33, supporting our interpretation that it is a very young and massive binary star.

A definitive determination of the progenitors of type Ia supernovae (SNIa) has been a conundrum for decades. The single degenerate scenario $-$ a white dwarf (WD) in a semi-detached binary system accreting mass from its secondary $-$ is a plausible path; however, no simulation to date has shown that such an outcome is possible. In this study, we allowed a WD with a near Chandrasekhar mass of $1.4M_\odot$ to evolve over tens of thousands of nova cycles, accumulating mass secularly while undergoing periodic nova eruptions. We present the mass accretion limits within which a SNIa can possibly occur. The results showed, for each parameter combination within the permitted limits, tens of thousands of virtually identical nova cycles where the accreted mass exceeded the ejected mass, i.e. the WD grew slowly but steadily in mass. Finally, the WD became unstable, the maximal temperature rose by nearly two orders of magnitude, heavy element production was enhanced by orders of magnitude and the nuclear and neutrino luminosities became enormous. We also found that this mechanism leading to WD collapse is robust, with WDs in the range $1.0 - 1.38M_\odot$, and an accretion rate of $5*10^{-7}M_\odot/yr$, all growing steadily in mass. These simulations of the onset of a SNIa event make observationally testable predictions about the light curves of pre-SN stars, and about the chemistry of SNIa ejecta.

Symbiotic binaries are putative progenitors of type Ia supernovae. The census of Galactic symbiotic binaries is so incomplete that we cannot reliably estimate the total population of these stars, and use it to check whether that number is consistent with the observed type Ia supernova rate in spiral galaxies. We have thus begun a survey of the nearest counterpart of our own Galaxy, namely M31, where a relatively complete census of symbiotic stars is achievable. We report the first detections and spectrographic characterizations of 35 symbiotic binaries in M31, and compare these stars with the symbiotic population in the Milky Way. These newly detected M31 symbiotic binaries are remarkably similar to galactic symbiotics, though we are clearly only sampling (in this feasibility study) the most luminous symbiotics in M31. We have also found, in M31, the symbiotic star (M31SyS J004233.17+412720.7) with the highest ionization level known amongst all symbiotics. An optical outburst of the M31 symbiotic star M31SyS J004322.50+413940.9 was probably a nova-like outburst, the first symbiotic outburst detected outside the Milky Way and Magellanic Clouds.

We present Hubble Space Telescope (HST) imaging observations of the site of the Type-Ia supernova SN2011fe in the nearby galaxy M101, obtained about one year prior to the event, in a narrow band centred on the HeII 4686 \AA{} emission line. In a "single-degenerate" progenitor scenario, the hard photon flux from an accreting white dwarf (WD), burning hydrogen on its surface over $\sim1$ Myr should, in principle, create a HeIII Str\"{o}mgren sphere or shell surrounding the WD. Depending on the WD luminosity, the interstellar density, and the velocity of an outflow from the WD, the HeIII region could appear unresolved, extended, or as a ring, with a range of possible surface brightnesses. We find no trace of HeII 4686 \AA{} line emission in the HST data. Using simulations, we set $2\sigma$ upper limits on the HeII 4686 \AA{} luminosity of $L_{\rm HeII} < 3.4 \times 10^{34}$ erg s$^{-1}$ for a point source, corresponding to an emission region of radius $r < 1.8$ pc. The upper limit for an extended source is $L_{\rm HeII} < 1.7 \times 10^{35}$ erg s$^{-1}$, corresponding to an extended region with $r\sim11$ pc. The largest detectable shell, given an interstellar-medium density of 1 cm$^{-3}$, has a radius of $\sim6$ pc. Our results argue against the presence, within the $\sim10^5$ yr prior to the explosion, of a supersoft X-ray source of luminosity $L_{\rm bol} \ge 3 \times 10^{37}$ erg s$^{-1}$, or of a super-Eddington accreting WD that produces an outflowing wind capable of producing cavities with radii of 2-6 pc.

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.

Using deep Hubble Space Telescope imaging, color-magnitude diagrams are constructed for the globular clusters 47 Tuc and NGC 6397. As expected, because of its lower metal abundance, the main sequence of NGC 6397 lies well to the blue of that of 47 Tuc. A comparison of the white dwarf cooling sequences of the two clusters, however, demonstrates that these sequences are indistinguishable over most of their loci - a consequence of the settling out of heavy elements in the dense white dwarf atmosphere and the near equality of their masses. Lower quality data on M4 continues this trend to a third cluster whose metallicity is intermediate between these two. While the path of the white dwarfs in the color-magnitude diagram is nearly identical in 47 Tuc and NGC 6397, the numbers of white dwarfs along the path are not. This results from the relatively rapid relaxation in NGC 6397 compared to 47 Tuc and provides a cautionary note that simply counting objects in star clusters in random locations as a method of testing stellar evolutionary theory is likely dangerous unless dynamical considerations are included.

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.

Globular clusters trace the formation history of the spheroidal components of both our Galaxy and others, which represent the bulk of star formation over the history of the universe. They also exhibit a range of metallicities, with metal-poor clusters dominating the stellar halo of the Galaxy, and higher metallicity clusters found within the inner Galaxy, associated with the stellar bulge, or the thick disk. Age differences between these clusters can indicate the sequence in which the components of the Galaxy formed, and in particular which clusters were formed outside the Galaxy and later swallowed along with their original host galaxies, and which were formed in situ. Here we present an age determination of the metal-rich globular cluster 47 Tucanae by fitting the properties of the cluster white dwarf population, which implies an absolute age of 9.9 (0.7) Gyr at 95% confidence. This is about 2.0 Gyr younger than inferred for the metal-poor cluster NGC 6397 from the same models, and provides quantitative evidence that metal-rich clusters like 47 Tucanae formed later than the metal-poor halo clusters like NGC 6397.

Based on the width of its main sequence, and an actual observed split when viewed through particular filters, it is widely accepted that 47 Tucanae contains multiple stellar populations. In this contribution, we divide the main-sequence of 47 Tuc into four color groups, which presumably represent stars of various chemical compositions. The kinematic properties of each of these groups is explored via proper-motions, and a strong signal emerges of differing proper-motion anisotropies with differing main-sequence color; the bluest main-sequence stars exhibit the largest proper-motion anisotropy which becomes undetectable for the reddest stars. In addition, the bluest stars are also the most centrally concentrated. A simiilar analysis for SMC stars, which are located in the background of 47 Tuc on our frames, yields none of the anisotropy exhibited by the 47 Tuc stars. We discuss implications of these results for possible formation scenarios of the various populations.

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.

We present a new measurement of the stellar initial mass function (IMF) based on ultra-deep, high-resolution photometry of >5,000 stars in the outskirts of the Small Magellanic Cloud (SMC) galaxy. The Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) observations reveal this rich, co-spatial population behind the foreground globular cluster 47 Tuc, which we targeted for 121 HST orbits. The stellar main sequence of the SMC is measured in the F606W, F814W color-magnitude diagram (CMD) down to ~30th magnitude, and is cleanly separated from the foreground star cluster population using proper motions. We simulate the SMC population by extracting stellar masses (single and unresolved binaries) from specific IMFs, and converting those masses to luminosities in our bandpasses. The corresponding photometry for these simulated stars is drawn directly from a rich cloud of 4 million artificial stars, thereby accounting for the real photometric scatter and completeness of the data. Over a continuous and well populated mass range of M = 0.37 - 0.93 Msun (i.e., down to a ~75% completeness limit at F606W = 28.7), we demonstrate that the IMF is well represented by a single power-law form with slope \alpha = -1.90 (^{+0.15}_{-0.10}) (3 sigma error) (i.e., dN/dM \propto M^{\alpha}). This is shallower than the Salpeter slope of \alpha = -2.35, which agrees with the observed stellar luminosity function at higher masses. Our results indicate that the IMF does {\it not} turn over to a more shallow power-law form within this mass range. We discuss implications of this result for the theory of star formation, the inferred masses of galaxies, and the (lack of a) variation of the IMF with metallicity.

Multi-epoch observations with ACS on HST provide a unique and comprehensive probe of stellar dynamics within NGC 6397. We are able to confront analytic models of the globular cluster with the observed stellar proper motions. The measured proper motions probe well along the main sequence from 0.8 to below 0.1 M$_\odot$ as well as white dwarfs younger than one gigayear. The observed field lies just beyond the half-light radius where standard models of globular cluster dynamics (e.g. based on a lowered Maxwellian phase-space distribution) make very robust predictions for the stellar proper motions as a function of mass. The observed proper motions show no evidence for anisotropy in the velocity distribution; furthermore, the observations agree in detail with a straightforward model of the stellar distribution function. We do not find any evidence that the young white dwarfs have received a natal kick in contradiction with earlier results. Using the observed proper motions of the main-sequence stars, we obtain a kinematic estimate of the distance to NGC 6397 of $2.2^{+0.5}_{-0.7}$ kpc and a mass of the cluster of $1.1 \pm 0.1 \times 10^5 \mathrm{M}_\odot$ at the photometric distance of 2.53 kpc. One of the main-sequence stars appears to travel on a trajectory that will escape the cluster, yielding an estimate of the evaporation timescale, over which the number of stars in the cluster decreases by a factor of e, of about 3 Gyr. The proper motions of the youngest white dwarfs appear to resemble those of the most massive main-sequence stars, providing the first direct constraint on the relaxation time of the stars in a globular cluster of greater than or about 0.7 Gyr.

We report on the results of a direct N-body simulation of a star cluster that started with N = 200 000, comprising 195 000 single stars and 5 000 primordial binaries. The code used for the simulation includes stellar evolution, binary evolution, an external tidal field and the effects of two-body relaxation. The model cluster is evolved to 12 Gyr, losing more than 80% of its stars in the process. It reaches the end of the main core-collapse phase at 10.5 Gyr and experiences core oscillations from that point onwards -- direct numerical confirmation of this phenomenon. However, we find that after a further 1 Gyr the core oscillations are halted by the ejection of a massive binary comprised of two black holes from the core, producing a core that shows no signature of the prior core-collapse. We also show that the results of previous studies with N ranging from 500 to 100 000 scale well to this new model with larger N. In particular, the timescale to core-collapse (in units of the relaxation timescale), mass segregation, velocity dispersion, and the energies of the binary population all show similar behaviour at different N.

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.

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.

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.

We present the results of a multi-wavelength (FUV to I-band) survey of the stellar populations of the globular cluster NGC 6752, using STIS, ACS and WFC3 on board the Hubble Space Telescope. We have confirmed that two previously identified CV candidates are, in fact, dwarf novae which underwent outbursts during our observations. We have also identified previously unknown optical counterparts to two X-ray sources. We estimate the position of the centre of the cluster, and show that the stellar density profile is not well described by a single King model, indicating that this cluster is in a core-collapsed or post-core collapse phase. The colour-magnitude diagram shows a well-populated horizontal branch, numerous blue stragglers and white dwarfs (WDs), as well as 87 sources in the gap region where we expect to find WD - main sequence binaries, including cataclysmic variables (CVs). The X-ray sources and WD binary systems are the most centrally concentrated populations, with dynamically estimated characteristic masses >1.1Msun and >0.8Msun, respectively.

Wolf-Rayet (WR) stars are the evolved descendants of massive O-type stars and are considered to be progenitor candidates for Type Ib/c core-collapse supernovae (SNe). Recent results of our HST/WFC3 survey of Wolf-Rayet stars in M101 are summarised based on the detection efficiency of narrow-band optical imaging compared to broad-band methods. Weshow that on average of 42% WR stars, increasing to ~85% in central regions, are only detected in the narrow-band imaging. Hence, the non-detection of a WR star at the location of ~10 Type Ib/c SNe in broad-band imaging is no longer strong evidence for a non-WR progenitor channel.

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.

We report parallax measurements for 70 ultracool dwarfs (UCDs). Using both literature values and our sample, we report new polynomial relations between spectral type and M$_{JHK}$. Including resolved L/T transition binaries in the relations, we find no reason to differentiate between a "bright" (unresolved binary) and "faint" (single source) sample across the L/T boundary. Isolating early T dwarfs, we find that the brightening of T0-T4 sources is prominent in M$_{J}$ where there is a [1.2 - 1.4] magnitude difference. A similar yet dampened brightening of [0.3 - 0.5] magnitude happens at M$_{H}$ and a plateau or dimming of [-0.2 - -0.3] magnitude is seen in M$_{K}$. Comparing with evolutionary models that vary gravity, metallicity, and cloud thickness we find that a near constant temperature of 1200 $\pm$100 K along a narrow spectral subtype of T0-T4 is required to account for the brightening and color magnitude diagram of the L-dwarf/T-dwarf transition. Furthermore, there is a significant population of both L and T dwarfs which are red or potentially "ultra-cloudy" compared to the models, many of which are known to be young indicating a correlation between enhanced photospheric dust and youth. For the low surface-gravity or young companion L dwarfs we find that 8 out of 10 are at least [0.2-1.0] magnitude underluminous in M$_{JH}$ and/or M$_{K}$ compared to equivalent spectral type objects. We speculate that this is a consequence of increased dust opacity and conclude that low-surface gravity L dwarfs require a completely new spectral-type/absolute magnitude polynomial for analysis.

On August 24 (UT) the Palomar Transient Factory (PTF) discovered PTF11kly (SN 2011fe), the youngest and most nearby type Ia supernova (SN Ia) in decades. We followed this event up in the radio (centimeter and millimeter bands) and X-ray bands, starting about a day after the estimated explosion time. We present our analysis of the radio and X-ray observations, yielding the tightest constraints yet placed on the pre-explosion mass-loss rate from the progenitor system of this supernova. We find a robust limit of dM/dt<10^-8 (w/100 km/s) [M_solar/yr] from sensitive X-ray non-detections, as well as a similar limit from radio data, which depends, however, on assumptions about microphysical parameters. We discuss our results in the context of single-degenerate models for SNe Ia and find that our observations modestly disfavor symbiotic progenitor models involving a red giant donor, but cannot constrain systems accreting from main-sequence or sub-giant stars, including the popular supersoft channel. In view of the proximity of PTF11kly and the sensitivity of our prompt observations we would have to wait for a long time (decade or longer) in order to more meaningfully probe the circumstellar matter of Ia supernovae.

We present a new distance determination to the Galactic globular cluster 47 Tucanae by fitting the spectral energy distributions of its white dwarfs to pure hydrogen atmosphere white dwarf models. Our photometric dataset is obtained from a 121 orbit Hubble Space Telescope program using the Wide Field Camera 3 UVIS/IR channels, capturing F390W, F606W, F110W, and F160W images. These images cover more than 60 square arcmins and extend over a radial range of 5-13.7 arcmin (6.5-17.9 pc) within the globular cluster. Using a likelihood analysis, we obtain a best fitting unreddened distance modulus of (m - M)o=13.36+/-0.02+/-0.06 corresponding to a distance of 4.70+/-0.04+/-0.13 kpc, where the first error is random and the second is systematic. We also search the white dwarf photometry for infrared excess in the F160W filter, indicative of debris disks or low mass companions, and find no convincing cases within our sample.

In HST Cycle 17, we imaged the well known globular star cluster 47 Tucanae for 121 orbits using ACS and both the UVIS and IR channels of WFC3 (GO-11677, PI - H. Richer). This unique data set was obtained to address many scientific questions that demand a very deep, panchromatic, and panoramic view of the cluster's stellar populations. In total, the program obtained over 0.75 Ms of imaging exposure time with the three HST cameras, over a time span of 9 months in 2010. The primary ACS field was imaged in the two broadband filters F606W and F814W filters, at 13 orientations, for all 121 orbits. The parallel WFC3 imaging provides a panchromatic (0.4 - 1.7 micron) and contiguous imaging swath over a 250 degree azimuthal range at impact radii of 6.5 -- 17.9 pc in 47 Tuc. This imaging totals over 60 arcmin^2 in area and utilizes the F390W and F606W broadband filters on WFC3/UVIS and the F110W and F160W broadband filters on WFC3/IR. In this paper, we describe the observational design of the new survey and one of the methods used to analyze all of the imaging data. This analysis combines over 700 full-frame images taken with the three HST cameras into a handful of ultra-deep, well-sampled combined images in each of the six filters. The results reveal unprecedented color-magnitude diagrams (CMDs) of the cluster extending to >30th magnitude in the optical, 29th magnitude in the UV, and 27th magnitude in the IR. The data set provides a characterization of the complete stellar populations of 47 Tuc, extending from the faintest hydrogen burning dwarfs through the main-sequence and giant branches, down to very cool white dwarf remnants in the cluster. The imaging also provides the deepest probe of the stellar populations of the background Small Magellanic Cloud (SMC) galaxy, resolving low mass main-sequence dwarfs with M < 0.2 Msun.

Type Ia supernovae (SNe) serve as a fundamental pillar of modern cosmology, owing to their large luminosity and a well-defined relationship between light-curve shape and peak brightness. The precision distance measurements enabled by SNe Ia first revealed the accelerating expansion of the universe, now widely believed (though hardly understood) to require the presence of a mysterious "dark" energy. General consensus holds that Type Ia SNe result from thermonuclear explosions of a white dwarf (WD) in a binary system; however, little is known of the precise nature of the companion star and the physical properties of the progenitor system. Here we make use of extensive historical imaging obtained at the location of SN 2011fe/PTF11kly, the closest SN Ia discovered in the digital imaging era, to constrain the visible-light luminosity of the progenitor to be 10-100 times fainter than previous limits on other SN Ia progenitors. This directly rules out luminous red giants and the vast majority of helium stars as the mass-donating companion to the exploding white dwarf. Any evolved red companion must have been born with mass less than 3.5 times the mass of the Sun. These observations favour a scenario where the exploding WD of SN 2011fe/PTF11kly, accreted matter either from another WD, or by Roche-lobe overflow from a subgiant or main-sequence companion star.

We report a nearby L7.5 dwarf discovered using the Preliminary Data Release of the Wide-field Infrared Survey Explorer (WISE) and the Two Micron All-Sky Survey (2MASS). WISEP J180026.60+013453.1 has a motion of 0.42 arcsec/yr and an estimated distance of 8.8 \pm 1.0 pc. With this distance, it currently ranks as the sixth closest known L dwarf, although a trigonometric parallax is needed to confirm this distance. It was previously overlooked because it lies near the Galactic Plane (b=12). As a relatively bright and nearby late L dwarf with normal near-infrared colors, W1800+0134 will serve as a benchmark for studies of cloud-related phenomena in cool substellar atmospheres.

The eruption of the recurrent nova U Scorpii on 28 January 2010 is now the all-time best observed nova event. We report 36,776 magnitudes throughout its 67 day eruption, for an average of one measure every 2.6 minutes. This unique and unprecedented coverage is the first time that a nova has any substantial amount of fast photometry. With this, two new phenomena have been discovered: the fast flares in the early light curve seen from days 9-15 (which have no proposed explanation) and the optical dips seen out of eclipse from days 41-61 (likely caused by raised rims of the accretion disk occulting the bright inner regions of the disk as seen over specific orbital phases). The expanding shell and wind cleared enough from days 12-15 so that the inner binary system became visible, resulting in the sudden onset of eclipses and the turn-on of the supersoft X-ray source. On day 15, a strong asymmetry in the out-of-eclipse light points to the existence of the accretion stream. The normal optical flickering restarts on day 24.5. For days 15-26, eclipse mapping shows that the optical source is spherically symmetric with a radius of 4.1 R_sun. For days 26-41, the optical light is coming from a rim-bright disk of radius 3.4 R_sun. For days 41-67, the optical source is a center-bright disk of radius 2.2 R_sun. Throughout the eruption, the colors remain essentially constant. We present 12 eclipse times during eruption plus five just after the eruption.

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.

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.

The total masses ejected during classical nova eruptions are needed to answer two questions with broad astrophysical implications: Can accreting white dwarfs be pushed towards the Chandrasekhar mass limit to yield type Ia supernovae? Are Ultra-luminous red variables a new kind of astrophysical phenomenon, or merely extreme classical novae? We review the methods used to determine nova ejecta masses. Except for the unique case of BT Mon (nova 1939), all nova ejecta mass determinations depend on untested assumptions and multi-parameter modeling. The remarkably simple assumption of equipartition between kinetic and radiated energy (E_kin and E_rad, respectively) in nova ejecta has been invoked as a way around this conundrum for the ultra-luminous red variable in M31. The deduced mass is far larger than that produced by any classical nova model. Our nova eruption simulations show that radiation and kinetic energy in nova ejecta are very far from being in energy equipartition, with variations of four orders of magnitude in the ratio E_kin/E_rad being commonplace. The assumption of equipartition must not be used to deduce nova ejecta masses; any such "determinations" can be overestimates by a factor of up to 10,000. We data-mined our extensive series of nova simulations to search for correlations that could yield nova ejecta masses. Remarkably, the mass ejected during a nova eruption is dependent only on (and is directly proportional to) E_rad. If we measure the distance to an erupting nova and its bolometric light curve then E_rad and hence the mass ejected can be directly measured.

T Pyxidis is the prototypical recurrent nova (RN) with a mysterious nova shell. We report new observations of the shell with HST. The knots in the shell are expanding with velocities 500-715 km/s, for a distance of 3500 pc. The fractional expansion of the knots is constant, and this implies no significant deceleration. Hence, the knots were ejected by an eruption close to the year 1866. Knots have turned on after 1995, and this demonstrates that the knots are powered by shocks from the collision of the 1866 ejecta with fast ejecta from later RN eruptions. The 1866 ejecta has a total mass of 10^-4.5 Msun, which with the low ejection velocity shows that the 1866 event was an ordinary nova eruption, not a RN eruption. The accretion rate before the ordinary nova event must have been low (around the 4x10^-11 Msun/yr expected for gravitational radiation alone) and the matter accumulated on the surface of the white dwarf for ~750,000 years. The current accretion rate (>10^-8 Msun/yr) is 1000X higher than expected for a system below the period gap, with the plausible reason being that the 1866 event started a continuing supersoft source that drives the accretion. A key fact about T Pyx is that its accretion rate has been secularly declining since before the 1890 eruption, with the current rate being only 3% of its earlier rate. The decline in the observed accretion rate shows that the supersoft source is not self-sustaining, and we calculate that the accretion in T Pyx will effectively stop in upcoming decades. With this, T Pyx will enter a state of hibernation, lasting for an estimated 2,600,000 years, before gravitational radiation brings the system into contact again. Thus, T Pyx has an evolutionary cycle going from an ordinary CV state, to its current RN state, to a future hibernation state, and then repeating this cycle.

We report on nine wide common proper motion systems containing late-type M, L, or T companions. We confirm six previously reported companions, and identify three new systems. The ages of these systems are determined using diagnostics for both stellar primaries and low--mass secondaries and masses for the secondaries are inferred using evolutionary models. Of our three new discoveries, the M3+T6.5 pair G 204-39 and SDSS J1758+4633 has an age constrained to 0.5-1.5 Gyr making the secondary a potentially useful brown dwarf benchmark. The G5+L4 pair G 200-28 and SDSS J1416+5006 has a projected separation of ~25,000 AU making it one of the widest and lowest binding energy systems known to date. The system containing NLTT 2274 and SDSS J0041+1341 is an older M4+L0 (>4.5 Gyr) pair which shows Halpha activity in the secondary but not the primary making it a useful tracer of age/mass/activity trends. We find a resolved binary frequency for widely-separated (> 100 AU) low--mass companions (i.e. at least a triple system) which is at least twice the frequency found for the field ultracool dwarf population. The ratio of triples to binaries and quadruples to binaries is also high for this sample: 3:5 and 1:4, respectively, compared to 8-parsec sample values of 1:4 and 1:26. The additional components in these wide companion systems indicates a formation mechanism that requires a third or fourth component to maintain gravitational stability or facilitate the exchange of angular momentum. The binding energies for the nine multiples discussed in this text are among the lowest known for wide low-mass systems, suggesting that weakly bound, low--to--intermediate mass (0.2M_sun < M_tot <1.0M_sun) multiples can form and survive to exist in the field (1-8 Gyr).

We demonstrate that stars beyond the virial radii of galaxies may be generated by the gravitational impulse received by a satellite as it passes through the pericenter of its orbit around its parent. These stars may become energetically unbound (escaped stars), or may travel to further than a few virial radii for longer than a few Gyr, but still remain energetically bound to the system (wandering stars). Larger satellites (10-100% the mass of the parent), and satellites on more radial orbits are responsible for the majority of this ejected population. Wandering stars could be observable on Mpc scales via classical novae, and on 100 Mpc scales via SNIa. The existence of such stars would imply a corresponding population of barely-bound, old, high velocity stars orbiting the Milky Way, generated by the same physical mechanism during the Galaxy's formation epoch. Sizes and properties of these combined populations should place some constraints on the orbits and masses of the progenitor objects from which they came, providing insight into the merging histories of galaxies in general and the Milky Way in particular.

The Palomar Transient Factory (PTF) is a fully-automated, wide-field survey aimed at a systematic exploration of the optical transient sky. The transient survey is performed using a new 8.1 square degree camera installed on the 48-inch Samuel Oschin telescope at Palomar Observatory; colors and light curves for detected transients are obtained with the automated Palomar 60-inch telescope. PTF uses eighty percent of the 1.2-m and fifty percent of the 1.5-m telescope time. With an exposure of 60-s the survey reaches a depth of approximately 21.3 in g' and 20.6 in R (5 sigma, median seeing). Four major experiments are planned for the five-year project: 1) a 5-day cadence supernova search; 2) a rapid transient search with cadences between 90 seconds and 1 day; 3) a search for eclipsing binaries and transiting planets in Orion; and 4) a 3-pi sr deep H-alpha survey. PTF provides automatic, realtime transient classification and follow up, as well as a database including every source detected in each frame. This paper summarizes the PTF project, including several months of on-sky performance tests of the new survey camera, the observing plans and the data reduction strategy. We conclude by detailing the first 51 PTF optical transient detections, found in commissioning data.

The Palomar Transient Factory (PTF) is a wide-field experiment designed to investigate the optical transient and variable sky on time scales from minutes to years. PTF uses the CFH12k mosaic camera, with a field of view of 7.9 deg^2 and a plate scale of 1 asec/pixel, mounted on the the Palomar Observatory 48-inch Samuel Oschin Telescope. The PTF operation strategy is devised to probe the existing gaps in the transient phase space and to search for theoretically predicted, but not yet detected, phenomena, such as fallback supernovae, macronovae, .Ia supernovae and the orphan afterglows of gamma-ray bursts. PTF will also discover many new members of known source classes, from cataclysmic variables in their various avatars to supernovae and active galactic nuclei, and will provide important insights into understanding galactic dynamics (through RR Lyrae stars) and the Solar system (asteroids and near-Earth objects). The lessons that can be learned from PTF will be essential for the preparation of future large synoptic sky surveys like the Large Synoptic Survey Telescope. In this paper we present the scientific motivation for PTF and describe in detail the goals and expectations for this experiment.

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.

We present astrometric parallaxes for 18 suspected nearby stars selected from the LSPM-north proper motion catalog. Sixteen objects are confirmed to be main sequence M dwarfs within 16 parsecs of the Sun, including three stars (LSPM J0011+5908, LSPM J0330+5413, LSPM J0510+2714) which lie just within the 10 parsec horizon. Two other targets (LSPM J1817+1328, LSPM J2325+1403) are confirmed to be nearby white dwarfs at distances of 14 and 22 parsecs, respectively. One of our targets, the common proper motion pair LSPM J0405+7116E + LSPM J0405+7116W, is revealed to be a triple system, with the western component resolved into a pair of 16th magnitude stars (LSPM J0405+7116W-A and LSPM J0405+7116W-B) with a 0.7"+/-0.1" angular separation. We find two stars (LSPM J1314+1320 and LSPM J1757+7042) to be significantly overluminous for their colors, and suggest that these may be unresolved doubles/multiples.

Blue hook (BHk) stars are a rare class of horizontal branch stars that so far have been found in only very few Galactic globular clusters (GCs). The dominant mechanism for producing these objects is currently still unclear. In order to test if the presence of BHk populations in a given GC is linked to specific physical or structural cluster properties, we have constructed a parent sample of GCs for which existing data is sufficient to establish the presence or absence of BHk populations with confidence. We then compare the properties of those clusters in our parent sample that do contain a BHk population to those that do not. We find that there is only one compelling difference between BHk and non-BHk clusters: all known BHk clusters are unusually massive. However, we also find that the BHk clusters are consistent with being uniformly distributed within the cumulative mass distribution of the parent sample. Thus, while it is attractive to suggest there is is a lower mass cut-off for clusters capable of forming BHk stars, the data do not require this. Instead, the apparent preference for massive clusters could still be a purely statistical effect: intrinsically rare objects can only be found by searching a sufficiently large number of stars.

A new all-sky catalog of stars with proper motions pm>0.15"/yr is presented. The catalog is largely a product of the SUPERBLINK survey, a data-mining initiative in which the entire Digitized Sky Surveys are searched for moving stellar sources. Findings from earlier proper motions surveys are also incorporated. The new all-sky catalog supersedes the great historic proper motion catalogs assembled by W. J. Luyten (LHS, NLTT), and provides a virtually complete (>98%) census of high proper motion stars down to magnitude R=19.

We report proper motion measurements for 427 late-type M, L and T dwarfs, 332 of which have been measured for the first time. Combining these new proper motions with previously published measurements yields a sample of 841 M7-T8 dwarfs. We combined parallax measurements or calculated spectrophotometric distances and computed tangential velocities for the entire sample. We find that kinematics for the full and volume-limited 20 pc samples are consistent with those expected for the Galactic thin disk, with no significant differences between late-type M, L, and T dwarfs. Applying an age-velocity relation we conclude that the average kinematic age of the 20 pc sample of ultracool dwarfs is older than recent kinematic estimates and more consistent with age results calculated with population synthesis models. There is a statistically distinct population of high tangential velocity sources whose kinematics suggest an even older population of ultracool dwarfs belonging to either the Galactic thick disk or halo. We isolate subsets of the entire sample, including low surface-gravity dwarfs, unusually blue L dwarfs, and photometric outliers in J-Ks color and investigate their kinematics. We find that the spectroscopically distinct class of unusually blue L dwarfs has kinematics clearly consistent with old age, implying that high surface-gravity and/or low metallicity may be relevant to their spectral properties. The low surface-gravity dwarfs are kinematically younger than the overall population, and the kinematics of the red and blue ultracool dwarfs suggest ages that are younger and older than the full sample, respectively. We also present a reduced proper motion diagram at 2MASS Ks for the entire population and find that a limit of H_Ks > 18 excludes M dwarfs from the L and T dwarf population regardless of near-infrared color.

We report parallax and distance estimates for twelve more cataclysmic binaries and related objects observed with the 2.4m Hiltner telescope at MDM Observatory. The final parallax accuracy is typically about 1 mas. For only one of the twelve objects, IR Gem, do we fail to detect a significant parallax. Notable results include distances for V396 Hya (CE 315), a helium double degenerate with a relatively long orbital period, and for MQ Dra (SDSSJ155331+551615), a magnetic system with a very low accretion rate. We find that the Z Cam star KT Persei is physically paired with a K main-sequence star lying 15 arcsec away. Several of the targets have distance estimates in the literature that are based on the white dwarf's effective temperature and flux; our measurements broadly corroborate these estimates, but tend to put the stars a bit closer, indicating that the white dwarfs may have rather larger masses than assumed. As a side note, we briefly describe radial velocity spectroscopy that refines the orbital period of V396 Hya to 65.07 +- 0.08 min.

Using new Chandra X-ray observations and existing XMM-Newton X-ray and Hubble far ultraviolet observations, we aim to detect and identify the faint X-ray sources belonging to the Galactic globular cluster NGC 2808 in order to understand their role in the evolution of globular clusters. We present a Chandra X-ray observation of the Galactic globular cluster NGC 2808. We classify the X-ray sources associated with the cluster by analysing their colours and variability. Previous observations with XMM-Newton and far ultraviolet observations with the Hubble Space Telescope are re-investigated to help identify the Chandra sources associated with the cluster. We compare our results to population synthesis models and observations of other Galactic globular clusters. We detect 113 sources, of which 16 fall inside the half-mass radius of NGC 2808 and are concentrated towards the cluster core. From statistical analysis, these 16 sources are very likely to be linked to the cluster. We detect short-term (1 day) variability in X-rays for 7 sources, of which 2 fall inside the half-mass radius, and long-term (28 months) variability for 10 further sources, of which 2 fall inside the half-mass radius. Ultraviolet counterparts are found for 8 Chandra sources in the core, of which 2 have good matching probabilities and have ultraviolet properties expected for cataclysmic variables. We find one likely neutron star-quiescent low-mass X-ray binary and 7 cataclysmic variable candidates in the core of NGC 2808. The other 8 sources are cataclysmic variable candidates, but some could possibly be active binaries or millisecond pulsars. We find a possible deficit of X-ray sources compared to 47 Tuc which could be related to the metallicity content and the complexity of the evolution of NGC 2808.

We present far-UV spectroscopy obtained with HST for 48 blue objects in the core of 47 Tuc. Based on their position in a FUV-optical colour-magnitude diagram, these were expected to include cataclysmic variables (CVs), blue stragglers (BSs), white dwarfs (WDs) and other exotic objects. For a subset of these sources, we also construct FUV-NIR SEDs. Based on our analysis of this extensive data set, we report the following main results. (1) We spectroscopically confirm 3 previously known or suspected CVs via the detection of emission lines and find new evidence for dwarf nova eruptions in two of these. (2) Only one other source in our spectroscopic sample exhibits marginal evidence for line emission, but predicted and observed CV numbers still agree to within a factor of about 2-3. (3) We have discovered a hot (T_eff = 8700 K), low-mass (M = 0.05 M_sun) secondary star in a previously known 0.8 day binary system. This exotic object is probably the remnant of a subgiant that has been stripped of its envelope and may represent the ``smoking gun'' of a recent dynamical encounter. (4) We have found a Helium WD, the second to be optically detected in 47 Tuc, and the first outside a millisecond-pulsar system. (5) We have discovered a BS-WD binary system, the first known in any globular cluster. (6) We have found two additional candidate WD binary systems with putative main sequence and subgiant companions. (7) We estimate the WD binary fraction in the core of 47 Tuc to be 15 +17/-9 (stat) +8/-7 (sys). (8) One BS in our sample may exceed twice the cluster turn-off mass, but the uncertainties are large. Taken as a whole, our study illustrates the wide range of stellar exotica that are lurking in the cores of GCs, most of which are likely to have undergone significant dynamical encounters. [abridged]

We present N-body models to complement deep imaging of the metal-poor core-collapsed cluster NGC6397 obtained with the Hubble Space Telescope. All simulations include stellar and binary evolution in-step with the stellar dynamics and account for the tidal field of the Galaxy. We focus on the results of a simulation that began with 100000 objects (stars and binaries), 5% primordial binaries and Population II metallicity. After 16 Gyr of evolution the model cluster has about 20% of the stars remaining and has reached core-collapse. We compare the color-magnitude diagrams of the model at this age for the central region and an outer region corresponding to the observed field of NGC6397 (about 2-3 half-light radii from the cluster centre). This demonstrates that the white dwarf population in the outer region has suffered little modification from dynamical processes - contamination of the luminosity function by binaries and white dwarfs with non-standard evolution histories is minimal and should not significantly affect measurement of the cluster age. We also show that the binary fraction of main-sequence stars observed in the NGC6397 field can be taken as representative of the primordial binary fraction of the cluster. For the mass function of the main-sequence stars we find that although this has been altered significantly by dynamics over the cluster lifetime, especially in the central and outer regions, that the position of the observed field is close to optimal for recovering the initial mass function of the cluster stars (below the current turn-off mass). More generally we look at how the mass function changes with radius in a dynamically evolved stellar cluster and suggest where the best radial position to observe the initial mass function is for clusters of any age.

We have obtained deep far- (FUV) and near-ultraviolet (NUV) images of the inner region of the dense globular cluster M15 with the Advanced Camera for Surveys on board the Hubble Space Telescope. The FUV-NUV colour-magnitude diagram shows a well defined track of horizontal branch stars, as well as a trail of blue stragglers and white dwarfs. The main sequence turn-off is clearly visible at FUV~23.5 mag and FUV-NUV~3 mag, and the main sequence stars form a prominent track that extends at least two magnitudes below the main sequence turn-off. As such, this is the deepest FUV-NUV colour-magnitude diagram of a globular cluster presented so far. Cataclysmic variable and blue straggler candidates are the most centrally concentrated stellar populations, which might either be an effect of mass segregation or reflect the preferred birthplace in the dense cluster core of such dynamically-formed objects. We find 41 FUV sources that exhibit significant variability. We classify the variables based on an analysis of their UV colours and variability properties. We find four previously known RR Lyrae and 13 further RR Lyrae candidates, one known Cepheid and six further candidates, six cataclysmic variables, one known and one probable SX Phoenicis star, and the well known low-mass X-ray binary AC211. Our analysis represents the first detection of SX Phoenicis pulsations in the FUV. We find that Cepheids, RR Lyraes and SX Phoenicis exhibit massive variability amplitudes in this waveband (several mags).

We present the CMD from deep HST imaging in the globular cluster NGC 6397. The ACS was used for 126 orbits to image a single field in two colors (F814W, F606W) 5 arcmin SE of the cluster center. The field observed overlaps that of archival WFPC2 data from 1994 and 1997 which were used to proper motion (PM) clean the data. Applying the PM corrections produces a remarkably clean CMD which reveals a number of features never seen before in a globular cluster CMD. In our field, the main sequence stars appeared to terminate close to the location in the CMD of the hydrogen-burning limit predicted by two independent sets of stellar evolution models. The faintest observed main sequence stars are about a magnitude fainter than the least luminous metal-poor field halo stars known, suggesting that the lowest luminosity halo stars still await discovery. At the bright end the data extend beyond the main sequence turnoff to well up the giant branch. A populous white dwarf cooling sequence is also seen in the cluster CMD. The most dramatic features of the cooling sequence are its turn to the blue at faint magnitudes as well as an apparent truncation near F814W = 28. The cluster luminosity and mass functions were derived, stretching from the turn off down to the hydrogen-burning limit. It was well modeled with either a very flat power-law or a lognormal function. In order to interpret these fits more fully we compared them with similar functions in the cluster core and with a full N-body model of NGC 6397 finding satisfactory agreement between the model predictions and the data. This exercise demonstrates the important role and the effect that dynamics has played in altering the cluster IMF.

The current classification system of M stars on the main sequence distinguishes three metallicity classes (dwarfs - dM, subdwarfs - sdM, and extreme subdwarfs - esdM). The spectroscopic definition of these classes is based on the relative strength of prominent CaH and TiO molecular absorption bands near 7000A, as quantified by three spectroscopic indices (CaH2, CaH3, and TiO5). We re-examine this classification system in light of our ongoing spectroscopic survey of stars with proper motion \mu > 0.45 "/yr, which has increased the census of spectroscopically identified metal-poor M stars to over 400 objects. Kinematic separation of disk dwarfs and halo subdwarfs suggest deficiencies in the current classification system. Observations of common proper motion doubles indicates that the current dM/sdM and sdM/esdM boundaries in the [TiO5,CaH2+CaH3] index plane do not follow iso-metallicity contours, leaving some binaries inappropriately classified as dM+sdM or sdM+esdM. We propose a revision of the classification system based on an empirical calibration of the TiO/CaH ratio for stars of near solar metallicity. We introduce the parameter \zeta_{TiO/CaH} which quantifies the weakening of the TiO bandstrength due to metallicity effect, with values ranging from \zeta_{TiO/CaH}=1 for stars of near-solar metallicity to \zeta_{TiO/CaH}~0 for the most metal-poor (and TiO depleted) subdwarfs. We redefine the metallicity classes based on the value of the parameter \zeta_{TiO/CaH}; and refine the scheme by introducing an additional class of ultra subdwarfs (usdM). We introduce sequences of sdM, esdM, and usdM stars to be used as formal classification standards.

We report the discovery of a companion to the high proper motion star LHS 1589, a nearby high-velocity, low-mass subdwarf. The companion (LHS 1589B) is located 0.224+/-0.004" to the southwest of the primary (LHS 1589A), and is 0.5 magnitude fainter than the primary in the K_s band. The pair was resolved with the IRCAL infrared camera at Lick Observatory, operating with the Laser Guide Star Adaptive Optics system. A low-resolution spectrum of the unresolved pair obtained at the MDM observatory shows the source to be consistent with a cool subdwarf of spectral subtype sdK7.5. A photometric distance estimate places the metal-poor system at a distance d=81+/-18 parsecs from the Sun. We also measure a radial velocity V_rad=67+/-8 km/s which, together with the proper motion and estimated distance, suggests that the pair is roaming the inner Galactic halo on a highly eccentric orbit. With a projected orbital separation s=18.1+/-4.8 AU, and a crude estimate of the system's total mass, we estimate the orbital period of the system to be in the range 75 yr < P < 500 yr. This suggests that the dynamical mass of the system could be derived astrometrically, after monitoring the orbital motion over a decade or so. The LHS 1589AB system could thus provide a much needed constraint to the mass-luminosity relationship of metal-poor, low-mass stars.

We investigate the evolution of binary fractions in star clusters using N-body models of up to 100000 stars. Primordial binary frequencies in these models range from 5% to 50%. Simulations are performed with the NBODY4 code and include a full mass spectrum of stars, stellar evolution, binary evolution and the tidal field of the Galaxy. We find that the overall binary fraction of a cluster almost always remains close to the primordial value, except at late times when a cluster is near dissolution. A critical exception occurs in the central regions where we observe a marked increase in binary fraction with time -- a simulation starting with 100000 stars and 5% binaries reached a core binary frequency as high as 40% at the end of the core-collapse phase (occurring at 16 Gyr with ~20000 stars remaining). Binaries are destroyed in the core by a variety of processes as a cluster evolves, but the combination of mass-segregation and creation of new binaries in exchange interactions produces the observed increase in relative number. We also find that binaries are cycled into and out of cluster cores in a manner that is analogous to convection in stars. For models of 100000 stars we show that the evolution of the core-radius up to the end of the initial phase of core-collapse is not affected by the exact value of the primordial binary frequency (for frequencies of 10% or less). We discuss the ramifications of our results for the likely primordial binary content of globular clusters.

We present the results of a deep Hubble Space Telescope (HST) exposure of the nearby globular cluster NGC6397, focussing attention on the cluster's white dwarf cooling sequence. This sequence is shown to extend over 5 magnitudes in depth, with an apparent cutoff at magnitude F814W=27.6. We demonstrate, using both artificial star tests and the detectability of background galaxies at fainter magnitudes, that the cutoff is real and represents the truncation of the white dwarf luminosity function in this cluster. We perform a detailed comparison between cooling models and the observed distribution of white dwarfs in colour and magnitude, taking into account uncertainties in distance, extinction, white dwarf mass, progenitor lifetimes, binarity and cooling model uncertainties. After marginalising over these variables, we obtain values for the cluster distance modulus and age of \mu_0 = 12.02 \pm 0.06 and T_c = 11.47 \pm 0.47Gyr (95% confidence limits). Our inferred distance and white dwarf initial-final mass relations are in good agreement with other independent determinations, and the cluster age is consistent with, but more precise than, prior determinations made using the main sequence turnoff method. In particular, within the context of the currently accepted \Lambda CDM cosmological model, this age places the formation of NGC6397 at a redshift z=3, at a time when the cosmological star formation rate was approaching its peak.

NGC 6397 is the second closest globular star cluster to the Sun. Using 5 days of time on the Hubble Space Telescope, we have constructed the deepest ever color-magnitude diagram for this cluster. We see a clear truncation in each of its two major stellar sequences. Faint red main sequence stars run out well above our observational limit and near to the theoretical prediction for the lowest mass stars capable of stable hydrogen-burning in their cores. We also see a truncation in the number counts of faint blue stars, namely white dwarfs. This reflects the limit to which the bulk of the white dwarfs can cool over the lifetime of the cluster. There is also a turn towards bluer colors in the least luminous of these objects. This was predicted for the very coolest white dwarfs with hydrogen-rich atmospheres as the formation of H2 causes their atmospheres to become largely opaque to infrared radiation due to collision-induced absorption.

As a by-product of high-precision, ultra-deep stellar photometry in the Galactic globular cluster NGC 6397 with the Hubble Space Telescope, we are able to measure a large population of background galaxies whose images are nearly point-like. These provide an extragalactic reference frame of unprecedented accuracy, relative to which we measure the most accurate absolute proper motion ever determined for a globular cluster. We find mu_alpha = 3.56 +/- 0.04 mas/yr and mu_delta = -17.34 +/- 0.04 mas/yr. We note that the formal statistical errors quoted for the proper motion of NGC 6397 do not include possible unavoidable sources of systematic errors, such as cluster rotation. These are very unlikely to exceed a few percent. We use this new proper motion to calculate NGC 6397's UVW space velocity and its orbit around the Milky Way, and find that the cluster has made frequent passages through the Galactic disk.

RBS (Rosat Bright Source) 0490 is a cataclysmic variable star (CV) with unusually strong emission lines. The strength of the emission lines led to a suggestion that the object is intrinsically faint and correspondingly nearby (approx. 33 pc), which, if true, would strongly affect estimates of the CV space density. Here we report astrometry, filter photometry, and time-series spectroscopy of this object. The astrometry gives an absolute parallax pi(abs) = 4.5 +- 1.5 mas and a relative proper motion of 102 mas/yr. A Bayesian procedure gives a very uncertain distance estimate of d approx. 300 pc, and the small parallax alone implies d > 133 pc (at two standard deviations). The mean V magnitude is 17.4, which implies M_V = 10.9 - 5 log (d/200 pc), neglecting extinction. At 200 pc, the space velocity would be over 90 km/s with respect to the LSR. The time-series spectroscopy shows a possible emission-line radial-velocity period near 46 min. This would be unusually short for an orbital period and it may represent some other clock in the system.

The Catalog and Atlas of Cataclysmic Variables has been a staple of the CV community for over 10 years. The catalog has grown from 751 CVs in 1993 to 1600 CVs at present. The catalog became a ``living'' edition in 2001, and its contents have been continually updated since that time. Effective 27 January 2006, the catalog will transition to an archival site, with no further updates to its contents. While it is antipicated that the site will remain active, we present the complete contents of the site as a precaution against a loss of the on-line data.

We show that the X-ray source W31 in the core of the globular cluster 47 Tucanae is physically associated with the bright blue straggler BSS-7. The two sources are astrometrically matched to 0.061\arcsec, with a chance coincidence probability of less than 1%. We then analyse optical time-series photometry obtained with the {\em Hubble Space Telescope} (HST) and find that BSS-7 displays a 1.56 day periodic signal in the I band. We also construct a broad-band (far-ultraviolet through far-red) spectral energy distribution for BSS-7 and fit this with single and binary models. The binary model is a better fit to the data, and we derive the corresponding stellar parameters. All of our findings are consistent with BSS-7 being a detached binary consisting of a blue straggler primary with an X-ray-active, upper-main-sequence companion. The formation of such a system would necessarily involve at least three stars, which is consistent with recent N-body models in which blue stragglers often form via multiple encounters that can involve both single and binary stars. However, we cannot yet entirely rule out the possibility that BSS-7 descended directly from a binary system via mass transfer. The system parameters needed to distinguish definitively between these scenarios may be obtainable from time-resolved spectroscopy.

We present a spectroscopic analysis of HST/STIS and FOS low- and intermediate-resolution spectroscopy of 55 stars (turn-off stars, horizontal branch stars and blue stragglers) in four globular clusters (47 Tucanae, M3, NGC6752, and NGC6397). Stars were analyzed with non-Local Thermodynamic Equilibrium model atmospheres, and values for their effective temperatures and gravities and some rotation rates were obtained. Using photometric fluxes, we also obtained radii, luminosities and spectroscopic masses.

We have used the Advanced Camera for Surveys on board the Hubble Space Telescope to image the core of the globular cluster M15 in the far-ultraviolet (FUV) waveband. Based on these observations, we identify the FUV counterpart of the recently discovered low-mass X-ray binary M15 X-2. Our time-resolved FUV photometry shows a modulation with 0.062+/-0.004 mag semi-amplitude and we clearly detect a period of 22.5806+/-0.0002 min. We have carried out extensive Monte Carlo simulations which show that the signal is consistent with being coherent over the entire observational time range of more than 3000 cycles. This strongly suggests that it represents the orbital period of the binary system. M15 X-2 is FUV bright (approx. 17 mag) and is characterized by an extremely blue spectral energy distribution (F_lambda ~ lambda^{-2.0}). We also find evidence for an excess of flux between 1500 and 1600 \AA and probably between 1600 and 2000 \AA, which might be due to CIV 1550 and HeII 1640 emission lines. We also show that M15 X-2's X-ray luminosity can be powered by accretion at the rate expected for gravitational-wave-driven mass transfer at this binary period. The observed FUV emission appears to be dominated by an irradiated accretion disk around the neutron star primary, and the variability can be explained by irradiation of the low-mass white dwarf donor if the inclination of the system is approx. 34 degree. We conclude that all observational characteristics of M15 X-2 are consistent with it being an ultracompact X-ray binary, only the third confirmed such object in a globular cluster.

We report the discovery of PM J13420-3415, a faint (V=17) white dwarf with a very high proper motion mu = 2.55 arcsec/yr. The star was found in the southern sky extension of the SUPERBLINK proper motion survey. A red spectrum shows the classical signature of a DA white dwarf, with a weak H-alpha line in absorption as the only prominent feature. The star is also found to have a large radial velocity V_rad = +212+/-15 km/s. At the adopted distance of 18 pc, the star has a very large space motion of 313 km/s relative to the Sun. An integration of the space motion shows that the star is on a nearly polar Galactic orbit, and is thus an unambiguous member of the Galactic halo. However, with an estimated effective temperature 5,000K < T_eff < 5,500K, the white dwarf appears to be much younger than expected for a denizen of the halo. The apparent paradox can be explained if the white dwarf is the relatively young (~2 Gyr) remnant of a longer-lived (10-14 Gyr) main sequence star, in which case the object is predicted to be a low-mass white dwarf with M ~ 0.45 M_sun.

Large populations of cataclysmic variables (CVs) in globular clusters have long been predicted, but the number of absolutely certain cluster CVs known in globulars is still less than 10. HST and Chandra observers have recently found dozens of cataclysmic variable candidates in several populous globular clusters. Confirmation and characterization of these candidates are extremely difficult, thus identification of unambiguous CVs remains important. We have searched all archival HST images of the dense globular cluster M80 for erupting dwarf novae (DN), and to check the outburst behaviors of two very blue objects first identified a decade ago. Two new erupting dwarf novae were found in 8 searched epochs, making M80 a record holder for erupting DN. The quiescent classical nova in M80 varies by no more than a few tenths of a magnitude on timescales of minutes to years, and a similar faint, blue object varies by a similar amount. Simulations and completeness tests indicate that there are at most 3 erupting DN like SS Cyg and at most 9 U Gem-like DN in M80. Either this very dense cluster contains about an order of magnitude fewer CVs than theory predicts, or most M80 CVs are extremely faint and/or erupt very infrequently like WZ Sge. We have detected a sequence of 54 objects running parallel to the main sequence and several tenths of a magnitude bluewards of it. These blue objects are significantly more centrally concentrated than the main sequence stars, but not as centrally concentrated as the blue stragglers. We suggest that these objects are white dwarf--red dwarf binaries, and that some are the faint CV population of M80.

NGC 6397 is the closest globular cluster, and hence the ideal place to search for faint stellar populations such as cataclysmic variables (CVs). HST and Chandra observers have identified nine certain and likely CVs in this nearby cluster, including several magnetic CV candidates. We have combined our recent UV imagery with archival HST images of NGC 6397 to search for new CV candidates and especially to look for dwarf nova-like eruptive events. We find remarkable and somewhat unexpected dwarf nova-like eruptions of the two well-known cataclysmic systems CV2 and CV3. These two objects have been claimed to be {\it magnetic} CVs, as indicated by their helium emission-line spectra. Magnetic fields in CVs are usually expected to prevent the disk instability that leads to dwarf nova eruptions. In fact, most field magnetic CVs are observed to not undergo eruptions. Our observations of the dwarf nova eruptions of CV2 and CV3 can be reconciled with these objects' HeII emission lines if both objects are infrequently-erupting intermediate polars, similar to EX Hya. If this is the case for most globular cluster CVs then we can reconcile the many X-ray and UV bright CV candidates seen by Chandra and HST with the very small numbers of erupting dwarf novae observed in cluster cores.

We present a reanalysis of far-ultraviolet (FUV) observations of the globular cluster NGC 2808 obtained with the Hubble Space Telescope. These data were first analyzed by Brown and coworkers, with an emphasis on the bright, blue horizontal branch (HB) stars in this cluster. Here, our focus is on the population of fainter FUV sources, which include white dwarfs (WDs), blue stragglers (BSs) and cataclysmic variables (CVs). We have therefore constructed the deepest FUV-NUV colour-magnitude diagram of NGC 2808 and also searched for variability among our FUV sources. Overall, we have found approx. 40 WD, approx. 60 BS and approx. 60 CV candidates; three of the BSs and two of the CV candidates are variable. We have also recovered a known RR Lyrae star in the core of NGC 2808, which exhibits massive (approx. 4 mag) FUV variability. We have investigated the radial distribution and found that our CV and BS candidates are more centrally concentrated than the HBs and WD candidates. This might be an effect of mass segregation, but could as well be due to the preferential formation of such dynamically-formed objects in the dense cluster core. For one of our CV candidates we found a counterpart in WFPC2 optical data published by Piotto and coworkers.

This paper is a sequel to an earlier paper devoted to multiple, multicycle nova evolution models (Prialnik & Kovetz 1995, first paper of the series), which showed that the different characteristics of nova outbursts can be reproduced by varying the values of three basic and independent parameters: the white dwarf mass-M_{WD}, the temperature of its isothermal core-T_{WD} and the mass transfer rate-Mdot. Apart from being the largest computational classical novae parameter-space survey, we show here that the parameter space is constrained by several analytical considerations and find its limiting surfaces. Consequently, we extend the grid of multicycle nova evolution models presented in Paper I almost to its limits, adding multicycle nova outburst calculations for a considerable number of new parameter combinations. In particular, the extended parameter pace that produces nova eruptions includes low mass transfer rates down to 5e-13 Msun/yr, and more models for low T_{WD}. Resulting characteristics of these runs are added to the former parameter combination results, to provide a full grid spanning the entire parameter space for Carbon-Oxygen white dwarfs. The full grid covers the entire range of observed nova characteristics, even those of peculiar objects, which have not been numerically reproduced until now. Most remarkably, runs for very low Mdot lead to very high values for some characteristics, such as outburst amplitude A>~20, high super-Eddington luminosities at maximum, heavy element abundance of the ejecta Z_{ej}~0.63 and high ejected masses m_{ej}~7e-4 Msun.

We are conducting a high angular resolution imaging survey of Galactic Wolf-Rayet stars using the Wide Field and Planetary Camera 2 aboard the Hubble Space Telescope. We have found a small stellar cluster associated with the faint, close pair WR 38 and WR 38a. We present astrometric measurements and photometry in the wide-band F336W (U), F439W (B), and F555W (V) filter system for these cluster and nearby stars. We compare their colors and magnitudes with calibrated model results for reddened stars to identify seven probable main sequence members of the cluster. A least-squares fit of the colors and magnitudes of this set yields a cluster reddening of E(B-V) = 1.45 +/- 0.14 mag and a distance of 3.7^{+3.8}_{-1.2} kpc. We discuss the relationship of this cluster to other objects along the line of sight, and we argue that the distance probably lies in the range 5 - 8 kpc (but is not as great as 14.5 kpc distance recently advocated by Shorlin, Turner, & Pedreros). At a distance of 8 kpc, the cluster would reside in a dense region of the Carina spiral arm, close to a giant molecular cloud and the starburst cluster NGC 3603.

We report the results of a preliminary nova survey of Local Group dwarf ellipticals. We used the Tenagra Observatory to observe M32, NGC 205, NGC 147, and NGC 185 in their entirety every clear night over a 4.5 month interval and discovered one nova in M32 and a candidate nova in NGC 205. The nova in M32 was verified spectroscopically. The nova candidate in NGC 205 had an unusually low peak luminosity (M_V = -5.1), and we were unable to obtain spectroscopic verification. We report a high bulk nova rate for M32 of 2(+2.4,-1.0) yr^-1 and, assuming the candidate nova is correctly identified, for NGC 205 of 2(+2.2,-1.0) yr^-1. If the NGC 205 variable is not a nova, we calculate an upper limit on the bulk nova rate for NGC 205 of 1.5 yr^-1. We report upper limits on the bulk nova rates in NGC 147 of 2 yr^-1 and NGC 185 of 1.8 yr^-1 and a combined bulk nova rate for the four galaxies of 4(+4.2,-1.4) yr^-1 (2{+3.9,-1.4} yr^-1 without the NGC 205 nova candidate). From the bulk rates, integrated and extinction corrected V-band photometry, and (V-K)_0 colors we derive a luminosity specific nova rate for M32 of 12.0(+14.4,-6.0) yr^-1 [10^10 L_Sun,K]^-1 and for NGC 205 of 29.3(+32.3,-14.7) yr^-1 [10^10 L_Sun,K]^-1 and for the combined 4 galaxies of 14.1(+14.8,-4.9) yr^-1 [10^10 L_Sun,K]^-1 (7.0{+13.7,-4.9} yr^-1[10^10 L_Sun,K]^-1 without the NGC 205 nova candidate). If the higher rate is confirmed by surveys in subsequent seasons, it would imply that either dwarf ellipticals have a higher interacting binary fraction than their higher mass counter parts, or that the completeness is higher for these less complex systems and the nova rates for larger, more distant systems are systematically underestimated.

The LSPM-north catalog is a comprehensive list of 61,977 stars north of the J2000 celestial equator that have proper motions larger than 0.15 seconds of arc per year (local-background-stars frame). The catalog has been generated primarily as a result of our systematic search for high proper motion stars in the Digitized Sky Surveys using our SUPERBLINK software. At brighter magnitudes, the catalog incorporates stars and data from the Tycho-2 Catalogue and also, to a lesser extent, from the All-sky Compiled Catalogue of 2.5 million stars. The LSPM catalog expands considerably over the old Luyten (LHS, NLTT) catalogs, superseding them for northern declinations. Corrections to the local-background-stars proper motions have been calculated, and absolute proper motions in the extragalactic frame are given. The catalog is estimated to be over 99% complete at high Galactic latitudes (|b>15|), and over 90% complete at low Galactic latitudes (|b<15|), down to a magnitude V=19.0, and has a limiting magnitude V=21.0. All the northern stars listed in the LHS and NLTT catalogs have been re-identified, and their positions, proper motions and magnitudes re-evaluated. The catalog also lists a large number of completely new objects, which promise to expand very significantly the census of red dwarfs, subdwarfs, and white dwarfs in the vicinity of the Sun.

We report the results of a survey for novae in and between the galaxies of the Fornax cluster. Our survey provides strong evidence that intracluster novae exist and that they provide a useful, independent measure of the intracluster light in Fornax. We discovered six strong nova candidates in six distinct epochs spanning eleven years from 1993 to 2004. The data were taken with the 4m and the 1.5m telescopes at CTIO. The spatial distribution of the nova candidates is consistent with $\sim$16-41% of the total light in the cluster being in the intracluster light, based on the ratio of the number of novae we discovered in intracluster space over the total number of novae discovered plus a simple completeness correction factor. This estimate is consistent with independent measures of intracluster light in Fornax and Virgo using intracluster planetary nebulae. The accuracy of the intracluster light measurement improves with each survey epoch as more novae are discovered.

Despite decades-old predictions of the expected presence of dozens of cataclysmic variables in the cores of globular clusters, the number of irrefutable, out-bursting candidates is still barely a handful. Using multi-wavelength, multi-epoch HST images we have produced outburst and quiescence light curves for the recently discovered large amplitude variable CV1 in the core of the post core-collapse globular cluster M15. The light curves and blue colors show that the object is a bona fide dwarf nova, with absolute magnitude at maximum light rivaling that of the most luminous known dwarf novae.

We present an analysis of optical HST/STIS and HST/FOS spectroscopy of 6 blue stragglers found in the globular clusters M3, NGC6752 and NGC6397. These stars are a subsample of a set of ~50 blue stragglers and stars above the main sequence turn-off in four globular clusters which will be presented in an forthcoming paper. All but the 6 stars presented here can be well fitted with non-LTE model atmospheres. The 6 misfits, on the other hand, possess Balmer jumps which are too large for the effective temperatures implied by their Paschen continua. We find that our data for these stars are consistent with models only if we account for extra absorption of stellar Balmer photons by an ionized circumstellar disk. Column densities of HI and CaII are derived as are the the disks' thicknesses. This is the first time that a circumstellar disk is detected around blue stragglers. The presence of magnetically-locked disks attached to the stars has been suggested as a mechanism to lose the large angular momentum imparted by the collision event at the birth of these stars. The disks implied by our study might not be massive enough to constitute such an angular momentum sink, but they could be the leftovers of once larger disks.

Far-ultraviolet (FUV) observations with the Hubble Space Telescope are an excellent way to find and study the hot, blue stellar populations in the cores of globular clusters. These populations include dynamically-formed blue stragglers and interacting binaries (such as cataclysmic variables), i.e. the products of stellar collisions and near misses. Using the cluster 47 Tuc as an example, we show how the combination of FUV imaging and slitless spectroscopy can be used to uncover and study these populations.

Only one certain classical nova eruption has ever been detected inside a globular cluster - nova 1860 A.D. (T Sco) in M80. During a survey of M87 we have detected an erupting star coincident (to within 0.08 pixels) with a globular cluster of that giant elliptical galaxy. We are able to discount variables in the foreground or background of M87. The light curve and color of the erupting star match those expected for a nova at the distance of M87. The chance superposition of an M87 field nova on the globular cluster is very unlikely but cannot be completely ruled out.Our detection hints at a globular cluster nova frequency $f \sim .004$ novae/cluster/year, much higher than previous observations have suggested.

We respond to the paper by De Marchi, Paresce, Straniero and Moroni (2003) on the white dwarf cooling age of M4. The authors question the data analysis and interpretation that led to the conclusions in Hansen et al. (2002). In their paper, De Marchi et al. are unable to obtain photometry as deep as ours from the same data set and therefore assert that only a lower limit to the white dwarf cooling age for this cluster of approximately 9 Gyr can be obtained. In this short contribution we show that shortcomings in the data analysis and reduction techniques of De Marchi et al. are responsible for their inability to reach the photometry limits that our study reports. In a forthcoming paper in which the complete techniques for age determination with white dwarfs are laid out, we demonstrate that their method of fitting the luminosity function gives a spuriously low white dwarf cooling age.

We investigate in detail the white dwarf cooling sequence of the globular cluster Messier 4. In particular we study the influence of various systematic uncertainties, both observational and theoretical, on the determination of the cluster age from the white dwarf cooling sequence. These include uncertainties in the distance to the cluster and the extinction along the line of sight, as well as the white dwarf mass, envelope and core compositions and the white dwarf --main sequence mass relation. We find that fitting to the full two-dimensional colour-magnitude diagram offers a more robust method for age determination than the traditional method of fitting the one-dimensional white dwarf luminosity function. After taking into account the various uncertainties, we find a best fit age of 12.1 Gyr, with a 95% lower limit of 10.3 Gyr. We also perform fits using two other sets of cooling models from the literature. The models of Chabrier et al (2000) yield an encouragingly similar result, although the models of Salaris et al (2000) do not provide as good a fit. Our results support our previous determination of a delay between the formation of the Galactic halo and the onset of star formation in the Galactic disk.

This paper presents new results from a photometric study of the main-sequence stars in M4, the closest globular cluster to the sun. Multi-field, multi-epoch observations at approximately 1, 2, and 6 core-radii were obtained with the WFPC2 camera on the HST. The multi-epoch observations allowed us to clean the data on the basis of proper motion. In all the fields the cluster main sequence can be traced to at least V = 27.0 but there remains a trail of stars to the limit of the data near V = 30 in the deepest outer field. There is no evidence that we have reached the end of the hydrogen burning main sequence in any of our fields, however, there is some indication that very few stars remain to be detected in the deepest data. The scatter about the cluster main sequence yields a surprisingly small and variable binary fraction; f_b = 2% in the inner parts of the cluster falling to the 1% range outside. For the currently visible main sequence stars, the cluster mass function is very flat (alpha =0.1) in the outer field and flattens further in the inner fields suggesting well developed mass segregation. The observed variation in the mass function is broadly consistent with isotropic, multi-mass Michie-King models. Because we have a large sample of white dwarfs in the outer field we are able to show that the cluster IMF above 0.8Msun was considerably steeper than the present day mass function for low mass stars. Two appendicies are included in this contribution; a detailed discussion of the techniques used to reduce the data and a direct comparison between the cluster stars and those belonging to the inner spheroid of the Galaxy yielding a Galactocentric distance of 7.5+/-0.6 kpc.

The position of the giant elliptical galaxy M87 at the center of the Virgo Cluster means that the inferred column density of dark matter associated with both the cluster halo and the galaxy halo is quite large. This system is thus an important laboratory for studying massive dark objects in elliptical galaxies and galaxy clusters by gravitational microlensing, strongly complementing the studies of spiral galaxy halos performed in the Local Group. We have performed a microlensing survey of M87 with the WFPC2 instrument on the Hubble Space Telescope. Over a period of thirty days, with images taken once daily, we discover seven variable sources. Four are variable stars of some sort, two are consistent with classical novae, and one exhibits an excellent microlensing lightcurve, though with a very blue color implying the somewhat disfavored possibility of a horizontal branch source being lensed. Based on sensitivity calculations from artificial stars and from artificial lightcurves, we estimate the expected microlensing rate. We find that the detection of one event is consistent with a dark halo with a 20% contribution of microlensing objects for both M87 and the Virgo Cluster, similar to the value found from observations in the Local Group. Further work is required to test the hypothesized microlensing component to the cluster.

We present the results of a preliminary H-alpha survey of M81 for novae conducted over a 5 month interval using the 5' field of view camera (WFCAM) on the Calypso Telescope at Kitt Peak, AZ. We observed M81 nearly every clear night during this interval, covering the entire galaxy, and discovered 12 novae. Our comprehensive time coverage allowed us to produce the most complete set of H-alpha light curves for novae in M81 to date. A raw nova rate for M81 gives 23 yr^-1 which, because of the nature of our survey, is a hard lower limit. An analysis of the completeness in our survey gives a corrected nova rate of 30 yr^-1. This agrees well with the rate of 33 (+13,-8) yr^-1, derived from Monte Carlo simulations using nova light curves and survey frame limits. The spatial distribution of the novae we discovered follows the bulge light much better than the disk or total light according to Kolmogorov - Smirnov tests of their radial distributions. The asymmetry in the distribution of novae across the major axis line of M81 implies a bulge-to-disk nova ratio of > 9 and supports the idea that novae originate primarily in older stellar populations.

We present time-resolved, far-ultraviolet (FUV) spectroscopy and photometry of the 1.1 day eclipsing binary system AKO 9 in the globular cluster 47 Tucanae. The FUV spectrum of AKO 9 is blue and exhibits prominent C IV and He II emission lines. The spectrum broadly resembles that of long-period, cataclysmic variables in the galactic field. Combining our time-resolved FUV data with archival optical photometry of 47 Tuc, we refine the orbital period of AKO 9 and define an accurate ephemeris for the system. We also place constraints on several other system parameters, using a variety of observational constraints. We find that all of the empirical evidence is consistent with AKO 9 being a long-period dwarf nova in which mass transfer is driven by the nuclear expansion of a sub-giant donor star. We therefore conclude that AKO 9 is the first spectroscopically confirmed cataclysmic variable in 47 Tuc. We also briefly consider AKO 9's likely formation and ultimate evolution. Regarding the former, we find that the system was almost certainly formed dynamically, either via tidal capture or in a 3-body encounter. Regarding the latter, we show that AKO 9 will probably end its CV phase by becoming a detached, double WD system or by exploding in a Type Ia supernova.

We establish an extragalactic, zero-motion frame of reference within the deepest optical image of a globular star cluster, an HST 123-orbit exposure of M4 (GO 8679, cycle 9). The line of sight beyond M4 (l,b (deg) = 351,16) intersects the inner halo (spheroid) of our Galaxy at a tangent-point distance of 7.6 kpc (for Ro = 8 kpc). We isolate these spheroid stars from the cluster based on their proper motions over the 6-year baseline between these and previous epoch HST data (GO 5461, cycle 4). Distant background galaxies are also found on the same sight line using image-morphology techniques. This fixed reference frame allows us to independently determine the fundamental Galactic constant, Vo/Ro = 25.3 +/- 2.6 km/s/kpc, thus providing a velocity of the Local Standard of Rest, v = 202.7 +/- 24.7 km/s for Ro = 8.0 +/- 0.5 kpc. Secondly, the galaxies allow a direct measurement of M4's absolute proper motion, mu_total = 22.57 +/- 0.76 mas/yr, in excellent agreement with recent studies. The clear separation of galaxies from stars in these deep data also allow us to search for inner-halo white dwarfs. We model the conventional Galactic contributions of white dwarfs along our line of sight and predict 7.9 (thin disk), 6.3 (thick disk) and 2.2 (spheroid) objects to the limiting magnitude at which we can clearly delineate stars from galaxies (V = 29). An additional 2.5 objects are expected from a 20% white dwarf dark halo consisting of 0.5 Mo objects, 70% of which are of the DA type. After considering the kinematics and morphology of the objects in our data set, we find the number of white dwarfs to be consistent with the predictions for each of the conventional populations. However, we do not find any evidence for dark halo white dwarfs.

Time-series data taken with the Hubble Space Telescope of a field six core radii (~5') from the center of the globular cluster Messier 4, and covering a period of about 10 weeks in early 2001, have been analyzed in search of variable objects. Various criteria were employed to select candidate variable stars. Period searches were performed on the selected candidates using phase dispersion minimization (PDM). The reliability of the PDM search results was tested using synthetic light curves of eclipsing binary stars and sinusoidal light curves of different periods. Results from this analysis showed that there are probably no eclipsing binary stars or periodic variables in our field with periods on the order of a few hours to a few days, down to limiting magnitudes of V~25 and I~24, which is consistent with the absence of contact binaries such as W Ursae Majoris systems. However, one candidate variable star does show an increase in brightness of ~0.1 magnitudes in both bandpasses, which seems to last for a few days. Possible explanations concerning the nature of this object include a binary system with a white dwarf primary and a low-mass main sequence secondary, or a BY Draconis variable star. We are able to set an upper limit to the observed fraction of photometric variability in this dataset of 0.05%.

We report the first likely detections of erupting Dwarf Novae (DN) in an external galaxy: the Large Magellanic Cloud. Six candidates were isolated from approximately a million stars observed every second night over 11 nights with the CTIO 8K x 8K Mosaic2 CCD imager. Artificial dwarf nova and completeness tests suggest that we are seeing only the brightest of the LMC DN, probably SS Cygni-like CVs, but possibly SU UMa-type cataclysmics undergoing superoutbursts. We derive crude but useful limits on the LMC DN surface density, and on the number of DN in the LMC. Many thousands of cataclysmic variables in the Magellanic Clouds can be discovered and characterized with 8 meter class telescopes.

In a continuation of our systematic search for high proper motion stars in the Digitized Sky Survey, we have completed the analysis of northern sky fields at galactic latitudes above 25 degrees. With the help of our SUPERBLINK software, a powerful automated blink comparator developed by us, we have identified 1146 stars in the magnitude range 8<r<20 with proper motions 0.500<mu<2.000 arcsec/yr. These include 1080 stars previously listed in Luyten's proper motion catalogs (LHS, NLTT), 9 stars not previously listed in the Luyten catalogs but reported elsewhere in the literature (including 1 previously reported by our team), and 57 new objects reported here for the first time. This paper includes a list of positions, proper motions, magnitudes, and finder charts for all the new high proper motion stars. Combined with our previous study of low galactic latitude fields (see Paper I), our survey now covers over 98% of the northern sky. We conclude that the Luyten catalogs were 90% complete in the northern sky for stars with 0.5<mu<2.0 arcsec/yr down to magnitude r=19. We discuss the incompleteness of the old Luyten proper motion survey, and estimate completeness limits for our new survey.

We report the discovery of LSR1610-0040, a previously unreported r=17.5 star with a very high proper motion mu=1.46"/yr. This very cool star (b-i=6.3) is found to have a peculiar spectrum that does not fit into the standard sequence of late-M dwarfs and L dwarfs. Rather, the spectrum is more typical of an ultra-cool subdwarf, with weak bands of TiO and no detectable VO. But because LSR1610-0040 is so much redder than any other sdM known, and because it does not appear to fit well into the sdM sequence, we propose that LSR1610-0040 be assigned a spectral type sdL, making it the first early-type L subdwarf known. Evidence suggests that LSR1610-0040 is an old, metal-poor star whose mass is just above the limit of hydrogen burning (M ~ 0.08 M_sun).

We report the discovery of LSR0602+3910, an L dwarf of class L1. The object was initially identified by Lepine et al. (2002) as a new high proper motion star lying close to the Galactic plane. Its 2MASS J-K_s=1.43 is consistent with an L dwarf, which we now confirm spectroscopically. In addition, we see a signature of Li I absorption, making LSR0602+3910 a brown dwarf, one of the brightest known (K_s=10.86). Among L dwarfs it is second in brightness to the combined light of 2MASS 0746+20, a close binary system. We see no indication that LSR0602+3910 is a binary, although high-resolution imaging will be required to confirm this. Spectroscopic and photometric distance estimates agree very well, placing LSR0602+3910 at d=10.6+-0.8 pc. LSR0602+3910 was most likely missed in previous searches because of its proximity to the plane, the region that most searches avoided. We estimate that some 40% of bright L dwarfs are missed because of this selection effect.

We use the results of realistic N-body simulations to investigate the appearance of the white dwarf population in dense star clusters. We show that the presence of a substantial binary population in a star cluster, and the interaction of this population with the cluster environment, has serious consequences for the morphology of the observed white dwarf sequence and the derived white dwarf cooling age of the cluster. We find that over time the dynamical evolution of the cluster -- mass-segregation, stellar interactions and tidal stripping -- hampers the use of white dwarfs as tracers of the initial mass function, and also leads to a significant enhancement of the white dwarf mass fraction. Future observations of star clusters should be conducted slightly interior to the half-mass radius of the cluster in order to best obtain information about the cluster age and initial mass function from the white dwarf luminosity function. The evolution of binary stars and the cluster environment must necessarily be accounted for when studying the white dwarf populations of dynamically evolved star clusters.

We report the discovery of the coolest subdwarf reported to date. The star LSR1425+7102 was discovered in our survey for faint high proper motion stars in the northern sky. Follow-up spectroscopy revealed the star to be a very red object with the characteristic signature of M subdwarfs: strong CaH bands but relatively weaker TiO bands. The CaH molecular band at $\approx\lambda$6900\AA\ is particularly strong. By extrapolating the empirical relationship between the strength of the CaH molecular band and the subdwarf subtype, we conclude that LSR1425+7102 is the first star to be discovered with spectral type sdM8.0.

We define an enhanced spectral classification scheme for M dwarf stars, and use it to derive spectral classification of 104 northern stars with proper motions larger than 0.5"/yr which we discovered in a survey of high proper motion stars at low galactic latitudes. The final tally is as follows: 54 M dwarfs, 25 sdK and sdM subdwarfs, 14 esdK and esdM extreme subdwarfs, and 11 DA and DC white dwarfs. Among the most interesting cases, we find one star to be the coolest subdwarf ever reported (LSR2036+5059, with spectral type sdM7.5), a new M9.0 dwarf only about 6pc distant (LSR1835+3259), and a new M6.5 dwarf only 7pc from the Sun (LSR2124+4003). Spectroscopic distances suggests that 27 of the M dwarfs, 3 of the white dwarfs, and one of the subdwarfs (LSR2036+5059) are within 25pc of the Sun, making them excellent candidates for inclusion in the solar neighborhood census. Estimated sky-projected velocities suggest that most of our subdwarfs and extreme subdwarfs have halo kinematics. We find that several white dwarfs and non metal-poor M dwarfs also have kinematics consistent with the halo, and we briefly discuss their possible origin.

We report the discovery of LSR1826+3014, a very faint (V=19.36) star with a very large proper motion (mu=2.38 arcsec/yr). A low resolution red spectrum reveals that LSR1826+3014 is an ultra-cool red dwarf with spectral type M8.5 V and with a radial velocity v_rad=+77+/-10 km/s. LSR1826+3014 is thus the faintest red dwarf ever discovered with a proper motion larger than 2 arcsec/yr. Optical and infrared photometry suggest that the star is at a distance d=13.9+/-3.5 pc from the Sun, which implies it is moving relative to the local standard of rest with a total velocity of 175+/-25 km/s. Numerical integration of its orbit suggests that LSR1826+3014 is on a halo-like galactic orbit.

We summarize the main results from MODEST-1, the first workshop on MOdeling DEnse STellar systems. Our goal is to go beyond traditional population synthesis models, by introducing dynamical interactions between single stars, binaries, and multiple systems. The challenge is to define and develop a software framework to enable us to combine in one simulation existing computer codes in stellar evolution, stellar dynamics, and stellar hydrodynamics. With this objective, the workshop brought together experts in these three fields, as well as other interested astrophysicists and computer scientists. We report here our main conclusions, questions and suggestions for further steps toward integrating stellar evolution and stellar (hydro)dynamics.

We have conducted a systematic search for high proper motion stars in the Digitized Sky Survey, in the area of the sky north of -2.8 degrees in declination and within 25 degrees of the galactic plane. Using the SUPERBLINK software, a powerful automated blink comparator developed by us, we have identified 601 stars in the magnitude range 9<r<20 with proper motions in the range 0.5 to 2.0 arcsec/yr which have not been observed with the HIPPARCOS satellite. Among those, we recovered 460 stars previously listed in Luyten's proper motion catalogs (LHS, NLTT), and discovered 141 new high proper motion stars, ranging in magnitude from R=13.0 to R=19.8. Only 9 stars from the Luyten catalogs that were not observed by HIPPARCOS could not be recovered with SUPERBLINK: 3 had proper motions larger than the search limit of the code (>2.0\arcsec/yr), and 5 were missed because they were either too bright for SUPERBLINK to handle or they are in the immediate proximity of very bright stars. Only one of Luyten's stars (LHS1657) could not be recovered at all, even by visual inspection of the POSS plates, and is now suspected to be bogus. The very high success rate in the recovery by SUPERBLINK of faint Luyten stars suggests that our new survey of high proper motion stars is at least 99 percent complete for stars with proper motions 0.5 to 2.0 arcsec/yr down to R=19. This paper includes a list of positions, proper motions, magnitudes, and finder charts for all the new high proper motion stars.

The deepest optical image ever in a globular star cluster, a Hubble Space Telescope (HST) 123 orbit exposure in a single field of Messier 4, was obtained in 2 filters (F606W, F814W) over a 10 week period in early 2001. A somewhat shallower image obtained in 1995 allowed us to select out cluster and field objects via their proper motion displacement resulting in remarkably clean color-magnitude diagrams that reach to V = 30, I = 28. The cluster main sequence luminosity function contains very few stars fainter than MV = 15.0, MI = 11.8 which, in both filters, is more than 2 magnitudes brighter than our limit. This is about the faintest luminosity seen among field Population II subdwarfs of the same metallicity. However, there remains a sprinkling of potential cluster stars to lower luminosity all the way down to our limiting magnitudes. These latter objects are significantly redder than any known metal-poor field subdwarf. Comparison with the current generation of theoretical stellar models implies that the masses of the lowest luminosity cluster stars observed are near 0.09 Msun. We derive the mass function of the cluster in our field and find that it is very slowly rising towards the lowest masses with no convincing evidence of a turnover even below 0.1 Msun. The formal slope between 0.65 and 0.09 Msun is alpha = 0.75 (Salpeter of 2.35) with a 99% confidence interval 0.55 - 1.05. A consistency check between these slopes and the number of observed cluster white dwarfs yields a range of possible conclusions, one of which is that we have indeed seen the termination of the white dwarf cooling sequence in M4.

We present results from the imaging portion of a far-ultraviolet (FUV) survey of the core of 47 Tucanae. We have detected 767 FUV sources, 527 of which have optical counterparts in archival HST/WFPC2 images of the same field. Most of our FUV sources are main-sequence (MS) turn-off stars near the detection limit of our survey. However, the FUV/optical color-magnitude diagram (CMD) also reveals 19 blue stragglers (BSs), 17 white dwarfs (WDs) and 16 cataclysmic variable (CV) candidates. The BSs lie on the extended cluster MS, and four of them are variable in the FUV data. The WDs occupy the top of the cluster cooling sequence, down to an effective temperature of T_{eff} \simeq 20,000 K. Our FUV source catalog probably contains many additional, cooler WDs without optical counterparts. Finally, the CV candidates are objects between the WD cooling track and the extended cluster MS. Four of the CV candidates are previously known or suspected cataclysmics. All of these are bright and variable in the FUV. Another CV candidate is associated with the semi-detached binary system V36 that was recently found by Albrow et al. (2001). V36 has an orbital period of 0.4 or 0.8 days, blue optical colors and is located within 1 arcsec of a Chandra x-ray source. A few of the remaining CV candidates may represent chance superpositions or SMC interlopers, but at least half are expected to be real cluster members with peculiar colors. However, only a few of these CV candidates are possible counterparts to Chandra x-ray sources. Thus it is not yet clear which, if any, of them are true CVs, rather than non-interacting MS/WD binaries or Helium WDs.

We present the white dwarf sequence of the globular cluster M4, based on a 123 orbit Hubble Space Telescope exposure, with limiting magnitude V = 30, I = 28. The white dwarf luminosity function rises sharply for I >25.5, consistent with the behaviour expected for a burst population. The white dwarfs of M4 extend to approximately 2.5 magnitudes fainter than the peak of the local Galactic disk white dwarf luminosity function. This demonstrates a clear and significant age difference between the Galactic disk and the halo globular cluster M4. Using the same standard white dwarf models (Hansen 1999) to fit each luminosity function yields ages of 7.3 +/- 1.5 Gyr for the disk and 12.7 +/- 0.7 Gyr for M4 (2-sigma statistical errors).

In an ongoing survey for high proper motion stars at low galactic latitudes, we have identified 5 previously uncatalogued common proper motion companions to stars listed in the Preliminary Version of the Third Catalogue of Nearby Stars (GJ stars). For each system, the relative proper motion between the components is less than 5% of the common proper motion of the pair. Spectra of the pairs have been obtained at the Lick 3-m Shane Telescope, confirming that the systems are indeed wide separation binaries. The systems are classified as follows: GJ 4047AB = K5 V + M5 V, GJ 718AB = K5 V + M4.5 V, GJ 747.2ABC = (K7 V + K7 V) + M4 Ve, GJ 4100AB = M1 V + M4.5 Ve, and GJ 4153AB = M0.5 V + M3.5 V. The total area surveyed contains 346 Gliese stars, which suggests that approximately 1.5% of the stars listed in the Preliminary Version of the Third Catalogue of Nearby Stars have unrecognized proper motion companions. We predict that 15-50 more new distant companions to GJ stars will be discovered in the Digitized Sky Survey.

We report the first orbital period determination for a Dwarf Nova (DN) in a glubular cluster: V101 in M5 has a period of 5.796 +- 0.036 hours. We derived this period from I-band photometry acquired with the Calypso Observatory High Resolution Camera operating with tip-tilt adaptive optics correction. Observations from the South African Astronomical Observatory in the V-band were also analyzed and exhibit a periodic signal of the same period. This orbital period suggests that V101 has a secondary of mid to late K spectral type with Mv = +8.2 +- 0.5. The predicted spectral type is consistent with previous spectral observations in quiescence which show a fairly red continuum. From the observed minimum brightness of V = 22.5, we derive a distance modulus of (m - M)v = 14.3 +- 0.5 to the DN which supports V101's membership in the globular cluster M5. Measurement of the ellipsoidality effect indicates that the orbital plane of the V101 system is moderately inclined, but not enough to exhibit eclipses.

We find a remarkably enhanced production rate in star clusters (relative to the field) of very short period, massive double-white-dwarf stars and of giant-white dwarf binaries. These results are based on N-body simulations performed with the new GRAPE-6 special purpose hardware and are important in identifying and characterizing the progenitors of type Ia supernovae. The high incidence of very close double-white-dwarf systems is the result of dynamical encounters between (mostly) primordial binaries and other cluster stars. Orbital hardening rapidly drives these degenerate binaries to periods under ~10 hours. Gravitational radiation emission and mergers producing supra-Chandrasekhar objects follow in less than a Hubble time. If most stars are born in clusters then estimates of the double white dwarf merger rates in galaxies (due to cluster dynamical interaction) must be increased more than tenfold. A majority of the Roche lobe overflow giant-white dwarf binaries are not primordial; they are produced in exchange reactions. Most cases resulted in a common-envelope and formation of a double-white-dwarf binary rather than Supersoft X-ray sources leading possibly to a type Ia supernova.

The recent availability of special purpose computers designed for calculating gravitational interactions of N-bodies at extremely high speed has provided the means to model globular clusters on a star-by-star basis for the first time. By endeavouring to make the N-body codes that operate on these machines as realistic as possible, the addition of stellar evolution being one example, much is being learnt about the interaction between the star cluster itself and the stars it contains. A fascinating aspect of this research is the ability to follow the orbits of individual stars in detail and to document the formation of observed exotic systems. This has revealed that many stars within a star cluster lead wildly promiscuous lives, interacting, often intimately and in rapid succession, with a variety of neighbours.

We present spatially resolved spectra of the visual WR+OB massive binaries WR86, WR146, and WR147, obtained with the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope. The systems are classified as follows: WR86 = WC7 + B0 III, WR146 = WC6 + O8 I-IIf, WR147 = WN8 + O5-7 I-II(f). Both WR146 and WR147 are known to have strong non-thermal radio emission arising in a wind-wind collision shock zone between the WR and OB components. We find that the spectra of their O companions show Halpha profiles in emission, indicative of large mass-loss rates, and consistent with the colliding-wind model. Our spectra indicate that the B component in WR86 has a low mass-loss rate, which possibly explains the fact that WR86, despite being a long period WR+OB binary, was not found to be a strong non-thermal radio emitter. Because of the small mass-loss rate of the B star component in WR86, the wind collision region must be closer to the B star and smaller in effective area, hence generating smaller amounts of non-thermal radio emission. Absolute magnitudes for all the stars are estimated based on the spectral types of the components (based on the tables by Schmidt-Kaler for OB stars, and van der Hucht for WR stars), and compared with actual, observed magnitude differences. While the derived luminosities for the WC7 and B0 III stars in WR86 are consistent with the observed magnitude difference, we find a discrepancy of at least 1.5 magnitudes between the observed luminosities of the components in each of WR146 and WR147 and the absolute magnitudes expected from their spectral types. In both cases, it looks as though either the WR components are about 2 magnitudes too bright for their spectral types, or that the O components are about 2 magnitudes too faint. We discuss possible explanations for this apparent discrepancy.

We find that free-floating planets can remain bound to a star cluster for much longer than was previously assumed: of the order of the cluster half-mass relaxation timescale as opposed to the crossing-time. This result is based on N-body simulations performed with the new GRAPE-6 special purpose hardware and is important in the context of the preliminary detection of a population of free-floating sub-stellar objects in the globular cluster M22. The planets in our N-body study are of Jupiter mass and are initially placed in circular orbits of between 0.05 and 50 AU about a parent star whose mass is chosen from a realistic initial mass function. The presence of the free-floating planets is the result of dynamical encounters between planetary systems and the cluster stars. Most planets are liberated from their parent star in, or near, the cluster core, and then drift outwards on a timescale of ~ 10^{8-9} year. This still implies the existence of many (~100) planets per star if the M22 result is confirmed.

We report new photometry for the cluster NGC 2121 in the Large Magellanic Cloud, which shows a prominent hydrogen core exhaustion gap at the turnoff, and a descending subgiant branch remniscent of Galactic open clusters. We achieve an excellent fit using the Girardi isochrones, finding an age of 3.2 +/- 0.5Gyr, with [Fe/H]=-0.6 +/- 0.2. The isochrones fit the color and shape of the turnoff and subgiant branch so precisely that we constrain the metallicity as well as the age. The smae isochrones also fit SL 663 and NGC 2155, although our photometry for these clusters has much larger errors. We find these clusters to be 0.8 Gyr younger, and 0.4 dex more metal rich, than recently reported in the literature. Consequently, we argue that NGC 2121, NGC 2155, and SL 663 are not properly assigned to the age gap in the LMC, but instead are among the first clusters to have formed int he relatively metal rich, younger group of LMC clusters. We propose a new definition of the LMC Age Gap as extending from 3.2 to 13 Gyr, with ESO121-SCO3 still the only remaining candidate for membership in the age gap.

N-body simulations of star cluster evolution have reached a high degree of realism, by incorporating more and more elements of stellar dynamics, stellar evolution, and hydrodynamics. At the end of this conference, six participants discussed how to present the increasingly realistic data from star cluster simulations in a way that is most useful for a direct comparison with observations.

The Catalog and Atlas of Cataclysmic Variables (Edition 1 - 1993) and Edition 2 - 1997) has been a valuable source of information for the cataclysmic variable (CV) community. However, the goal of having a central location for all objects is slowly being lost as each new edition is generated. There can also be a long time delay between new information becoming available on an object and its publication in the catalog. To eliminate these concerns, as well as to make the catalog more accessible, we have created a web site which will contain a ``living'' edition of the catalog. We have also added orbital period information, as well as finding charts for novae, to the catalog.

In dense clusters a bewildering variety of interactions between stars can be observed, ranging from simple encounters to collisions and other mass-transfer encounters. With faster and special-purpose computers like GRAPE, the amount of data per simulation is now exceeding 1TB. Visualization of such data has now become a complex 4D data-mining problem, combining space and time, and finding interesting events in these large datasets. We have recently starting using the virtual reality simulator, installed in the Hayden Planetarium in the American Museum for Natural History, to tackle some of these problem. This work (http://www.astro.umd.edu/nemo/amnh/) reports on our first ``observations'', modifications needed for our specific experiments, and perhaps field ideas for other fields in science which can benefit from such immersion. We also discuss how our normal analysis programs can be interfaced with this kind of visualization.

We present preliminary results from a FUV spectroscopic and photometric survey of the globular cluster 47 Tuc. Our main goal is to either confirm or rule out the existence of a large population of CVs in the core of this cluster. We have so far identified approximately 425 FUV sources, most of which are probably hot, young WDs. However, we have also found approximately 30 FUV sources whose position in a FUV-optical color-magnitude diagram makes them strong CV candidates. If most or all of these objects are eventually confirmed as CVs, the large CV population predicted by tidal capture theory has finally been found. Our data has also allowed us to resolve the long-standing puzzle surrounding AKO 9, a UV-bright and highly variable 1.1 day binary system in 47 Tuc. AKO 9 is the brightest FUV source in our data and presents a blue FUV spectrum with strong CIV and HeII emission lines. Its spectrum is similar to that of the long-period, DN-type field CV GK Per. Based on this similarity and other evidence, we suggest that AKO 9 is a long-period CV in which mass-transfer is driven by the nuclear evolution of a sub-giant donor star.

We report new photometry of populous intermediate-age clusters in the SMC using the {\sl Hubble Space Telescope}. In contrast to the accepted picture, these clusters appear to have formed in two brief intervals, one $8\pm 2$ Gyr ago, and a more recent burst $2\pm 0.5$ Gyr ago. When the ridgelines of the four clusters (NGC 339, 361, 416, and Kron 3) in the 8 Gyr burst are aligned, the dispersion in turnoff luminosities is $<0.2$ mag, corresponding to a maximum age spread of $\pm 0.7$ Gyr. When the ridgelines of three clusters (NGC 152, 411, and 419) in the 2 Gyr burst are aligned, the maximum dispersion of 0.2 mag in turnoff luminosity corresponds to a permitted age spread of $\pm 0.2$ Gyr. Within each group of clusters, the entire cluster loci (including red giant branches and clumps) are nearly identical, consistent with indistinguishable metallicities and ages. In contrast to the wide dispersion in ages previously reported in the literature, our sample with more precise photometry and age measurements supports a burst-punctuated rather than a continuous cluster formation history for the 2 Gyr and 8 Gyr SMC clusters.

Data obtained in the 1950-1955 Palomar campaign for the discovery of classical novae in M81 are set out in detail. Positions and apparent B magnitudes are listed for the 23 novae that were found. There is modest evidence that the spatial distribution of the novae does not track the B brightness distribution of either the total light or the light beyond an isophotal radius that is $70\arcsec$ from the center of M81. The nova distribution is more extended than the aforementioned light, with a significant fraction of the sample appearing in the outer disk/spiral arm region. We suggest that many (perhaps a majority) of the M81 novae that are observed at any given epoch (compared with say $10^{10}$ years ago) are daughters of Population I interacting binaries. The conclusion that the present day novae are drawn from two population groups, one from low mass white dwarf secondaries of close binaries identified with the bulge/thick disk population, and the other from massive white dwarf secondaries identified with the outer thin disk/spiral arm population, is discussed. We conclude that the M81 data are consistent with the two population division as argued previously from (1) the observational studies on other grounds by Della Valle et al. (1992, 1994), Della Valle & Livio (1998), and Shafter et al. (1996) of nearby galaxies, (2) the Hatano et al. (1997a,b) Monte Carlo simulations of novae in M31 and in the Galaxy, and (3) the Yungelson et al. (1997) population synthesis modeling of nova binaries. Two different methods of using M81 novae as distance indicators give a nova distance modulus for M81 as $(m-M)_0 = 27.75$, consistent with the Cepheid modulus that is the same value.

We present an analysis of the deepest pure-UV observations with the highest angular resolution ever performed, a set of 12 exposures with the HST WFPC2 and F160BW filter obtained in parallel observing mode, which cover $\sim$12 square arcminutes in the LMC, North of the bar and in the ``general field'' regime of the LMC. The 341 independent measurements of 198 objects represent an accumulated exposure of $\geq 2 10^4$ sec and reveal stars as faint as m$_{UV}\simeq$22 mag. The observations show that $\sim 2/3$ of the UV emission from the LMC is emitted by our HST-detected UV stars in the field, i.e., not in clusters or associations. We identified optical counterparts in the ROE/NRL photometric catalog for $\sim 1/3$ of the objects. The results are used to discuss the nature of these UV sources, to estimate the diffuse UV emission from the LMC as a prototype of dwarf galaxies, and to evaluate the contamination by field stars of UV observations of globular and open clusters in the LMC. We find that the projected density of UV stars in the general field of the LMC is a few times higher than in the Galactic disk close to the Sun. Combining our data with observations by UIT allows us to define the stellar UV luminosity function from $m_{UV}=8$ to 18 mag, and to confirm that the field regions in the LMC have been forming stars at a steady rate during the last 1 Gyr, with an IMF close to the Salpeter law.

We report the discovery of a variable blue straggler in the core of the globular cluster M5, based on a 12-hour long series of images obtained with the Planetary Camera aboard the Hubble Space Telescope. In addition, we present the light curves of 28 previously unknown or poorly studied large-amplitude variable stars (all but one are RR Lyrae) in the cluster core. A (V, U-I) color-magnitude diagram shows 24 blue stragglers within 2 core radii of the cluster center. The blue straggler population is significantly more centrally concentrated than the horizontal branch and red giant stars.

We present our analysis of Hubble Space Telescope Wide Field Planetary Camera 2 observations in F555W (broadband V) and F450W (broadband B) of the globular cluster Hodge 11 in the Large Magellanic Cloud galaxy. The resulting V vs. (B-V) color-magnitude diagram reaches 2.4 mag below the main-sequence turnoff (which is at V_TO = 22.65 +- 0.10 mag or M_V^TO = 4.00 +- 0.16 mag). Comparing the fiducial sequence of Hodge 11 with that of the Galactic globular cluster M92, we conclude that, within the accuracy of our photometry, the age of Hodge 11 is identical to that of M92 with a relative age-difference uncertainty ranging from 10% to 21%. Provided that Hodge 11 has always been a part of the Large Magellanic Cloud and was not stripped from the halo of the Milky Way or absorbed from a cannibalized dwarf spheroidal galaxy, then the oldest stars in the Large Magellanic Clouds and the Milky Way appear to have the same age.