A dozen short-period detached binaries are known to host transiting circumbinary planets. In all circumbinary systems so far, the planetary and binary orbits are aligned within a couple of degrees. However, the obliquity of the primary star, which is an important tracer of their formation, evolution, and tidal history, has only been measured in one circumbinary system until now. EBLM J0608-59/TOI-1338 is a low-mass eclipsing binary system with a recently discovered circumbinary planet identified by TESS. Here, we perform high-resolution spectroscopy during primary eclipse to measure the projected stellar obliquity of the primary component. The obliquity is low, and thus the primary star is aligned with the binary and planetary orbits with a projected spin-orbit angle $\beta = 2.8 \pm 17.1$ deg. The rotation period of $18.1 \pm 1.6$ days implied by our measurement of $v\sin{i_\star}$ suggests that the primary has not yet pseudo-synchronized with the binary orbit, but is consistent with gyrochronology and weak tidal interaction with the binary companion. Our result, combined with the known coplanarity of the binary and planet orbits, is suggestive of formation from a single disc. Finally, we considered whether the spectrum of the faint secondary star could affect our measurements. We show through simulations that the effect is negligible for our system, but can lead to strong biases in $v\sin{i_\star}$ and $\beta$ for higher flux ratios. We encourage future studies in eclipse spectroscopy test the assumption of a dark secondary for flux ratios $\gtrsim 1$ ppt.

Accurate masses and radii for normal stars derived from observations of detached eclipsing binary stars are of fundamental importance for testing stellar models and may be useful for calibrating free parameters in these model if the masses and radii are sufficiently precise and accurate. We aim to measure precise masses and radii for the stars in the bright eclipsing binary AI Phe, and to quantify the level of systematic error in these estimates. We use several different methods to model the TESS light curve of AI Phe combined with spectroscopic orbits from multiple sources to estimate precisely the stellar masses and radii together with robust error estimates. We find that the agreement between different methods for the light curve analysis is very good but some methods underestimate the errors on the model parameters. The semi-amplitudes of the spectroscopic orbits derived from spectra obtained with modern echelle spectrographs are consistent to within 0.1%. The masses of the stars in AI Phe are $M_1 = 1.1938 \pm 0.0008 M_{\odot}$ and $M_2 = 1.2438 \pm 0.0008M_{\odot}$, and the radii are $R_1 = 1.8050 \pm 0.0022 R_{\odot}$ and $R_2 = 2.9332 \pm 0.0023 R_{\odot}$. We conclude that it is possible to measure accurate masses and radii for stars in bright eclipsing binary stars to a precision of 0.2% or better using photometry from TESS and spectroscopy obtained with modern echelle spectrographs. We provide recommendations for publishing masses and radii of eclipsing binary stars at this level of precision.

We report the detection of the first circumbinary planet found by TESS. The target, a known eclipsing binary, was observed in sectors 1 through 12 at 30-minute cadence and in sectors 4 through 12 at two-minute cadence. It consists of two stars with masses of 1.1 MSun and 0.3 MSun on a slightly eccentric (0.16), 14.6-day orbit, producing prominent primary eclipses and shallow secondary eclipses. The planet has a radius of ~6.9 REarth and was observed to make three transits across the primary star of roughly equal depths (~0.2%) but different durations -- a common signature of transiting circumbinary planets. Its orbit is nearly circular (e ~ 0.09) with an orbital period of 95.2 days. The orbital planes of the binary and the planet are aligned to within ~1 degree. To obtain a complete solution for the system, we combined the TESS photometry with existing ground-based radial-velocity observations in a numerical photometric-dynamical model. The system demonstrates the discovery potential of TESS for circumbinary planets, and provides further understanding of the formation and evolution of planets orbiting close binary stars.

The black hole MAXI J1820+070 was discovered during its 2018 outburst and was extensively monitored across the electromagnetic spectrum. Following the detection of relativistic radio jets, we obtained four Chandra X-ray observations taken between 2018 November and 2019 May, along with radio observations conducted with the VLA and MeerKAT arrays. We report the discovery of X-ray sources associated with the radio jets moving at relativistic velocities with a possible deceleration at late times. The broadband spectra of the jets are consistent with synchrotron radiation from particles accelerated up to very high energies (>10 TeV) by shocks produced by the jets interacting with the interstellar medium. The minimal internal energy estimated from the X-ray observations for the jets is $\sim 10^{41}$ erg, significantly larger than the energy calculated from the radio flare alone, suggesting most of the energy is possibly not radiated at small scales but released through late-time interactions.

We present the discovery of a likely new redback millisecond pulsar binary associated with the \emph{Fermi} $\gamma$-ray source 4FGL J2333.1--5527. Using optical photometric and spectroscopic observations from the SOAR telescope, we identify a low-mass, main sequence-like companion in a 6.9-hr, highly inclined orbit around a suspected massive neutron star primary. Archival XMM-Newton X-ray observations show this system has a hard power-law spectrum $\Gamma = 1.6\pm0.3$ and $L_X \sim 5 \times 10^{31}$ erg s$^{-1}$, consistent with redback millisecond pulsar binaries. Our data suggest that for secondary masses typical of redbacks, the mass of the neutron star is likely well in excess of $\sim1.4\,M_{\odot}$, but future timing of the radio pulsar is necessary to bolster this tentative conclusion. This work shows that a bevy of nearby compact binaries still await discovery, and that unusually massive neutron stars continue to be common in redbacks.

We report the discovery of a Neptune-size (R_p = 3.87 +/- 0.06 R_Earth) transiting circumbinary planet, Kepler-1661 b, found in the Kepler photometry. The planet has a period of ~175 days and its orbit precesses with a period of only 35 years. The precession causes the alignment of the orbital planes to vary, and the planet is in a transiting configuration only ~7% of the time as seen from Earth. As with several other Kepler circumbinary planets, Kepler-1661 b orbits close to the stability radius, and is near the (hot) edge of habitable zone. The planet orbits a single-lined, grazing eclipsing binary, containing a 0.84 M_Sun and 0.26 M_Sun pair of stars in a mildly eccentric (e=0.11), 28.2-day orbit. The system is fairly young, with an estimated age of ~1-3 Gyrs, and exhibits significant starspot modulations. The grazing-eclipse configuration means the system is very sensitive to changes in the binary inclination, which manifests itself as a change in the eclipse depth. The starspots contaminate the eclipse photometry, but not in the usual way of inducing spurious eclipse timing variations. Rather, the starspots alter the normalization of the light curve, and hence the eclipse depths. This can lead to spurious eclipse depth variations, which are then incorrectly ascribed to binary orbital precession.

There is a well-known stellar parameter discrepancy for late K and M dwarfs, in that the observed radii and temperatures are often respectively larger and cooler than predicted by theory by several percent. In an on-going effort to elucidate this issue, we examine the double-lined Kepler eclipsing binary star system KIC 8736245. We supplement the near-continuous 4-year Kepler light curve with ground-based multicolor photometry from Mount Laguna Observatory and spectroscopy from the Hobby-Eberly Telescope. The binary has an edge-on, circular 5.07 day orbit with stellar masses equal to $0.987\pm0.009$ and $0.782\pm0.009\,\text{M}_\odot$ and radii of $1.311 \pm 0.006$ and $0.804 \pm 0.004\,\text{R}_\odot$, respectively, and an estimated age of 7-9 Gyr. We find that the stellar radii are consistent with theoretical models within the uncertainties, whereas the temperature of the secondary star is $\sim$6% cooler than predicted. An important aspect of this work is that the uncertainties derived from a single epoch (individual night of observations) underestimates the overall system parameter uncertainties due to the effect of the 1-4% fluctuations caused by stellar activity. Our error estimates come from the spread in parameters measured at 8 epochs. From the periodicities in the light curve and from the eclipse times, we measure candidate spin periods to be approximately 4.98 and 5.87 days for the primary and secondary star. Surprisingly, these imply super- and sub-synchronous rotation compared to the orbital period. Thus KIC 8736245 serves as an interesting case study for the exchange of angular momentum and general stellar astrophysics as stars in binaries evolve off the main sequence.

To date a dozen transiting "Tatooines" or circumbinary planets (CBPs) have been discovered, by eye, in the data from the Kepler mission; by contrast, thousands of confirmed circumstellar planets orbiting around single stars have been detected using automated algorithms. Automated detection of CBPs is challenging because their transits are strongly aperiodic with irregular profiles. Here, we describe an efficient and automated technique for detecting circumbinary planets that transit their binary hosts in Kepler light curves. Our method accounts for large transit timing and duration variations (TTVs and TDVs), induced by binary reflex motion, in two ways: 1) We directly correct for large-scale TTVs and TDVs in the light curves by using Keplerian models to approximate binary and CBP orbits; and 2) We allow additional aperiodicities on the corrected light curves by employing the Quasi-periodic Automated Transit Search algorithm (QATS). We demonstrate that our method dramatically improves detection significance using simulated data and two previously identified CBP systems, Kepler-35 and Kepler-64.

Recently there has been a renewed interest in the power-2 limb darkening law for modeling exoplanet transits. This law provides a better match to the intensities generated by spherical stellar atmosphere models than other 2-parameter laws. To help facilitate a wider use of the power-2 law we correct a minor error and, expanding on previous work, suggest a parametrization that can improve the sampling required by some numerical methods such as MCMC.

Kepler-1647b is the most recently discovered planet that transits two stars, i.e., a circumbinary planet (CBP). Due to its large orbital separation, Kepler-1647b stands out from the rest of the Kepler CBPs, which mostly reside on much tighter orbits near the stability limit. The large separation of Kepler-1647b challenges inward disk migration as a dominant formation pathway, suggested by the other Kepler CBPs. In this paper, we consider the possibility of an undetected planet near the stability limit by examining observational consequences of such a planet. We calculate the transit probability of the putative planet, transit timing variations (TTVs) of the known planet, and eclipsing timing variations (ETVs) of the host binary caused by the putative planet. We find the presence of a $\gtrsim$30M$_{\oplus}$ inner planet to be highly unlikely near the stability limit. In addition, we provide future TTV observation windows, which will further constrain possible undetected planets with lower masses.

Black holes in binary star systems are vital for understanding the process of pr oducing gravitational wave sources, understanding how supernovae work, and for p roviding fossil evidence for the high mass stars from earlier in the Universe. At the present time, sample sizes of these objects, and especially of black hole s in binaries, are quite limited. Furthermore, more precise measurements of the binary parameters are needed, as well. With improvements primarily in X-ray an d radio astronomy capabilities, it should be possible to build much larger sampl es of much better measured black hole binaries.

Of the nine confirmed transiting circumbinary planet systems, only Kepler-47 is known to contain more than one planet. Kepler-47 b (the "inner planet") has an orbital period of 49.5 days and a radius of about $3\,R_{\oplus}$. Kepler-47 c (the "outer planet") has an orbital period of 303.2 days and a radius of about $4.7\,R_{\oplus}$. Here we report the discovery of a third planet, Kepler-47 d (the "middle planet"), which has an orbital period of 187.4 days and a radius of about $7\,R_{\oplus}$. The presence of the middle planet allows us to place much better constraints on the masses of all three planets, where the $1\sigma$ ranges are less than $26\,M_{\oplus}$, between $7-43\,M_{\oplus}$, and between $2-5\,M_{\oplus}$ for the inner, middle, and outer planets, respectively. The middle and outer planets have low bulk densities, with $\rho_{\rm middle} < 0.68$ g cm$^{-3}$ and $\rho_{\rm outer} < 0.26$ g cm$^{-3}$ at the $1\sigma$ level. The two outer planets are "tightly packed," assuming the nominal masses, meaning no other planet could stably orbit between them. All of the orbits have low eccentricities and are nearly coplanar, disfavoring violent scattering scenarios and suggesting gentle migration in the protoplanetary disk.

We present the first optical spectroscopy of five confirmed (or strong candidate) redback millisecond pulsar binaries, obtaining complete radial velocity curves for each companion star. The properties of these millisecond pulsar binaries with low-mass, hydrogen-rich companions are discussed in the context of the 14 confirmed and 10 candidate field redbacks. We find that the neutron stars in redbacks have a median mass of 1.78 +/- 0.09 M_sun with a dispersion of sigma = 0.21 +/- 0.09. Neutron stars with masses in excess of 2 M_sun are consistent with, but not firmly demanded by, current observations. Redback companions have median masses of 0.36 +/- 0.04 M_sun with a scatter of sigma = 0.15 +/- 0.04, and a tail possibly extending up to 0.7-0.9 M_sun. Candidate redbacks tend to have higher companion masses than confirmed redbacks, suggesting a possible selection bias against the detection of radio pulsations in these more massive candidate systems. The distribution of companion masses between redbacks and the less massive black widows continues to be strongly bimodal, which is an important constraint on evolutionary models for these systems. Among redbacks, the median efficiency of converting the pulsar spindown energy to gamma-ray luminosity is ~10%.

We present here an efficient method for computing the visible flux for each body during a multi-body eclipsing event for all commonly used limb darkening laws. Our approach follows the idea put forth by Pal (2012) to apply Green's Theorem on the limb darkening integral, thus transforming the two-dimensional flux integral over the visible disk into a one-dimensional integral over the visible boundary. We implement this idea through an iterative process which combines a fast method for describing the visible boundary of each body with a fast numerical integration scheme to compute the integrals. For the two-body case, our method compares well in speed with both that of Mandel & Agol (2002) and that of Gimenez (2006a). The strength of the method is that it works for any number of spherical bodies, with a computational accuracy that is adjustable through the use of a tolerance parameter. Most significantly, the method offers two main advantages over previously used techniques: (i) it can employ a multitude of limb darkening laws, including all of the commonly used ones; (ii) it can compute the Rossiter-McLaughlin effect for rigid body rotation with an arbitrary orientation of the rotation axis, using any of these limb darkening laws. In addition, we can compute the Rossiter-McLaughlin effect for stars exhibiting differential rotation, using the quadratic limb darkening law. We provide the mathematical background for the method and explain in detail how to implement the technique with the help of several examples and codes which we make available.

KIC 9832227 is a contact binary whose 11 hr orbital period is rapidly changing. Based on the apparent exponential decay of its period, the two stars were predicted to merge in early 2022 resulting in a rare red nova outburst. Fortunately KIC 832227 was observed in 2003 as part of the NASA Ames pre-Kepler Vulcan Project to search for transiting exoplanets. We find that the Vulcan timing measurement does not agree with the previous exponential decay model. This led us to re-evaluate the other early epoch non-Kepler data sets, the Northern Sky Variability Survey (NSVS) and Wide Angle Search for Planets (WASP) survey. We find that the WASP times are in good agreement with the previous prediction, but the NSVS eclipse time differs by nearly an hour. The very large disagreement of the Vulcan and NSVS eclipse times with an exponentially decaying model forces us to reject the merger hypothesis. Although period variations are common in contact binaries, the physical cause of the period changes in KIC 9832227 remains unexplained; a third star scenario is unlikely. This study shows the data collected by the Vulcan photometer to be extremely valuable for extending the baseline for measurements of variable stars in the Kepler field.

We determine the properties of the binary star V106 in the old open cluster NGC6791. We identify the system to be a blue straggler cluster member by using a combination of ground-based and Kepler photometry and multi-epoch spectroscopy. The properties of the primary component are found to be $M_{\rm p}\sim1.67 \rm M_{\odot}$, more massive than the cluster turn-off, with $R_{\rm p}\sim1.91 \rm R_{\odot}$ and $T_{\rm eff}=7110\pm100$ K. The secondary component is highly oversized and overluminous for its low mass with $M_{\rm s}\sim0.182 \rm M_{\odot}$, $R_{\rm s}\sim0.864 \rm R_{\odot}$ and $T_{\rm eff}=6875\pm200$ K. We identify this secondary star as a bloated (proto) extremely low-mass helium white dwarf. These properties of V106 suggest that it represents a typical Algol-paradox system and that it evolved through a mass-transfer phase which provides insight into its past evolution. We present a detailed binary stellar evolution model for the formation of V106 using the MESA code and find that the mass-transfer phase only ceased about 40 Myr ago. Due to the short orbital period (P=1.4463 d) another mass-transfer phase is unavoidable once the current primary star evolves towards the red giant phase. We argue that V106 will evolve through a common-envelope phase within the next 100 Myr and merge to become a single over-massive giant. The high mass will make it appear young for its true age, which is revealed by the cluster properties. Therefore, V106 is potentially a prototype progenitor of old field giants masquerading as young.

We present an analysis of a slightly eccentric ($e=0.05$), partially eclipsing long-period ($P = 69.73$ d) main sequence binary system (WOCS 12009, Sanders 1247) in the benchmark old open cluster M67. Using Kepler K2 and ground-based photometry along with a large set of new and reanalyzed spectra, we derived highly precise masses ($1.111\pm0.015$ and $0.748\pm0.005 M_\odot$) and radii ($1.071\pm0.008\pm0.003$ and $0.713\pm0.019\pm0.026 R_\odot$, with statistical and systematic error estimates) for the stars. The radius of the secondary star is in agreement with theory. The primary, however, is approximately $15\%$ smaller than reasonable isochrones for the cluster predict. Our best explanation is that the primary star was produced from the merger of two stars, as this can also account for the non-detection of photospheric lithium and its higher temperature relative to other cluster main sequence stars at the same $V$ magnitude. To understand the dynamical characteristics (low measured rotational line broadening of the primary star and the low eccentricity of the current binary orbit), we believe that the most probable (but not the only) explanation is the tidal evolution of a close binary within a primordial triple system (possibly after a period of Kozai-Lidov oscillations), leading to merger approximately 1Gyr ago. This star appears to be a future blue straggler that is being revealed as the cluster ages and the most massive main sequence stars die out.

