Source author record

Michael M. Shara

Michael M. Shara appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

49works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

49 published item(s)

preprint2020arXiv

A unified theory of cataclysmic variables from self-consistent numerical simulations

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.

preprint2020arXiv

Hot Jupiter and ultra-cold Saturn formation in dense star clusters

The discovery of high incidence of hot Jupiters in dense clusters challenges the field-based hot Jupiter formation theory. In dense clusters, interactions between planetary systems and flyby stars are relatively common. This has a significant impact on planetary systems, dominating hot Jupiter formation. In this paper, we perform high precision, few-body simulations of stellar flybys and subsequent planet migration in clusters. A large parameter space exploration demonstrates that close flybys that change the architecture of the planetary system can activate high eccentricity migration mechanisms: Lidov-Kozai and planet-planet scattering, leading to high hot Jupiter formation rate in dense clusters. Our simulations predict that many of the hot Jupiters are accompanied by "ultra-cold Saturns", expelled to apastra of thousands of AU. This increase is particularly remarkable for planetary systems originally hosting two giant planets with semi-major axis ratios $\sim$ 4 and the flyby star approaching nearly perpendicular to the planetary orbital plane. The estimated lower limit to the hot Jupiter formation rate of a virialized cluster is $\sim 1.6\times10^{-4}(σ/{\rm 1kms^{-1}})^5({a_{\rm p}}/{\rm 20 AU})({M_{\rm c}}/{\rm 1000M_\odot})^{-2}$Gyr$^{-1}$ per star, where $σ$ is the cluster velocity dispersion, $a_{\rm p}$ is the size of the planetary system and $M_{\rm c}$ is the mass of the cluster. Our simulations yield a hot Jupiter abundance which is $\sim$ 50 times smaller than that observed in the old open cluster M67. We expect that interactions involving binary stars, as well as a third or more giant planets, will close the discrepancy.

preprint2020arXiv

The spin rates of O stars in WR + O Magellanic Cloud binaries

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.

preprint2019arXiv

Chemical Properties of the Local Galactic Disk and Halo. I. Fundamental Properties of 1,544 Nearby, High Proper-Motion M dwarfs and subdwarfs

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.

preprint2017arXiv

A comparison between SALT/SAAO observations and kilonova models for AT 2017gfo: the first electromagnetic counterpart of a gravitational wave transient - GW170817

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.

preprint2016arXiv

A Hubble Space Telescope Survey for Novae in M87. I. Light and Color Curves, Spatial Distributions and the Nova Rate

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.

preprint2016arXiv

A Survey of the Local Group of Galaxies for Symbiotic Binary Stars. I. First detection of symbiotic stars in M33

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.

preprint2016arXiv

Late-Time Photometry of Type Ia Supernova SN 2012cg Reveals the Radioactive Decay of $^{57}$Co

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.

preprint2016arXiv

On the origins of enigmatic stellar populations in Local Group galactic nuclei

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.

preprint2016arXiv

The chaotic four-body problem in Newtonian gravity I: Identical point-particles

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.

preprint2016arXiv

The spin rates of O stars in WR + O binaries. I. Motivation, methodology and first results from SALT

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.

preprint2016arXiv

When does a star cluster become a multiple star system? I. Lifetimes of equal-mass small-N systems

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.

preprint2016arXiv

When does an old nova become a dwarf nova? Kinematics and age of the nova shell of the dwarf nova AT Cnc

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.

preprint2015arXiv

A Near-Infrared Survey of the Inner Galactic Plane for Wolf-Rayet Stars III. New Methods: Faintest WR Stars

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.

preprint2015arXiv

A Survey for hot Central Stars of Planetary Nebulae I. Methods and First Results

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.

preprint2015arXiv

Characterization of the Most Luminous Star in M33: A Super Symbiotic Binary

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.

preprint2015arXiv

Growing White Dwarfs to the Chandrasekhar Limit: The Parameter Space of the Single Degenerate SNIa Channel

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.

preprint2015arXiv

HST Images Flash Ionization of Old Ejecta by the 2011 Eruption of Recurrent Nova T Pyxidis

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.

preprint2015arXiv

The first transition Wolf-Rayet WN/C star in M31

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.

