We present a serendipitous multiwavelength campaign of optical photometry simultaneous with Integral X-ray monitoring of the 2015 outburst of the black hole V404 Cyg. Large amplitude optical variability is generally correlated with X-rays, with lags of order a minute or less compatible with binary light travel timescales or jet ejections. Rapid optical flaring on time-scales of seconds or less is incompatible with binary light-travel timescales and has instead been associated with synchrotron emission from a jet. Both this rapid jet response and the lagged and smeared one can be present simultaneously. The optical brightness is not uniquely determined by the X-ray brightness, but the X-ray/optical relationship is bounded by a lower-envelope such that at any given optical brightness there is a maximum X-ray brightness seen.} This lower-envelope traces out a Fopt proportional to Fx^0.54 relation which can be approximately extrapolated back to quiescence. Rapid optical variability is only seen near this envelope, and these periods correspond to the hardest hard X-ray colours. This correlation between hard X-ray colour and optical variability (and anti-correlation with optical brightness) is a novel finding of this campaign, and apparently a facet of the outburst behaviour in V404 Cyg. It is likely that these correlations are driven by changes in the central accretion rate and geometry.

We use Hubble Space Telescope multicolour photometry of the globular cluster 47 Tucanae to uncover a population of 24 objects with no previous classification that are outliers from the single-star model tracks in the colour-magnitude diagram and yet are likely cluster members. By comparing those sources with evolutionary models and X-ray source catalogues, we were able to show that the majority of those sources are likely binary systems that do not have any X-ray source detected nearby, most possibly formed by a white dwarf and a main-sequence star and a small number of possible double-degenerate systems.

We used the light curve code XRBinary to model the quiescent K2 light curves of three low-inclination cataclysmic variables (CVs): 1RXS\,J0632+2536 (J0632+2536), RZ\,Leo, TW\,Vir and the pre-CV WD\,1144+011. Optimized light curve models were obtained using a nonlinear fitting code NMfit and visualized by Phoebe 2.0. The disk model of J0632+2536 shows that one hotspot at the edge of the disk is enough to describe its light curve, while the other two dwarf nova (DN): RZ\,Leo and TW\,Vir require two hotspots. A typical pre-CV model with a weak irradiation effect for WD\,1144+011 can explain its single-hump modulation, and the newly observed spectrum confirms its previous classification. The synthetic analyses for the DN clearly indicate that phase zero of the double-hump modulations occurs around the secondary minimum and the primary hump is mainly caused by the hotspot at the edge of the disk. The quiescent disk has a flat temperature distribution with a power index of $\sim0.11$. The disk model of RZ\,Leo implies a truncated disk, supporting its previously speculated classification as an intermediate polar (IP). Except for the IP model of RZ\,Leo, which lacks a component related to the inferred accretion curtain, the models of J0632+2536, TW\,Vir and WD\,1144+011 are consistent with results from the Gaia mission. The derived masses and radii of the secondaries of the three DN are consistent with the semi-empirical relations for CV donor stars, while their effective temperatures are higher than the predictions. Irradiation of the donor stars is investigated to explain this discrepancy.

The initial-final mass relation (IFMR) of white dwarfs (WDs) plays an important role in stellar evolution. To derive precise estimates of IFMRs and explore how they may vary among star clusters, we propose a Bayesian hierarchical model that pools photo- metric data from multiple star clusters. After performing a simulation study to show the benefits of the Bayesian hierarchical model, we apply this model to five star clus- ters: the Hyades, M67, NGC 188, NGC 2168, and NGC 2477, leading to reasonable and consistent estimates of IFMRs for these clusters. We illustrate how a cluster-specific analysis of NGC 188 using its own photometric data can produce an unreasonable IFMR since its WDs have a narrow range of zero-age main sequence (ZAMS) masses. However, the Bayesian hierarchical model corrects the cluster-specific analysis by bor- rowing strength from other clusters, thus generating more reliable estimates of IFMR parameters. The data analysis presents the benefits of Bayesian hierarchical modelling over conventional cluster-specific methods, which motivates us to elaborate the pow- erful statistical techniques in this article.