We present an in-depth study of the large-scale, western jet of the microquasar XTE J1550-564, based on X-ray and radio observations performed in 2002-2003. The jet is spatially resolved in both observing windows. The X-ray jet is expanding in time along the axis of the jet's propagation: we observe the formation of a tail (~2.25"), which appears to extend backwards with an apparent velocity ~-0.10c. The origin of this feature is discussed in the framework of scenarios of energy dissipation. A single power-law adequately describes the broadband spectra, supporting a synchrotron origin of the X-ray emission. However, a spectral break at ~10^{15} Hz is necessary in coincidence with a re-flare at 8.64 GHz in September 2002. This finding may be indicative of emission from newly accelerated low-energy particles. The first detection of the jet is in February 2001 (F_{8.64GHz}=0.25+/-0.09 mJy) in the flux rising phase. A phase of stable emission is followed by a rapid decay (t_{decay}=167+/-5 days). The decay at radio frequencies is significantly shorter than in X-rays (t_{decay}=338+/-14 days). We detected a high fraction (up to ~9%) of linearly polarized radiation at 4.8 GHz and 8.6 GHz. The orientation of the electric vector is consistent with the picture of a shock-compressed magnetic field, and there are hints of variations on month-timescales, possibly connected with the evolution of the jet structure.

The body of work presented here revolves around the investigation of the existence and nature of extra-solar planetary systems. The fitting of stellar radial velocity time series data is attempted by constructing a model to quantify the orbital properties of a star-planetary system. This is achieved with the Planetary Orbit Fitting Process (POFP). Though specific to the investigated problem, the POFP is founded on two separate, more general ideas. One is a Solver producing the gravitational dynamics of a Three-Body system by integrating its Newtonian equations of motion. The other is an independent optimisation scheme. Both have been devised using MATLAB. Applying the optimisation to the Solver results in a realistic Three-Body dynamics that best describes the radial velocity data under the model-specific orbital-observational constraints. Combining these aspects also allows for the study of dynamical instability derived from interaction, which is reaffirmed as a necessary criterion for evaluating the fit. The validity of POFP solutions with respect to the observations and other models is discussed in this context. The underlying generality and fundamental principles demonstrate a larger frame of operation where problems in Physics and Mathematics can be solved with a multitude of techniques.

In this letter we report the possible existence of a second planet in the transiting extrasolar planet system HD 17156 and its interactive dynamics with the previously known planet. The analysis is achieved through the \POFP\ optimization software which is based on a full integration of the system's multiple-body Newtonian equations of motion. The two-planet solution yields a significantly improved fit to the previously published radial velocities. The two planets are strongly interacting and exchange angular momentum in a 5:1 mean motion resonance, yet remain stable as they mutually excite orbital eccentricities and periastron advances.

We present an extensive set of photometry and radial velocities for the detached eclipsing binary KV 29 in the intermediate-aged open cluster M11 (NGC 6705). Spectroscopy shows that the system is double-lined and all available evidence (proper motion, photometry, and position on the CMD) indicates it is a member of the cluster. We find the period of the binary to be $4.64276\pm0.00001$ days. We find masses $3.604^{+0.002}_{-0.011}M_\odot$ and $1.837^{+0.001}_{-0.006}M_\odot$, and radii $5.392^{+0.018}_{-0.035}R_\odot$ and $1.656^{+0.007}_{-0.044}R_\odot$ for the primary and secondary stars, respectively. Because the primary star in the binary is rapidly evolving and is brighter than the cluster turnoff in a color-magnitude diagram, the measurement of its radius leads to a strong constraint on the cluster age. We find the age of M11 to be $222^{+2}_{-3}\pm15$ Myr, where the quoted uncertainties come from statistical errors in the calculated masses and radii, and systematic uncertainties due to the ambiguity of the metallicity of the open cluster and variations within the isochrone models concerning heavy elements and convective overshooting.

We present the analysis of an eccentric, partially eclipsing long-period ($P=19.23$ d) binary system KIC 9777062 that contains main sequence stars near the turnoff of the intermediate age open cluster NGC 6811. The primary is a metal-lined Am star with a possible convective blueshift to its radial velocities, and one star (probably the secondary) is likely to be a $\gamma$ Dor pulsator. The component masses are $1.603\pm0.006$(stat.)$\pm0.016$(sys.) and $1.419\pm0.003\pm0.008 M_\odot$, and the radii are $1.744\pm0.004\pm0.002$ and $1.544\pm0.002\pm0.002 R_\odot$. The isochrone ages of the stars are mildly inconsistent: the age from the mass-radius combination for the primary ($1.05\pm0.05\pm0.09$ Gyr, where the last quote was systematic uncertainty from models and metallicity) is smaller than that from the secondary ($1.21\pm0.05\pm0.15$ Gyr) and is consistent with the inference from the color-magnitude diagram ($1.00\pm0.05$ Gyr). We have improved the measurements of the asteroseismic parameters $\Delta \nu$ and $\nu_{\rm max}$ for helium-burning stars in the cluster. The masses of the stars appear to be larger (or alternately, the radii appear to be smaller) than predicted from isochrones using the ages derived from the eclipsing stars. The majority of stars near the cluster turnoff are pulsating stars: we identify a sample of 28 $\delta$ Sct, 15 $\gamma$ Dor, and 5 hybrid types. We used the period-luminosity relation for high-amplitude $\delta$ Sct stars to fit the ensemble of the strongest frequencies for the cluster members, finding $(m-M)_V = 10.37\pm0.03$. This is larger than most previous determinations, but smaller than values derived from the eclipsing binary ($10.47\pm0.05$).

We report the discovery of a new Kepler transiting circumbinary planet (CBP). This latest addition to the still-small family of CBPs defies the current trend of known short-period planets orbiting near the stability limit of binary stars. Unlike the previous discoveries, the planet revolving around the eclipsing binary system Kepler-1647 has a very long orbital period (~1100 days) and was at conjunction only twice during the Kepler mission lifetime. Due to the singular configuration of the system, Kepler-1647b is not only the longest-period transiting CBP at the time of writing, but also one of the longest-period transiting planets. With a radius of 1.06+/-0.01 RJup it is also the largest CBP to date. The planet produced three transits in the light-curve of Kepler-1647 (one of them during an eclipse, creating a syzygy) and measurably perturbed the times of the stellar eclipses, allowing us to measure its mass to be 1.52+/-0.65 MJup. The planet revolves around an 11-day period eclipsing binary consisting of two Solar-mass stars on a slightly inclined, mildly eccentric (e_bin = 0.16), spin-synchronized orbit. Despite having an orbital period three times longer than Earth's, Kepler-1647b is in the conservative habitable zone of the binary star throughout its orbit.

Eclipsing binaries serve as a valuable source of stellar masses and radii that inform stellar evolutionary models and provide insight into additional astrophysical processes. The exquisite light curves generated by space-based missions such as Kepler offer the most stringent tests to date. We use the Kepler light curve of the 4.8-day eclipsing binary KIC 5739896 with ground based optical spectra to derive fundamental parameters for the system. We reconstruct the component spectra to determine the individual atmospheric parameters, and model the Kepler photometry with the binary synthesis code ELC to obtain accurate masses and radii. The two components of KIC 5738698 are F-type stars with M1 = 1.39+/-0.04M, M2 = 1.34+/-0.06M, and R1 = 1.84+/-0.03R, R2 = 1.72+/-0.03R. We also report a small eccentricity (e < 0.0017) and unusual albedo values that are required to match the detailed shape of the Kepler light curve. Comparisons with evolutionary models indicate an approximate age of 2.3 Gyr for the system.

As part of our study of the old (~2.5 Gyr) open cluster NGC 6819 in the Kepler field, we present photometric (Kepler and ground-based BVRcIc) and spectroscopic observations of the detached eclipsing binary WOCS 24009 (Auner 665; KIC 5023948) with a short orbital period of 3.6 days. WOCS 24009 is a triple-lined system, and we verify that the brightest star is physically orbiting the eclipsing binary using radial velocities and eclipse timing variations. The eclipsing binary components have masses M_B =1.090+/-0.010 Msun and M_C =1.075+/-0.013 Msun, and radii R_B =1.095+/-0.007 Rsun and R_C =1.057+/-0.008 Rsun. The bright non-eclipsing star resides at the cluster turnoff, and ultimately its mass will directly constrain the turnoff mass: our preliminary determination is M_A =1.251+/-0.057 Msun. A careful examination of the light curves indicates that the fainter star in the eclipsing binary undergoes a very brief period of total eclipse, which enables us to precisely decompose the light of the three stars and place them in the color-magnitude diagram. We also present improved analysis of two previously discussed detached eclipsing stars in NGC 6819 (WOCS 40007 and WOCS 23009) en route to a combined determination of the cluster's distance modulus (m-M)_V=12.38+/-0.04. Because this paper significantly increases the number of measured stars in the cluster, we can better constrain the age of the color-magnitude diagram to be 2.21+/-0.10+/-0.20 Gyr. Additionally, using all measured eclipsing binary star masses and radii, we constrain the age to 2.38+/-0.05+/-0.22 Gyr. The quoted uncertainties are estimates of measurement and systematic uncertainties (due to model physics differences and metal content), respectively.

We combine Kepler photometry with ground-based spectra to present a comprehensive dynamical model of the double red giant eclipsing binary KIC 9246715. While the two stars are very similar in mass (M1 = 2.171 [+0.006 / -0.008], M2 = 2.149 [+0.006 / -0.008] Msun) and radius (R1 = 8.37 [+0.03 / -0.07], R2 = 8.30 [+0.04 / -0.03] Rsun), an asteroseismic analysis finds one main set of solar-like oscillations with unusually low-amplitude, wide modes. A second set of oscillations from the other star may exist, but this marginal detection is extremely faint. Because the two stars are nearly twins, KIC 9246715 is a difficult target for a precise test of the asteroseismic scaling relations, which yield M = 2.17 +/- 0.14 Msun and R = 8.26 +/- 0.18 Rsun. Both stars are consistent with the inferred asteroseismic properties, but we suspect the main oscillator is Star 2 because it is less active than Star 1. We find evidence for stellar activity and modest tidal forces acting over the 171-day eccentric orbit, which are likely responsible for the essential lack of solar-like oscillations in one star and weak oscillations in the other. Mixed modes indicate the main oscillating star is on the secondary red clump (a core-He-burning star), and stellar evolution modeling supports this with a coeval history for a pair of red clump stars. This system is a useful case study and paves the way for a detailed analysis of more red giants in eclipsing binaries, an important benchmark for asteroseismology.

We present a dynamical study of the Galactic black hole binary system Nova Muscae 1991 (GS/GRS 1124-683). We utilize 72 high resolution Magellan Echellette (MagE) spectra and 72 strictly simultaneous V-band photometric observations; the simultaneity is a unique and crucial feature of this dynamical study. The data were taken on two consecutive nights and cover the full 10.4-hour orbital cycle. The radial velocities of the secondary star are determined by cross-correlating the object spectra with the best-match template spectrum obtained using the same instrument configuration. Based on our independent analysis of five orders of the echellette spectrum, the semi-amplitude of the radial velocity of the secondary is measured to be K_2 = 406.8+/-2.7 km/s, which is consistent with previous work, while the uncertainty is reduced by a factor of 3. The corresponding mass function is f(M) = 3.02+/-0.06 M_\odot. We have also obtained an accurate measurement of the rotational broadening of the stellar absorption lines (v sin i = 85.0+/-2.6 km/s) and hence the mass ratio of the system q = 0.079+/-0.007. Finally, we have measured the spectrum of the non-stellar component of emission that veils the spectrum of the secondary. In a future paper, we will use our veiling-corrected spectrum of the secondary and accurate values of K_2 and q to model multi-color light curves and determine the systemic inclination and the mass of the black hole.

The primary Kepler Mission provided nearly continuous monitoring of ~200,000 objects with unprecedented photometric precision. We present the final catalog of eclipsing binary systems within the 105 square degree Kepler field of view. This release incorporates the full extent of the data from the primary mission (Q0-Q17 Data Release). As a result, new systems have been added, additional false positives have been removed, ephemerides and principal parameters have been recomputed, classifications have been revised to rely on analytical models, and eclipse timing variations have been computed for each system. We identify several classes of systems including those that exhibit tertiary eclipse events, systems that show clear evidence of additional bodies, heartbeat systems, systems with changing eclipse depths, and systems exhibiting only one eclipse event over the duration of the mission. We have updated the period and galactic latitude distribution diagrams and included a catalog completeness evaluation. The total number of identified eclipsing and ellipsoidal binary systems in the Kepler field of view has increased to 2878, 1.3% of all observed Kepler targets. An online version of this catalog with downloadable content and visualization tools is maintained at http://keplerEBs.villanova.edu.

The bright soft X-ray transient Nova Muscae 1991 was intensively observed during its entire 8-month outburst using the Large Area Counter (LAC) onboard the Ginga satellite. Recently, we obtained accurate estimates of the mass of the black hole primary, the orbital inclination angle of the system, and the distance. Using these crucial input data and Ginga X-ray spectra, we have measured the spin of the black hole using the continuum-fitting method. For four X-ray spectra of extraordinary quality we have determined the dimensionless spin parameter of the black hole to be a/M = 0.63 (-0.19, +0.16) (1 sigma confidence level), a result that we confirm using eleven additional spectra of lower quality. Our spin estimate challenges two published results: It is somewhat higher than the value predicted by a proposed relationship between jet power and spin; and we find that the spin of the black hole is decidedly prograde, not retrograde as has been claimed.

The optical counterpart of the black-hole soft X-ray transient Nova Muscae 1991 has brightened by $\Delta{V}\approx0.8$ mag since its return to quiescence 23 years ago. We present the first clear evidence that the brightening of soft X-ray transients in quiescence occurs at a nearly linear rate. This discovery, and our precise determination of the disk component of emission obtained using our $simultaneous$ photometric and spectroscopic data, have allowed us to identify and accurately model archival ellipsoidal light curves of the highest quality. The simultaneity, and the strong constraint it provides on the component of disk emission, is a key element of our work. Based on our analysis of the light curves, and our earlier measurements of the mass function and mass ratio, we have obtained for Nova Muscae 1991 the first accurate estimates of its systemic inclination $i=43.2^{+2.1}_{-2.7}$ deg, and black hole mass $M=11.0^{+2.1}_{-1.4}\ M_\odot$. Based on our determination of the radius of the secondary, we estimate the distance to be $D=4.95^{+0.69}_{-0.65}$ kpc. We discuss the implications of our work for future dynamical studies of black-hole soft X-ray transients.

Many short-period binary stars have distant orbiting companions that have played a role in driving the binary components into close separation. Indirect detection of a tertiary star is possible by measuring apparent changes in eclipse times of eclipsing binaries as the binary orbits the common center of mass. Here we present an analysis of the eclipse timings of 41 eclipsing binaries observed throughout the NASA Kepler mission of long duration and precise photometry. This subset of binaries is characterized by relatively deep and frequent eclipses of both stellar components. We present preliminary orbital elements for seven probable triple stars among this sample, and we discuss apparent period changes in seven additional eclipsing binaries that may be related to motion about a tertiary in a long period orbit. The results will be used in ongoing investigations of the spectra and light curves of these binaries for further evidence of the presence of third stars.

Mazeh, Holczer, and Shporer (2015) have presented an approach that can, in principle, use the derived transit timing variation (TTV) of some transiting planets observed by the $Kepler$ mission to distinguish between prograde and retrograde motion of their orbits with respect to their parent stars' rotation. The approach utilizes TTVs induced by spot-crossing events that occur when the planet moves across a spot on the stellar surface, looking for a correlation between the derived TTVs and the stellar brightness derivatives at the corresponding transits. This can work even in data that cannot temporally resolve the spot-crossing events themselves. Here we apply this approach to the $Kepler$ KOIs, identifying nine systems where the photometric spot modulation is large enough and the transit timing accurate enough to allow detection of a TTV-brightness-derivatives correlation. Of those systems five show highly significant prograde motion (Kepler-17b, Kepler-71b, KOI-883.01, KOI-895.01, and KOI-1074.01), while no system displays retrograde motion, consistent with the suggestion that planets orbiting cool stars have prograde motion. All five systems have impact parameter $0.2\lesssim b\lesssim0.5$, and all systems within that impact parameter range show significant correlation, except HAT-P-11b where the lack of a correlation follows its large stellar obliquity. Our search suffers from an observational bias against detection of high impact parameter cases, and the detected sample is extremely small. Nevertheless, our findings may suggest that stellar spots, or at least the larger ones, tend to be located at a low stellar latitude, but not along the stellar equator, similar to the Sun.

KIC 2835289 is a triple stellar system that consists of an inner ellipsoidal variable with an orbital period of ~0.86 days and an outer star that eclipses the inner pair every ~750 days. Two eclipse events were observed by the Kepler mission, but they do not fully constrain our photodynamical models. The next eclipse event will occur on May 14 UT, and we solicit community involvement to follow it up from the ground. All details are available in the attached call. Please contact the authors to join the follow-up campaign.

The Kepler mission has provided high quality light curves for more than 2000 eclipsing binaries. Tertiary companions to these binaries can be detected if they transit one or both stars in the binary or if they perturb the binary enough to cause deviations in the observed times of the primary and secondary eclipses (in a few cases both effects are observed in the same eclipsing binary). From the study of eclipse timing variations, it is estimated that 15 to 20% of the Kepler eclipsing binaries have close-in tertiary companions. I will give an overview of recent results and discuss some specific systems of interest.

