Near-Infrared Observations of Compact Binary Systems

Low mass X-ray binaries (LMXBs) are a subset of compact binary systems in which a main-sequence or slightly evolved star fills its Roche lobe and donates mass to a neutron star or a black hole (BH) via an accretion disk. Robust estimates of compact object masses in these systems are required to enhance our current understanding of the physics of compact object formation, accretion disks and jets. Compact object masses are typically determined at near-infrared (NIR) wavelengths when the system is in quiescence and the donor star is the dominant source of flux. Previous studies have assumed that any non-stellar contribution at these wavelengths is minimal. However, this assumption is rarely true. By performing NIR spectroscopy, we determined the fractional donor star contribution to the NIR flux and the compact object masses in two LMXBs: V404 Cyg and Cen X-4. In our analysis, it was assumed that the light curve morphology remains consistent throughout quiescence. It has now been shown in several systems that veiling measurements from non-stellar sources are meaningful only if acquired contemporaneously with light curve measurements. We accounted for this in the measurement of the BH mass in the LMXB, XTE J1118+480. LMXBs are also considered to be the most likely candidates responsible for the formation of milli-second pulsars (MSP). Here, I present the unique case of PSR J1903+0327 that challenges this currently accepted theory of MSP formation and is a potential candidate for testing General Relativity. Observations in the NIR come with their own set of challenges. NIR detector arrays used in these observations generally have high dark current and readout noise. In an effort to lower the read noise in NICFPS at APO, we present a study done on the Hawaii--1RG engineering grade chip that served as a test bed for reducing the read noise in NICFPS

The Mass of the Black Hole in XTE J1118+480

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

Near-Infrared Spectroscopy of Low Mass X-ray Binaries : Accretion Disk Contamination and Compact Object Mass Determination in V404 Cyg and Cen X-4

We present near-infrared (NIR) broadband (0.80--2.42 $\mu$m) spectroscopy of two low mass X-ray binaries: V404 Cyg and Cen X-4. One important parameter required in the determination of the mass of the compact objects in these systems is the binary inclination. We can determine the inclination by modeling the ellipsoidal modulations of the Roche-lobe filling donor star, but the contamination of the donor star light from other components of the binary, particularly the accretion disk, must be taken into account. To this end, we determined the donor star contribution to the infrared flux by comparing the spectra of V404 Cyg and Cen X-4 to those of various field K-stars of known spectral type. For V404 Cyg, we determined that the donor star has a spectral type of K3 III. We determined the fractional donor contribution to the NIR flux in the H- and K-bands as $0.98 \pm .05$ and $0.97 \pm .09$, respectively. We remodeled the H-band light curve from \citet{sanwal1996} after correcting for the donor star contribution to obtain a new value for the binary inclination. From this, we determined the mass of the black hole in V404 Cyg to be $M_{BH}= 9.0^{+.2}_{-.6}M_{\odot}$. We performed the same spectral analysis for Cen X-4 and found the spectral type of the donor star to be in the range K5 -- M1V. The donor star contribution in Cen X-4 is $0.94\pm.14$ in the H-band while in the K-band, the accretion disk can contribute up to 10% of the infrared flux. We remodeled the H-band light curve from \citet{shahbaz1993}, again correcting for the fractional contribution of the donor star to obtain the inclination. From this, we determined the mass of the neutron star as $M_{NS}= 1.5^{+.1}_{-.4}M_{\odot}$. However, the masses obtained for both systems should be viewed with some caution since contemporaneous light curve and spectral data are required to obtain definitive masses

The Mass Distribution of Stellar-Mass Black Holes

We perform a Bayesian analysis of the mass distribution of stellar-mass black holes using the observed masses of 15 low-mass X-ray binary systems undergoing Roche lobe overflow and five high-mass, wind-fed X-ray binary systems. Using Markov Chain Monte Carlo calculations, we model the mass distribution both parametrically---as a power law, exponential, gaussian, combination of two gaussians, or log-normal distribution---and non-parametrically---as histograms with varying numbers of bins. We provide confidence bounds on the shape of the mass distribution in the context of each model and compare the models with each other by calculating their relative Bayesian evidence as supported by the measurements, taking into account the number of degrees of freedom of each model. The mass distribution of the low-mass systems is best fit by a power-law, while the distribution of the combined sample is best fit by the exponential model. We examine the existence of a "gap" between the most massive neutron stars and the least massive black holes by considering the value, M_1%, of the 1% quantile from each black hole mass distribution as the lower bound of black hole masses. The best model (the power law) fitted to the low-mass systems has a distribution of lower-bounds with M_1% > 4.3 Msun with 90% confidence, while the best model (the exponential) fitted to all 20 systems has M_1% > 4.5 Msun with 90% confidence. We conclude that our sample of black hole masses provides strong evidence of a gap between the maximum neutron star mass and the lower bound on black hole masses. Our results on the low-mass sample are in qualitative agreement with those of Ozel, et al (2010).Comment: 56 pages, 22 figures, 9 tables, as accepted by Ap

Correlated optical, X-ray, and $-ray flaring activity seen with INTEGRAL during the 2015 outburst of V404 Cygni

