High-energy particle transport in 3D hydrodynamic models of colliding-wind binaries

Massive stars in binary systems (as WR140, WR147 or $\eta$ Carinae) have long been regarded as potential sources of high-energy $\gamma$-rays. The emission is thought to arise in the region where the stellar winds collide and produce relativistic particles which subsequently might be able to emit $\gamma$-rays. Detailed numerical hydrodynamic simulations have already offered insight in the complex dynamics of the wind collision region (WCR), while independent analytical studies, albeit with simplified descriptions of the WCR, have shed light on the spectra of charged particles. In this paper, we describe a combination of these two approaches. We present a 3D-hydrodynamical model for colliding stellar winds and compute spectral energy distributions of relativistic particles for the resulting structure of the WCR. The hydrodynamic part of our model incorporates the line-driven acceleration of the winds, gravity, orbital motion and the radiative cooling of the shocked plasma. In our treatment of charged particles we consider diffusive shock acceleration in the WCR and the subsequent cooling via inverse Compton losses (including Klein-Nishina effects), bremsstrahlung, collisions and other energy loss mechanisms.Comment: 28 pages, 9 figures / accepted for publication in The Astrophysical Journa

High-energy gamma-ray observations of the accreting black hole V404 Cygni during its June 2015 outburst

We report on Fermi/Large Area Telescope observations of the accreting black hole low-mass X-ray binary V404 Cygni during its outburst in June-July 2015. Detailed analyses reveal a possible excess of $\gamma$-ray emission on 26 June 2015, with a very soft spectrum above $100$ MeV, at a position consistent with the direction of V404 Cyg (within the $95\%$ confidence region and a chance probability of $4 \times 10^{-4}$). This emission cannot be associated with any previously-known Fermi source. Its temporal coincidence with the brightest radio and hard X-ray flare in the lightcurve of V404 Cyg, at the end of the main active phase of its outburst, strengthens the association with V404 Cyg. If the $\gamma$-ray emission is associated with V404 Cyg, the simultaneous detection of $511\,$keV annihilation emission by INTEGRAL requires that the high-energy $\gamma$ rays originate away from the corona, possibly in a Blandford-Znajek jet. The data give support to models involving a magnetically-arrested disk where a bright $\gamma$-ray jet can re-form after the occurrence of a major transient ejection seen in the radio.Comment: 5 pages, 3 figures, accepted for publication in MNRA

The Brightest Black Holes

I suggest that there are two classes of ultraluminous X-ray sources (ULXs), corresponding to super-Eddington mass inflow in two situations: (a) thermal-timescale mass transfer in high-mass X-ray binaries, and (b) long-lasting transient outbursts in low-mass X-ray binaries. These two classes are exemplified by SS433 and microquasars like GRS 1915+105 respectively. The observed ULX population is a varying mixture of the two, depending on the star formation history of the host galaxy. ULXs in galaxies with vigorous star formation (such as the Antennae) are generally SS433--like, while ULXs in elliptical galaxies must be of the microquasar type. The latter probably have significantly anisotropic radiation patterns. They should also be variable, but demonstrating this may require observations over decades. The close analogy between models of X-ray binaries and active galactic nuclei (AGN) suggests that there should exist an apparently super-Eddington class of the latter, which may be the ultrasoft AGN, and a set of X-ray binaries with Doppler--boosted X-ray emission. These are presumably a subset of the ULXs, but remain as yet unidentified.Comment: 4 pages, no figures; accepted for MNRAS Letter

INTEGRAL observation of hard X-ray variability of the TeV binary LS5039 / RX J1826.2-1450

LS 5039/RX J1826.2-1450 is one of the few High Mass X-ray binary systems from which radio and high energy TeV emission has been observed. Moreover, variability of the TeV emission with orbital period was detected. We investigate the hard X-ray (25 - 200keV) spectral and timing properties of the source with the monitoring IBIS/ISGRI instrument on-board the INTEGRAL satellite. We present the analysis of INTEGRAL observations for a total of about 3 Msec exposure time, including both public data and data from the Key Programme. We search for flux and spectral variability related to the orbital phase. The source is observed to emit from 25 up to 200 keV and the emission is concentrated around inferior conjunction. Orbital variability in the hard X-ray band is detected and established to be in phase with the orbitally modulated TeV emission observed with H.E.S.S. For this energy range we determine an average flux for the inferior conjunction phase interval of $(3.54 \pm 2.30) \times 10^{-11}$ erg cm$^{-2}$ s$^{-1}$, and a flux upper limit for the superior conjunction phase interval of $1.45 \times 10^{-11}$ erg cm$^{-2}$ s$^{-1}$ (90% conf. level respectively). The spectrum for the inferior conjunction phase interval follows a power law with an index $\Gamma = 2.0^{+0.2}_{-0.2}$ (90% conf. level).Comment: 4 pages, 4 figures, accepted by A&