We present a spectroscopic analysis of the black hole binary Nova Muscae 1991 in quiescence using data obtained in 2009 with MagE on the Magellan Clay telescope and in 2010 with IMACS on the Magellan Baade telescope at the Las Campanas Observatory. Emission from the disc is observed in H alpha, H beta and Ca II (8662 A). A prominent hotspot is observed in the Doppler maps of all three emission lines. The existence of this spot establishes ongoing mass transfer from the donor star in 2009-2010 and, given its absence in the 1993-1995 observations, demonstrates the presence of a variable hotspot in the system. We find the radial distance to the hotspot from the black hole to be consistent with the circularization radius. Our tomograms are suggestive of stream-disc overflow in the system. We also detect possible Ca II (8662 A) absorption from the donor star.

We present the discovery of KIC 9632895b, a 6.2 Earth-radius planet in a low-eccentricity, 240.5-day orbit about an eclipsing binary. The binary itself consists of a 0.93 and 0.194 solar mass pair of stars with an orbital period of 27.3 days. The plane of the planet's orbit is rapidly precessing, and its inclination only becomes sufficiently aligned with the primary star in the latter portion of the Kepler data. Thus three transits are present in the latter half of the light curve, but none of the three conjunctions that occurred during the first half of the light curve produced transits. The precession period is ~103 years, and during that cycle, transits are visible only ~8% of the time. This has the important implication that for every system like KIC 9632895 that we detect, there are ~12 circumbinary systems that exist but are not currently exhibiting transits. The planet's mass is too small to noticeably perturb the binary, consequently its mass is not measurable with these data; but our photodynamical model places a 1-sigma upper limit of 16 Earth masses. With a period 8.8 times that of the binary, the planet is well outside the dynamical instability zone. It does, however, lie within the habitable zone of the binary, and making it the third of ten Kepler circumbinary planets to do so.

We analyze a large set of new and archival photometric and spectroscopic observations of LMC X-3 to arrive at a self-consistent dynamical model for the system. Using echelle spectra obtained with the MIKE instrument on the 6.5m Magellan Clay telescope and the UVES instrument on the second 8.2m Very Large Telescope we find a velocity semiamplitude for the secondary star of $K_2=241.1\pm 6.2$ km s$^{-1}$, where the uncertainty includes an estimate of the systematic error caused by X-ray heating. Using the spectra, we also find a projected rotational velocity of $V_{\rm rot}\sin i=118.5\pm 6.6$ km s$^{-1}$. From an analysis of archival $B$ and $V$ light curves as well as new $B$ and $V$ light curves from the SMARTS 1.3m telescope, we find an inclination of $i=69.84\pm 0.37^{\circ}$ for models that do not include X-ray heating and an inclination of $i=69.24\pm 0.72^{\circ}$ for models that incorporate X-ray heating. Adopting the latter inclination measurement, we find masses of $3.63\pm 0.57\,M_{\odot}$ and $6.98\pm 0.56\,M_{\odot}$ for the companion star and the black hole, respectively. We briefly compare our results with earlier work and discuss some of their implications.

In Gou et al. (2011), we reported that the black hole primary in the X-ray binary Cygnus X-1 is a near-extreme Kerr black hole with a spin parameter a*>0.95(3{\sigma}). We confirm this result while setting a new and more stringent limit: a*>0.983 at the 3{\sigma}(99.7%) level of confidence. The earlier work, which was based on an analysis of all three useful spectra that were then available, was possibly biased by the presence in these spectra of a relatively strong Compton power-law component: The fraction of the thermal seed photons scattered into the power law was f_s=23-31%, while the upper limit for reliable application of the continuum-fitting method is f_s<25%. We have subsequently obtained six additional spectra of Cygnus X-1 suitable for the measurement of spin. Five of these spectra are of high quality with f_s in the range 10% to 19%, a regime where the continuum-fitting method has been shown to deliver reliable results. Individually, the six spectra give lower limits on the spin parameter that range from a*>0.95 to a*>0.98, allowing us to conservatively conclude that the spin of the black hole is a*>0.983 (3{\sigma}).

Building upon a new dynamical model for the X-ray binary LMC X-3, we measure the spin of its black hole primary via the continuum-fitting method. We consider over one thousand thermal-state RXTE X-ray spectra of LMC X-3. Using a large subset of these spectra, we constrain the spin parameter of the black hole to be spin = 0.21(+0.18,-0.22), 90% confidence. Our estimate of the uncertainty in spin takes into account a wide range of systematic errors. We discuss evidence for a correlation between a black hole's spin and the complexity of its X-ray spectrum.

We report the discovery of a transiting, Rp = 4.347+/-0.099REarth, circumbinary planet (CBP) orbiting the Kepler K+M Eclipsing Binary (EB) system KIC 12351927 (Kepler-413) every ~66 days on an eccentric orbit with ap = 0.355+/-0.002AU, ep = 0.118+/-0.002. The two stars, with MA = 0.820+/-0.015MSun, RA = 0.776+/-0.009RSun and MB = 0.542+/-0.008MSun, RB = 0.484+/-0.024RSun respectively revolve around each other every 10.11615+/-0.00001 days on a nearly circular (eEB = 0.037+/-0.002) orbit. The orbital plane of the EB is slightly inclined to the line of sight (iEB = 87.33+/-0.06 degrees) while that of the planet is inclined by ~2.5 degrees to the binary plane at the reference epoch. Orbital precession with a period of ~11 years causes the inclination of the latter to the sky plane to continuously change. As a result, the planet often fails to transit the primary star at inferior conjunction, causing stretches of hundreds of days with no transits (corresponding to multiple planetary orbital periods). We predict that the next transit will not occur until 2020. The orbital configuration of the system places the planet slightly closer to its host stars than the inner edge of the extended habitable zone. Additionally, the orbital configuration of the system is such that the CBP may experience Cassini-States dynamics under the influence of the EB, in which the planet's obliquity precesses with a rate comparable to its orbital precession. Depending on the angular precession frequency of the CBP, it could potentially undergo obliquity fluctuations of dozens of degrees (and complex seasonal cycles) on precession timescales.

The X-ray persistence and characteristically soft spectrum of the black hole X-ray binary LMC X-3 make this source a touchstone for penetrating studies of accretion physics. We analyze a rich, 10-year collection of optical/infrared (OIR) time-series data in conjunction with all available contemporaneous X-ray data collected by the ASM and PCA detectors aboard the Rossi X-ray Timing Explorer. A cross-correlation analysis reveals an X-ray lag of ~2 weeks. Motivated by this result, we develop a model that reproduces the complex OIR light curves of LMC X-3. The model is comprised of three components of emission: stellar light; accretion luminosity from the outer disk inferred from the time-lagged X-ray emission; and light from the X-ray-heated star and outer disk. Using the model, we filter a strong noise component out of the ellipsoidal light curves and derive an improved orbital period for the system. Concerning accretion physics, we find that the local viscous timescale in the disk increases with the local mass accretion rate; this in turn implies that the viscosity parameter alpha decreases with increasing luminosity. Finally, we find that X-ray heating is a strong function of X-ray luminosity below ~50% of the Eddington limit, while above this limit X-ray heating is heavily suppressed. We ascribe this behavior to the strong dependence of the flaring in the disk upon X-ray luminosity, concluding that for luminosities above ~50% of Eddington, the star lies fully in the shadow of the disk.

We report on the masses, sizes, and orbits of the planets orbiting 22 Kepler stars. There are 49 planet candidates around these stars, including 42 detected through transits and 7 revealed by precise Doppler measurements of the host stars. Based on an analysis of the Kepler brightness measurements, along with high-resolution imaging and spectroscopy, Doppler spectroscopy, and (for 11 stars) asteroseismology, we establish low false-positive probabilities for all of the transiting planets (41 of 42 have a false-positive probability under 1%), and we constrain their sizes and masses. Most of the transiting planets are smaller than 3X the size of Earth. For 16 planets, the Doppler signal was securely detected, providing a direct measurement of the planet's mass. For the other 26 planets we provide either marginal mass measurements or upper limits to their masses and densities; in many cases we can rule out a rocky composition. We identify 6 planets with densities above 5 g/cc, suggesting a mostly rocky interior for them. Indeed, the only planets that are compatible with a purely rocky composition are smaller than ~2 R_earth. Larger planets evidently contain a larger fraction of low-density material (H, He, and H2O).

We examine ~10 years of photometric data and find that the black hole X-ray binary V4641 Sgr has two optical states, passive and active, during X-ray quiescence. The passive state is dominated by ellipsoidal variations and is stable in the shape and variability of the light curve. The active state is brighter and more variable. Emission during the active state varies over the course of the orbital period and is redder than the companion star. These optical/infrared states last for weeks or months. V4641 Sgr spends approximately 85% of X-ray quiescence in the passive state and 15% in the active. We analyze passive colors and spectroscopy of V4641 Sgr and show that they are consistent with a reddened B9III star (with E(B-V) = 0.37 +/- 0.19) with little or no contribution from the accretion disk. We use X-ray observations with an updated ephemeris to place an upper limit on the duration of an X-ray eclipse of <8.3 deg in phase (~1.6 hours). High resolution spectroscopy yields a greatly improved measurement of the rotational velocity of the companion star of V_rot_sin(i) = 100.9 +/- 0.8 km s^-1. We fit ellipsoidal models to the passive state data and find an inclination angle of i = 72.3 +/- 4.1 deg, a mass ratio of Q = 2.2 +/- 0.2, and component masses for the system of M_BH = 6.4 +/- 0.6 M_sun and M_2 = 2.9 +/- 0.4 M_sun. Using these values we calculate an updated distance to V4641 Sgr of 6.2 +/- 0.7 kpc.

We present a catalog of precise eclipse times and analysis of third body signals among 1279 close binaries in the latest Kepler Eclipsing Binary Catalog. For these short period binaries, Kepler's 30 minute exposure time causes significant smearing of light curves. In addition, common astrophysical phenomena such as chromospheric activity, as well as imperfections in the light curve detrending process, can create systematic artifacts that may produce fictitious signals in the eclipse timings. We present a method to measure precise eclipse times in the presence of distorted light curves, such as in contact and near-contact binaries which exhibit continuously changing light levels in and out of eclipse. 236 systems for which we find a timing variation signal compatible with the presence of a third body are identified. These are modeled for the light time travel effect and the basic properties of the third body are derived. This study complements Orosz et al. (2013; in prep), which focuses on eclipse timing variations of longer period binaries with flat out-of-eclipse regions. Together, these two papers provide comprehensive eclipse timings for all binaries in the Kepler Eclipsing Binary Catalog, as an ongoing resource freely accessible online to the community.

We show that the Kepler spacecraft in two-reaction wheel mode of operation is very well suited for the study of eclipsing binary star systems. Continued observations of the Kepler field will provide the most enduring and long-term valuable science. It will enable the discovery and characterization of eclipsing binaries with periods greater than 1 year - these are the most important, yet least understood binaries for habitable-zone planet background considerations. The continued mission will also enable the investigation of hierarchical multiple systems (discovered through eclipse timing variations), and provide drastically improved orbital parameters for circumbinary planetary systems.

Ranked near the top of the long list of exciting discoveries made with NASA's Kepler photometer is the detection of transiting circumbinary planets. In just over a year the number of such planets went from zero to seven, including a multi-planet system with one of the planets in the habitable zone (Kepler-47). We are quickly learning to better detect and characterize these planets, including the recognition of their transit timing and duration variation "smoking gun" signature. Even with only a handful of such planets, some exciting trends are emerging.

We analyze extensive BVR_cI_c time-series photometry and radial-velocity measurements for WOCS 40007 (Auner 259; KIC 5113053), a double-lined detached eclipsing binary and a member of the open cluster NGC 6819. Utilizing photometric observations from the 1-meter telescope at Mount Laguna Observatory and spectra from the WIYN 3.5-meter telescope, we measure precise and accurate masses (~1.6% uncertainty) and radii (~0.5%) for the binary components. In addition, we discover a third star orbiting the binary with a period greater than 3000 days using radial velocities and Kepler eclipse timings. Because the stars in the eclipsing binary are near the cluster turnoff, they are evolving rapidly in size and are sensitive to age. With a metallicity of [Fe/H]=+0.09+/-0.03, we find the age of NGC 6819 to be about 2.4 Gyr from CMD isochrone fitting and 3.1+/-0.4 Gyr by analyzing the mass-radius (M-R) data for this binary. The M-R age is above previous determinations for this cluster, but consistent within 1 sigma uncertainties. When the M-R data for the primary star of the additional cluster binary WOCS 23009 is included, the weighted age estimate drops to 2.5+/-0.2 Gyr, with a systematic uncertainty of at least 0.2 Gyr. The age difference between our CMD and M-R findings may be the result of systematic error in the metallicity or helium abundance used in models, or due to slight radius inflation of one or both stars in the WOCS 40007 binary.

We analyze extensive BVRI photometry and radial velocity measurements for three double-lined deeply-eclipsing binary stars in the field of the old open cluster NGC 7142. The short period (P = 1.9096825 d) detached binary V375 Cep is a high probability cluster member, and has a total eclipse of the secondary star. The characteristics of the primary star (M=1.288+/-0.017 Sun) at the cluster turnoff indicate an age of 3.6 Gyr (with a random uncertainty of 0.25 Gyr), consistent with earlier analysis of the color-magnitude diagram. The secondary star (M = 0.871+/-0.008 Msun) is not expected to have evolved significantly, but its radius is more than 10% larger than predicted by models. Because this binary system has a known age, it is useful for testing the idea that radius inflation can occur in short period binaries for stars with significant convective envelopes due to the inhibition of energy transport by magnetic fields. The brighter star in the binary also produces a precision estimate of the distance modulus, independent of reddening estimates: (m-M)_V=12.86+/-0.07. The other two eclipsing binary systems are not cluster members, although one of the systems (V2) could only be conclusively ruled out as a present or former member once the stellar characteristics were determined. That binary is within 0.5 degrees of edge-on, is in a fairly long-period eccentric binary, and contains two almost indistiguishable stars. The other binary (V1) has a small but nonzero eccentricity (e=0.038) in spite of having an orbital period under 5 d.

Eight years of optical photometry from OGLE-III are presented of the optical counterpart to the High Mass X-ray Binary system, SMC X-1. The optical data provide the best view to date of the modulation of the light from the system at the binary period of 3.9d. In addition, it is shown for the first time that the light is also modulated at the superorbital period - a period associated with the precessing of a warped accretion disk around the neutron star partner. The implications for the sources of the various components of the optical light in this system are discussed.

We present the detection of five planets -- Kepler-62b, c, d, e, and f -- of size 1.31, 0.54, 1.95, 1.61 and 1.41 Earth radii, orbiting a K2V star at periods of 5.7, 12.4, 18.2, 122.4 and 267.3 days, respectively. The outermost planets (Kepler-62e & -62f) are super-Earth-size (1.25 < planet radius/earth radius < 2.0) planets in the habitable zone (HZ) of their host star, receiving 1.2 +- 0.2 and 0.41 +- 0.05 times the solar flux at Earth's orbit. Theoretical models of Kepler-62e and -62f for a stellar age of ~7 Gyr suggest that both planets could be solid: either with a rocky composition or composed of mostly solid water in their bulk.

We report the discovery and confirmation of a transiting circumbinary planet (PH1b) around KIC 4862625, an eclipsing binary in the Kepler field. The planet was discovered by volunteers searching the first six Quarters of publicly available Kepler data as part of the Planet Hunters citizen science project. Transits of the planet across the larger and brighter of the eclipsing stars are detectable by visual inspection every ~137 days, with seven transits identified in Quarters 1-11. The physical and orbital parameters of both the host stars and planet were obtained via a photometric-dynamical model, simultaneously fitting both the measured radial velocities and the Kepler light curve of KIC 4862625. The 6.18 +/- 0.17 Earth radii planet orbits outside the 20-day orbit of an eclipsing binary consisting of an F dwarf (1.734 +/- 0.044 Solar radii, 1.528 +/- 0.087 Solar masses) and M dwarf (0.378+/- 0.023 Solar radii, 0.408 +/- 0.024 Solar masses). For the planet, we find an upper mass limit of 169 Earth masses (0.531 Jupiter masses) at the 99.7% confidence level. With a radius and mass less than that of Jupiter, PH1b is well within the planetary regime. Outside the planet's orbit, at ~1000 AU,a previously unknown visual binary has been identified that is likely bound to the planetary system, making this the first known case of a quadruple star system with a transiting planet.

We present the discovery of the totally eclipsing long-period (P = 771.8 d) binary system WOCS 23009 in the old open cluster NGC 6819 that contains both an evolved star near central hydrogen exhaustion and a low-mass (0.45 Msun) star. This system was previously known to be a single-lined spectroscopic binary, but the discovery of an eclipse near apastron using data from the Kepler space telescope makes it clear that the system has an inclination that is very close to 90 degrees. Although the secondary star has not been identified in spectra, the mass of the primary star can be constrained using other eclipsing binaries in the cluster. The combination of total eclipses and a mass constraint for the primary star allows us to determine a reliable mass for the secondary star and radii for both stars, and to constrain the cluster age. Unlike well-measured stars of similar mass in field binaries, the low-mass secondary is not significantly inflated in radius compared to model predictions. The primary star characteristics, in combination with cluster photometry and masses from other cluster binaries, indicates a best age of 2.62+/-0.25 Gyr, although stellar model physics may introduce systematic uncertainties at the ~10% level. We find preliminary evidence that the asteroseismic predictions for red giant masses in this cluster are systematically too high by as much as 8%.