In astrophysics, we often aim to estimate one or more parameters for each member object in a population and study the distribution of the fitted parameters across the population. In this paper, we develop novel methods that allow us to take advantage of existing software designed for such case-by-case analyses to simultaneously fit parameters of both the individual objects and the parameters that quantify their distribution across the population. Our methods are based on Bayesian hierarchical modelling which is known to produce parameter estimators for the individual objects that are on average closer to their true values than estimators based on case-by-case analyses. We verify this in the context of estimating ages of Galactic halo white dwarfs (WDs) via a series of simulation studies. Finally, we deploy our new techniques on optical and near-infrared photometry of ten candidate halo WDs to obtain estimates of their ages along with an estimate of the mean age of Galactic halo WDs of [11.25, 12.96] Gyr. Although this sample is small, our technique lays the ground work for large-scale studies using data from the Gaia mission.

The mass distribution of stellar mass black holes can provide important clues to supernova modeling, but observationally it is still ill constrained. Therefore it is of importance to make black hole mass measurements as accurate as possible. The X-ray transient 1A 0620-00 is well studied, with a published black hole mass of $6.61\pm0.25\,$M$_{\odot}$, based on an orbital inclination $i$ of $51.0\pm0.9$ degrees. This was obtained by Cantrell et al. (2010), as an average of independent fits to $V$-, $I$- and $H$-band light curves. In this work we perform an independent check on the value of $i$ by re-analyzing existing YALO/SMARTS $V$-, $I$- and $H$-band photometry, using different modeling software and fitting strategy. Performing a fit to the three light curves simultaneously, we obtain a value for $i$ of $54.1\pm1.1$ degrees, resulting in a black hole mass of $5.86\pm0.24\,$M$_{\odot}$. Applying the same model to the light curves individually, we obtain $58.2\pm1.9$, $53.6\pm1.6$ and $50.5\pm2.2$ degrees for $V$-, $I$- and $H$-band, respectively, where the differences in best-fitting $i$ are caused by the contribution of the residual accretion disc light in the three different bands. We conclude that the mass determination of this black hole may still be subject to systematic effects exceeding the statistical uncertainty. Obtaining more accurate masses would be greatly helped by continuous phase-resolved spectroscopic observations simultaneous with photometry.

SS~433 is an X-ray binary and the source of sub-relativistic, precessing, baryonic jets. We present high-resolution spectrograms of SS 433 in the infrared H and K bands. The spectrum is dominated by hydrogen and helium emission lines. The precession phase of the emission lines from the jet continues to be described by a constant period, P_jet= 162.375 d. The limit on any secularly changing period is $|\dot P| \lesssim 10^{-5}$. The He I 2.0587 micron line has complex and variable P Cygni absorption features produced by an inhomogeneous wind with a maximum outflow velocity near 900 km/s. The He II emission lines in the spectrum also arise in this wind. The higher members of the hydrogen Brackett lines show a double-peaked profile with symmetric wings extending more than +/-1500 km/s from the line center. The lines display radial velocity variations in phase with the radial velocity variation expected of the compact star, and they show a distortion during disk eclipse that we interpret as a rotational distortion. We fit the line profiles with a model in which the emission comes from the surface of a symmetric, Keplerian accretion disk around the compact object. The outer edge of the disk has velocities that vary from 110 to 190 km/s. These comparatively low velocities place an important constraint on the mass of the compact star: Its mass must be less than 2.2 M_solar and is probably less than 1.6 M_solar.

We use Hubble Space Telescope (HST) imaging from the ACS Treasury Survey to determine fits for single population isochrones of 69 Galactic globular clusters. Using robust Bayesian analysis techniques, we simultaneously determine ages, distances, absorptions, and helium values for each cluster under the scenario of a "single" stellar population on model grids with solar ratio heavy element abundances. The set of cluster parameters is determined in a consistent and reproducible manner for all clusters using the Bayesian analysis suite BASE-9. Our results are used to re-visit the age-metallicity relation. We find correlations with helium and several other parameters such as metallicity, binary fraction, and proxies for cluster mass. The helium abundances of the clusters are also considered in the context of CNO abundances and the multiple population scenario.