Optical and Infrared Light Curves of the Eclipsing X-ray Binary V395 Car = 2S 0921-630

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

Long time-scale variability in GRS1915+105

We present very high resolution hydrodynamical simulations of accretion discs in black hole X-ray binaries accreting near the Eddington limit. The results show that mass loss, irradiation and tidal interactions all have a profound effect on the observed behaviour of long period X-ray transients. In particular, the interplay of all of these effects in the outer regions of the accretion disc is able to drive long time-scale (weeks to years) variability is these objects, and is a possible origin for some of the extreme variability of GRS1915+105.Comment: 7 pages, 9 figures (2 in colour), accepted for publication in MNRA

Broad and Luminous [OIII] and [NII] in Globular Cluster ULXs

We consider an accretion-disc origin for the broad and luminous forbidden-line emission observed in ultraluminous X-ray (ULX) sources CXOJ033831.8-352604 and XMMU 122939.7+075333 in globular clusters hosted by elliptical galaxies NGC 1399 and NGC 4472, respectively. We will refer to the latter by the globular cluster name RZ2109. The first has strong [OIII] and [NII], the second only [OIII]. Both H$\alpha$ and H$\beta$ are very weak or undetected in both objects. We assume that the large line widths are due to Keplerian rotation around a compact object and derive expressions for maximum line luminosities. These idealized models require central masses $\gtrsim100$ and \gtrsim30000\Msun for CXOJ033831.8-352604 and RZ2109, respectively. An independent, bootstrap argument for the total disc mass yields, for both systems, M_{\mathrm{disc}}\gtrsim10^{-4}\Msun for a purely metallic disc (and two orders of magnitude larger for solar metallicities). If Roche-lobe overflow is implicated, viscous time-scales are $\gtrsim300$ yr. Standard disc theory then offers another limit on the central masses. Lobe radii for a \sim1\Msun donor are $\gtrsim10^{13}$ cm. We therefore rule out Roche-lobe overflow of a white dwarf in both systems. Red giants could fill the necessary lobes. Whether they are too metal-poor to produce the strong forbidden lines without strong hydrogen emission is unclear.Comment: Accepted to MNRAS Letters, 5 pages, 6 figure

The origin and fate of short-period low-mass black-hole binaries

We present results of a population synthesis study for semidetached short orbital period binaries which contain low-mass(<1.5 Msun) donors and black hole (>4 Msun) accretors. Evolution of these binaries is determined by nuclear evolution of the donors and/or orbital angular momentum loss due to magnetic braking by the stellar wind of the donors and gravitational wave radiation. According to our model, the estimated total number of this type of black-hole binaries in the Galaxy is about 10000. If the magnetic braking is described by the Verbunt & Zwaan formula, the model predicts around 3000 transient systems with periods >2 hours and around 300 luminous stable systems with periods between 3 and 8 hours. Several dozens of these bright systems should be above the RXTE ASM sensitivity limit. The absence of such systems implies that angular momentum losses are reduced by a factor more than 2 with respect to the Verbunt & Zwaan prescription. We show that it is unlikely that the transient behaviour of black-hole short-period X-ray binaries is explained by the evolved nature of the stellar companion. A substantial fraction of black-hole binaries with periods >3 hours could be faint with truncated, stable cold accretion discs as proposed by Menou et al. Most of the semidetached black-hole binaries are expected to have periods shorter than ~2 hours. Properties of such, still to be observed, very small mass-ratio (q<0.02) binaries are different from those of their longer period cousins.Comment: 13 pages, 6 figures, accepted for publication in A&

Revisiting a fundamental test of the disc instability model for X-ray binaries

We revisit a core prediction of the disc instability model (DIM) applied to X-ray binaries. The model predicts the existence of a critical mass transfer rate, which depends on disc size, separating transient and persistent systems. We therefore selected a sample of 52 persistent and transient neutron star and black hole X-ray binaries and verified if observed persistent (transient) systems do lie in the appropriate stable (unstable) region of parameter space predicted by the model. We find that, despite the significant uncertainties inherent to these kinds of studies, the data are in very good agreement with the theoretical expectations. We then discuss some individual cases that do not clearly fit into this main conclusion. Finally, we introduce the transientness parameter as a measure of the activity of a source and show a clear trend of the average outburst recurrence time to decrease with transientness in agreement with the DIM predictions. We therefore conclude that, despite difficulties in reproducing the complex details of the lightcurves, the DIM succeeds to explain the global behaviour of X-ray binaries averaged over a long enough period of time.Comment: 12 pages, 4 figures. Accepted for publication in MNRAS. Version 2: some typos corrected and references adde