KIC 6131659 is a long-period (17.5 days) eclipsing binary discovered by the Kepler mission. We analyzed six quarters of Kepler data along with supporting ground-based photometric and spectroscopic data to obtain accurate values for the mass and radius of both stars, namely M_1=0.922 +/- 0.007 M_sun, R_1=0.8800 +/- 0.0028 R_sun, and M_2=0.685 +/- 0.005 M_sun, R_2=0.6395 +/- 0.0061 R_sun. There is a well-known issue with low mass (M <<0.8 M_sun) stars (in cases where the mass and radius measurement uncertainties are smaller than two or three percent) where the measured radii are almost always 5 to 15 percent larger than expected from evolutionary models, i.e. the measured radii are all above the model isochrones in a mass-radius plane. In contrast, the two stars in KIC 6131659 were found to sit on the same theoretical isochrone in the mass-radius plane. Until recently, all of the well-studied eclipsing binaries with low-mass stars had periods less than about three days. The stars in such systems may have been inflated by high levels of stellar activity induced by tidal effects in these close binaries. KIC 6131659 shows essentially no evidence of enhanced stellar activity, and our measurements support the hypothesis that the unusual mass-radius relationship observed in most low-mass stars is influenced by strong magnetic activity created by the rapid rotation of the stars in tidally-locked, short-period systems. Finally, using short cadence data, we show that KIC 6131657 has one of the smallest measured non-zero eccentricities of a binary with two main sequence stars, where e cos omega=(4.57 +/-0.02)*10^-5.

We report the detection of Kepler-47, a system consisting of two planets orbiting around an eclipsing pair of stars. The inner and outer planets have radii 3.0 and 4.6 times that of the Earth, respectively. The binary star consists of a Sun-like star and a companion roughly one-third its size, orbiting each other every 7.45 days. With an orbital period of 49.5 days, eighteen transits of the inner planet have been observed, allowing a detailed characterization of its orbit and those of the stars. The outer planet's orbital period is 303.2 days, and although the planet is not Earth-like, it resides within the classical "habitable zone", where liquid water could exist on an Earth-like planet. With its two known planets, Kepler-47 establishes that close binary stars can host complete planetary systems.

We discuss the discovery and characterization of the circumbinary planet Kepler-38b. The stellar binary is single-lined, with a period of 18.8 days, and consists of a moderately evolved main-sequence star (M_A = 0.949 +/- 0.059 solar masses and R_A = 1.757 +/- 0.034 solar radii) paired with a low-mass star (M_B = 0.249 +/- 0.010 solar masses and R_B = 0.2724 +/- 0.0053 solar radii) in a mildly eccentric (e=0.103) orbit. A total of eight transits due to a circumbinary planet crossing the primary star were identified in the Kepler light curve (using Kepler Quarters 1 through 11), from which a planetary period of 105.595 +/- 0.053 days can be established. A photometric dynamical model fit to the radial velocity curve and Kepler light curve yields a planetary radius of 4.35 +/- 0.11 Earth radii, or 1.12 +/- 0.03 Neptune radii. Since the planet is not sufficiently massive to observably alter the orbit of the binary from Keplerian motion, we can only place an upper limit on the mass of the planet of 122 Earth masses (7.11 Neptune masses or 0.384 Jupiter masses) at 95% confidence. This upper limit should decrease as more Kepler data become available.