Reproduced with permission from Astronomy & Astrophysics. © 2015 ESO.After 25 years of quiescence, the microquasar V404 Cyg entered a new period of activity in June 2015. This X-ray source is known to undergo extremely bright and variable outbursts seen at all wavelengths. It is therefore an object of prime interest to understand the accretion-ejection connections. These can, however, only be probed through simultaneous observations at several wavelengths. We made use of the INTEGRAL instruments to obtain long, almost uninterrupted observations from 2015 June 20th, 15:50 UTC to June 25th, 4:05 UTC, from the optical V-band, up to the soft γ-rays. V404 Cyg was extremely variable in all bands, with the detection of 18 flares with fluxes exceeding 6 Crab (20--40 keV) within 3 days. The flare recurrence can be as short as ∼ 20~min from peak to peak. A model-independent analysis shows that the >6 Crab flares have a hard spectrum. A simple 10--400 keV spectral analysis of the off-flare and flare periods shows that the variation in intensity is likely to be due to variations of a cut-off power law component only. The optical flares seem to be at least of two different types: one occurring in simultaneity with the X-ray flares, the other showing a delay greater than 10 min. The former could be associated with X-ray reprocessing by either an accretion disk or the companion star. We suggest that the latter are associated with plasma ejections that have also been seen in radio.Peer reviewe

BlackCAT: A catalogue of stellar-mass black holes in X-ray transients

During the last ~50 years, the population of black hole candidates in X-ray binaries has increased considerably with 59 Galactic objects detected in transient low-mass X-ray binaries, plus a few in persistent systems (including ~5 extragalactic binaries). We collect near-infrared, optical and X-ray information spread over hundreds of references in order to study the population of black holes in X-ray transients as a whole. We present the most updated catalogue of black hole transients, which contains X-ray, optical and near-infrared observations together with their astrometric and dynamical properties. It provides new useful information in both statistical and observational parameters providing a thorough and complete overview of the black hole population in the Milky Way. Analysing the distances and spatial distribution of the observed systems, we estimate a total population of ~1300 Galactic black hole transients. This means that we have already discovered less than ~5% of the total Galactic distribution. The complete version of this catalogue will be continuously updated online and in the Virtual Observatory, including finding charts and data in other wavelengths.Comment: http://www.astro.puc.cl/BlackCAT - Accepted for publication in Astronomy & Astrophysics. 20 pages, 8 figures, 5 Table

The radio/X-ray domain of black hole X-ray binaries at the lowest radio luminosities

We report on deep, coordinated radio and X-ray observations of the black hole X-ray binary XTE J1118+480 in quiescence. The source was observed with the Karl G. Jansky Very Large Array for a total of 17.5 h at 5.3 GHz, yielding a 4.8 ± 1.4 μJy radio source at a position consistent with the binary system. At a distance of 1.7 kpc, this corresponds to an integrated radio luminosity between 4 and 8 × 1025 erg s−1, depending on the spectral index. This is the lowest radio luminosity measured for any accreting black hole to date. Simultaneous observations with the Chandra X-ray Telescope detected XTE J1118+480 at 1.2 × 10−14 erg s−1 cm−2 (1–10 keV), corresponding to an Eddington ratio of ~4 × 10−9 for a 7.5 M☉ black hole. Combining these new measurements with data from the 2005 and 2000 outbursts available in the literature, we find evidence for a relationship of the form lr = α+βlX (where l denotes logarithmic luminosities), with β = 0.72 ± 0.09. XTE J1118+480 is thus the third system – together with GX339-4 and V404 Cyg – for which a tight, non-linear radio/X-ray correlation has been reported over more than 5 dex in lX. Confirming previous results, we find no evidence for a dependence of the correlation normalization of an individual system on orbital parameters, relativistic boosting, reported black hole spin and/or black hole mass. We then perform a clustering and linear regression analysis on what is arguably the most up-to-date collection of coordinated radio and X-ray luminosity measurements from quiescent and hard-state black hole X-ray binaries, including 24 systems. At variance with previous results, a two-cluster description is statistically preferred only for random errors <~0.3 dex in both lr and lX, a level which we argue can be easily reached when the known spectral shape/distance uncertainties and intrinsic variability are accounted for. A linear regression analysis performed on the whole data set returns a best-fitting slope β = 0.61 ± 0.03 and intrinsic scatter σ0 = 0.31 ± 0.03 dex

Positron annihilation signatures associated with the outburst of the microquasar V404 Cygni

This document is the Accepted Manuscript version of the following article: Thomas Siegert, et al, ‘Positron annihilation signatures associated with the outburst of the microquasar V404 Cygni’, Nature: International Journal of Science, Vol. 531: 341-343, March 2016, DOI: https://doi.org/10.1038/nature16978. Content in the UH Research Archive is made available for personal research, educational, and non-commercial purposes only. Unless otherwise stated, all content is protected by copyright, and in the absence of an open license, permissions for further re-use should be sought from the publisher, the author, or other copyright holder.Microquasars1, 2, 3, 4 are stellar-mass black holes accreting matter from a companion star5 and ejecting plasma jets at almost the speed of light. They are analogues of quasars that contain supermassive black holes of 106 to 1010 solar masses. Accretion in microquasars varies on much shorter timescales than in quasars and occasionally produces exceptionally bright X-ray flares6. How the flares are produced is unclear, as is the mechanism for launching the relativistic jets and their composition. An emission line near 511 kiloelectronvolts has long been sought in the emission spectrum of microquasars as evidence for the expected electron–positron plasma. Transient high-energy spectral features have been reported in two objects7, 8, but their positron interpretation9 remains contentious. Here we report observations of γ-ray emission from the microquasar V404 Cygni during a recent period of strong flaring activity10. The emission spectrum around 511 kiloelectronvolts shows clear signatures of variable positron annihilation, which implies a high rate of positron production. This supports the earlier conjecture that microquasars may be the main sources of the electron–positron plasma responsible for the bright diffuse emission of annihilation γ-rays in the bulge region of our Galaxy11. Additionally, microquasars could be the origin of the observed megaelectronvolt continuum excess in the inner Galaxy.Peer reviewe

Multiwavelength Observations of A0620-00 in Quiescence

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