In this paper we apply a Bayesian technique to determine the best fit of stellar evolution models to find the main sequence turn off age and other cluster parameters of four intermediate-age open clusters: NGC 2360, NGC 2477, NGC 2660, and NGC 3960. Our algorithm utilizes a Markov chain Monte Carlo technique to fit these various parameters, objectively finding the best-fit isochrone for each cluster. The result is a high-precision isochrone fit. We compare these results with the those of traditional "by-eye" isochrone fitting methods. By applying this Bayesian technique to NGC 2360, NGC 2477, NGC 2660, and NGC 3960, we determine the ages of these clusters to be 1.35 +/- 0.05, 1.02 +/- 0.02, 1.64 +/- 0.04, and 0.860 +/- 0.04 Gyr, respectively. The results of this paper continue our effort to determine cluster ages to higher precision than that offered by these traditional methods of isochrone fitting.

We use Cycle 21 Hubble Space Telescope (HST) observations and HST archival ACS Treasury observations of 30 Galactic Globular Clusters to characterize two distinct stellar populations. A sophisticated Bayesian technique is employed to simultaneously sample the joint posterior distribution of age, distance, and extinction for each cluster, as well as unique helium values for two populations within each cluster and the relative proportion of those populations. We find the helium differences among the two populations in the clusters fall in the range of ~0.04 to 0.11. Because adequate models varying in CNO are not presently available, we view these spreads as upper limits and present them with statistical rather than observational uncertainties. Evidence supports previous studies suggesting an increase in helium content concurrent with increasing mass of the cluster and also find that the proportion of the first population of stars increases with mass as well. Our results are examined in the context of proposed globular cluster formation scenarios. Additionally, we leverage our Bayesian technique to shed light on inconsistencies between the theoretical models and the observed data.

It takes years of effort employing the best telescopes and instruments to obtain high-quality stellar photometry, astrometry, and spectroscopy. Stellar evolution models contain the experience of lifetimes of theoretical calculations and testing. Yet most astronomers fit these valuable models to these precious datasets by eye. We show that a principled Bayesian approach to fitting models to stellar data yields substantially more information over a range of stellar astrophysics. We highlight advances in determining the ages of star clusters, mass ratios of binary stars, limitations in the accuracy of stellar models, post-main-sequence mass loss, and the ages of individual white dwarfs. We also outline a number of unsolved problems that would benefit from principled Bayesian analyses.