Most Sun-like stars in the Galaxy reside in gravitationally-bound pairs of stars called "binary stars". While long anticipated, the existence of a "circumbinary planet" orbiting such a pair of normal stars was not definitively established until the discovery of Kepler-16. Incontrovertible evidence was provided by the miniature eclipses ("transits") of the stars by the planet. However, questions remain about the prevalence of circumbinary planets and their range of orbital and physical properties. Here we present two additional transiting circumbinary planets, Kepler-34 and Kepler-35. Each is a low-density gas giant planet on an orbit closely aligned with that of its parent stars. Kepler-34 orbits two Sun-like stars every 289 days, while Kepler-35 orbits a pair of smaller stars (89% and 81% of the Sun's mass) every 131 days. Due to the orbital motion of the stars, the planets experience large multi-periodic variations in incident stellar radiation. The observed rate of circumbinary planets implies > ~1% of close binary stars have giant planets in nearly coplanar orbits, yielding a Galactic population of at least several million.

We present an automated classification of 2165 \textit{Kepler} eclipsing binary (EB) light curves that accompanied the second \textit{Kepler} data release. The light curves are classified using Locally Linear Embedding, a general nonlinear dimensionality reduction tool, into morphology types (detached, semi-detached, overcontact, ellipsoidal). The method, related to a more widely used Principal Component Analysis, produces a lower-dimensional representation of the input data while preserving local geometry and, consequently, the similarity between neighboring data points. We use this property to reduce the dimensionality in a series of steps to a one-dimensional manifold and classify light curves with a single parameter that is a measure of "detachedness" of the system. This fully automated classification correlates well with the manual determination of morphology from the data release, and also efficiently highlights any misclassified objects. Once a lower-dimensional projection space is defined, the classification of additional light curves runs in a negligible time and the method can therefore be used as a fully automated classifier in pipeline structures. The classifier forms a tier of the \textit{Kepler} EB pipeline that pre-processes light curves for the artificial intelligence based parameter estimator.

We present the first results of a Kepler survey of 41 eclipsing binaries that we undertook to search for third star companions. Such tertiaries will periodically alter the eclipse timings through light travel time and dynamical effects. We discuss the prevalence of starspots and pulsation among these binaries and how these phenomena influence the eclipse times. There is no evidence of short period companions (P < 700 d) among this sample, but we do find evidence for long term timing variations in 14 targets (34%). We argue that this finding is consistent with the presence of tertiary companions among a significant fraction of the targets, especially if many have orbits measured in decades. This result supports the idea that the formation of close binaries involves the deposition of angular momentum into the orbital motion of a third star.

We present photometric and spectroscopic observations of a peculiar variable (designated DASCH J075731.1+201735 or J0757) discovered from our DASCH project using the digitized Harvard College Observatory archival photographic plates. It brightened by about 1.5 magnitudes in B within a year starting in 1942, and then slowly faded back to its pre-outburst brightness from 1943 to the 1950s. The mean brightness level was stable before and after the outburst, and ellipsoidal variations with a period of $P=119.18\pm0.07$ days are seen, suggesting that the star is tidally distorted. Radial-velocity measurements indicate that the orbit is nearly circular ($e=0.02\pm0.01$) with a spectroscopic period that is the same as the photometric period. The binary consists of a $1.1\pm0.3 M_\odot$ M0III star, and a $0.6\pm0.2 M_\odot$ companion, very likely a white dwarf (WD). Unlike other symbiotic binaries, there is no sign of emission lines or a stellar wind in the spectra. With an outburst timescale of ~10 years and estimated B band peak luminosity M_B~0.7, J0757 is different from any other known classic or symbiotic novae. The most probable explanation of the outburst is Hydrogen shell-burning on the WD, although an accretion-powered flare cannot be ruled out.

Dynamical mass measurements hold the key to answering whether ultraluminous X-ray sources (ULXs) are intermediate mass black holes (IMBHs) or stellar mass black holes with special radiation mechanisms. NGC 1313 X-2 is so far the only ULX with HST light curves, the orbital period, and the black hole's radial velocity amplitude based on the He II $\lambda4686$\AA\ disk emission line shift of $\sim200$ km/s. We constrain its black hole mass and other parameters by fitting observations to a binary light curve code with accommodations for X-ray heating of the accretion disk and the secondary. Given the dynamical constraints from the observed light curves and the black hole radial motion and the observed stellar environment age, the only acceptable models are those with 40-50 Myrs old intermediate mass secondaries in their helium core and hydrogen shell burning phase filling 40%-80% of their Roche lobes. The black hole can be a massive black hole of a few tens of $M_\odot$ that can be produced from stellar evolution of low metalicity stars, or an IMBH of a few hundred to above 1000$M_\odot$ if its true radial velocity $2K^\prime<40$ km/s. Further observations are required to better measure the black hole radial motion and the light curves in order to determine whether NGC 1313 X-2 is a stellar black hole or an IMBH.

We report a direct and accurate measurement of the distance to the X-ray binary Cygnus X-1, which contains the first black hole to be discovered. The distance of 1.86 (-0.11,+0.12) kpc was obtained from a trigonometric parallax measurement using the Very Long Baseline Array. The position measurements are also sensitive to the 5.6 d binary orbit and we determine the orbit to be clockwise on the sky. We also measured the proper motion of Cygnus X-1 which, when coupled to the distance and Doppler shift, gives the three-dimensional space motion of the system. When corrected for differential Galactic rotation, the non-circular (peculiar) motion of the binary is only about 21 km/s, indicating that the binary did not experience a large "kick" at formation.

Cygnus X-1 is a binary star system that is comprised of a black hole and a massive giant companion star in a tight orbit. Building on our accurate distance measurement reported in the preceding paper, we first determine the radius of the companion star, thereby constraining the scale of the binary system. To obtain a full dynamical model of the binary, we use an extensive collection of optical photometric and spectroscopic data taken from the literature. By using all of the available observational constraints, we show that the orbit is slightly eccentric (both the radial velocity and photometric data independently confirm this result) and that the companion star rotates roughly 1.4 times its pseudosynchronous value. We find a black hole mass of M =14.8\pm1.0 M_{\sun}, a companion mass of M_{opt}=19.2\pm1.9 M_{\sun}, and the angle of inclination of the orbital plane to our line of sight of i=27.1\pm0.8 deg.

The compact primary in the X-ray binary Cygnus X-1 was the first black hole to be established via dynamical observations. We have recently determined accurate values for its mass and distance, and for the orbital inclination angle of the binary. Building on these results, which are based on our favored (asynchronous) dynamical model, we have measured the radius of the inner edge of the black hole's accretion disk by fitting its thermal continuum spectrum to a fully relativistic model of a thin accretion disk. Assuming that the spin axis of the black hole is aligned with the orbital angular momentum vector, we have determined that Cygnus X-1 contains a near-extreme Kerr black hole with a spin parameter a/M>0.95 (3\sigma). For a less probable (synchronous) dynamical model, we find a/M>0.92 (3\sigma). In our analysis, we include the uncertainties in black hole mass, orbital inclination angle and distance, and we also include the uncertainty in the calibration of the absolute flux via the Crab. These four sources of uncertainty totally dominate the error budget. The uncertainties introduced by the thin-disk model we employ are particularly small in this case given the extreme spin of the black hole and the disk's low luminosity.

Kepler-16 is an eccentric low-mass eclipsing binary with a circumbinary transiting planet. Here we investigate the angular momentum of the primary star, based on Kepler photometry and Keck spectroscopy. The primary star's rotation period is 35.1 +/- 1.0 days, and its projected obliquity with respect to the stellar binary orbit is 1.6 +/- 2.4 degrees. Therefore the three largest sources of angular momentum---the stellar orbit, the planetary orbit, and the primary's rotation---are all closely aligned. This finding supports a formation scenario involving accretion from a single disk. Alternatively, tides may have realigned the stars despite their relatively wide separation (0.2 AU), a hypothesis that is supported by the agreement between the measured rotation period and the "pseudosynchronous" period of tidal evolution theory. The rotation period, chromospheric activity level, and fractional light variations suggest a main-sequence age of 2-4 Gyr. Evolutionary models of low-mass stars can match the observed masses and radii of the primary and secondary stars to within about 3%.

We report the detection of a planet whose orbit surrounds a pair of low-mass stars. Data from the Kepler spacecraft reveal transits of the planet across both stars, in addition to the mutual eclipses of the stars, giving precise constraints on the absolute dimensions of all three bodies. The planet is comparable to Saturn in mass and size, and is on a nearly circular 229-day orbit around its two parent stars. The eclipsing stars are 20% and 69% as massive as the sun, and have an eccentric 41-day orbit. The motions of all three bodies are confined to within 0.5 degree of a single plane, suggesting that the planet formed within a circumbinary disk.

Kepler observations of the star HD 187091 (KID 8112039, hereafter KOI-54) revealed a remarkable light curve exhibiting sharp periodic brightening events every 41.8 days with a superimposed set of oscillations forming a beating pattern in phase with the brightenings. Spectroscopic observations revealed that this is a binary star with a highly eccentric orbit, e=0.83. We are able to match the Kepler light curve and radial velocities with a nearly face-on (i=5.5 degree) binary star model in which the brightening events are caused by tidal distortion and irradiation of nearly identical A stars during their close periastron passage. The two dominant oscillations in the light curve, responsible for the beating pattern, have frequencies that are the 91st and 90th harmonic of the orbital frequency. The power spectrum of the light curve, after removing the binary star brightening component, reveals a large number of pulsations, 30 of which have a signal-to-noise ratio > 7. Nearly all of these pulsations have frequencies that are either integer multiples of the orbital frequency or are tidally-split multiples of the orbital frequency. This pattern of frequencies unambiguously establishes the pulsations as resonances between the dynamic tides at periastron and the free oscillation modes of one of the stars. KOI-54 is only the 4th star to show such a phenomenon, and is by far the richest in terms of excited modes.

The Kepler Mission provides nearly continuous monitoring of ~156 000 objects with unprecedented photometric precision. Coincident with the first data release, we presented a catalog of 1879 eclipsing binary systems identified within the 115 square degree Kepler FOV. Here, we provide an updated catalog augmented with the second Kepler data release which increases the baseline nearly 4-fold to 125 days. 386 new systems have been added, ephemerides and principle parameters have been recomputed. We have removed 42 previously cataloged systems that are now clearly recognized as short-period pulsating variables and another 58 blended systems where we have determined that the Kepler target object is not itself the eclipsing binary. A number of interesting objects are identified. We present several exemplary cases: 4 EBs that exhibit extra (tertiary) eclipse events; and 8 systems that show clear eclipse timing variations indicative of the presence of additional bodies bound in the system. We have updated the period and galactic latitude distribution diagrams. With these changes, the total number of identified eclipsing binary systems in the Kepler field-of-view has increased to 2165, 1.4% of the Kepler target stars.

The Kepler spacecraft has been monitoring the light from 150,000 stars in its primary quest to detect transiting exoplanets. Here we report on the detection of an eclipsing stellar hierarchical triple, identified in the Kepler photometry. KOI-126 (A,(B, C)), is composed of a low-mass binary (masses M_B = 0.2413+/-0.0030 M_Sun, M_C = 0.2127+/-0.0026 M_Sun; radii R_B = 0.2543+/-0.0014 R_Sun, R_C = 0.2318+/-0.0013 R_Sun; orbital period P_1 = 1.76713+/-0.00019 days) on an eccentric orbit about a third star (mass M_A = 1.347+/-0.032 M_Sun; radius R_A = 2.0254+/-0.0098 R_Sun; period of orbit around the low-mass binary P_2 = 33.9214+/-0.0013 days; eccentricity of that orbit e_2 = 0.3043+/-0.0024). The low-mass pair probe the poorly sampled fully-convective stellar domain offering a crucial benchmark for theoretical stellar models.

When an extrasolar planet passes in front of its star (transits), its radius can be measured from the decrease in starlight and its orbital period from the time between transits. Multiple planets transiting the same star reveal more: period ratios determine stability and dynamics, mutual gravitational interactions reflect planet masses and orbital shapes, and the fraction of transiting planets observed as multiples has implications for the planarity of planetary systems. But few stars have more than one known transiting planet, and none has more than three. Here we report Kepler spacecraft observations of a single Sun-like star that reveal six transiting planets, five with orbital periods between 10 and 47 days plus a sixth one with a longer period. The five inner planets are among the smallest whose masses and sizes have both been measured, and these measurements imply substantial envelopes of light gases. The degree of coplanarity and proximity of the planetary orbits imply energy dissipation near the end of planet formation.

A typical galaxy is thought to contain tens of millions of stellar-mass black holes, the collapsed remnants of once massive stars, and a single nuclear supermassive black hole. Both classes of black holes accrete gas from their environments. The accreting gas forms a flattened orbiting structure known as an accretion disk. During the past several years, it has become possible to obtain measurements of the spins of the two classes of black holes by modeling the X-ray emission from their accretion disks. Two methods are employed, both of which depend upon identifying the inner radius of the accretion disk with the innermost stable circular orbit (ISCO), whose radius depends only on the mass and spin of the black hole. In the Fe K method, which applies to both classes of black holes, one models the profile of the relativistically-broadened iron line with a special focus on the gravitationally redshifted red wing of the line. In the continuum-fitting method, which has so far only been applied to stellar-mass black holes, one models the thermal X-ray continuum spectrum of the accretion disk. We discuss both methods, with a strong emphasis on the continuum-fitting method and its application to stellar-mass black holes. Spin results for eight stellar-mass black holes are summarized. These data are used to argue that the high spins of at least some of these black holes are natal, and that the presence or absence of relativistic jets in accreting black holes is not entirely determined by the spin of the black hole.

The Kepler space mission is devoted to finding Earth-size planets in habitable zones orbiting other stars. Its large, 105-deg field-of-view features over 156,000 stars that are observed continuously to detect and characterize planet transits. Yet this high-precision instrument holds great promise for other types of objects as well. Here we present a comprehensive catalog of eclipsing binary stars observed by Kepler in the first 44 days of operation, the data which are publicly available through MAST as of 6/15/2010. The catalog contains 1879 unique objects. For each object we provide its Kepler ID (KID), ephemeris (BJD0, P0), morphology type, physical parameters (Teff, log g, E(B-V), crowding), and principal parameters (T2/T1, q, fillout factor and sin i for overcontacts, and T2/T1, (R1+R2)/a, e sin(w), e cos(w), and sin i for detached binaries). We present statistics based on the determined periods and measure an average occurence rate of eclipsing binaries to be ~1.2% across the Kepler field. We further discuss the distribution of binaries as function of galactic latitude, and thoroughly explain the application of artificial intelligence to obtain principal parameters in a matter of seconds for the whole sample. The catalog was envisioned to serve as a bridge between the now public Kepler data and the scientific community interested in eclipsing binary stars.

We present an improved dynamical model of the X-ray binary and microquasar XTE J1550-564 based on new moderate-resolution optical spectroscopy and near-infrared photometry. By combining our new radial velocity measurements with previous measurements obtained 2001 May at the 8.2m VLT and with light curves, we find an orbital period of P=1.5420333 +/- 0.0000024 days and a radial velocity semiamplitude of K_2=363.14 +/- 5.97$ km/sec, which together imply an optical mass function of f(M)=7.65 +/- 0.38 solar masses. We find that the projected rotational velocity of the secondary star is 55 +/- 5 km/sec, which implies a very extreme mass ratio of Q=M/M_2=30. Using a model of a Roche lobe-filling star and an azimuthally symmetric accretion disk, we fit simultaneously optical light curves from 2001, near-infrared light curves from 2008 and all of the radial velocity measurements to derive system parameters. We find an inclination of 74.7 +/- 3.8 deg and component masses of M_2=0.30 +/- 0.07 solar masses and M=9.10 +/- 0.61 solar masses for the secondary star and black hole, respectively. The radius of the secondary star for the adopted model is 1.75 +/- 0.12 solar radii. Using this radius, the average K_S magnitude, and an extinction of A_K=0.507 +/- 0.050 mag, we find a distance of 4.38^{+0.58}_{-0.41} kpc, which is in good agreement with a recent distance estimate based on HI absorption lines (abstract shortened).

We present an improved method for determining the mass of neutron stars in eclipsing X-ray pulsar binaries and apply the method to six systems, namely Vela X-1, 4U 1538-52, SMC X-1, LMC X-4, Cen X-3, and Her X-1. In previous studies to determine neutron star mass, the X-ray eclipse duration has been approximated analytically by assuming the companion star is spherical with an effective Roche lobe radius. We use a numerical code based on Roche geometry with various optimizers to analyze the published data for these systems, which we supplement with new spectroscopic and photometric data for 4U 1538-52. This allows us to model the eclipse duration more accurately and thus calculate an improved value for the neutron star mass. The derived neutron star mass also depends on the assumed Roche lobe filling factor beta of the companion star, where beta = 1 indicates a completely filled Roche lobe. In previous work a range of beta between 0.9 and 1.0 was usually adopted. We use optical ellipsoidal lightcurve data to constrain beta. We find neutron star masses of 1.77 +/- 0.08 M_{sun} for Vela X-1, 0.87 +/- 0.07 M_{sun} for 4U 1538-52 (eccentric orbit), 1.00 +/- 0.10 M_{sun} for 4U 1538-52 (circular orbit), 1.04 +/- 0.09 M_{sun} for SMC X-1, 1.29 +/- 0.05 M_{sun} for LMC X-4, 1.49 +/- 0.08 M_{sun} for Cen X-3, and 1.07 +/- 0.36 M_{sun} for Her X-1. We discuss the limits of the approximations that were used to derive the earlier mass determinations, and we comment on the implications our new masses have for observationally refining the upper and lower bounds of the neutron star mass distribution.

M33 X-7 is among the most massive X-Ray binary stellar systems known, hosting a rapidly spinning 15.65 Msun black hole orbiting an underluminous 70 Msun Main Sequence companion in a slightly eccentric 3.45 day orbit. Although post-main-sequence mass transfer explains the masses and tight orbit, it leaves unexplained the observed X-Ray luminosity, star's underluminosity, black hole's spin, and eccentricity. A common envelope phase, or rotational mixing, could explain the orbit, but the former would lead to a merger and the latter to an overluminous companion. A merger would also ensue if mass transfer to the black hole were invoked for its spin-up. Here we report that, if M33 X-7 started as a primary of 85-99 Msun and a secondary of 28-32 Msun, in a 2.8-3.1 day orbit, its observed properties can be consistently explained. In this model, the Main Sequence primary transferred part of its envelope to the secondary and lost the rest in a wind; it ended its life as a ~16 Msun He star with a Fe-Ni core which collapsed to a black hole (with or without an accompanying supernova). The release of binding energy and, possibly, collapse asymmetries "kicked" the nascent black hole into an eccentric orbit. Wind accretion explains the X-Ray luminosity, while the black hole spin can be natal.

Black hole (BH) X-ray binaries (XRBs) are X-ray luminous binary systems comprising a BH accreting matter from a companion star. Understanding their origins sheds light on the still not well understood physics of BH formation. M33 X-7 hosts one of the most massive stellar-mass BH among all XRBs known to date, a 15.65 Msun BH orbiting a 70 Msun companion star in a 3.45 day orbit. The high masses of the two components and the tight orbit relative to the large H-rich stellar component challenge our understanding of the typically invoked BH-XRBs formation channels. The measured underluminosity of the optical component further complicates the picture. A solution to the evolutionary history of this system that can account for all its observed properties has yet to be presented, and here we propose the first scenario that is consistent with the complete set of current observational constraints. In our model, M33 X-7 started its life hosting a 85-99 Msun primary and a 28-32 Msun companion in a Keplerian orbit of 2.8-3.1 days. In order to form a BH of 15.65 Msun, the initially most massive component transferred part of its envelope to the companion star and lost the rest in a strong stellar wind. During this dynamically stable mass transfer phase the companion accreted matter, to become the presently underluminous 70 Msun star.

We measure the spin of XTE J1550-564 in two ways: by modelling the thermal continuum spectrum of the accretion disc, and independently by modeling the broad red wing of the reflection fluorescence Fe-K line. We find that the spin measurements conducted independently using both leading methods are in agreement with one another. For the continuum-fitting analysis, we use a data sample consisting of several dozen RXTE spectra, and for the Fe-K analysis, we use a pair of ASCA spectra from a single epoch. Our spin estimate for the black hole primary using the continuum-fitting method is -0.11 < a* < 0.71 (90 per cent confidence), with a most likely spin of a* = 0.34. In obtaining this result, we have thoroughly explored the dependence of the spin value on a wide range of model-dependent systematic errors and observational errors; our precision is limited by uncertainties in the distance and orbital inclination of the system. For the Fe-line method, our estimate of spin is a* = 0.55(+0.15,-0.22). Combining these results, we conclude that the spin of this black hole is moderate, a* = 0.49(+0.13,-0.20), which suggests that the jet of this microquasar is powered largely by its accretion disc rather than by the spin energy of the black hole.

Based on high-resolution ultraviolet spectroscopy obtained with FUSE and COS, we present new detections of O VI and N V emission from the black-hole X-ray binary (XRB) system LMC X-3. We also update the ephemeris of the XRB using recent radial velocity measurements obtained with the echelle spectrograph on the Magellan-Clay telescope. We observe significant velocity variability of the UV emission, and we find that the O VI and N V emission velocities follow the optical velocity curve of the XRB. Moreover, the O VI and N V intensities regularly decrease between binary phase = 0.5 and 1.0, which suggests that the source of the UV emission is increasingly occulted as the B star in the XRB moves from superior to inferior conjunction. These trends suggest that illumination of the B-star atmosphere by the intense X-ray emission from the accreting black hole creates a hot spot on one side of the B star, and this hot spot is the origin of the O VI and N V emission. However, the velocity semiamplitude of the ultraviolet emission, K_{UV}~180 km/s, is lower than the optical semiamplitude; this difference could be due to rotation of the B star. If our hypothesis about the origin of the highly ionized emission is correct, then careful analysis of the emission occultation could, in principle, constrain the inclination of the XRB and the mass of the black hole.

During its yearlong outburst in 1975--76, the transient source A0620--00 reached an intensity of 50 Crab, an all-time record for any X-ray binary. The source has been quiescent since. We recently determined accurate values for the black hole mass, orbital inclination angle and distance. Building on these results, we have measured the radius of the inner edge of the accretion disk around the black hole primary by fitting its thermal continuum spectrum to our version of the relativistic Novikov-Thorne thin-disk model. We have thereby estimated the spin of the black hole. Although our spin estimate depends on a single high quality spectrum, which was obtained in 1975 by OSO-8, we are confident of our result because of the consistent values of the inner-disk radius that we have obtained for hundreds of observations of other sources: H1743-322, XTE J1550-564, and notably LMC X-3. We have determined the dimensionless spin parameter of the black hole to be a*=0.12+/-0.19, with a*<0.49 and a*>-0.59 at the 3 sigma level of confidence. This result takes into account all sources of observational and model-parameter uncertainties. Despite the low spin, the intensity and properties of the radio counterpart, both in outburst and quiescence, attest to the presence of a strong jet. If jets are driven by black hole spin, then current models indicate that jet power should be a steeply increasing function of a*. Consequently, the low spin of A0620--00 suggests that its jet may be disk-driven.

In prior work, {\it Chandra} and Gemini-North observations of the eclipsing X-ray binary M33 X-7 have yielded measurements of the mass of its black hole primary and the system's orbital inclination angle of unprecedented accuracy. Likewise, the distance to the binary is known to a few percent. In an analysis based on these precise results, fifteen {\it Chandra} and {\it XMM-Newton} X-ray spectra, and our fully relativistic accretion disk model, we find that the dimensionless spin parameter of the black hole primary is $a_* = 0.77 \pm 0.05$. The quoted 1-$\sigma$ error includes all sources of observational uncertainty. Four {\it Chandra} spectra of the highest quality, which were obtained over a span of several years, all lead to the same estimate of spin to within statistical errors (2%), and this estimate is confirmed by 11 spectra of lower quality. There are two remaining uncertainties: (1) the validity of the relativistic model used to analyze the observations, which is being addressed in ongoing theoretical work; and (2) our assumption that the black hole spin is approximately aligned with the angular momentum vector of the binary, which can be addressed by a future X-ray polarimetry mission.

GU Boo is one of only a relatively small number of well studied double-lined eclipsing binaries that contain low-mass stars. L'opez-Morales & Ribas (2005, hereafter LR05) present a comprehensive analysis of multi-color light and radial velocity curves for this system. The GU Boo light curves presented in LR05 had substantial asymmetries, which were attributed to large spots. In spite of the asymmetry LR05 derived masses and radii accurate to about 2%. We obtained additional photometry of GU Boo using both a CCD and a single-channel photometer and modeled the light curves with the ELC software to determine if the large spots in the light curves give rise to systematic errors at the few percent level. We also modeled the original light curves from LR05 using models with and without spots. We derived a radius of the primary of 0.6329 +- 0.0026 Solar-Radii, 0.6413 +- 0.0049 SR, and 0.6373 +- 0.0029 SR from the CCD, photoelectric, and LR05 data, respectively. Each of these measurements agrees with the value reported by LR05 (R1=0.623 +- 0.016 SR) at the level of about 2%. In addition, the spread in these values is about 1-2% from the mean. For the secondary, we derive radii of 0.6074 +- 0.0035 SR, 0.5944 +- 0.0069 SR, and 0.5976 +- 0.0059 SR from the three respective data sets. The LR05 value is R2=0.620 +- 0.020 SR, which is about 2-3% larger than each of the three values we found. The spread in these values is about 2% from the mean. The systematic difference between our three determinations of the secondary radius and that of LR05 might be attributed to differences in the modelling process and codes used. Our own fits suggest that, for GU Boo at least, using accurate spot modelling of a single set of multicolor light curves results in radii determinations accurate at the 2% level.

We present an analysis of the early Kepler observations of the previously discovered transiting planet HAT-P-7b. The light curve shows the transit of the star, the occultation of the planet, and the orbit phase-dependent light from the planet. In addition, phase-dependent light from the star is present, known as "ellipsoidal variations". The very nearby planet (only 4 stellar radii away) gravitationally distorts the star and results in a flux modulation twice per orbit. The ellipsoidal variations can confuse interpretation of the planetary phase curve if not self-consistently included in the modeling. We fit the light curve using the Roche potential approximation and derive improved planet and orbit parameters.

We performed an optical/infrared study of the counterpart of the low-mass X-ray binary KS1731-260 to test its identification and obtain information about the donor. Optical and infrared images of the counterpart of KS1731-260 were taken in two different epochs (2001 and 2007) after the source returned to quiescence in X-rays. We compared those observations with obtained when KS 1731-260 was still active. We confirm the identification of KS1731-260 with the previously proposed counterpart and improve its position to RA=17:34:13.46 and DEC=-26:05:18.60. The H-band magnitude of this candidate showed a decline of ~1.7 mags from outburst to quiescence. In 2007 April we obtained R=22.8+-0.1 and I=20.9+-0.1 for KS1731-260. Similar optical brightness was measured in June 2001 and July 2007. The intrinsic optical color R-I is consistent with spectral types from F to G for the secondary although there is a large excess over that from the secondary at the infrared wavelengths. This may be due to emission from the cooler outer regions of the accretion disk. We cannot rule out a brown dwarf as a donor star, although it would require that the distance to the source is significantly lower than the 7 kpc reported by Muno et al. 2000.

We analyze photometry of the Soft X-ray Transient A0620-00 spanning nearly 30 years, including previously published and previously unpublished data. Previous attempts to determine the inclination of A0620 using subsets of these data have yielded a wide range of measured values of i. Differences in the measured value of i have been due to changes in the shape of the light curve and uncertainty regarding the contamination from the disk. We give a new technique for estimating the disk fraction and find that disk light is significant in all light curves, even in the infrared. We also find that all changes in the shape and normalization of the light curve originate in a variable disk component. After accounting for this disk component, we find that all the data, including light curves of significantly different shapes, point to a consistent value of i. Combining results from many separate data sets, we find i=51 plus or minus 0.9 degrees, implying M=6.6 plus or minus 0.25 solar masses. Using our dynamical model and zero-disk stellar VIH magnitudes, we find d=1.06 plus or minus 0.12 kpc. Understanding the disk origin of non-ellipsoidal variability may assist with making reliable determinations of i in other systems, and the fluctuations in disk light may provide a new observational tool for understanding the three-dimensional structure of the accretion disk.

The first extragalactic X-ray binary, LMC X-1, was discovered in 1969. In the 1980s, its compact primary was established as the fourth dynamical black-hole candidate. Recently, we published accurate values for the mass of the black hole and the orbital inclination angle of the binary system. Building on these results, we have analyzed 53 X-ray spectra obtained by RXTE and, using a selected sample of 18 of these spectra, we have determined the dimensionless spin parameter of the black hole to be a* = 0.92(-0.07,+0.05). This result takes into account all sources of observational and model-parameter uncertainties. The standard deviation around the mean value of a* for these 18 X-ray spectra, which were obtained over a span of several years, is only 0.02. When we consider our complete sample of 53 RXTE spectra, we find a somewhat higher value of the spin parameter and a larger standard deviation. Finally, we show that our results based on RXTE data are confirmed by our analyses of selected X-ray spectra obtained by the XMM-Newton, BeppoSAX and Ginga missions.

The bright X-ray transient H 1743-322 was observed daily by the Rossi X-ray Timing Explorer (RXTE) during most of its 8-month outburst in 2003. We present a detailed spectral analysis and a supporting timing analysis of all of these data, and we discuss the behavior and evolution of the source in terms of the three principal X-ray states defined by Remillard and McClintock. These X-ray results are complemented by Very Large Array (VLA) data obtained at six frequencies that provide quite complete coverage of the entire outburst cycle at 4.860 GHz and 8.460 GHz. We also present photometric data and finding charts for the optical counterpart in both outburst and quiescence. We closely compare H 1743-322 to the well-studied black-hole X-ray transient XTE J1550-564 and find the behaviors of these systems to be very similar. As reported elsewhere, both H 1743-322 and XTE J1550-564 are relativistic jet sources and both exhibit a pair of high-frequency QPO oscillations with a 3:2 frequency ratio. The many striking similarities between these two sources argue strongly that H 1743-322 is a black hole binary, although presently no dynamical data exist to support this conclusion.

We present a dynamical model of the high mass X-ray binary LMC X-1 based on high-resolution optical spectroscopy and extensive optical and near-infrared photometry. From our new optical data we find an orbital period of P=3.90917 +/- 0.00005 days. We present a refined analysis of the All Sky Monitor data from RXTE and find an X-ray period of P=3.9094 +/- 0.0008 days, which is consistent with the optical period. A simple model of Thomson scattering in the stellar wind can account for the modulation seen in the X-ray light curves. The V-K color of the star (1.17 +/- 0.05) implies A_V = 2.28 +/- 0.06, which is much larger than previously assumed. For the secondary star, we measure a radius of R_2 = 17.0 +/- 0.8 solar radii and a projected rotational velocity of V_rot*sin(i) = 129.9 +/- 2.2 km/s. Using these measured properties to constrain the dynamical model, we find an inclination of i = 36.38 +/- 1.92 deg, a secondary star mass of M_2 = 31.79 +/- 3.48 solar masses, and a black hole mass of 10.91 +/- 1.41 solar masses. The present location of the secondary star in a temperature-luminosity diagram is consistent with that of a star with an initial mass of 35 solar masses that is 5 Myr past the zero-age main sequence. The star nearly fills its Roche lobe (~90% or more), and owing to the rapid change in radius with time in its present evolutionary state, it will encounter its Roche lobe and begin rapid and possibly unstable mass transfer on a timescale of a few hundred thousand years.

We present results from an optical spectroscopic investigation of the massive binary system [L72] LH~54-425 in the LH 54 OB association in the Large Magellanic Cloud. We revise the ephemeris of [L72] LH 54-425 and find an orbital period of 2.247409 +/- 0.000010 days. We find spectral types of O3 V for the primary and O5 V for the secondary. We made a combined solution of the radial velocities and previously published V-band photometry to determine the inclination for two system configurations, i = 52 degrees for the configuration of the secondary star being more tidally distorted and i = 55 degrees for the primary as the more tidally distorted star. We argue that the latter case is more probable, and this solution yields masses and radii of M_1 = 47 +/- 2 M_Sun and R_1 = 11.4 +/- 0.1 R_Sun for the primary, and M_2 = 28 +/- 1 M_Sun and R_2 = 8.1 +/- 0.1 R_Sun for the secondary. Our analysis places LH 54-425 amongst the most massive stars known. Based on the position of the two stars plotted on a theoretical HR diagram, we find the age of the system to be about 1.5 Myr.

In this letter we present several new three and four-planet solutions based on the most current available radial velocity data for HD160691 ($\mu$Ara). These solutions are optimised using the Planetary Orbit Fitting Process (POFP) which is programmed and executed in MATLAB. POFP is based on a full integration of the system's multiple-body Newtonian equations of motion and on a multi level optimisation utilizing a variety of algorithms. The POFP solutions are presented in the context of the Keplerian-based solutions already appearing in the literature which we have reproduced here for comparison. The different solutions and their properties are compared over all data sets separately and combined. The new solutions do not seem to exhibit instabilities and are both co-planar and three-dimensional. We also provide a comparative prediction of the published and new solutions showing their diversion in the near future. In the short term, this projection will allow to choose between the variety of solutions as further observations are made.

This Supplementary Information provides details about the spectral extraction (crowding issues and the removal of nebular lines), a discussion about the distance to M33, a model for the O-star wind and the measurement of the true photospheric X-ray eclipse width Theta, and details about ellipsoidal modelling. It also contains five related figures, one related table, and additional references.

Stellar-mass black holes are discovered in X-ray emitting binary systems, where their mass can be determined from the dynamics of their companion stars. Models of stellar evolution have difficulty producing black holes in close binaries with masses >10 solar masses, which is consistent with the fact that the most massive stellar black holes known so all have masses within 1 sigma of 10 solar masses. Here we report a mass of 15.65 +/- 1.45 solar masses for the black hole in the recently discovered system M33 X-7, which is located in the nearby galaxy Messier 33 (M33) and is the only known black hole that is in an eclipsing binary. In order to produce such a massive black hole, the progenitor star must have retained much of its outer envelope until after helium fusion in the core was completed. On the other hand, in order for the black hole to be in its present 3.45 day orbit about its 70.0 +/- 6.9 solar mass companion, there must have been a ``common envelope'' phase of evolution in which a significant amount of mass was lost from the system. We find the common envelope phase could not have occured in M33 X-7 unless the amount of mass lost from the progenitor during its evolution was an order of magnitude less than what is usually assumed in evolutionary models of massive stars.

We explore the possibility that the transit signature of an Earth-size planet can be detected in spectroscopic velocity shifts via the Rossiter effect. Under optimistic but not unrealistic conditions, it should be possible to detect a large terrestrial-size planet. While not suitable for discovering planets, this method can be used to confirm suspected planets.

We report the discovery of a new detached eclipsing binary known as NSVS01031772 with component masses M_1= 0.5428 +/- 0.0027Msun, M_2= 0.4982 +/- 0.0025Msun, and radii R_1= 0.5260 +/- 0.0028Rsun, R_2= 0.5088 +/- 0.0030Rsun. The system has an orbital period of 0.3681414(3) days and an apparent magnitude V ~ 12.6. The estimated effective temperatures of the stars are T_eff1 = 3615 +/- 72K and T_eff2 = 3513 +/- 31K. The space velocities of the system suggest that it is a main sequence binary and with a metallicity that is approximately solar. The two stars in this binary are located in a region of the Mass-Radius relation where no accurate observational data were previously available. Similarly to other low-mass binaries recently studied, the radius of each star in NSVS01031772 exceeds the best evolutionary model predictions by about 8.5% on average.

We show that the star known as 2MASS J05162881+2607387 (hereafter J0516) is a double-lined eclipsing binary with nearly identical low-mass components. The spectroscopic elements derived from 18 spectra obtained with the High Resolution Spectrograph on the Hobby-Eberly Telescope during the Fall of 2005 are K_1=88.45 +/- 0.48 km/s and K_2=90.43 +/- 0.60 km/s, resulting in a mass ratio of$q=K_1/K_2 = 0.978 +/- 0.018 and minimum masses of M_1 sin^{3}i=0.775 +/- 0.016 solar masses and M_2 sin^{3}i=0.759 +/- 0.012 solar masses, respectively. We have extensive differential photometry of J0516 obtained over several nights between 2004 January-March (epoch 1) and 2004 October-2005 January plus 2006 January (epoch 2) using the 1m telescope at the Mount Laguna Observatory. The source was roughly 0.1 mag brighter in all three bandpasses during epoch 1 when compared to epoch 2. Also, phased light curves from epoch 1 show considerable out-of-eclipse variability, presumably due to bright spots on one or both stars. In contrast, the phased light curves from epoch 2 show little out-of-eclipse variability. The light curves from epoch 2 and the radial velocity curves were analyzed using our ELC code with updated model atmospheres for low-mass stars. We find the following: M_1=0.787 +/- 0.012 solar masses, R_1=0.788 +/- 0.015 solar radii, M_2=0.770 +/- 0.009 solar masses, and R_2=0.817 +/- 0.010 solar radii. The stars in J0516 have radii that are significantly larger than model predictions for their masses, similar to what is seen in a handful of other well-studied low-mass double-lined eclipsing binaries. We compiled all recent mass and radius determinations from low-mass binaries and determine an empirical mass-radius relation of the form R = 0.0324 + 0.9343M + 0.0374M^2, where the quantities are in solar units.

We derive improved system parameters for the HD 209458 system using a model that simultaneously fits both photometric transit and radial velocity observations. The photometry consists of previous Hubble Space Telescope STIS and FGS observations, twelve I-band transits observed between 2001-2003 with the Mt. Laguna Observatory 1m telescope, and six Stromgren b+y transits observed between 2001-2004 with two of the Automatic Photometric Telescopes at Fairborn Observatory. The radial velocities were derived from Keck/HIRES observations. The model properly treats the orbital dynamics of the system, and thus yields robust and physically self-consistent solutions. Our set of system parameters agrees with previously published results though with improved accuracy. For example, applying robust limits on the stellar mass of 0.93-1.20Msun, we find 1.26 < Rplanet < 1.42 Rjup and 0.59 < Mplanet < 0.70 Mjup. We can reduce the uncertainty on these estimates by including a stellar mass-radius relation constraint, yielding Rplanet = 1.35 +/- 0.07 Rjup and Mplanet = 0.66 +/- 0.04 Mjup. Our results verify that the planetary radius is 10-20% larger than predicted by planet evolution models, confirming the need for an additional mechanism to slow the evolutionary contraction of the planet. A revised ephemeris is derived, T0=2452854.82545 + 3.52474554E (HJD), which now contains an uncertainty in the period of 0.016s and should facilitate future searches for planetary satellites and other bodies in the HD 209458 system.

We report on the formation and evolution of two large-scale, synchrotron-emitting jets from the black hole candidate H 1743-322 following its reactivation in 2003. In November 2003 after the end of its 2003 outburst, we noticed, in observations with the Australia Telescope Compact Array, the presence of a new and variable radio source about 4.6" to the East of H 1743-322, that was later found to move away from H 1743-322. In February 2004, we detected a radio source to the West of H 1743-322, symmetrically placed relative to the Eastern jet. In 2004, follow-up X-ray observations with {\em Chandra} led to the discovery of X-ray emission associated with the two radio sources. This likely indicates that we are witnessing the interaction of relativistic jets from H 1743-322 with the interstellar medium causing in-situ particle acceleration. The spectral energy distribution of the jets during the decay phase is consistent with a classical synchrotron spectrum of a single electron distribution from radio up to X-rays, implying the production of very high energy ($>$ 10 TeV) particles in those jets. We discuss the jet kinematics, highlighting the presence of a significantly relativistic flow in H 1743-322 almost a year after the ejection event.

(Shortened) We present R-band photometry of the X-ray transient and candidate black hole binary XTE J1650-500 obtained between 2003 May and August with the 6.5m Clay Telescope. A timing analysis of these data reveals a photometric period of 0.3205 +/- 0.0007 days (i.e. 7.63 hr) with a possible alias at 0.3785 days (9.12 hr). Our photometry completely rules out the previously published spectroscopic period of 0.212 days (5.09 hr). Consequently, we reanalyzed the 15 archival ESO/VLT spectra (obtained 2002 June by Sanchez-Fernandez et al.) that were the basis of the previously published spectroscopic period. We used a ``restframe search'' technique that is well suited for cases when the signal-to-noise ratio of individual spectra is low. The results confirmed the photometric period of 0.3205 days, and rule out the alias period near 0.38 days. The best value for the velocity semiamplitude of the companion star is K_2 = 435 +/- 30 km/sec, and the corresponding optical mass function is f(M) = 2.73 +/- 0.56 solar masses. The amplitude of the phased R-band light curve is 0.2 magnitudes, which gives a lower limit to the inclination of 50 +/- 3 degrees in the limiting case of no contribution to the R-band light curve from the accretion disk. If the mass ratio of XTE J1650-500 is similar to the mass ratios of other black hole binaries like A0620-00 or GRS 1124-683 (e.g. Q >~ 10), then our lower limit to the inclination gives an upper limit to the mass of the black hole in XTE J1650-500 of M_1 <~ 7.3 solar masses. However, the mass can be considerably lower if the R-band flux is dominated by the accretion disk. For example, if the accretion disk does contribute 80% of the flux, as our preliminary results suggest, then the black hole mass would be only about 4 solar masses.

We observed the 2003 X-ray outburst of H1743-322 in a series of 130 pointed observation with RXTE. We searched individual observations for high-frequency QPOs (HFQPOs) and found only weak or marginal detections near 240 and 160 Hz. We next grouped the observations in several different ways and computed the average power-density spectra (PDS) in a search for further evidence of HFQPOs. This effort yielded two significant results for those observations defined by the presence of low-frequency QPOs (0.1-20 Hz) and an absence of ``band-limited'' power continua: (1) The 9 time intervals with the highest 7-35 keV count rates yielded an average PDS with a QPO at $166 \pm 5$ Hz. ($4.1 \sigma$; 3--35 keV); and (2) a second group with lower 7-35 keV count rates (26 intervals) produced an average PDS with a QPO at $242 \pm 3$ Hz ($6.0 \sigma$; 7--35 keV). The ratio of these two frequencies is $1.46 \pm 0.05$. This finding is consistent with results obtained for three other black hole systems that exhibit commensurate HFQPOs in a 3:2 ratio. Furthermore, the occurrence of H1743-322's slower HFQPO at times of higher X-ray luminosity closely resembles the behavior of XTE J1550-564 and GRO J1655-40. We discuss our results in terms of a resonance model that invokes frequencies set by general relativity for orbital motions near a black-hole event horizon.

We present an X-ray spectral and timing analysis of 4U 1543-47 during its 2002 outburst based on 49 pointed observations obtained using the Rossi X-ray Timing Explorer (RXTE). The outburst reached a peak intensity of 4.2 Crab in the 2-12 keV band and declined by a factor of 32 throughout the month-long observation. A 21.9 +- 0.6 mJy radio flare was detected at 1026.75 MHz two days before the X-ray maximum; the radio source was also detected late in the outburst, after the X-ray source entered the low hard state. The X-ray light curve exhibits the classic shape of a rapid rise and an exponential decay. The spectrum is soft and dominated by emission from the accretion disk. The continuum is fit with a multicolor disk blackbody (kT_{max} = 1.04 keV) and a power-law (Gamma ~ 2.7). Midway through the decay phase, a strong low-frequency QPO (nu = 7.3-8.1 Hz) was present for several days. The spectra feature a broad Fe K alpha line that is asymmetric, suggesting that the line is due to relativistic broadening rather than Comptonization. Relativistic Laor models provide much better fits to the line than non-relativistic Gaussian models, particularly near the beginning and end of our observations. The line fits yield estimates for the inner disk radius that are within 6 R_g; this result and additional evidence indicates that this black hole may have a non-zero angular momentum.

(Abridged) We present densely sampled BVRI light curves of the optical transient associated with the gamma-ray burst GRB 030329, the result of a coordinated observing campaign conducted at five observatories. Augmented with published observations of this GRB, the compiled optical dataset contains 2687 photometric measurements, obtained between 78 minutes and 79 days after the burst. We show that the underlying supernova 2003dh evolved faster than, and was probably somewhat fainter than the type Ic SN 1998bw, associated with GRB 980425. We find that our data can be described by a broken power-law decay perturbed by a complex variable component. The early- and late-time decay slopes are determined to be ~1.1 and ~2, respectively. Assuming this single power-law model, we constrain the break to lie between ~3 and ~8 days after the burst. This simple, singly-broken power-law model, derived only from the analysis of our optical observations, may also account for available multi-band data, provided that the break happened ~8 days after the burst. The more complex double-jet model of Berger et al. provides a comparable fit to the optical, X-ray, mm and radio observations of this event. We detect a significant change in optical colors during the first day. Our color analysis is consistent with a cooling break frequency sweeping through the optical band during the first day. The light curves of GRB 030329 reveal a rich array of variations, superposed over the mean power-law decay. We find that the early variations are asymmetric, with a steep rise followed by a relatively slower (by a factor of about two) decline. The variations maintain a similar time scale during the first four days, and then get significantly longer.

Using the Chandra X-ray Observatory, we have detected the black hole transients V4641 Sgr and XTE J1859+226 in their low luminosity, quiescent states. The 0.3-8 keV luminosities are (4.0^(+3.3)_(-2.4))E31 (d/7 kpc)^2 erg/s and (4.2^(+4.8)_(-2.2))E31 (d/11 kpc)^2 erg/s for V4641 Sgr and XTE J1859+226, respectively. With the addition of these 2 systems, 14 out of the 15 transients with confirmed black holes (via compact object mass measurements) now have measured quiescent luminosities or sensitive upper limits. The only exception is GRS 1915+105, which has not been in quiescence since its discovery in 1992. The luminosities for V4641 Sgr and XTE J1859+226 are consistent with the median luminosity of 2E31 erg/s for the systems with previous detections. Our analysis suggests that the quiescent X-ray spectrum of V4641 Sgr is harder than for the other systems in this group, but, due to the low statistical quality of the spectrum, it is not clear if V4641 Sgr is intrinsically hard or if the column density is higher than the interstellar value. Focusing on V4641 Sgr, we compare our results to theoretical models for X-ray emission from black holes in quiescence. Also, we obtain precise X-ray positions for V4641 Sgr and XTE J1859+226 via cross-correlation of the X-ray sources detected near our targets with IR sources in the 2 Micron All-Sky Survey catalog.

An extensive grid of synthetic mid- and far-ultraviolet spectra for accretion disks in cataclysmic variables has been presented by Wade and Hubeny (1998). In those models, the disk was assumed to be in steady-state, that is T_eff(r) is specified completely by the mass M_WD and radius R_WD of the accreting white dwarf star and the mass transfer rate M_dot which is constant throughout the disk. In these models, T_eff(r) is proportional to r^{-3/4} except as modified by a cutoff term near the white dwarf. Actual disks may vary from the steady-state prescription for T_eff(r), however, e.g. owing to outburst cycles in dwarf novae M_dot not constant with radius) or irradiation (in which case T_eff in the outer disk is raised above T_steady). To show how the spectra of such disks might differ from the steady case, we present a study of the ultraviolet (UV) spectra of models in which power-law temperature profiles T_eff(r) is proportional to r^{-gamma} with gamma < 3/4 are specified. Otherwise, the construction of the models is the same as in the Wade & Hubeny grid, to allow comparison. We discuss both the UV spectral energy distributions and the appearance of the UV line spectra. We also briefly discuss the eclipse light curves of the non-standard models. Comparison of these models with UV observations of novalike variables suggests that better agreement may be possible with such modified T_eff(r) profiles.