preprint2014arXiv

Characterizing Wolf-Rayet Stars in the Near and Mid Infrared

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.

preprint2014arXiv

First detection and characterization of symbiotic stars in M31

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.

preprint2014arXiv

Observational Signatures of SNIa Progenitors, as Predicted by Models

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.

preprint2014arXiv

Progenitor constraints on the Type-Ia supernova SN2011fe from pre-explosion Hubble Space Telescope HeII narrow-band observations

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 Å 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ö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 Å line emission in the HST data. Using simulations, we set $2σ$ upper limits on the HeII 4686 Å 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.

preprint2013arXiv

A Dynamical Signature of Multiple Stellar Populations in 47 Tucanae

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.

preprint2013arXiv

Comparing the White Dwarf Cooling Sequences in 47 Tuc and NGC 6397

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.

preprint2013arXiv

Dynamical Fragmentation of the T Pyxidis Nova Shell During Recurrent Eruptions

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.

preprint2013arXiv

The Vast Population of Wolf-Rayet and Red Supergiant Stars in M101: I. Motivation and First Results

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.

preprint2012arXiv

A direct N-body model of core-collapse and core oscillations

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.

preprint2012arXiv

A Near-Infrared Survey of the Inner Galactic Plane for Wolf-Rayet Stars II. Going Fainter: 71 More New WR Stars

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.

preprint2012arXiv

AT Cnc: A Second Dwarf Nova with a Classical Nova Shell

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.

preprint2012arXiv

GK Per (Nova Persei 1901): HST Imagery and Spectroscopy of the Ejecta, and First Spectrum of the Jet-Like Feature

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.

preprint2012arXiv

The Brown Dwarf Kinematics Project (BDKP). III. Parallaxes for 70 Ultracool Dwarfs

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.

preprint2012arXiv

The Inter-Eruption Timescale of Classical Novae from Expansion of the Z Camelopardalis Shell

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.

preprint2011arXiv

A Deep, Wide-Field, and Panchromatic View of 47 Tuc and the SMC with HST: Observations and Data Analysis Methods

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.

preprint2011arXiv

Constraints on the Progenitor System of the Type Ia Supernova SN 2011fe/PTF11kly

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.

preprint2011arXiv

WISEP J180026.60+013453.1: A Nearby Late L Dwarf Near the Galactic Plane

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.

preprint2010arXiv

An Extended Grid of Nova Models. III. Very Luminous, Red Novae

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.

preprint2010arXiv

Non-Equipartition of Energy, Masses of Nova Ejecta, and Type Ia Supernovae

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.

preprint2010arXiv

The Red Nova-like Variable in M31 - A Blue Candidate in Quiescence

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.

preprint2009arXiv

A Near-Infrared Survey of the Inner Galactic Plane for Wolf-Rayet Stars I. Methods and First Results: 41 New WR Stars

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.

preprint2009arXiv

The Brown Dwarf Kinematics Project (BDKP). II. Details on Nine Wide Common Proper Motion Very Low-Mass Companions to Nearby Stars

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).

preprint2008arXiv

The Brown Dwarf Kinematics Project (BDKP) I. Proper Motions and Tangential Velocities for a Large Sample of Late-type M, L and T Dwarfs

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.

preprint2005arXiv

A Possible High Nova Rate for Two Local Group Dwarf Galaxies: M32 and NGC 205

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.

preprint2004arXiv

A Catalog of Northern Stars With Annual Proper Motions Larger Than 0.15 Seconds of Arc (LSPM catalog - North)

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.

preprint2003arXiv

The H-alpha Light Curves and Spatial Distribution of Novae in M81

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.

preprint2002arXiv

New High Proper Motion Stars from the Digitized Sky Survey. I. Northern Stars with proper motion between 0.5 and 2.0 arcsec/yr at Low Galactic Latitudes

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.

preprint2002arXiv

The First Orbital Period for a Dwarf Nova in a Globular Cluster: V101 in M5

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.

preprint2002arXiv

The White Dwarf Cooling Sequence of the Globular Cluster Messier 4

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).

preprint2001arXiv

Free-Floating Planets: Not So Surprising

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.