We develop a Bayesian model for globular clusters composed of multiple stellar populations, extending earlier statistical models for open clusters composed of simple (single) stellar populations (vanDyk et al. 2009, Stein et al. 2013). Specifically, we model globular clusters with two populations that differ in helium abundance. Our model assumes a hierarchical structuring of the parameters in which physical properties---age, metallicity, helium abundance, distance, absorption, and initial mass---are common to (i) the cluster as a whole or to (ii) individual populations within a cluster, or are unique to (iii) individual stars. An adaptive Markov chain Monte Carlo (MCMC) algorithm is devised for model fitting that greatly improves convergence relative to its precursor non-adaptive MCMC algorithm. Our model and computational tools are incorporated into an open-source software suite known as BASE-9. We use numerical studies to demonstrate that our method can recover parameters of two-population clusters, and also show model misspecification can potentially be identified. As a proof of concept, we analyze the two stellar populations of globular cluster NGC 5272 using our model and methods. (BASE-9 is available from GitHub: https://github.com/argiopetech/base/releases).

We use Cycle 21 Hubble Space Telescope (HST) observations and HST archival ACS Treasury observations of Galactic Globular Clusters to find and characterize two stellar populations in NGC 5024 (M53), NGC 5272 (M3), and NGC 6352. For these three clusters, both single and double-population analyses are used to determine a best fit isochrone(s). We employ a sophisticated Bayesian analysis technique to simultaneously fit the cluster parameters (age, distance, absorption, and metallicity) that characterize each cluster. For the two-population analysis, unique population level helium values are also fit to each distinct population of the cluster and the relative proportions of the populations are determined. We find differences in helium ranging from $\sim$0.05 to 0.11 for these three clusters. Model grids with solar $\alpha$-element abundances ([$\alpha$/Fe] =0.0) and enhanced $\alpha$-elements ([$\alpha$/Fe]=0.4) are adopted.

There are very few confirmed black holes with a mass that could be $\sim\! 4\, M_\odot$ and no neutron stars with masses greater than $\sim\! 2\, M_\odot$, creating a gap in the observed distribution of compact star masses. Some black holes with masses between 2 and $4\, M_\odot$ might be hiding among other X-ray sources, whose masses are difficult to measure. We present new high-speed optical photometry of the low-mass X-ray binary V1408 Aql (= 4U 1957+115), which is a persistent X-ray source thought to contain a black hole. The optical light curve of V1408~Aql shows a nearly sinusoidal modulation at the orbital period of the system superimposed on large night-to-night variations in mean intensity. We combined the new photometry with previously-published photometry to derive a more precise orbital period, $P = 0.388893(3)$\ d, and to better define the orbital light curve and night-to-night variations. The orbital light curve agrees well with a model in which the modulation is caused entirely by the changing aspect of the heated face of the secondary star. The lack of eclipses rules out orbital inclinations greater than $65^{\circ}$. Our best models for the orbital light curve favor inclinations near $13^{\circ}$ and black hole masses near $3\, M_\odot$ with a 90\% upper bound of $6.2\, M_\odot$, and a lower bound of $2.0\, M_\odot$ imposed solely by the maximum mass of neutron stars. We favor a black hole primary over a neutron star primary based on evidence from the X-ray spectra, the high spin of the compact object, and the fact that a type I X-ray burst has not been observed for this system. Although uncertainties in the data and the models allow higher masses, possibly much higher masses, the compact star in V1408~Aql is a viable candidate for a black hole lying in the mass gap.

BASE-9 is a Bayesian software suite that recovers star cluster and stellar parameters from photometry. BASE-9 is useful for analyzing single-age, single-metallicity star clusters, binaries, or single stars, and for simulating such systems. BASE-9 uses Markov chain Monte Carlo and brute-force numerical integration techniques to estimate the posterior probability distributions for the age, metallicity, helium abundance, distance modulus, and line-of-sight absorption for a cluster, and the mass, binary mass ratio, and cluster membership probability for every stellar object. BASE-9 is provided as open source code on a version-controlled web server. The executables are also available as Amazon Elastic Compute Cloud images. This manual provides potential users with an overview of BASE-9, including instructions for installation and use.

We present contemporaneous X-ray, ultraviolet, optical and near-infrared observations of the black hole binary system, Swift J1753.5-0127, acquired in 2012 October. The UV observations, obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope, are the first UV spectra of this system. The dereddened UV spectrum is characterized by a smooth, blue continuum and broad emission lines of CIV and HeII. The system was stable in the UV to <10% during our observations. We estimated the interstellar reddening by fitting the 2175 A absorption feature and fit the interstellar absorption profile of Ly$\alpha$ to directly measure the neutral hydrogen column density along the line of sight. By comparing the UV continuum flux to steady-state thin accretion disk models, we determined upper limits on the distance to the system as a function of black hole mass. The continuum is well fit with disk models dominated by viscous heating rather than irradiation. The broadband spectral energy distribution shows the system has declined at all wavelengths since previous broadband observations in 2005 and 2007. If we assume that the UV emission is dominated by the accretion disk the inner radius of the disk must be truncated at radii above the ISCO to be consistent with the X-ray flux, requiring significant mass loss from outflows and/or energy loss via advection into the black hole to maintain energy balance.

Compact binary systems with neutron stars or black holes are one of the most promising sources for ground-based gravitational wave detectors. Gravitational radiation encodes rich information about source physics; thus parameter estimation and model selection are crucial analysis steps for any detection candidate events. Detailed models of the anticipated waveforms enable inference on several parameters, such as component masses, spins, sky location and distance that are essential for new astrophysical studies of these sources. However, accurate measurements of these parameters and discrimination of models describing the underlying physics are complicated by artifacts in the data, uncertainties in the waveform models and in the calibration of the detectors. Here we report such measurements on a selection of simulated signals added either in hardware or software to the data collected by the two LIGO instruments and the Virgo detector during their most recent joint science run, including a "blind injection" where the signal was not initially revealed to the collaboration. We exemplify the ability to extract information about the source physics on signals that cover the neutron star and black hole parameter space over the individual mass range 1 Msun - 25 Msun and the full range of spin parameters. The cases reported in this study provide a snap-shot of the status of parameter estimation in preparation for the operation of advanced detectors.

We report a search for gravitational waves from the inspiral, merger and ringdown of binary black holes (BBH) with total mass between 25 and 100 solar masses, in data taken at the LIGO and Virgo observatories between July 7, 2009 and October 20, 2010. The maximum sensitive distance of the detectors over this period for a (20,20) Msun coalescence was 300 Mpc. No gravitational wave signals were found. We thus report upper limits on the astrophysical coalescence rates of BBH as a function of the component masses for non-spinning components, and also evaluate the dependence of the search sensitivity on component spins aligned with the orbital angular momentum. We find an upper limit at 90% confidence on the coalescence rate of BBH with non-spinning components of mass between 19 and 28 Msun of 3.3 \times 10^-7 mergers /Mpc^3 /yr.

We present the results of the first search for gravitational wave bursts associated with high energy neutrinos. Together, these messengers could reveal new, hidden sources that are not observed by conventional photon astronomy, particularly at high energy. Our search uses neutrinos detected by the underwater neutrino telescope ANTARES in its 5 line configuration during the period January - September 2007, which coincided with the fifth and first science runs of LIGO and Virgo, respectively. The LIGO-Virgo data were analysed for candidate gravitational-wave signals coincident in time and direction with the neutrino events. No significant coincident events were observed. We place limits on the density of joint high energy neutrino - gravitational wave emission events in the local universe, and compare them with densities of merger and core-collapse events.

We present the first multi-wavelength follow-up observations of two candidate gravitational-wave (GW) transient events recorded by LIGO and Virgo in their 2009-2010 science run. The events were selected with low latency by the network of GW detectors and their candidate sky locations were observed by the Swift observatory. Image transient detection was used to analyze the collected electromagnetic data, which were found to be consistent with background. Off-line analysis of the GW data alone has also established that the selected GW events show no evidence of an astrophysical origin; one of them is consistent with background and the other one was a test, part of a "blind injection challenge". With this work we demonstrate the feasibility of rapid follow-ups of GW transients and establish the sensitivity improvement joint electromagnetic and GW observations could bring. This is a first step toward an electromagnetic follow-up program in the regime of routine detections with the advanced GW instruments expected within this decade. In that regime multi-wavelength observations will play a significant role in completing the astrophysical identification of GW sources. We present the methods and results from this first combined analysis and discuss its implications in terms of sensitivity for the present and future instruments.

We present contemporaneous, broadband, near-infrared spectroscopy (0.9-2.45 micron) and H-band photometry of the black hole X-ray binary, XTE J1118+480. We determined the fractional dilution of the NIR ellipsoidal light curves of the donor star from other emission sources in the system by comparing the absorption features in the spectrum with field stars of known spectral type. We constrained the donor star spectral type to K7 V - M1 V and determined that the donor star contributed 54+/-27% of the H-band flux at the epoch of our observations. This result underscores the conclusion that the donor star cannot be assumed to be the only NIR emission source in quiescent X-ray binaries. The H-band light curve shows a double-humped asymmetric modulation with extra flux at orbital phase 0.75. The light curve was fit with a donor star model light curve, taking into account a constant second flux component based on the dilution analysis. We also fit models that included emission from the donor star, a constant component from the accretion disk, and a phase-variable component from the bright spot where the mass accretion stream impacts the disk. These simple models with reasonable estimates for the component physical parameters can fully account for the observed light curve, including the extra emission at phase 0.75. From our fits, we constrained the binary inclination to 68 <= i <= 79 deg. This leads to a black hole mass of 6.9 <= M_BH <= 8.2 solar masses. Long-term variations in the NIR light curve shape in XTE J1118+480 are similar to those seen in other X-ray binaries and demonstrate the presence of continued activity and variability in these systems even when in full quiescence.

We present the results of a search for gravitational waves associated with 154 gamma-ray bursts (GRBs) that were detected by satellite-based gamma-ray experiments in 2009-2010, during the sixth LIGO science run and the second and third Virgo science runs. We perform two distinct searches: a modeled search for coalescences of either two neutron stars or a neutron star and black hole; and a search for generic, unmodeled gravitational-wave bursts. We find no evidence for gravitational-wave counterparts, either with any individual GRB in this sample or with the population as a whole. For all GRBs we place lower bounds on the distance to the progenitor, under the optimistic assumption of a gravitational-wave emission energy of 10^-2 M c^2 at 150 Hz, with a median limit of 17 Mpc. For short hard GRBs we place exclusion distances on binary neutron star and neutron star-black hole progenitors, using astrophysically motivated priors on the source parameters, with median values of 16 Mpc and 28 Mpc respectively. These distance limits, while significantly larger than for a search that is not aided by GRB satellite observations, are not large enough to expect a coincidence with a GRB. However, projecting these exclusions to the sensitivities of Advanced LIGO and Virgo, which should begin operation in 2015, we find that the detection of gravitational waves associated with GRBs will become quite possible.

This paper presents results of an all-sky searches for periodic gravitational waves in the frequency range [50, 1190] Hz and with frequency derivative ranges of [-2 x 10^-9, 1.1 x 10^-10] Hz/s for the fifth LIGO science run (S5). The novelty of the search lies in the use of a non-coherent technique based on the Hough-transform to combine the information from coherent searches on timescales of about one day. Because these searches are very computationally intensive, they have been deployed on the Einstein@Home distributed computing project infrastructure. The search presented here is about a factor 3 more sensitive than the previous Einstein@Home search in early S5 LIGO data. The post-processing has left us with eight surviving candidates. We show that deeper follow-up studies rule each of them out. Hence, since no statistically significant gravitational wave signals have been detected, we report upper limits on the intrinsic gravitational wave amplitude h0. For example, in the 0.5 Hz-wide band at 152.5 Hz, we can exclude the presence of signals with h0 greater than 7.6 x 10^-25 with a 90% confidence level.

Uncertainty in the calibration of gravitational-wave (GW) detector data leads to systematic errors which must be accounted for in setting limits on the strength of GW signals. When cross-correlation measurements are made using data from a pair of instruments, as in searches for a stochastic GW background, the calibration uncertainties of the individual instruments can be combined into an uncertainty associated with the pair. With the advent of multi-baseline GW observation (e.g., networks consisting of multiple detectors such as the LIGO observatories and Virgo), a more sophisticated treatment is called for. We describe how the correlations between calibration factors associated with different pairs can be taken into account by marginalizing over the uncertainty associated with each instrument.

Many population synthesis and stellar evolution studies have addressed the evolution of close binary systems in which the primary is a compact remnant and the secondary is filling its Roche lobe, thus triggering mass transfer. Although tidal locking is expected in such systems, most studies have neglected the rotationally-induced mixing that may occur. Here we study the possible effects of mixing in the mass-losing stars for a range in secondary star masses and metallicities. We find that tidal locking can induce rotational mixing prior to contact and thus affect the evolution of the secondary star if the effects of the Spruit-Tayler dynamo are included both for angular momentum and chemical transport. Once contact is made, the effect of mass transfer tends to be more rapid than the evolutionary time scale, so the effects of mixing are no longer directly important, but the mass transfer strips matter to inner layers that may have been affected by the mixing. These effects are enhanced for secondaries of 1-1.2 Msun and for lower metallicities. We discuss the possible implications for the paucity of carbon in the secondaries of the cataclysmic variable SS Cyg and the black hole candidate XTE J1118+480 and for the progenitor evolution of Type Ia supernovae. We also address the issue of the origin of blue straggler stars in globular and open clusters. We find that for models that include rotation consistent with that observed for some blue straggler stars, evolution is chemically homogeneous. This leads to tracks in the HR diagram that are brighter and bluer than the non-rotating main-sequence turn-off point. Rotational mixing could thus be one of the factors that contribute to the formation of blue stragglers.

We present results of optical and infrared photometric monitoring of the eclipsing low-mass X-ray binary V395 Car (2S 0921-630). Our observations reveal a clear, repeating orbital modulation with an amplitude of about one magnitude in B, and V and a little less in J. Combining our data with archival observations spanning about 20 years, we derive an updated ephemeris with orbital period 9.0026+/-0.0001d. We attribute the modulation to a combination of the changing aspect of the irradiated face of the companion star and eclipses of the accretion disk around the neutron star. Both appear to be necessary as a secondary eclipse of the companion star is clearly seen. We model the B, V, and J lightcurves using a simple model of an accretion disk and companion star and find a good fit is possible for binary inclinations of 82.2+/-1.0 degrees. We estimate the irradiating luminosity to be about 8x10^35 erg/s, in good agreement with X-ray constraints.

We present the results of a LIGO search for gravitational waves (GWs) associated with GRB 051103, a short-duration hard-spectrum gamma-ray burst (GRB) whose electromagnetically determined sky position is coincident with the spiral galaxy M81, which is 3.6 Mpc from Earth. Possible progenitors for short-hard GRBs include compact object mergers and soft gamma repeater (SGR) giant flares. A merger progenitor would produce a characteristic GW signal that should be detectable at the distance of M81, while GW emission from an SGR is not expected to be detectable at that distance. We found no evidence of a GW signal associated with GRB 051103. Assuming weakly beamed gamma-ray emission with a jet semi-angle of 30 deg we exclude a binary neutron star merger in M81 as the progenitor with a confidence of 98%. Neutron star-black hole mergers are excluded with > 99% confidence. If the event occurred in M81 our findings support the the hypothesis that GRB 051103 was due to an SGR giant flare, making it the most distant extragalactic magnetar observed to date.

We present HST/ACS ultraviolet photometry of three quiescent black hole X-ray transients: X-ray Nova Muscae 1991 (GU Mus), GRO J0422+32 (V518 Per), and X-ray Nova Vel 1993 (MM Vel), and one neutron star system, Aql X-1. These are the first quiescent UV detections of these objects. All are detected at a much higher level than expected from their companion stars alone and are significant detections of the accretion flow. Three of the four UV excesses can be characterized by a black body of temperature 5000-13,000K, hotter than expected for the quiescent outer disk. A good fit could not be found for MM Vel. The source of the black-body-like emission is most likely a heated region of the inner disk. Contrary to initial indications from spectroscopy there does not appear to be a systematic difference in the UV luminosity or spectral shape between black holes and neutron star systems. However combining our new data with earlier spectroscopy and published X-ray luminosities there is a significant difference in the X-ray to UV flux ratios with the neutron stars exhibiting Lx/Luv about 10x higher than the black hole systems. Since both bandpasses are expected to be dominated by accretion light this suggests the difference in X-ray luminosities cannot simply reflect differences in quiescent accretion rates and so is a more robust discriminator between the black hole and neutron star populations than the comparison of X-ray luminosities alone.

Aims. A transient astrophysical event observed in both gravitational wave (GW) and electromagnetic (EM) channels would yield rich scientific rewards. A first program initiating EM follow-ups to possible transient GW events has been developed and exercised by the LIGO and Virgo community in association with several partners. In this paper, we describe and evaluate the methods used to promptly identify and localize GW event candidates and to request images of targeted sky locations. Methods. During two observing periods (Dec 17 2009 to Jan 8 2010 and Sep 2 to Oct 20 2010), a low-latency analysis pipeline was used to identify GW event candidates and to reconstruct maps of possible sky locations. A catalog of nearby galaxies and Milky Way globular clusters was used to select the most promising sky positions to be imaged, and this directional information was delivered to EM observatories with time lags of about thirty minutes. A Monte Carlo simulation has been used to evaluate the low-latency GW pipeline's ability to reconstruct source positions correctly. Results. For signals near the detection threshold, our low-latency algorithms often localized simulated GW burst signals to tens of square degrees, while neutron star/neutron star inspirals and neutron star/black hole inspirals were localized to a few hundred square degrees. Localization precision improves for moderately stronger signals. The correct sky location of signals well above threshold and originating from nearby galaxies may be observed with ~50% or better probability with a few pointings of wide-field telescopes.

We report the discovery of a 1mHz optical quasi-periodic oscillation (QPO) in the candidate ultracompact low-mass X-ray binary 4U 0614+091. The ultra-low frequency QPO has no X-ray counterpart in contemporaneous RXTE/PCA data and is likely a signature of structure in the accretion disk. The QPO can be reasonably fitted with a single sine wave but with a phase jump part way through the observation, indicating that it is not coherent.We also identify a 48 min modulation, approximately consistent with the suggested orbital period of O'Brien (2005) and Shahbaz et al. (2008). If this is indeed orbital, it supports an identification of 4U 0614+091 as an ultra-compact source.

The gravitational-wave (GW) sky may include nearby pointlike sources as well as astrophysical and cosmological stochastic backgrounds. Since the relative strength and angular distribution of the many possible sources of GWs are not well constrained, searches for GW signals must be performed in a model-independent way. To that end we perform two directional searches for persistent GWs using data from the LIGO S5 science run: one optimized for pointlike sources and one for arbitrary extended sources. The latter result is the first of its kind. Finding no evidence to support the detection of GWs, we present 90% confidence level (CL) upper-limit maps of GW strain power with typical values between 2-20x10^-50 strain^2 Hz^-1 and 5-35x10^-49 strain^2 Hz^-1 sr^-1 for pointlike and extended sources respectively. The limits on pointlike sources constitute a factor of 30 improvement over the previous best limits. We also set 90% CL limits on the narrow-band root-mean-square GW strain from interesting targets including Sco X-1, SN1987A and the Galactic Center as low as ~7x10^-25 in the most sensitive frequency range near 160 Hz. These limits are the most constraining to date and constitute a factor of 5 improvement over the previous best limits.

[Abridged.] We present multiwavelength observations of the black hole binary system, A0620-00. Using the Cosmic Origins Spectrograph on the Hubble Space Telescope, we have obtained the first FUV spectrum of A0620-00. The observed spectrum is flat in the FUV and very faint (with continuum fluxes \simeq 1e - 17 ergs/cm^2/s/A). We compiled the dereddened, broadband spectral energy distribution of A0620-00 and compared it to previous SEDs as well as theoretical models. The SEDs show that the source varies at all wavelengths for which we have multiple samples. Contrary to previous observations, the optical-UV spectrum does not continue to drop to shorter wavelengths, but instead shows a recovery and an increasingly blue spectrum in the FUV. We created an optical-UV spectrum of A0620-00 with the donor star contribution removed. The non-stellar spectrum peaks at \simeq3000 {\deg}A. The peak can be fit with a T=10,000 K blackbody with a small emitting area, probably originating in the hot spot where the accretion stream impacts the outer disk. However, one or more components in addition to the blackbody are needed to fit the FUV upturn and the red optical fluxes in the optical-UV spectrum. By comparing the mass accretion rate determined from the hot spot luminosity to the mean accretion rate inferred from the outburst history, we find that the latter is an order of magnitude smaller than the former, indicating that \sim90% of the accreted mass must be lost from the system if the predictions of the disk instability model and the estimated interoutburst interval are correct. The mass accretion rate at the hot spot is 10^5 the accretion rate at the black hole inferred from the X-ray luminosity. To reconcile these requires that outflows carry away virtually all of the accreted mass, a very low rate of mass transfer from the outer cold disk into the inner hot region, and/or radiatively inefficient accretion.