We present an X-ray/UV/optical spectrum of the black hole primary in the X-ray nova XTE J1118+480 in quiescence at L_x = 4 x 10^{-9} of the Eddington luminosity. The Chandra, HST and MMT spectroscopic observations were performed simultaneously on 2002 January 12 UT. Because this 4.1-hr binary is located at b = 62 deg, the transmission of the ISM is very high (e.g., 70% at 0.3 keV). We present many new results for the quiescent state, such as the first far-UV spectrum and evidence for an 0.35 mag orbital modulation in the near-UV flux. However, the centerpiece of our work is the multiwavelength spectrum of XTE J1118+480, which we argue represents the canonical spectrum of a stellar-mass black hole radiating at L_x = 4 x 10^{-8.5} of the Eddington luminosity. This spectrum is comprised of two apparently disjoint components: a hard X-ray spectrum with a photon index Gamma = 2.02 +/- 0.16, and an optical/UV continuum that resembles a 13,000 K disk blackbody spectrum punctuated by several strong emission lines. We present a model of the source in which the accretion flow has two components: (1) an X-ray-emitting interior region where the flow is advection-dominated, and (2) a thin, exterior accretion disk with a truncated inner edge (R_tr ~ 10^4 Schwarzschild radii) that is responsible for the optical/UV spectrum. For D = 1.8 kpc, the luminosity of the X-ray component is L_x = 3.5 x 10^{30} erg/s (0.3-7 keV); the bolometric luminosity of the optical/UV component is 20 times greater.

In this brief note we point out that the four spatially resolved relativistic jets among the 14 dynamically confirmed black hole X-ray novae are all in systems with long orbital periods. Many shorter period systems show transient radio outbursts which are attributed to jets, but these jets have not been spatially resolved. Super-Eddington accretion has been suggested as a requirement for jet formation and may be consistent with our compilation of luminosities, but some super-Eddington outbursts did not form spatially resolved jets. We speculate that some as yet unknown process (or combination of processes) favors formation of substantially larger jets in long period systems. Two short period systems show evidence for extent, but have not been resolved into multiple components as the long period systems have.

