Variability in high-mass X-ray binaries

Strongly magnetized, accreting neutron stars show periodic and aperiodic variability over a wide range of time scales. By obtaining spectral and timing information on these different time scales, we can have a closer look into the physics of accretion close to the neutron star and the properties of the accreted material. One of the most prominent time scales is the strong pulsation, i.e., the rotation period of the neutron star itself. Over one rotation, our view of the accretion column and the X-ray producing region changes significantly. This allows us to sample different physical conditions within the column but at the same time requires that we have viewing-angle-resolved models to properly describe them. In wind-fed high-mass X-ray binaries, the main source of aperiodic variability is the clumpy stellar wind, which leads to changes in the accretion rate (i.e., luminosity) as well as absorption column. This variability allows us to study the behavior of the accretion column as a function of luminosity, as well as to investigate the structure and physical properties of the wind, which we can compare to winds in isolated stars.Comment: 6 pages, 4 figures, accepted for publication in Astronomische Nachrichten (proceedings of the XMM-Newton Workshop 2019

GeV detection of HESS J0632+057

HESS J0632+057 is the only gamma-ray binary that has been detected at TeV energies, but not at GeV energies yet. Based on nearly nine years of Fermi Large Area Telescope (LAT) Pass 8 data, we report here on a deep search for the gamma-ray emission from HESS J0632+057 in the 0.1-300 GeV energy range. We find a previously unknown gamma-ray source, Fermi J0632.6+0548, spatially coincident with HESS J0632+057. The measured flux of Fermi J0632.6+0548 is consistent with the previous flux upper limit on HESS J0632+057 and shows variability that can be related to the HESS J0632+057 orbital phase. We propose that Fermi J0632.6+0548 is the GeV counterpart of HESS J0632+057. Considering the Very High Energy (VHE) spectrum of HESS J0632+057, a possible spectral turnover above 10 GeV may exist in Fermi J0632.6+0548, as appears to be common in other established gamma-ray binaries.Comment: 17 pages, 4 figures, 1 table; Accepted for publication in Ap

Investigation of the energy dependence of the orbital light curve in LS 5039

LS 5039 is so far the best studied $\gamma$-ray binary system at multi-wavelength energies. A time resolved study of its spectral energy distribution (SED) shows that above 1 keV its power output is changing along its binary orbit as well as being a function of energy. To disentangle the energy dependence of the power output as a function of orbital phase, we investigated in detail the orbital light curves as derived with different telescopes at different energy bands. We analysed the data from all existing \textit{INTEGRAL}/IBIS/ISGRI observations of the source and generated the most up-to-date orbital light curves at hard X-ray energies. In the $\gamma$-ray band, we carried out orbital phase-resolved analysis of \textit{Fermi}-LAT data between 30 MeV and 10 GeV in 5 different energy bands. We found that, at $\lesssim$100 MeV and $\gtrsim$1 TeV the peak of the $\gamma$-ray emission is near orbital phase 0.7, while between $\sim$100 MeV and $\sim$1 GeV it moves close to orbital phase 1.0 in an orbital anti-clockwise manner. This result suggests that the transition region in the SED at soft $\gamma$-rays (below a hundred MeV) is related to the orbital phase interval of 0.5--1.0 but not to the one of 0.0--0.5, when the compact object is "behind" its companion. Another interesting result is that between 3 and 20 GeV no orbital modulation is found, although \textit{Fermi}-LAT significantly ($\sim$18$\sigma$) detects LS 5039. This is consistent with the fact that at these energies, the contributions to the overall emission from the inferior conjunction phase region (INFC, orbital phase 0.45 to 0.9) and from the superior conjunction phase region (SUPC, orbital phase 0.9 to 0.45) are equal in strength. At TeV energies the power output is again dominant in the INFC region and the flux peak occurs at phase $\sim$0.7.Comment: 7 pages, 6 figures, accepted for publication in MNRA

Multiple cyclotron line-forming regions in GX 301-2

We present two observations of the high-mass X-ray binary GX 301-2 with NuSTAR, taken at different orbital phases and different luminosities. We find that the continuum is well described by typical phenomenological models, like a very strongly absorbed NPEX model. However, for a statistically acceptable description of the hard X-ray spectrum we require two cyclotron resonant scattering features (CRSF), one at ~35 keV and the other at ~50 keV. Even though both features strongly overlap, the good resolution and sensitivity of NuSTAR allows us to disentangle them at >=99.9% significance. This is the first time that two CRSFs are seen in GX 301-2. We find that the CRSFs are very likely independently formed, as their energies are not harmonically related and, if it were a single line, the deviation from a Gaussian shape would be very large. We compare our results to archival Suzaku data and find that our model also provides a good fit to those data. We study the behavior of the continuum as well as the CRSF parameters as function of pulse phase in seven phase bins. We find that the energy of the 35 keV CRSF varies smoothly as function of phase, between 30-38 keV. To explain this variation, we apply a simple model of the accretion column, taking the altitude of the line-forming region, the velocity of the in-falling material, and the resulting relativistic effects into account. We find that in this model the observed energy variation can be explained simply due to a variation of the projected velocity and beaming factor of the line forming region towards us.Comment: 18 pages, 10 figures, accepted for publication in A&

A large spin-up rate measured with INTEGRAL in the High Mass X-ray Binary Pulsar SAXJ2103.5+4545

The High Mass X-ray Binary Pulsar SAXJ2103.5+4545 has been observed with INTEGRAL several times during the last outburst in 2002-2004. We report a comprehensive study of all INTEGRAL observations, allowing a study of the pulse period evolution during the recent outburst. We measured a very rapid spin-up episode, lasting 130days, which decreased the pulse period by 1.8s. The spin-up rate, pdot=-1.5e-7 s/s, is the largest ever measured for SAXJ2103.5+4545, and it is among the fastest for an accreting pulsar. The pulse profile shows evidence for temporal variability, apparently not related to the source flux or to the orbital phase. The X-ray spectrum is hard and there is significant emission up to 150keV. A new derivation of the orbital period, based on RXTE data, is also reported.Comment: 8 pages, 7 figures, accepted for publication in A&

Multi-wavelength INTEGRAL NEtwork (MINE) observations of the microquasar GRS 1915+105

We present the international collaboration MINE (Multi-lambda Integral NEtwork) aimed at conducting multi-wavelength observations of X-ray binaries and microquasars simultaneously with the INTEGRAL gamma-ray satellite. We will focus on the 2003 March-April campaign of observations of the peculiar microquasar GRS 1915+105 gathering radio, IR and X-ray data. The source was observed 3 times in the plateau state, before and after a major radio and X-ray flare. It showed strong steady optically thick radio emission corresponding to powerful compact jets resolved in the radio images, bright near-infrared emission, a strong QPO at 2.5 Hz in the X-rays and a power law dominated spectrum without cutoff in the 3-300 keV range. We compare the different observations, their multi-wavelength light curves, including JEM-X, ISGRI and SPI, and the parameters deduced from fitting the spectra obtained with these instruments on board INTEGRAL.Comment: 4 pages, 9 fig., Proc. of the 5th INTEGRAL Workshop (Feb. 16-20 2004), to be published by ES