We report on X-ray observations of the the large-scale jets recently discovered in the radio and detected in X-rays from the black hole candidate X-ray transient and microquasar XTE J1550-564. On 11 March 2002, X-ray emission was detected 23 arcsec to the West of the black hole candidate and was extended along the jet axis with a full width at half maximum of 1.2 arcsec and a full width at 10% of maximum intensity of 5 arcsec. The morphology of the X-ray emission matched well to that of the radio emission at the same epoch. The jet moved by 0.52 +/- 0.13 arcsec between 11 March and 19 June 2002. The apparent speed during that interval was 5.2 +/- 1.3 mas/day. This is significantly less than the average apparent speed of 18.1 +/- 0.4 mas/day from 1998 to 2002, assuming that the jet was ejected in September 1998, and indicates that the jet has decelerated. The X-ray spectrum is adequately described by a powerlaw with a photon index near 1.8 subject to interstellar absorption. The unabsorbed X-ray flux was 3.4 x 10^-13 erg cm^-2 s^-1 in the 0.3-8 keV band in March 2002, and decreased to 2.9 x 10^-13 erg cm^-2 s^-1 in June. We also detect X-rays from the eastern jet in March 2002 and show that it has decelerated and dimmed since the previous detections in 2000.

We have discovered an X-ray jet due to material ejected from the black hole X-ray transient XTE J1550-564. The discovery was first reported by Corbel et al. (Science, 298, 196 and astro-ph/0210224), and here, we present an analysis of the three Chandra observations made between 2000 June and 2000 September. For these observations, a source is present that moves in an eastward direction away from the point source associated with the compact object. The separation between the new source and the compact object changes from 21.3 arcseconds in June to 23.4 arcseconds in September, implying a proper motion of 21.2 +/- 7.2 mas/day, a projected separation of 0.31-0.85 pc and an apparent jet velocity between 0.34 +/- 0.12 and 0.93 +/- 0.32 times the speed of light for a source distance range of d = 2.8-7.6 kpc. These observations represent the first time that an X-ray jet proper motion measurement has been obtained for any accretion powered Galactic or extra-galactic source. While this work deals with the jet to the east of the compact object, the western jet has also been detected in the X-ray and radio bands. The most likely scenario is that the eastern jet is the approaching jet and that the jet material was ejected from the black hole in 1998. Along with a 1998 VLBI proper motion measurement, the Chandra proper motion indicates that the eastern jet decelerated between 1998 and 2000. We present results on the morphology and energy spectrum of the jet. We cannot definitively determine the X-ray emission mechanism, but a synchrotron origin is viable and may provide the simplest explanation for the observations.

We have discovered at x-ray and radio wavelengths large-scale moving jets from the microquasar XTE J1550-654. Plasma ejected from near the black hole traveled at relativistic velocities for at least four years. We present direct evidence for gradual deceleration in a relativistic jet. The broadband spectrum of the jets is consistent with synchrotron emission from high energy (up to 10 TeV) particles accelerated in shock waves formed within the relativistic ejecta or by the interaction of the jets with the interstellar medium. XTE J1550-654 offers a unique opportunity to study the dynamical evolution of relativistic jets on time scales inaccessible for active galactic nuclei jets, with implications for our understanding of relativistic jets from Galactic x-ray binaries and active galactic nuclei.

We present a high-resolution echelle spectrum taken with the Very Large Telescope and analyse archival Hubble Space Telescope photometry of the recently identified companion of the eclipsing millisecond radio pulsar PSR J1740-5340 in the globular cluster NGC 6397. From the spectrum, we show that the companion is metal poor, as expected for a member of NGC 6397. Using synthetic photometry, and assuming a true distance modulus of 12.13 +/- 0.15 mag and a colour excess of E(B-V) = 0.179, we derive a radius of 1.60 +/- 0.17 solar radii and an effective temperature of 5410 +/- 50 K, implying a luminosity of 2.0 +/- 0.4 solar luminosities (1 sigma errors). These properties make it similar to the so-called `sub-subgiants' in the old open cluster M67 and `red stragglers' in 47 Tuc, which have luminosities comparable to those of turn-off stars, but cooler temperatures and larger radii. The light curve of the companion is well described by ellipsoidal variations, and despite the incomplete (~60%) phase coverage, we are able to derive good constraints on a number of the system parameters. In particular, for the inclination we find a 2 sigma lower limit of i > 48 deg (i > 46 deg at 3 sigma). Assuming a pulsar mass of 1.2 M_MSP < 2.4 solar masses, this implies a companion mass in the range 0.14 < M_comp < 0.38 solar masses. Combined with the photometric constraint, we find a best fit for i ~ 50 deg and M_comp ~ 0.3 solar masses. We infer a Roche lobe filling factor by radius of ~97%. [Abstract shortened]

A small group of X-ray binaries currently provides the best evidence for the existence of stellar-mass black holes. These objects are interacting binary systems where the X-rays arise from accretion of material onto a compact object (i.e. an object with a radius of less than a few hundred km). In some favourable cases, optical studies of the companion star lead to dynamical mass estimates for both components. In 17 cases, the mass of the compact object an X-ray binary has been shown to exceed the maximum mass of a stable neutron star (about 3 solar masses), which leads to the conclusion that these objects are black holes. In this contribution I will review the basic properties of these black hole binaries.

We briefly review the properties and physical consequences of quasiperiodic oscillations (QPOs) seen on many occasions in the X-ray emission from black-hole binary systems. We also report on the discovery of QPOs at 164 and 328 Hz in GRS1915+105. These appear to be analogous to the harmonic systems of high-frequency QPOs seen in XTEJ1550-564 and GROJ1655-40.

We continue to investigate the X-ray timing and spectral properties of the black hole binary, XTE J1550-564. For three different groups of observations, we show evidence that two high-frequency QPOs (HFQPOs) occur simultaneously near 184 and 276 Hz. In one group, there is also evidence of a broad feature at 92 Hz. In each case, we can model the QPO profiles while assuming that the central frequencies are related as integral harmonics of a single frequency. We next investigate the X-ray energy spectra, and we find a systematic increase in the strength of the power-law component as the stronger of the two HFQPOs shifts from 276 to 184 Hz. A strikingly similar result is seen for GRO J1655-40 when the stronger HFQPO shifts from 450 to 300 Hz. The fundamental HFQPO frequencies for the two X-ray sources scale inversely with black hole mass, which is consistent with the hypotheses that these HFQPOs represent some oscillation rooted in general relativity (GR) and that the two black holes have similar values of the dimensionless spin parameter. We discuss physical mechanisms that may explain these HFQPOs. In particular, a resonance between the orbital and radial coordinate frequencies in the innner accretion disk, as proposed by Abramowicz & Kluzniak, would imply moderate values for the dimensionless spin parameter (0.1 < a < 0.6) for both black holes.

Optical spectroscopic observations of the companion star (type G8IV to K4III) in the microquasar system XTE J1550-564 reveal a radial velocity curve with a best fitting spectroscopic period of P_sp = 1.552 +/- 0.010 days and a semiamplitude of K_2 = 349 +/- 12 km/s. The optical mass function is f(M) = 6.86 +/- 0.71 solar masses (1 sigma). We tentatively measure the rotational velocity of the companion star to be V_rot*sin(i) = 90 +/- 10 km/s, which when taken at face value implies a mass ratio of Q = M_1/M_2 = 6.6 (+2.5 / -1.6) (1 sigma), using the above value of K_2. We derive constraints on the binary parameters from simultaneous modelling of the ellipsoidal light and radial velocity curves. We find the most likely value of the mass of the compact object is 9.41 solar masses with a 1 sigma range of 8.36 < M_1 < 10.76 solar masses. If we apply our tentative value of V_rot*sin(i) = 90 +/- 10 km/s as an additional constraint in the ellipsoidal modelling, we find the most likely value of the mass of the compact object is 10.56 solar masses with a 1 sigma range of 9.68 < M_1 < 11.58 solar masses. In either case the mass of the compact object is well above the maximum mass of a stable neutron star, and we therefore conclude XTE J1550-564 contains a black hole.

We have obtained optical and infrared photometry of the soft x-ray transient prototype V616 Mon (A0620-00). From this photometry, we find a spectral type of K4 for the secondary star in the system, which is consistent with spectroscopic observations. We present J-, H-, and K-band light curves modeled with WD98 and ELC. Combining detailed, independently run models for ellipsoidal variations due to a spotted, non-spherical secondary star, and the observed ultraviolet to infrared spectral energy distribution of the system, we show that the most likely value for the orbital inclination is 40.75 +/- 3 deg. This inclination angle implies a primary black hole mass of 11.0 +/- 1.9 solar masses.

We present extensive spectroscopy and photometry of the dwarf nova V630 Cassiopeiae. A late-type (K4-5) absorption spectrum is easily detectable, from which we derive the orbital parameters. We find a spectroscopic period of P=2.56387 +/- (4 times 10^{-5}) days and a semiamplitude of K_2=132.9 +/- 4.0 km/s. The resulting mass function, which is a firm lower limit on the mass of the white dwarf, is then f(M)=0.624 +/- 0.056 solar masses. The secondary star is a ``stripped giant'', and using relations between the core mass and the luminosity and the core mass and the radius we derive a lower limit of M_2 > 0.165 solar masses for the secondary star. The rotational velocity of the secondary star is not resolved in our spectra and we place a limit of V_rot*sin(i) < 40 km/s. The long-term light curve shows variations of up to 0.4 mag on short (1-5 days) time scales, and variations of 0.2-0.4 mag on longer (3-9 months) time scales. In spite of these variations, the ellipsoidal light curve of the secondary star is easily seen when the data are folded on the spectroscopic ephemeris. Ellipsoidal models fit to the mean light curve give an inclination in the range 66.96 < i < 78.08 degrees (90 per cent confidence). This inclination range, and the requirement that M_2 > 0.165 solar masses and V_rot*sin(i) < 40 km/s yields a white dwarf mass of M_1=0.977^{+0.168}_{-0.098} solar masses and a secondary star mass of M_2=0.172^{+0.029}_{-0.012} solar masses (90 per cent confidence limits). Our findings confirm the suggestion of Warner (1994), namely that V630 Cas is rare example of a dwarf nova with a long orbital period.

We report optical, infrared, and X-ray light curves for the outburst, in 2000, of the black hole candidate XTE J1550-564. We find that the start of the outburst in the H and V bands precedes that seen in the RXTE All Sky Monitor by 11.5 +/- 0.9 and 8.8 +/- 0.6 days, respectively; a similar delay has been observed in two other systems. About 50 days after the primary maxima in the VIH light curves, we find secondary maxima, most prominently in H. This secondary peak is absent in the X-ray light curve, but coincides with a transition to the low/hard state. We suggest that this secondary peak may be due to non-thermal emission associated with the formation of a jet.

We present a photometric and spectroscopic study of the blue straggler S1082 in the open cluster M67. Our observations confirm the previously reported 1.07 day eclipse light curve and the absence of large radial-velocity variations of the narrow-lined star. However, we find two more spectral components which do vary on the 1.07 day period. We conclude that the system is triple. We fit the light and radial-velocity curves and find that the total mass of the inner binary is more than twice the turnoff mass and that the outer companion to the binary is a blue straggler on its own account. We briefly discuss formation scenarios for this multiple system.

(shortened) Spectroscopic observations of the fast X-ray transient and superluminal jet source SAX J1819.3-2525 (V4641 Sgr) reveal a best fitting period of P_spect=2.81678 +/- 0.00056 days and a semiamplitude of K_2=211.0 +/- 3.1 km/sec. The optical mass function is f(M)=2.74 +/- 0.12 solar masses. We find a photometric period of P_photo=2.81730 +/- 0.00001 days using a light curve measured from photographic plates. The folded light curve resembles an ellipsoidal light curve with two maxima of roughly equal height and two minima of unequal depth per orbital cycle. The secondary star is a late B-type star which has evolved off the main sequence. Using a moderate resolution spectrum (R=7000) we measure T_eff=10500 +/- 200K, log(g)=3.5 +/- 0.1, and V_rot*sin(i)=123 +/- 4 km/sec (1 sigma errors). Assuming synchronous rotation, our measured value of the projected rotational velocity implies a mass ratio of Q=M_1/M_2=1.50 +/- 0.08 (1sigma). The lack of X-ray eclipses implies an upper limit to the inclination of i<70.7 deg. On the other hand, the large amplitude of the folded light curve (about 0.5 mag) implies a large inclination (i>60 deg). Using the above mass function, mass ratio, and inclination range, the mass of the compact object is in the range 8.73 < M_1 < 11.70 solar masses and the mass of the secondary star is in the range 5.49 < M_2 < 8.14 solar masses (90% confidence). The mass of the compact object is well above the maximum mass of a stable neutron star and we conclude that V4641 Sgr contains a black hole.

We present BVIJK photometry of the black-hole candidate GRO J1655-40 in full quiescence. We report a refined orbital period of 2.62191 +/- 0.00020 days. The light curves are dominated by ellipsoidal variations from the secondary star. We model the light curves with an upgraded code which includes a more accurate treatment of limb darkening. Previous models containing a large cool disk are ruled out, and indeed our data can be fit with a pure ellipsoidal light curve without any disk contribution. In general agreement with previous results, we derive a confidence region of the correlated quantities of inclination and mass ratio, centered on an inclination of 70.2 +/ 1.9 degrees, and mass ratio 2.6 +/- 0.3, resulting in a primary mass M = 6.3 +/- 0.5 Mo (all 95% confidence). The complex limits and errors on these values, and on the possible disk contribution to the light curve, are discussed.

We have written a light curve synthesis code that makes direct use of model atmosphere specific intensities, in particular the NEXTGEN model atmosphere grid for cool giants (T_eff <= 6800 K and log(g) <= 3.5, Hauschildt et al. 1999). We point out that these models (computed using spherical geometry) predict a limb darkening behaviour that deviates significantly from a simple linear or two-parameter law (there is less intensity at the limb of the star). The presence of a significantly nonlinear limb darkening law has two main consequences. First, the ellipsoidal light curve computed for a tidally distorted giant using the NEXTGEN intensities is in general different from the light curve computed using the same geometry but with the black body approximation and a one- or two-parameter limb darkening law. In most cases the light curves computed with the NEXTGEN intensities have deeper minima than their black body counterparts. Thus the light curve solutions for binaries with a giant component obtained with models with near linear limb darkening (either black body or plane-parallel model atmosphere intensities) are biased. Observations over a wide wavelength range (i.e. both the optical and infrared) are particularly useful in discriminating between models with nearly linear limb darkening and the NEXTGEN models. Second, we show that rotational broadening kernels for Roche lobe filling (or nearly filling) giants can be significantly different from analytic kernels due to a combination of the nonspherical shape of the star and the radical departure from a simple limb darkening law. As a result, geometrical information inferred from V_rot*sin(i) measurements of cool giants in binary systems are likewise biased.

We report optical monitoring of the soft X-ray transient XTE J1550-564 during the 1999 season (4 January 1999 to 24 August 1999). The first optical observations available in 1999 show that the peak ``re-flare'' brightness had exceeded the peak brightness of the initial optical flare in September 1998 by over half a magnitude. We compare the optical re-flare light curves with the total X-ray flux, the power-law flux and disk flux light curves constructed from the spectral fits to RXTE/PCA data made by Sobczak et al. (1999, 2000). During the first 60 days of the observed optical re-flare, we find no correspondence between the thermal component of the X-rays often associated with a disk and the optical flux -- the former remains essentially flat whereas the latter declines exponentially and exhibits three substantial dips. However, the power law flux is anti-correlated with the optical dips, suggesting that the optical flux may by up-scattered into the X-ray by the hot corona. Periodic modulations were discovered during the final stage of the outburst (May to June), with P=1.546+/-0.038 days, and during quiescence (July and August), with P=1.540+/-0.008 days. The analysis of the combined data set reveals a strong signal for a unique period at P=1.541+/-0.009 days, which we believe to be the orbital period.

We report on the X-ray spectral behavior of the exceptionally bright X-ray nova XTE J1550-564 during its 1998-99 outburst. Our study is based on 209 pointed observations using the PCA and HEXTE instruments onboard the Rossi X-ray Timing Explorer spanning 250 days and covering the entire double-peaked eruption that occurred from 1998 September until 1999 May. The spectra are fit to a model including multicolor blackbody disk and power-law components. The source is observed in the very high and high/soft outburst states of black hole X-ray novae. During the very high state, when the power-law component dominated the spectrum, the inner disk radius is observed to vary by more than an order of magnitude; the radius decreased by a factor of 16 in one day during a 6.8 Crab flare. If the larger of these observed radii is taken to be the last stable orbit, then the smaller observed radius would imply that the inner edge of the disk is inside the event horizon! However, we conclude that the apparent variations of the inner disk radius observed during periods of increased power-law emission are probably caused by the failure of the multicolor disk/power-law model; the actual physical radius of the inner disk may remain fairly constant. This interpretation is supported by the fact that the observed inner disk radius remains approximately constant over 120 days in the high state, when the power-law component is weak, even though the disk flux and total flux vary by an order of magnitude. The mass of the black hole inferred by equating the approximately constant inner disk radius observed in the high/soft state with the last stable orbit for a Schwarzschild black hole is M_BH = 7.4 M_sun (D/6 kpc) (cos i)^{-1/2}.

Utilizing observations obtained with the Rossi X-ray Timing Explorer, we examine correlations between the properties of 0.08-22 Hz variable frequency QPOs and the X-ray spectral parameters for the black hole candidates XTE J1550-564 and GRO J1655-40. The spectra were fitted to a model including a multi-temperature blackbody disk and a power-law component. We find that the QPO frequency and amplitude are well correlated with the spectral parameters, although the correlations found for XTE J1550-564 are generally opposite to those for GRO J1655-40. There is one exeption: Both sources exhibit a general increase in the QPO frequency as the disk flux increases (or as the mass accretion rate through the disk increases). In addition, these QPOs are observed only when the power-law component contributes more than 20% of the 2-20 keV flux, which indicates that both the disk and the power-law components are linked to the QPO phenomenon.

This paper presents analysis and interpretation of OGLE photometric data of four X-ray binary pulsar systems in the Small Magellanic Cloud: 1WGA J0054.9-7226, RX J0050.7-7316, RX J0049.1-7250, and 1SAX J0103.2-7209. In each case, the probable optical counterpart is identified on the basis of its optical colours. In the case of RX J0050.7-7316 the regular modulation of its optical light curve appears to reveal an ellipsoidal modulation with a period of 1.416 days. Using reasonable masses for the neutron star and the B star, we show that the amplitude and relative depths of the minima of the I-band light curve of RX J0050.7-7316 can be matched with an ellipsoidal model where the B star nearly fills its Roche lobe. For mass ratios in the range of 0.12 to 0.20, the corresponding best-fitting inclinations are about 55 degrees or larger. The neutron star would be eclipsed by the B star at inclinations larger than 60 degrees for this particular mass ratio range. Thus RX J0050.7-7316 is a good candidate system for further study. In particular, we would need additional photometry in several colours, and most importantly, radial velocity data for the B star before we could draw more quantitative conclusions about the component masses.

We report additional photometric CCD observations of KPD 0422+5421, a binary with an orbital period of 2.16 hours which contains a subdwarf B star (sdB) and a white dwarf. There are two main results of this work. First, the light curve of KPD 0422+5421 contains two distinct periodic signals, the 2.16 hour ellipsoidal modulation discovered by Koen, Orosz, & Wade (1998) and an additional modulation at 7.8 hours. This 7.8 hour modulation is clearly not sinusoidal: the rise time is about 0.25 in phase, whereas the decay time is 0.75 in phase. Its amplitude is roughly half of the amplitude of the ellipsoidal modulation. Second, after the 7.8 hour modulation is removed, the light curve folded on the orbital period clearly shows the signature of the transit of the white dwarf across the face of the sdB star and the signature of the occultation of the white dwarf by the sdB star. We used the Wilson-Devinney code to model the light curve to obtain the inclination, the mass ratio, and the Omega potentials, and a Monte Carlo code to compute confidence limits on interesting system parameters. We find component masses of M_sdB = 0.36 +/- 0.16 solar masses and M_WD = 0.47 +/- 0.16 solar masses (M_total = 0.86 +/- 0.35 solar masses, 68 per cent confidence limits). If we impose an additional constraint and require the computed mass and radius of the white dwarf to be consistent with a theoretical mass-radius relation, we find M_sdB = 0.511 +0.047 -0.050 solar masses and M_WD = 0.526 +0.033 -0.030 solar masses (68 per cent confidence limits). In this case the total mass of the system is less than 1.4 solar masses at the 99.99 per cent confidence level. We briefly discuss possible interpretations of the 7.8 hour modulation and the importance of KPD 0422+5421 as a member of a rare class of evolved binaries.

Excellent coverage of the 1998 outburst of the X-ray Nova XTE J1550-564 was provided by the Rossi X-ray Timing Explorer. XTE J1550-564 exhibited an intense (6.8 Crab) flare on 1998 September 19 (UT), making it the brightest new X-ray source observed with RXTE. We present a spectral analysis utilizing 60 Proportional Counter Array spectra from 2.5-20 keV spanning 71 days, and a nearly continuous All Sky Monitor light curve. The spectra were fit to a model including multicolor blackbody disk and power-law components. XTE J1550-564 is observed in the very high, high/soft, and intermediate canonical outburst states of Black Hole X-ray Novae.

We report the identification of the optical counterpart of the X-ray transient XTE J1550-564 described in two companion papers by Sobczak et al (1999) and Remillard et al (1999). We find that the optical source brightened by approximately 4 magnitudes over the quiescent counterpart seen at B~22 on a SERC survey plate, and then decayed by approximately 1.5 magnitudes over the 7 week long observation period. There was an optical response to the large X-ray flare described by Sobczak et al (1999), but it was much smaller and delayed by roughly 1 day.

We have investigated the X-ray timing properties of XTE J1550-564 during 60 RXTE PCA observations made between 1998 September 18 and November 28. We detect quasi-periodic oscillations (QPOs) near 185 Hz during four time intervals. The QPO widths (FWHM) are near 50 Hz, and the rms amplitudes are about 1% of the mean flux at 2-30 keV. This is the third Galactic black hole candidate known to exhibit a transient X-ray timing signature above 50 Hz, following the 67 Hz QPO in GRS1915+105 and the 300 Hz QPO in GRO J1655-40. However, unlike the previous cases, which appear to show stationary frequencies, the QPO frequency in XTE J1550-564 must vary by at least 10% to be consistent with observations. The occurrences and properties of the QPO were insensitive to large changes in the X-ray intensity (1.5 to 6.8 Crab). However, the QPO appearance was accompanied by changes in the energy spectrum, namely, an increase in the temperature and a decrease in the normalization of the thermal component. The QPO is also closely related to the hard X-ray power-law component of the energy spectrum since the fractional amplitude of the QPO increases with photon energy. The fast QPOs in accreting black hole binaries are thought to be effects of general relativity; however, the relevance of the specific physical models that have been proposed remains largely uncertain. Low frequency QPOs in the range 3-13 Hz were often observed. Occasionally at high luminosity the rms QPO amplitude was 15% of the flux, a level previously reached only by GRS1915+105.

Excellent coverage of the entire 16-month 1996-97 outburst cycle of GRO J1655-40 was provided by RXTE. We present a full spectral analysis of these data, which includes 52 PCA spectra from 2.5-20 keV and HEXTE spectra above 20 keV. We also include a nearly continuous ASM light curve with several intensity measurements per day. The data are interpreted in the context of the multicolor blackbody disk/power-law model. The source is observed in the very high, high/soft, and low/hard outburst states. During the very high state, the source exhibits intense hard flares on time scales of hours to days which are correlated with changes in both the fitted temperature and radius of the inner accretion disk. During the high/soft state, the spectrum is dominated by the soft thermal emission from the accretion disk with spectral parameters that suggest approximately constant inner disk radius and temperature. We find that a tight relationship exists between the observed inner radius of the disk and the flux in the power-law component. During intense hard flares, the inner disk radius is observed to decrease by as much as a factor of three on a time scale of days. The apparent decrease of the inner disk radius observed during the flares may be due to the failure of the multicolor disk model caused by a steepening of the radial temperature profile in the disk coupled with increased spectral hardening and not physical changes of the inner disk radius. Assuming that our spectral model is valid during periods of weak power-law emission, our most likely value for the inner disk radius implies a* < 0.7. Such a low value for the black hole angular momentum is inconsistent with the relativistic frame dragging and the `diskoseismic' models as interpretations for the 300 Hz X-ray QPO seen during some of these RXTE observations.

We present U, B and V light curves (taken from the literature) of the low mass X-ray binary Cygnus X-2. The ``lower envelope'' of the light curves folded on the orbital period are ellipsoidal. We fit an ellipsoidal model to the lower envelopes of the B and V light curves to derive inclination constraints. If we assume the accretion disc is steady-state where its radial temperature profile goes as T(r) \propto r^{-3/4}, we find an inclination of i = 62.5 +/- 4 deg. However, the predicted ratio of the disc flux to the total flux in B (the ``disc fraction'') is larger than what is observed (about 0.55 compared to < 0.3). If we use a flatter radial temperature profile of the disc expected for strongly irradiated discs (T(r) \propto r^{-3/7}), then we find an inclination of i = 54.6 deg and a disc fraction in B of approximately 0.30. However, in this case the value of chi^2 is much larger (48.4 with 36 degrees of freedom compared to 40.9 for the steady-state case). Adopting i = 62.5 +/- 4 deg and using a previous determination of the mass ratio (q = M_c/M_x = 0.34 +/- 0.04) and the optical mass function (f(M) = 0.69 +/- 0.03 solar masses), we find that the mass of the neutron star is M_x = 1.78 +/- 0.23 solar masses and the mass of the secondary star is M_c = 0.60 +/- 0.13 solar masses. We derive a distance of d = 7.2 +/- 1.1 kpc, which is significantly smaller than a recent distance measurement of d = 11.6 +/- 0.3 kpc derived from an observation of a type I radius-expansion X-ray burst, but consistent with earlier distance estimates. (abridged)

We have made extensive spectroscopic and photometric observations of PG 1224+309, a close binary containing a DA white dwarf primary and an M4+ secondary. The H alpha line is in emission due to irradiation of the M-star by the hot white dwarf and is seen to vary around the orbit. From the radial velocities of the H alpha line we derive a period of P = 0.258689 +/- 0.000004 days and a semi-amplitude of K_Halpha = 160 +/- 8 km/s. We estimate a correction Delta_K = 21 +/- 2 km/s, where K_M = K_Halpha + Delta_K. Radial velocity variations of the white dwarf reveal a semi-amplitude of K_WD = 112 +/- 14 km/s. The blue spectrum of the white dwarf is well fit by a synthetic spectrum having T_eff = 29,300 K and log(g) = 7.38. The white dwarf contributes 97% of the light at 4500 Angstroms and virtually all of the light blueward of 3800 Angstroms. No eclipses are observed. The mass inferred for the white dwarf depends on the assumed mass of the thin residual hydrogen envelope: 0.40 < M_WD < 0.45 solar masses for hydrogen envelope masses of 0 < M_H < 4.0E-4 solar masses. We argue that the mass of the white dwarf is closer to 0.45 solar masses, hence it appears that the white dwarf has a relatively large residual hydrogen envelope. The mass of the M-star is then M_M = 0.28 +/- 0.05 solar masses, and the inclination is i = 77 +/- 7 degrees. We discuss briefly how PG 1224+309 may be used to constrain theories of close binary star evolution, and the past and future histories of PG 1224+309 itself. The star is both a ``post-common envelope'' star and a ``pre-cataclysmic binary'' star. Mass transfer by Roche-lobe overflow should commence in about 10 Gyr.

The sdB star KPD 0422+5421 was discovered to be a single-lined spectroscopic binary with a period of P=0.0901795 +/- (3\times 10^{-7}) days (2 hours, 10 minutes). The U and B light curves display an ellipsoidal modulation with amplitudes of about 0.02 magnitudes. The sdB star contributes nearly all of the observed flux. This and the absence of any reflection effect suggest that the unseen companion star is small (i.e. R_comp ~ 0.01 solar radii) and therefore degenerate. We modeled the U and B light curves and derived i = 78.05 +/- 0.50 degrees and a mass ratio of q = M_comp/M_sdB = 0.87 +/- 0.15. The sdB star fills 69% of its Roche lobe. These quantities may be combined with the mass function of the companion (f(M) = 0.126 +/- 0.028 solar masses) to derive M_sdB = 0.72 +/- 0.26 solar masses and M_comp = 0.62 +/- 0.18 solar masses. We used model spectra to derive the effective temperature, surface gravity, and helium abundance of the sdB star. We found T_eff = 25,000 +/- 1500K, log g = 5.4 +/- 0.1, and [He/H] = -1.0. With a period of 2 hours and 10 minutes, KPD 0422+5421 has one of the shortest known orbital periods of a detached binary. This system is also one of only a few known binaries which contain a subdwarf B star and a white dwarf. Thus KPD 0422+5421 represents a relatively unobserved, and short-lived, stage of binary star evolution.

We have investigated 52 RXTE pointed observations of GRO J1655-40 spanning the X-ray outburst that commenced on 1996 April 25 and lasted for 16 months. Our X-ray timing analyses reveal four types of QPOs: three with relatively stable central frequencies at 300 Hz, 9 Hz, and 0.1 Hz, and a fourth that varied over the range 14-28 Hz. The 300 Hz and 0.1 Hz QPOs appear only at the highest observed luminosities (Lx > 0.15 Ledd), where the power-law component dominates the X-ray spectrum. At lower luminosity, the disk flux exceeds the power-law flux and only two of the QPOs are observed: the spectrally "soft" 9 Hz QPO, and the narrow, "hard" QPO that varies from 14-28 Hz as the hard flux decreases. The 300 Hz QPO is likely to be analogous to the stationary QPO at 67 Hz seen in the microquasar GRS1915+105. We discuss models of these high-frequency QPOs which depend on effects due to general relativity. The 9 Hz QPO displays a spectrum consistent with a thermal origin, but this frequency does not appear to be consistent with any of the natural time scales associated with the disk, or with the inferred values of the mass and rapid spin of the black hole. The mechanism for the 14-28 Hz QPOs appears to be linked to the power-law component, as do the 1-10 Hz QPOs in GRS1915+105. Finally, we show data for GRO J1655-40 and GRS1915+105 as each source teeters between relative stability and a state of intense oscillations at 0.1 Hz. A comparison of the sources' spectral parameters allows us to speculate that the black hole mass in GRS1915+105 is very large, possibly in the range 39-70 Msun.

(shortened) Spectroscopic observations of the soft X-ray transient 4U 1543-47 reveal a radial velocity curve with a period of P=1.123 +/- 0.008 days and a semi-amplitude of K_2 = 124 +/- 4 km/sec. The mass function is f(M) = 0.22 +/- 0.02 solar masses. We derive a distance of d = 9.1 +/-1.1 kpc if the secondary is on the main sequence. The V and I light curves exhibit two waves per orbital cycle with amplitudes of about 0.08 mag. We modeled the light curves as ellipsoidal variations in the secondary star and derive extreme inclination limits of 20 <= i <= 40 deg and formal 3 sigma limits of 24 <= i <= 36 deg for a mass ratio Q = M_1/M_2 > 1. However, there are systematic effects in the data that the model does not account for, so the above constraints should be treated with caution. We argue that the secondary star is still on the main sequence and if the secondary star has a mass near the main sequence values for early A-stars (2.3 <= M_2 <= 2.6 solar masses), then the best fits for the 3 sigma inclination range (24 <= i <= 36 deg) and the 3 sigma mass function range (0.16 <= f(M) <= 0.28 solar masses) imply a primary mass in the range 2.7 <= M_1 <= 7.5 solar masses. Thus the mass of the compact object in 4U 1543-47 is likely to be in excess of approximately 3 solar masses and we conclude 4U 1543-47 most likely contains a black hole.

We examine the distribution of masses of black holes in transient low mass X-ray binary systems. A Bayesian analysis suggests that it is probable that six of the seven systems with measured mass functions have black hole masses clustered near seven solar masses. There appears to be a significant gap between the masses of these systems and those of the observed neutron stars. The remaining source, V404 Cyg, has a mass significantly larger than the others, and our analysis suggests that it is probably drawn from a different distribution. Selection effects do not appear to play a role in producing the observed mass distribution, which may be explained by currently unknown details of the supernova explosions and of binary evolution prior to the supernova.

The two galactic `microquasars' with superluminal radio jets have been quite active during 1996, generatinga variety of studies involving both NASA and ground-based observatories. GRS~1915+105 has displayed dramatic accretion instability in observations with RXTE, revealing X-ray light curves and emission states unlike anything previously seen. Variable QPOs in the range of 0.07--10 Hz have been monitored with the capability to track the individual oscillations. The QPO amplitude is as high as 40% of the mean flux, while both amplitude and phase lag increase with photon energy. The results imply a direct link between the QPO mechanism and the origin of the energetic electrons believed to radiate the X-ray power-law component. GRS1915+105 also displays a transient yet stationary QPO at 67 Hz. The other source, GRO J1655-40, is an optically established black hole binary. Recent optical reports include an excellent model for the binary inclination and masses, while an optical precursor to the April 1996 X-ray outburst has been measured. We report new results from recent RXTE observations. GROJ1655-40 has displayed both the canonical ``soft/high'' and ``very high'' X-ray states, with QPOs at 8--22 Hz during the latter state. In addition, there is a high-frequency QPO at 300 Hz. The rapid oscillations in these sources are suspected of providing a measure of the mass and rotation of the accreting black holes, although several competing models may be applied when evaluating the results.

The All Sky Monitor on the Rossi X-ray Timing Explorer detected an X-ray (2-12 keV) outburst from the black hole binary GRO J1655-40 beginning near April 25, 1996. Optical photometry obtained April 20-24, 1996 shows a steady brightening of the source in B, V, R, and I beginning about six days before the start of the X-ray outburst. The onset of the optical brightening was earliest in I and latest in B. However, the rate of the optical brightening was fastest in B and slowest in I. The order of the increases in the different optical filters suggests that the event was an "outside-in" disturbance of the accretion disk. The substantial delay between the optical rise and the rise of the X-rays may provide indirect support for the advection-dominated accretion flow model of the inner regions of the accretion disk.

We report photometric and spectroscopic observations of the black hole binary GRO J1655-40 in complete quiescence. In contrast to the 1995 photometry, the light curves from 1996 are almost completely dominated by ellipsoidal modulations from the secondary star. Model fits to the light curves, which take into account the temperature profile of the accretion disk and eclipse effects, yield an inclination of i=69.50 +/- 0.08 degrees and a mass ratio of Q=2.99 +/- 0.08. The precision of our determinations of i and Q allow us to determine the black hole mass to an accuracy of approximately 4% (7.02 +/- 0.22 solar masses). The secondary star's mass is 2.34 +/- 0.12 solar masses. The position of the secondary on the Hertzsprung-Russell diagram is consistent with that of a 2.3 solar mass star which has evolved off the main sequence and is halfway to the start of the giant branch. Using the new spectra we present an improved value of the spectroscopic period (2.62157 +/- 0.00015 days), radial velocity semiamplitude (228.2 +/- 2.2 km/sec), and mass function (3.24 +/- 0.09 solar masses). Evolutionary models suggest an average mass transfer rate for such a system of 3.4E-9 solar masses per year (2.16E+17 grams per second), which is much larger than the average mass transfer rates implied in the other six transient black hole systems, but still barely below the critical mass transfer rate required for stability.