The iron line and the continuum-fitting methods are currently the two leading techniques for measuring black hole spins with electromagnetic radiation. They can be naturally extended for probing the spacetime geometry around black holes and testing general relativity in the strong field regime. In the past couple of years, there has been significant work to use the iron line method to test the nature of black holes. In this Letter, we use the continuum-fitting method and we show its capability of constraining the spacetime geometry around black holes by analyzing 17 RXTE data of the X-ray binary LMC X-1.

In high mass X-ray binaries (HMXBs), an accreting compact object orbits a high mass star which loses mass through a dense and inhomogeneous wind. Using the compact object as an X-ray backlight, the time variability of the absorbing column density in the wind can be exploited in order to shed light on the micro-structure of the wind and obtain unbiased stellar mass loss rates for high mass stars. We explore the impact of clumpiness on the variability of the column density with a simplified wind model. In particular, we focus on the standard deviation of the column density and the characteristic duration of enhanced absorption episodes, and compare them with analytical predictions based on the porosity length. We identified the favorable systems and orbital phases to determine the wind micro-structure. The coherence time scale of the column density is shown to be the self-crossing time of a clump in front of the compact object. We provide a recipe to get accurate measurements of the size and of the mass of the clumps, purely based on the observable time variability of the column density. The coherence time scale grants direct access to the size of the clumps while their mass can be deduced separately from the amplitude of the variability. If it is due to unaccreted passing-by clumps, the high column density variations in some HMXBs requires high mass clumps to reproduce the observed peak-to-peak amplitude and coherence time scales. These clump properties are hardly compatible with the ones derived from first principles. Alternatively, other components could contribute to the variability of the column density: larger orbital scale structures produced by a mechanism still to be identified, or a dense environment in the immediate vicinity of the accretor such as an accretion disk, an outflow or a spherical shell around the magnetosphere of the accreting neutron star.

The analysis of the thermal spectrum of geometrically thin and optically thick accretion disks of black holes, the so-called continuum-fitting method, is one of the leading techniques for measuring black hole spins. Current models normally approximate the disk as infinitesimally thin, while in reality the disk thickness is finite and increases as the black hole mass accretion rate increases. Here we present an XSPEC model to calculate the multi-temperature blackbody spectrum of a thin accretion disk of finite thickness around a Kerr black hole. We test our new model with an RXTE observation of the black hole binary GRS 1915+105. We find that the spin value inferred with the new model is slightly higher than the spin value obtained with a model with an infinitesimally thin disk, but the difference is small and the effect is currently subdominant with respect to other sources of uncertainties in the final spin measurement.

We explore the accretion properties of the black hole X-ray binary \j1550\ during its outbursts in 1998/99 and 2000. We model the disk, corona, and reflection components of X-ray spectra taken with the {\it Rossi X-ray Timing Explorer} (\rxte), using the {\tt relxill} suite of reflection models. The key result of our modeling is that the reflection spectrum in the very soft state is best explained by disk self-irradiation, i.e., photons from the inner disk are bent by the strong gravity of the black hole, and reflected off the disk surface. This is the first known detection of thermal disk radiation reflecting off the inner disk. There is also an apparent absorption line at $\sim6.9$ keV which may be evidence of an ionized disk wind. The coronal electron temperature ($kT_{\rm e}$) is, as expected, lower in the brighter outburst of 1998/99, explained qualitatively by more efficient coronal cooling due to irradiating disk photons. The disk inner radius is consistent with being within a few times the innermost stable circular orbit (ISCO) throughout the bright-hard-to-soft states (10s of $r_{\rm g}$ in gravitational units). The disk inclination is low during the hard state, disagreeing with the binary inclination value, and very close to $90^{\circ}$ in the soft state, recovering to a lower value when adopting a blackbody spectrum as the irradiating continuum.

We present the analysis of several observations of the black hole binary GX 339--4 during its bright intermediate states from two different outbursts (2002 and 2004), as observed by RXTE/PCA. We perform a consistent study of its reflection spectrum by employing the relxill family of relativistic reflection models to probe the evolutionary properties of the accretion disk including the inner disk radius ($R_{\rm in}$), ionization parameter ($\xi$), temperatures of the inner disk ($T_{\rm in}$), corona ($kT_{\rm e}$), and its optical depth ($\tau$). Our analysis indicates that the disk inner edge approaches the inner-most stable circular orbit (ISCO) during the early onset of bright hard state, and that the truncation radius of the disk remains low ($< 9 R_{\rm g}$) throughout the transition from hard to soft state. This suggests that the changes observed in the accretion disk properties during the state transition are driven by variation in accretion rate, and not necessarily due to changes in the inner disk's radius. We compare the aforementioned disk properties in two different outbursts, with state transitions occurring at dissimilar luminosities, and find identical evolutionary trends in the disk properties, with differences only seen in corona's $kT_{\rm e}$ and $\tau$. We also perform an analysis by employing a self-consistent Comptonized accretion disk model accounting for the scatter of disk photons by the corona, and measure low inner disk truncation radius across the bright intermediate states, using the temperature dependent values of spectral hardening factor, thereby independently confirming our results from the reflection spectrum analysis.

We present the spectral and timing evolution of the persistent black hole X-ray binary GRS 1758-258 based on almost 12 years of observations using the Rossi X-ray Timing Explorer Proportional Counter Array. While the source was predominantly found in the hard state during this time, it entered the thermally dominated soft state seven times. In the soft state GRS 1758-258 shows a strong decline in flux above 3 keV rather than the pivoting flux around 10 keV more commonly shown by black hole transients. In its 3-20 keV hardness intensity diagram, GRS 1758-258 shows a hysteresis of hard and soft state fluxes typical for transient sources in outburst. The RXTE-PCA and RXTE-ASM long-term light curves do not show any orbital modulations in the range of 2 to 30 d. However, in the dynamic power spectra significant peaks drift between 18.47d and 18.04d for the PCA data, while less significant signatures between 19d and 20d are seen for the ASM data as well as for the Swift/BAT data. We discuss different models for the hysteresis behavior during state transitions as well as possibilities for the origin of the long term variation in the context of a warped accretion disk.

We present an analysis of X-ray spectra of the high mass X-ray binary 4U 0115+634 as observed with Suzaku and RXTE in 2011 July, during the fading phase of a giant X-ray outburst. We used a continuum model consisting of an absorbed cutoff power-law and an ad-hoc Gaussian emission feature centered around 8.5 keV, which we discuss to be due to cyclotron emission. Our results are consistent with a fundamental cyclotron absorption line centered at ${\sim}10.2$ keV for all observed flux ranges. At the same time we rule out significant influence of the 8.5 kev Gaussian on the CRSF parameters, which are not consistent with the cyclotron line energies and depths of previously reported flux-dependent descriptions. We also show that some continuum models can lead to artificial line-like residuals in the analyzed spectra, which are then misinterpreted as unphysically strong cyclotron lines. Specifically, our results do not support the existence of a previously claimed additional cyclotron feature at ${\sim}15$ keV. Apart from these features, we find for the first time evidence for a He-like Fe XXV emission line at ${\sim}6.7$ keV and weak H-like Fe XXVI emission close to ${\sim}7.0$ keV.

With the opening of the X and gamma--ray windows in the sixties, thanks to to sounding rockets and satellite-borne instruments, extremely energetic and violent phenomena were discovered and subsequently found to be ubiquitous in the Universe. Observations in the high energy domain are fundamental for understanding how matter is organized and behaves around black holes; unravelling how these extreme objects influence their environments on a very large scale; and finding the still elusive obscured massive objects in the centre of galaxies. Other major problems in contemporary astrophysics, such as the understanding of acceleration processes at shocks of all sizes (those of pulsar wind nebulae, supernova remnants, but also at larger scales those of Active Galactic Nuclei radio lobes) in relation to the origin of cosmic-rays, or the definitive characterization of the debated non-thermal X-ray energy content of clusters of galaxies, also requires observations at very high energies. An observatory type medium mission operating from around 1 keV to about 600 keV can provide direct insights into these major questions. The essential characteristics will be coverage of the full energy range by telescopes featuring a large throughput and arc-second resolution optics, coupled to a compact focal plane assembly, with excellent imaging resolution and spectroscopy. In addition, the mission will provide unique polarimetry measurements in the hard X-ray domain, an important new diagnostic tool at energies for which the non-thermal processes dominate. The Polarimetric High-Energy Modular Telescope Observatory (PHEMTO) is designed to have performance several orders of magnitude better than the present hard X-ray instruments. This gives to PHEMTO the improvements in scientific performance needed for a mission in the 2050 era.

Accreting black holes are unique tools to understand the physics under extreme gravity. While black hole X-ray binaries differ vastly in mass from AGN, their accretion and ejection flows are assumed to be essentially similar. Black hole X-ray binaries or microquasars are, however, quasars for the impatient as variability timescales scale directly with mass. State changes, i.e., strong variations in emission properties, in black hole X-ray binaries can happen within hours and whole outburst cycles within months to years. But our understanding of the drivers of such changes and the contributions of individual accretion and ejection components to the overall emission is still lacking. Here, I highlight some of the INTEGRAL's unique contributions to the understanding of black hole X-ray binaries through its coverage of the energies above the spectral cutoff, its long uninterrupted monitoring observations and the measurements of hard X-ray / soft $\gamma$-ray polarization.

We propose the development of X-ray interferometry (XRI), to reveal the universe at high energies with ultra-high spatial resolution. With baselines which can be accommodated on a single spacecraft, XRI can reach 100 $\mu$as resolution at 10 \AA (1.2 keV) and 20 $\mu$as at 2 \AA (6 keV), enabling imaging and imaging-spectroscopy of (for example) X-ray coronae of nearby accreting supermassive black holes (SMBH) and the SMBH `shadow'; SMBH accretion flows and outflows; X-ray binary winds and orbits; stellar coronae within ~100 pc and many exoplanets which transit across them. For sufficiently luminous sources XRI will resolve sub-pc scales across the entire observable universe, revealing accreting binary SMBHs and enabling trigonometric measurements of the Hubble constant with X-ray light echoes from quasars or explosive transients. A multi-spacecraft `constellation' interferometer would resolve well below 1 $\mu$as, enabling SMBH event horizons to be resolved in many active galaxies and the detailed study of the effects of strong field gravity on the dynamics and emission from accreting gas close to the black hole.

We report on the spectroscopic analysis of the black hole binary GX 339-4 during its recent 2017-2018 outburst, observed simultaneously by the Swift and NuSTAR observatories. Although during this particular outburst the source failed to make state transitions, and despite Sun constraints during the peak luminosity, we were able to trigger four different observations sampling the evolution of the source in the hard state. We show that even for the lowest luminosity observations the NuSTAR spectra show clear signatures of X-ray reprocessing (reflection) in an accretion disk. Detailed analysis of the highest signal-to-noise spectra with our family of relativistic reflection models RELXILL indicates the presence of both broad and narrow reflection components. We find that a dual-lamppost model provides a superior fit when compared to the standard single lamppost plus distant neutral reflection. In the dual lamppost model two sources at different heights are placed on the rotational axis of the black hole, suggesting that the narrow component of the Fe K emission is likely to originate in regions far away in the disk, but still significantly affected by its rotational motions. Regardless of the geometry assumed, we find that the inner edge of the accretion disk reaches a few gravitational radii in all our fits, consistent with previous determinations at similar luminosity levels. This confirms a very low degree of disk truncation for this source at luminosities above ~1% Eddington. Our estimates of Rin reinforces the suggested behavior for an inner disk that approaches the inner-most regions as the luminosity increases in the hard state.

XTE J1550-564 is a black hole X-ray binary for which the dynamical characteristics are well established, and the broadband spectral evolution of the source has been well studied. Its orbital inclination is known to be high, at $\sim75^{\circ}$, with the jet estimated to align well with the orbital axis. We explore simultaneous observations made with ASCA and RXTE, covering the $1$--$200$~keV band, during the early stages of the first outburst of XTE J1550-564 in its hard-intermediate state, on 1998-09-23/24. We show that the most up-to-date reflection models, applied to these data, yield an inclination estimate much lower than found in previous studies, at $\sim40^{\circ}$, grossly disagreeing with the dynamically estimated orbital inclination. We discuss the possible explanations for this disagreement and its implications for reflection models, including possible physical scenarios in which either the inner disk is misaligned both with binary orbit and the outer jet, or either the inner accretion flow, corona, and/or jet have vertical structure which leads to lower inferred disk inclination through various physical means.

The X-ray Polarization Probe (XPP) is a second generation X-ray polarimeter following up on the Imaging X-ray Polarimetry Explorer (IXPE). The XPP will offer true broadband polarimetery over the wide 0.2-60 keV bandpass in addition to imaging polarimetry from 2-8 keV. The extended energy bandpass and improvements in sensitivity will enable the simultaneous measurement of the polarization of several emission components. These measurements will give qualitatively new information about how compact objects work, and will probe fundamental physics, i.e. strong-field quantum electrodynamics and strong gravity.

Vela X-1 is an eclipsing high mass X-ray binary (HMXB) consisting of a 283s accreting X-ray pulsar in a close orbit of 8.964 days around the B0.5Ib supergiant HD77581 at a distance of just 2.4 kpc. The system is considered a prototype of wind-accreting HMXB and it has been used as a baseline in different theoretical or modelling studies. We discuss the observational properties of the system and the use of the observational data as laboratory to test recent developments in modelling the accretion process in High-Mass X-ray Binaries (e.g., Sander et al. 2018; El Mellah et al. 2018), which range from detailed descriptions of the wind acceleration to modelling of the structure of the flow of matter close to the neutron star and its variations.

We present an analysis of three Chandra High Energy Transmission Gratings observations of the black hole binary Cyg X-1/HDE 226868 at different orbital phases. The stellar wind that is powering the accretion in this system is characterized by temperature and density inhomogeneities including structures, or "clumps", of colder, more dense material embedded in the photoionized gas. As these clumps pass our line of sight, absorption dips appear in the light curve. We characterize the properties of the clumps through spectral changes during various dip stages. Comparing the silicon and sulfur absorption line regions (1.6-2.7 keV $\equiv$ 7.7-4.6 {\AA}) in four levels of varying column depth reveals the presence of lower ionization stages, i.e., colder or denser material, in the deeper dip phases. The Doppler velocities of the lines are roughly consistent within each observation, varying with the respective orbital phase. This is consistent with the picture of a structure that consists of differently ionized material, in which shells of material facing the black hole shield the inner and back shells from the ionizing radiation. The variation of the Doppler velocities compared to a toy model of the stellar wind, however, does not allow us to pin down an exact location of the clump region in the system. This result, as well as the asymmetric shape of the observed lines, point at a picture of a complex wind structure.

We utilize the joint capabilities of IBIS and SPI to perform a state-resolved 20--2000\,keV analysis of the microquasar Cygnus X-1. In both LHS and HSS, the spectral analysis reveals the presence of a high-energy tail above 400\,keV in addition to a standard, Compton-like continuum in the 20--400\,keV range. We study the polarisation properties of the hard X-ray radiation, as well as of the polarisation of this high-energy tail. Polarisation is detected in the LHS in agreement with previous work. We find potential variations of the LHS polarisation with time and no polarisation detection in the HSS.

High-resolution X-ray spectroscopy has advanced our understanding of the hot Universe by revealing physical properties like kinematics, temperature, and abundances of the astrophysical plasmas. Despite the technical and scientific achievements, the lack of scientific products at a level higher than count spectra is hampering full scientific exploitation of high-quality data. This paper introduces the Catalogue of Ionized Emission Lines Observed by the Reflection Grating Spectrometer (CIELO-RGS) onboard the XMM-Newton space observatory. The CIELO-RGS catalogue aims to facilitate the exploitation of emission features in the public RGS spectra archive, in particular, to perform the correlation between X-ray spectral diagnostics parameters with measurements at other wavelengths. This paper focuses on the methodology of catalogue generation, describing the automated line detection algorithm. A moderate sample (~2400 observations) of high-quality RGS spectra available at XMM-Newton Science Archive is used as our starting point. A list of potential emission lines is selected based on a multi-scale peak detection algorithm in a uniform and automated way without prior assumption on the underlying astrophysical model. The candidate line list is validated via spectral fitting with simple continuum and line profile models. We also compare the catalogue content with published literature results on a small number of exemplary sources. We generate a catalogue of emission lines ~12000 detected in ~1600 observations toward stars, X-ray binaries, supernovae remnants, active galactic nuclei, and groups and clusters of galaxies. For each line, we report the observed wavelength, broadening, energy and photon flux, and equivalent width, etc.

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.

Thanks to high-resolution and non-dispersive spectrometers onboard future X-ray missions such as XRISM and Athena, we are finally poised to answer important questions about the formation and evolution of galaxies and large-scale structure. However, we currently lack an adequate understanding of many atomic processes behind the spectral features we will soon observe. Large error bars on parameters as critical as transition energies and atomic cross sections can lead to unacceptable uncertainties in the calculations of e.g., elemental abundance, velocity, and temperature. Unless we address these issues, we risk limiting the full scientific potential of these missions. Laboratory astrophysics, which comprises theoretical and experimental studies of the underlying physics behind observable astrophysical processes, is therefore central to the success of these missions.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

GX~339$-$4 is a black hole X-ray binary that is a key focus of accretion studies since it goes into outburst roughly every two-to-three years. Tracking of its radio, IR and X-ray flux during multiple outbursts reveals tight broadband correlations. The radio emission originates in a compact, self-absorbed jet, however the origin of the X-ray emission is still debated: jet base or corona? We fit 20 quasi-simultaneous radio, IR, optical and X-ray observations of GX~339$-$4 covering three separate outbursts in 2005, 2007, 2010--2011, with a composite corona + jet model, where inverse Compton emission from both regions contributes to the X-ray emission. Using a recently-proposed identifier of the X-ray variability properties known as power-spectral hue, we attempt to explain both the spectral and evolving timing characteristics, with the model. We find the X-ray spectra are best fit by inverse Compton scattering in a dominant hot corona ($kT_{\rm e}\sim$ hundreds of keV). However, radio and IR-optical constraints imply a non-negligible contribution from inverse Compton scattering off hotter electrons ($kT_{\rm e} \ge 511$~keV) in the base of the jets, ranging from a few up to $\sim50$\% of the integrated 3--100~keV flux. We also find that the physical properties of the jet show interesting correlations with the shape of the broadband X-ray variability of the source, posing intriguing suggestions for the connection between the jet and corona.

In this paper we describe the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to accretion flows in the strong field gravity regime around both stellar-mass and supermassive black-holes. eXTP has the unique capability of using advanced 'spectral-timing-polarimetry' techniques to analyze the rapid variations with three orthogonal diagnostics of the flow and its geometry, yielding unprecedented insight into the inner accreting regions, the effects of strong field gravity on the material within them and the powerful outflows which are driven by the accretion process.

We present the analysis of an extended \textit{INTEGRAL} dataset of the high-mass microquasar Cygnus X-1. We first classify, in a model-independent way, all the \textit{INTEGRAL} individual pointings taken between 2003 and 2016 in three basic spectral states. This, in particular, allows us to triple the exposure time of the soft state in comparison with previous publication. We then study the spectral properties of the 5--400 keV stacked spectra of the soft and hard states and provide the parameters obtained with our modelling. Using a refined alternative method of extracting the Compton double events of the IBIS telescope, we then extract high-energy ($>$400 keV) spectra in the two states. We do detect an hard tail in both states. Our refined analysis allows us to obtain a hard state (count) spectrum at a flux lower than previously published by our team. Although a full estimate of the calibration property of this improved software is still needed, this seems to be more inline with the hard state hard tail seen with other instruments.

We present a detailed spectral analysis of the Black Hole Binary XTE J1752-223 in the hard state of its 2009 outburst. Regular monitoring of this source by RXTE provided high signal-to-noise spectra along the outburst rise and decay. During one full month this source stalled at $\sim$30\% of its peak count rate at a constant hardness and intensity. By combining all the data in this exceptionally-stable hard state, we obtained an aggregate PCA spectrum (3-45 keV) with 100 million counts, and a corresponding HEXTE spectrum (20-140 keV) with 5.8 million counts. Implementing a version of our reflection code with a physical model for Comptonization, we obtain tight constraints on important physical parameters for this system. In particular, the inner accretion disk is measured very close in, at $R_\mathrm{in}=1.7\pm0.4$ $R_g$. Assuming $R_\mathrm{in}=R_\mathrm{ISCO}$, we find a relatively high black hole spin ($a_*=0.92\pm0.06$). Imposing a lamppost geometry, we obtain a low inclination ($i=35\pm4$ deg), which agrees with the upper limit found in the radio ($i<49$ deg). However, we note that this model cannot be statistically distinguished from a non-lamppost model with free emissivity index, for which the inclination is markedly higher. Additionally, we find a relatively cool corona ($57-70$ keV), and large iron abundance ($3.3-3.7$ solar). We further find that properly accounting for Comptonization of the reflection emission improves the fit significantly and causes an otherwise low reflection fraction ($\sim 0.2-0.3$) to increase by an order of magnitude, in line with geometrical expectations for a lamppost corona. We compare these results with similar investigations reported for GX 339-4 in its bright hard state.

Here we study a 1-200 keV energy spectrum of the black hole binary Cygnus X-1 taken with NuSTAR and Suzaku. This is the first report of a NuSTAR observation of Cyg X-1 in the intermediate state, and the observation was taken during the part of the binary orbit where absorption due to the companion's stellar wind is minimal. The spectrum includes a multi-temperature thermal disk component, a cutoff power-law component, and relativistic and non-relativistic reflection components. Our initial fits with publicly available constant density reflection models (relxill and reflionx) lead to extremely high iron abundances (>9.96 and 10.6(+1.6)(-0.9) times solar, respectively). Although supersolar iron abundances have been reported previously for Cyg X-1, our measurements are much higher and such variability is almost certainly unphysical. Using a new version of reflionx that we modified to make the electron density a free parameter, we obtain better fits to the spectrum even with solar iron abundances. We report on how the higher density (n_e = (3.98(+0.12)(-0.25))E20 cm-3) impacts other parameters such as the inner radius and inclination of the disk.

We present results of multiband optical photometry of the black hole X-ray binary system V404 Cygni obtained using Wheaton College Observatory's 0.3m telescope, along with strictly simultaneous INTEGRAL and Swift observations during 2015 June 25.15--26.33 UT, and 2015 June 27.10--27.34 UT. These observations were made during the 2015 June outburst of the source when it was going through an epoch of violent activity in all wavelengths ranging from radio to $\gamma$-rays. The multiwavelength variability timescale favors a compact emission region, most likely originating in a jet outflow, for both observing epochs presented in this work. The simultaneous INTEGRAL/Imager on Board the Integral Satellite (IBIS) 20--40 keV light curve obtained during the June 27 observing run correlates very strongly with the optical light curve, with no detectable delay between the optical bands as well as between the optical and hard X-rays. The average slope of the dereddened spectral energy distribution was roughly flat between the $I_C$- and $V$-bands during the June 27 run, even though the optical and X-ray flux varied by $>$25$\times$ during the run, ruling out an irradiation origin for the optical and suggesting that the optically thick to optically thin jet synchrotron break during the observations was at a frequency larger than that of $V$-band, which is quite extreme for X-ray binaries. These observations suggest that the optical emission originated very close to the base of the jet. A strong H$\alpha$ emission line, probably originating in a quasi-spherical nebula around the source, also contributes significantly in the $R_C$-band. Our data, in conjunction with contemporaneous data at other wavelengths presented by other groups, strongly suggest that the jet-base was extremely compact and energetic during this phase of the outburst.

Bright and eclipsing, the high-mass X-ray binary Vela X-1 offers a unique opportunity to study accretion onto a neutron star from clumpy winds of O/B stars and to disentangle the complex accretion geometry of these systems. In Chandra-HETGS spectroscopy at orbital phase ~0.25, when our line of sight towards the source does not pass through the large-scale accretion structure such as the accretion wake, we observe changes in overall spectral shape on timescales of a few kiloseconds. This spectral variability is, at least in part, caused by changes in overall absorption and we show that such strongly variable absorption cannot be caused by unperturbed clumpy winds of O/B stars. We detect line features from high and low ionization species of silicon, magnesium and neon whose strengths and presence depend on the overall level of absorption. They imply a co-existence of cool and hot gas phases in the system that we interpret as a highly variable, structured accretion flow close to the compact object such as has been recently seen in simulations of wind accretion in high-mass X-ray binaries.

We present a spectral and timing study of the NuSTAR and Swift observations of the black hole candidate IGR J17091-3624 in the hard state during its outburst in 2016. Disk reflection is detected in each of the NuSTAR spectra taken in three epochs. Fitting with relativistic reflection models reveals that the accretion disk is truncated during all epochs with $R_{\rm in}>10~r_{\rm g}$, with the data favoring a low disk inclination of $\sim 30^{\circ}-40^{\circ}$. The steepening of the continuum spectra between epochs is accompanied by a decrease in the high energy cut-off: the electron temperature $kT_{\rm e}$ drops from $\sim 64$ keV to $\sim 26$ keV, changing systematically with the source flux. We detect type-C QPOs in the power spectra with frequency varying between 0.131 Hz and 0.327 Hz. In addition, a secondary peak is found in the power spectra centered at about 2.3 times the QPO frequency during all three epochs. The nature of this secondary frequency is uncertain, however a non-harmonic origin is favored. We investigate the evolution of the timing and spectral properties during the rising phase of the outburst and discuss their physical implications.

In September 2016, the microquasar Cygnus X-3 underwent a giant radio flare, which was monitored for 6 days with the Medicina Radio Astronomical Station and the Sardinia Radio Telescope. Long observations were performed in order to follow the evolution of the flare on a hourly scale, covering six frequency ranges from 1.5 GHz to 25.6 GHz. The radio emission reached a maximum of 13.2 +/- 0.7 Jy at 7.2 GHz and 10 +/- 1 Jy at 18.6 GHz. Rapid flux variations were observed at high radio frequencies at the peak of the flare, together with rapid evolution of the spectral index: alpha steepened from 0.3 to 0.6 within 5 hours. This is the first time that such fast variations are observed, giving support to the evolution from optically thick to optically thin plasmons in expansion moving outward from the core. Based on the Italian network (Noto, Medicina and SRT) and extended to the European antennas (Torun, Yebes, Onsala), VLBI observations were triggered at 22 GHz on five different occasions, four times prior to the giant flare, and once during its decay phase. Flux variations of 2-hour duration were recorded during the first session. They correspond to a mini-flare that occurred close to the core ten days before the onset of the giant flare. From the latest VLBI observation we infer that four days after the flare peak the jet emission was extended over 30 mas.

The extragalactic VHE gamma-ray sky is rich in blazars. These are jetted active galactic nuclei viewed at a small angle to the line-of-sight. Only a handful of objects viewed at a larger angle are known so far to emit above 100 GeV. Multi-wavelength studies of such objects up to the highest energies provide new insights into the particle and radiation processes of active galactic nuclei. We report the results from the first multi-wavelength campaign observing the TeV detected nucleus of the active galaxy IC 310, whose jet is observed at a moderate viewing angle of 10 deg - 20 deg. The multi-instrument campaign was conducted between 2012 Nov. and 2013 Jan., and involved observations with MAGIC, Fermi, INTEGRAL, Swift, OVRO, MOJAVE and EVN. These observations were complemented with archival data from the AllWISE and 2MASS catalogs. A one-zone synchrotron self-Compton model was applied to describe the broad-band spectral energy distribution. IC 310 showed an extraordinary TeV flare at the beginning of the campaign, followed by a low, but still detectable TeV flux. Compared to previous measurements, the spectral shape was found to be steeper during the low emission state. Simultaneous observations in the soft X-ray band showed an enhanced energy flux state and a harder-when-brighter spectral shape behaviour. No strong correlated flux variability was found in other frequency regimes. The broad-band spectral energy distribution obtained from these observations supports the hypothesis of a double-hump structure. The harder-when-brighter trend in the X-ray and VHE emission is consistent with the behaviour expected from a synchrotron self-Compton scenario. The contemporaneous broad-band spectral energy distribution is well described with a one-zone synchrotron self-Compton model using parameters that are comparable to those found for other gamma-ray-emitting misaligned blazars.

We have measured the energies of the strongest 1s-2ell (ell=s,p) transitions in He- through Ne-like silicon and sulfur ions to an accuracy of better than 1eV using Lawrence Livermore National Laboratory's electron beam ion traps, EBIT-I and SuperEBIT, and the NASA/GSFC EBIT Calorimeter Spectrometer (ECS). We identify and measure the energies of 18 and 21 X-ray features from silicon and sulfur, respectively. The results are compared to new Flexible Atomic Code calculations and to semi-relativistic Hartree Fock calculations by Palmeri et al. (2008). These results will be especially useful for wind diagnostics in high mass X-ray binaries, such as Vela X-1 and Cygnus X-1, where high-resolution spectral measurements using Chandra's high energy transmission grating has made it possible to measure Doppler shifts of 100km/s. The accuracy of our measurements is consistent with that needed to analyze Chandra observations, exceeding Chandra's 100km/s limit. Hence, the results presented here not only provide benchmarks for theory, but also accurate rest energies that can be used to determine the bulk motion of material in astrophysical sources. We show the usefulness of our results by applying them to redetermine Doppler shifts from Chandra observations of Vela X-1.

We present an analysis of NuSTAR observations of a hard intermediate state of the transient black hole GX 339-4 taken in January 2015. As the source softened significantly over the course of the 1.3 d-long observation we split the data into 21 sub-sets and find that the spectrum of all of them can be well described by a power-law continuum with an additional relativistically blurred reflection component. The photon index increases from ~1.69 to ~1.77 over the course of the observation. The accretion disk is truncated at around 9 gravitational radii in all spectra. We also perform timing analysis on the same 21 individual data sets, and find a strong type-C quasi-periodic oscillation (QPO), which increase in frequency from ~0.68 to ~1.05 Hz with time. The frequency change is well correlated with the softening of the spectrum. We discuss possible scenarios for the production of the QPO and calculate predicted inner radii in the relativistic precession model as well as the global disk mode oscillations model. We find discrepancies with respect to the observed values in both models unless we allow for a black hole mass of ~100 M_sun , which is highly unlikely. We discuss possible systematic uncertainties, in particular with the measurement of the inner accretion disk radius in the relativistic reflection model. We conclude that the combination of observed QPO frequencies and inner accretion disk radii, as obtained from spectral fitting, is difficult to reconcile with current models.

We present a multi-epoch hard X-ray analysis of Cygnus X-1 in its soft state based on four observations with NuSTAR. Despite the basic similarity of the observed spectra, there is clear spectral variability between epochs. To investigate this variability, we construct a model incorporating both the standard disk-corona continuum and relativistic reflection from the accretion disk, based on prior NuSTAR work on Cygnus X-1, and apply this model to each epoch independently. We find excellent consistency for the black hole spin, and the iron abundance of the accretion disk, which are expected to remain constant on observational timescales. In particular, we confirm that Cygnus X-1 hosts a rapidly rotating black hole, 0.93<a*<0.96, in broad agreement with the majority of prior studies of the relativistic disk reflection and constraints on the spin obtained through studies of the thermal accretion disk continuum. Our work also confirms the apparent misalignment between the inner disk and the orbital plane of the binary system reported previously, finding the magnitude of this warp to be ~10-15deg. This level of misalignment does not significantly change (and may even improve) the agreement between our reflection results and the thermal continuum results regarding the black hole spin. The spectral variability observed by NuSTAR is dominated by the primary continuum, implying variability in the temperature of the scattering electron plasma. Finally, we consistently observe absorption from ionized iron at ~6.7 keV, which varies in strength as a function of orbital phase in a manner consistent with the absorbing material being an ionized phase of the focused stellar wind from the supergiant companion star.

We present a spectral analysis of a brief Chandra/HETG observation of the neutron star low-mass X-ray binary GX~340+0. The high-resolution spectrum reveals evidence of ionized absorption in the Fe K band. The strongest feature, an absorption line at approximately 6.9 keV, is required at the 5 sigma level of confidence via an F-test. Photoionization modeling with XSTAR grids suggests that the line is the most prominent part of a disk wind with an apparent outflow speed of v = 0.04c. This interpretation is preferred at the 4 sigma level over a scenario in which the line is H-like Fe XXVI at a modest red-shift. The wind may achieve this speed owing to its relatively low ionization, enabling driving by radiation pressure on lines; in this sense, the wind in GX 340+0 may be the stellar-mass equivalent of the flows in broad absorption line quasars (BALQSOs). If the gas has a unity volume filling factor, the mass ouflow rate in the wind is over 10^-5 Msun/year, and the kinetic power is nearly 10^39 erg/s (or, 5-6 times the radiative Eddington limit for a neutron star). However, geometrical considerations - including a small volume filling factor and low covering factor - likely greatly reduce these values.

The accretion onto the black hole in the system HDE 226868/Cygnus X-1 is powered by the strong line driven stellar wind of the O-type donor star. We study the X-ray properties of the stellar wind in the hard state of Cyg X-1 as determined with data from the Chandra High Energy Transmission Gratings. Large density and temperature inhomogeneities are present in the wind, with a fraction of the wind consisting of clumps of matter with higher density and lower temperature embedded in a photoionized gas. Absorption dips observed in the light curve are believed to be caused by these clumps. This work concentrates on the non-dip spectra as a function of orbital phase. The spectra show lines of H-like and He-like ions of S, Si, Na, Mg, Al and highly ionized Fe (Fe xvii-Fe xxiv). We measure velocity shifts, column densities, and thermal broadening of the line series. The excellent quality of these five observations allows us to investigate the orbital phase dependence of these parameters. We show that the absorber is located close to the black hole. Doppler shifted lines point at a complex wind structure in this region, while emission lines seen in some observations are from a denser medium than the absorber. The observed line profiles are phase dependent. Their shapes vary from pure, symmetric absorption at the superior conjunction to P Cygni profiles at the inferior conjunction of the black hole.

We present results from spectral fitting of the very high state of GX~339-4 with NuSTAR and Swift. We use relativistic reflection modelling to measure the spin of the black hole and inclination of the inner disk, and find a spin of $a=0.95^{+0.02}_{-0.08}$ and inclination of $30${\deg}$\pm1$ (statistical errors). These values agree well with previous results from reflection modelling. With the exceptional sensitivity of NuSTAR at the high-energy side of the disk spectrum, we are able to constrain multiple physical parameters simultaneously using continuum fitting. By using the constraints from reflection as input for the continuum fitting method, we invert the conventional fitting procedure to estimate the mass and distance of GX~339-4 using just the X-ray spectrum, finding a mass of $9.0^{+1.6}_{-1.2}M_\odot$ and distance of $8.4\pm0.9$ kpc (statistical errors).

We have performed a full time- and luminosity-resolved spectral analysis of the high-mass X-ray binary 4U 1538-522 using the available RXTE, INTEGRAL, and Suzaku data, examining both phase-averaged and pulse-phase-constrained datasets and focusing on the behavior of the cyclotron resonance scattering feature (CRSF). No statistically significant trend between the energy of the CRSF and luminosity is observed in the combined dataset. However, the CRSF energy appears to have increased by ~1.5 keV in the ~8.5 years between the RXTE and Suzaku measurements, with Monte Carlo simulations finding the Suzaku measurement 4.6$\sigma$ above the RXTE points. Interestingly, the increased Suzaku CRSF energy is much more significant and robust in the pulse-phase-constrained spectra from the peak of the main pulse, suggesting a change that is limited to a single magnetic pole. The 7 years of RXTE measurements do not show any strongly-significant evolution with time on their own. We discuss the significance of the CRSF's behavior with respect to luminosity and time in the context of historical observations of this source as well as recent observational and theoretical work concerning the neutron star accretion column, and suggest some mechanisms by which the observed change over time could occur.

We report on the timing and spectral analysis of two Suzaku observations with different flux levels of the high-mass X-ray binary KS 1947+300 during its 2013 outburst. In agreement with simultaneous NuSTAR observations, the continuum is well described by an absorbed power law with a cut-off and an additional black body component. In addition we find fluorescent emission from neutral, He-like, and even H-like iron. We determine a pulse period of ~18.8 s with the source showing a spin-up between the two observations. Both Suzaku observations show a very similar behavior of the pulse profile, which is strongly energy dependent, with an evolution from a profile with one peak at low energies to a profile with two peaks of different widths towards higher energies seen in both, the Suzaku and NuSTAR data. Such an evolution to a more complex profile at higher energies is rarely seen in X-ray pulsars, most cases show the opposite behavior. Pulse phase-resolved spectral analysis shows a variation in the absorbing column density, NH , over pulse phase. Spectra taken during the pulse profile minima are intrinsically softer compared to the pulse phase-averaged spectrum.

We report on the analysis of the broad Fe Kalpha line feature of Cygnus X-1 in the spectra of four simultaneous hard intermediate state observations made with the X-ray Multiple Mirror mission (XMM-Newton), the Rossi X-ray Timing Explorer (RXTE), and the International Gamma-Ray Astrophysics Laboratory (INTEGRAL). The high quality of the XMM-Newton data taken in the Modified Timing Mode of the EPIC-pn camera provides a great opportunity to investigate the broadened Fe Kalpha reflection line at 6.4keV with a very high signal to noise ratio. The 4-500keV energy range is used to constrain the underlying continuum and the reflection at higher energies. We first investigate the data by applying a phenomenological model that consists of the sum of an exponentially cutoff power law and relativistically smeared reflection. Additionally, we apply a more physical approach and model the irradiation of the accretion disk directly from the lamp post geometry. All four observations show consistent values for the black hole parameters with a spin of $a\sim 0.9$, in agreement with recent measurements from reflection and disk continuum fitting. The inclination is found to be $i\sim30^\circ$, consistent with the orbital inclination and different from inclination measurements made during the soft state, which show a higher inclination. We speculate that the difference between the inclination measurements is due to changes in the inner region of the accretion disk.

We have developed a correction tool to improve the quality of RXTE HEXTE spectra by employing the same method we used earlier to improve the quality of RXTE PCA spectra. We fit all of the hundreds of HEXTE spectra of the Crab individually to a simple power-law model, some 37 million counts in total for Cluster A and 39 million counts for Cluster B, and we create for each cluster a combined spectrum of residuals. We find that the residual spectrum of Cluster A is free of instrumental artifacts while that of Cluster B contains significant features with amplitudes ~1%; the most prominent is in the energy range 30-50 keV, which coincides with the iodine K edge. Starting with the residual spectrum for Cluster B, via an iterative procedure we created the calibration tool hexBcorr for correcting any Cluster B spectrum of interest. We demonstrate the efficacy of the tool by applying it to Cluster B spectra of two bright black holes, which contain several million counts apiece. For these spectra, application of the tool significantly improves the goodness of fit, while affecting only slightly the broadband fit parameters. The tool may be important for the study of spectral features, such as cyclotron lines, a topic that is beyond the scope of this paper.

We studied time variability and spectral evolution of the Galactic black hole transient Swift J174510.8-262411 during the first phase of its outburst. INTEGRAL and Swift observations collected from 2012 September 16 until October 30 have been used. The total squared fractional rms values did not drop below 5% and QPOs, when present, were type-C, indicating that the source never made the transition to the soft-intermediate state. Even though the source was very bright (up to 1 Crab in hard X-rays), it showed a so called failed outburst as it never reached the soft state. XRT and IBIS broad band spectra, well represented by a hybrid thermal/non-thermalComptonisationmodel, showed physical parameters characteristic of the hard and intermediate states. In particular, the derived temperature of the geometrically thin disc black body was about 0.6 keV at maximum.We found a clear decline of the optical depth of the corona electrons (close to values of 0.1), as well as of the total compactness ratio lh/ls. The hard-to-hard/intermediate state spectral transition is mainly driven by the increase in the soft photon flux in the corona, rather than small variations of the electron heating. This, associated with the increasing of the disc temperature, is consistent with a disc moving towards the compact object scenario, i.e. the truncated-disc model. Moreover, this scenario is consistent with the decreasing fractional squared rms and increasing of the noise and QPO frequency. In our final group of observations, we found that the contribution from the non-thermal Comptonisation to the total power supplied to the plasma is 0.59+0.02/-0.05 and that the thermal electrons cool to kTe<26 keV.

We present a timing and spectral analysis of the X-ray pulsar XTE J1946+274 observed with Suzaku during an outburst decline in 2010 October and compare with previous results. XTE J1946+274 is a transient X-ray binary consisting of a Be-type star and a neutron star with a 15.75 s pulse period in a 172 d orbit with 2-3 outbursts per orbit during phases of activity. We improve the orbital solution using data from multiple instruments. The X-ray spectrum can be described by an absorbed Fermi-Dirac cutoff power law model along with a narrow Fe K line at 6.4 keV and a weak Cyclotron Resonance Scattering Feature (CRSF) at ~35 keV. The Suzaku data are consistent with the previously observed continuum flux versus iron line flux correlation expected from fluorescence emission along the line of sight. However, the observed iron line flux is slightly higher, indicating the possibility of a higher iron abundance or the presence of non-uniform material. We argue that the source most likely has only been observed in the subcritical (non-radiation dominated) state since its pulse profile is stable over all observed luminosities and the energy of the CRSF is approximately the same at the highest (~$5 \times 10^{37} $erg s$^{-1}$) and lowest (~$5 \times 10^{36} $erg s$^{-1}$) observed 3-60 keV luminosities.

We analyze {\it simultaneously} six composite {\it RXTE} spectra of GX 339--4 in the hard state comprising 77 million counts collected over 196 ks. The source spectra are ordered by luminosity and spanthe range 1.6\% to 17\% of the Eddington luminosity. Crucially, using our new tool {\tt pcacorr}, we re-calibrate the data to a precision of 0.1\%, an order of magnitude improvement over all earlier work. Using our advanced reflection model {\tt relxill}, we target the strong features in the component of emission reflected from the disk, namely, the relativistically-broadened Fe K emission line, the Fe K edge and the Compton hump. We report results for two joint fits to the six spectra: For the first fit, we fix the spin parameter to its maximal value ($a_*=0.998$) and allow the inner disk radius $R_{\rm in}$ to vary. Results include (i) precise measurements of $R_{\rm in}$, with evidence that the disk becomes slightly truncated at a few percent of Eddington; and (ii) an order-of-magnitude swing with luminosity in the high energy cutoff, which reaches $>890$ keV at our lowest luminosity. For the second fit, we make the standard assumption in estimating spin that the inner edge of the accretion disk is located at the innermost stable circular orbit ($R_\mathrm{in} = R_\mathrm{ISCO}$) and find $a_* = 0.95^{+0.03}_{-0.05}$ (90\% confidence, statistical). For both fits, and at the same level of statistical confidence, we estimate that the disk inclination is $i = 48\pm 1$ deg and that the Fe abundance is super-solar, $A_\mathrm{Fe} = 5\pm1$.

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 20$^{\mathrm{th}}$, 15:50 UTC to June 25$^{\mathrm{th}}$, 4:05 UTC, from the optical V-band, up to the soft $\gamma$-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 $\sim$ 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.

We present spectral analysis of five NuSTAR and Swift observations of GX 339-4 taken during a failed outburst in summer 2013. These observations cover Eddington luminosity fractions in the range ~0.9-6%. Throughout this outburst, GX 339-4 stayed in the hard state, and all five observations show similar X-ray spectra with a hard power-law with a photon index near 1.6 and significant contribution from reflection. Using simple reflection models we find unrealistically high iron abundances. Allowing for different photon indices for the continuum incident on the reflector relative to the underlying observed continuum results in a statistically better fit and reduced iron abundances. With a photon index around 1.3, the input power-law on the reflector is significantly harder than that which is directly observed. We study the influence of different emissivity profiles and geometries and consistently find an improvement when using separate photon indices. The inferred inner accretion disk radius is strongly model dependent, but we do not find evidence for a truncation radius larger than 100 r_g in any model. The data do not allow independent spin constraints but the results are consistent with the literature (i.e., a>0). Our best-fit models indicate an inclination angle in the range 40-60 degrees, consistent with limits on the orbital inclination but higher than reported in the literature using standard reflection models. The iron line around 6.4 keV is clearly broadened, and we detect a superimposed narrow core as well. This core originates from a fluorescence region outside the influence of the strong gravity of the black hole and we discuss possible geometries.

We present simultaneous Nuclear Spectroscopic Telescope Array (NuSTAR ) and Suzaku observations of the X-ray binary Cygnus X-1 in the hard state. This is the first time this state has been observed in Cyg X-1 with NuSTAR, which enables us to study the reflection and broad-band spectra in unprecedented detail. We confirm that the iron line cannot be fit with a combination of narrow lines and absorption features, and instead requires a relativistically blurred profile in combination with a narrow line and absorption from the companion wind. We use the reflection models of Garcia et al. (2014) to simultaneously measure the black hole spin, disk inner radius, and coronal height in a self-consistent manner. Detailed fits to the iron line profile indicate a high level of relativistic blurring, indicative of reflection from the inner accretion disk. We find a high spin, a small inner disk radius, and a low source height, and rule out truncation to greater than three gravitational radii at the 3{\sigma} confidence level. In addition, we find that the line profile has not changed greatly in the switch from soft to hard states, and that the differences are consistent with changes in the underlying reflection spectrum rather than the relativistic blurring. We find that the blurring parameters are consistent when fitting either just the iron line or the entire broad-band spectrum, which is well modelled with a Comptonized continuum plus reflection model.

We present spectral analysis of NuSTAR and Swift observations of Cep X-4 during its outburst in 2014. We observed the source once during the peak of the outburst and once during the decay, finding good agreement in the spectral shape between the observations. We describe the continuum using a powerlaw with a Fermi-Dirac cutoff at high energies. Cep X-4 has a very strong cyclotron resonant scattering feature (CRSF) around 30 keV. A simple absorption-like line with a Gaussian optical depth or a pseudo-Lorentzian profile both fail to describe the shape of the CRSF accurately, leaving significant deviations at the red side of the line. We characterize this asymmetry with a second absorption feature around 19 keV. The line energy of the CRSF, which is not influenced by the addition of this feature, shows a small but significant positive luminosity dependence. With luminosities between (1-6)e36 erg/s, Cep X-4 is below the theoretical limit where such a correlation is expected. This behavior is similar to Vela X-1 and we discuss parallels between the two systems.

Polarization of the >~400 keV hard tail of the microquasar Cygnus X-1 has been independently reported by INTEGRAL/IBIS, and INTEGRAL/SPI and interpreted as emission from a compact jet. These conclusions were, however, based on the accumulation of all INTEGRAL data regardless of the spectral state. We utilize additional INTEGRAL exposure accumulated until December 2012, and include the AMI/Ryle (15 GHz) radio data in our study. We separate the observations into hard, soft, and intermediate/transitional states and detect radio emission from a compact jet in hard and intermediate states, but not in the soft. The 10-400 keV INTEGRAL (JEM-X and IBIS) state resolved spectra are well modeled with thermal Comptonization and reflection components. We detect a hard tail in the 0.4-2 MeV range for the hard state only. We extract the state dependent polarigrams of Cyg X-1, which all are compatible to no or undetectable level of polarization except in 400-2000 keV range in the hard state where the polarization fraction is 75$\pm$32 % and the polarization angle 40.0 +-14 deg. An upper limit on the 0.4-2 MeV soft state polarization fraction is 70%. Due to the short exposure, we obtain no meaningful constraint for the intermediate state. The likely detection of a >400 keV polarized tail in the hard state, together with the simultaneous presence of a radio jet, reinforce the notion of a compact jet origin of the 400 keV emission.

Binary systems with an accreting compact object are a unique chance to investigate the strong, clumpy, line-driven winds of early type supergiants by using the compact object's X-rays to probe the wind structure. We analyze the two-component wind of HDE 226868, the O9.7Iab giant companion of the black hole Cyg X-1 using 4.77 Ms of RXTE observations of the system taken over the course of 16 years. Absorption changes strongly over the 5.6 d binary orbit, but also shows a large scatter at a given orbital phase, especially at superior conjunction. The orbital variability is most prominent when the black hole is in the hard X-ray state. Our data are poorer for the intermediate and soft state, but show signs for orbital variability of the absorption column in the intermediate state. We quantitatively compare the data in the hard state to a toy model of a focussed Castor-Abbott-Klein-wind: as it does not incorporate clumping, the model does not describe the observations well. A qualitative comparison to a simplified simulation of clumpy winds with spherical clumps shows good agreement in the distribution of the equivalent hydrogen column density for models with a porosity length on the order of the stellar radius at inferior conjunction; we conjecture that the deviations between data and model at superior conjunction could be either due to lack of a focussed wind component in the model or a more complicated clump structure.

We fitted all of the several hundred {\it RXTE} PCA spectra of the Crab individually to a simple power-law model; the total number of counts in the composite spectrum is $>10^9$. We then used the spectrum of residuals to derive a calibration tool, called {\tt pcacorr}, that we apply to large samples of spectra for GX~339--4, H1743--322, and XTE J1550--564. Application of the tool improved the quality of all the fits, and the improvement is dramatic for spectra with $\gtrsim10^7$ counts. The Crab residual spectrum is somewhat different for each of the five PCA detectors, but it was relatively stable over the course of the mission. We recommend that {\tt pcacorr} be routinely applied to spectra with $\gtrsim10^6$ counts and that one include a systematic error of 0.1\%, rather than the 0.5--1\% value that has customarily been used. We expect that application of the tool will result in an increase in sensitivity of the PCA to faint spectral features by up to an order of magnitude.

We present the most extensive analysis of Fourier-based X-ray timing properties of the black hole binary Cygnus X-1 to date, based on 12 years of bi-weekly monitoring with RXTE from 1999 to 2011. Our aim is a comprehensive study of timing behavior across all spectral states, including the elusive transitions and extreme hard and soft states. We discuss the dependence of the timing properties on spectral shape and photon energy, and study correlations between Fourier-frequency dependent coherence and time lags with features in the power spectra. Our main results are: (a) The fractional rms in the 0.125-256 Hz range in different spectral states shows complex behavior that depends on the energy range considered. It reaches its maximum not in the hard state, but in the soft state in the Comptonized tail above 10 keV. (b) The shape of power spectra in hard and intermediate states and the normalization in the soft state are strongly energy dependent in the 2.1-15 keV range. This emphasizes the need for an energy-dependent treatment of power spectra and a careful consideration of energy- and mass-scaling when comparing the variability of different source types, e.g., black hole binaries and AGN. PSDs during extremely hard and extremely soft states can be easily confused for energies above ~5 keV in the 0.125-256 Hz range. (c) The coherence between energy bands drops during transitions from the intermediate into the soft state but recovers in the soft state. (d) The time lag spectra in soft and intermediate states show distinct features at frequencies related to the frequencies of the main variability components seen in the power spectra and show the same shift to higher frequencies as the source softens. [...abridged] In particular, we discuss how the timing properties of Cyg X-1 can be used to assess the evolution of variability with spectral shape in other black hole binaries. [...abridged]

A study of archival RXTE, Swift, and Suzaku pointed observations of the transient high mass X-ray binary GRO J1008-57 is presented. A new orbital ephemeris based on pulse arrival timing shows times of maximum luminosities during outbursts of GRO J1008-57 to be close to periastron at orbital phase -0.03. This makes the source one of a few for which outburst dates can be predicted with very high precision. Spectra of the source in 2005, 2007, and 2011 can be well described by a simple power law with high energy cutoff and an additional black body at lower energies. The photon index of the power law and the black body flux only depend on the 15-50 keV source flux. No apparent hysteresis effects are seen. These correlations allow to predict the evolution of the pulsar's X-ray spectral shape over all outbursts as a function of just one parameter, the source's flux. If modified by an additional soft component, this prediction even holds during GRO J1008-57's 2012 type II outburst.

The knowledge of the spectral state of a black hole is essential for the interpretation of data from black holes in terms of their emission models. Based on pointed observations of Cyg X-1 with the Rossi X-ray timing Explorer (RXTE) that are used to classify simultaneous RXTE-ASM observations, we develop a scheme based on RXTE -ASM colors and count rates that can be used to classify all observations of this canonical black hole that were performed between 1996 and 2011. We show that a simple count rate criterion, as used previously, leads to a significantly higher fraction of misclassified observations. This scheme enables us to classify single INTEGRAL-IBIS science windows and to obtain summed spectra for the soft, intermediate and hard state with low contamination by other states.

We present a scheme to determine the spectral state of the canonical black hole Cyg X-1 using data from previous and current X-ray all sky monitors (RXTE-ASM, Swift-BAT, MAXI, and Fermi-GBM). State determinations of the hard/intermediate and soft state agree to better than 10% between different monitors, facilitating the determination of the state and its context for any observation of the source, potentially over the lifetimes of different individual monitors. A separation of the hard and the intermediate state, which strongly differ in their spectral shape and short-term timing behavior, is only possible when monitor data in the soft X-rays (<5 keV) are available. A statistical analysis of the states confirms the different activity patterns of the source (e.g., months to years long hard state periods or phases during which numerous transitions occurs). It also shows the hard and soft states to be stable, with the probability of Cyg X-1 remaining in a given state for at least one week to be larger than 85% for the hard state and larger than 75%, for the soft state. Intermediate states are short lived, with a probability of 50% that the source leaves the intermediate state within three days. A reliable detection of these potentially short-lived events is only possible with monitor data with a time resolution of better than 1 d.

We report on an outburst of the high mass X-ray binary 4U 0115+63 with a pulse period of 3.6s in 2008 March/April as observed with RXTE and INTEGRAL. During the outburst the neutron star's luminosity varied by a factor of 10 in the 3--50\,keV band. In agreement with earlier work we find evidence for five cyclotron resonance scattering features at ~10.7, 21.8, 35.5, 46.7, and 59.7keV. Previous work had found an anti-correlation between the fundamental cyclotron line energy and the X-ray flux. We show that this apparent anti-correlation is probably due to the unphysical interplay of parameters of the cyclotron line with the continuum models used previously, e.g., the negative and positive exponent power law (NPEX). For this model, we show that cyclotron line modeling erroneously leads to describing part of the exponential cutoff and the continuum variability, and not the cyclotron lines. When the X-ray continuum is modeled with a simple exponentially cutoff power law modified by a Gaussian emission feature around 10keV, the correlation between the line energy and the flux vanishes and the line parameters remain virtually constant over the outburst. We therefore conclude that the previously reported anti-correlation is an artifact of the assumptions adopted in the modeling of the continuum.

IGR J18179-1621 is an obscured accreting X-ray pulsar discovered by INTEGRAL on 2012 February 29. We report on our 20 ksec Chandra-High Energy Transmission Gratings Spectrometer observation of the source performed on 2012 March 17, on two short contemporaneous Swift observations, and on our two near-infrared (K_s, H_n, and J_n) observations performed on 2012 March 13 and March 26. We determine the most accurate X-ray position of IGR J18179-1621, alpha_{J2000}=18 h 17m 52.18 s, delta_{J2000}=-16^\circ 21', 31.68" (90% uncertainty of 0.6"). A strong periodic variability at 11.82 s is clearly detected in the Chandra data, confirming the pulsating nature of the source, with the lightcurve softening at the pulse peak. The quasi-simultaneous Chandra-Swift spectra of IGR J18179-1621 can be well fit by a heavily absorbed hard power-law N_H =2.2+/-0.3 \times 10^{23} cm^{-2}, and photon index \Gamma = 0.4+/-0.1) with an average absorbed 2-8 keV flux of 1.4\times 10^{-11} erg cm^{-2} s^{-1}. At the Chandra-based position, a source is detected in our near infrared (NIR) maps with K_s = 13.14+/-0.04 mag, H_n = 16+/-0.1 mag, and no J_n band counterpart down to ~18 mag. The NIR source, compatible with 2MASS J18175218-1621316, shows no variability between 2012 March 13 and March 26. Searches of the UKIDSS database show similar NIR flux levels at epochs six months prior to and after a 2007 February 11 archival Chandra observation where the source's X-ray flux was at least 87 times fainter. In many ways IGR J18179-1621 is unusual: its combination of a several week long outburst (without evidence of repeated outbursts in the historical record), high absorption column (a large fraction of which is likely local to the system), and 11.82 s period does not fit neatly into existing X-ray binary categories.

The most peculiar radio characteristics of the TeV emitting high-mass X-ray binary LS I +61{\deg}303 are two periodicities: A large periodic outburst which exhibits the same period as the orbit (phase \Phi) and a second periodicity of 1667 days (phase \Theta) which modulates the orbital phase and amplitude of the large outburst. Recent analysis of the radio spectral index present strong evidence for the presence of the critical transition from optically thick emission (related to a steady jet) to an optically thin outburst (related to a transient jet) as in other microquasars. In parallel, a switch from a low/hard X-ray state to a transitional state (e.g. steep power law state) would be expected. We show how the critical transition from optically thick emission to an optically thin outburst is modulated by \Theta. Folding over a too large \Theta-interval mixes up important information about the outbursts and can yield a false picture of the emission behaviour of the source along the orbit. We therefore analyse the implications of this long period for treatment of hard X-ray/high energy data obtained from LS I +61{\deg}303, e.g. with Fermi-LAT or INTEGRAL, taking into account this long-term periodicity.

We present an analysis of the highly variable accreting X-ray pulsar 3A 1954+319 using 2005-2009 monitoring data obtained with INTEGRAL and Swift. This considerably extends the pulse period history and covers flaring episodes in 2005 and 2008. In 2006 the source was identified as one of only a few known symbiotic X-ray binaries (SyXBs), i.e., systems composed of a neutron star accreting from the inhomogeneous medium around an M-giant star. The extremely long pulse period of 5.3 hr is directly visible in the 2008 INTEGRAL-ISGRI outburst light curve. The pulse profile is double peaked and generally not significantly energy dependent although there is an indication of possible softening during the main pulse. During the outburst a strong spin-up of -1.8 10^(-4) hr hr^(-1) occurred. Between 2005 and 2008 a long-term spin-down trend of 2.1 10^-5 hr hr^(-1) was observed for the first time for this source. The 3-80 keV pulse peak spectrum of 3A 1954+319 during the 2008 flare could be well described by a thermal Comptonization model. We interpret the results within the framework of a recently developed quasi-spherical accretion model for SyXBs.

We present the analysis of two long, quasi-uninterrupted RXTE observations of Cygnus X-1 that span several days within a 10 d interval. The spectral characteristics during this observation cover the region where previous observations have shown the source to be most dynamic. Despite that the source behavior on time scales of hours and days is remarkably similar to that on year time scales. This includes a variety of spectral/temporal correlations that previously had only been observed over Cyg X-1's long-term evolution. Furthermore, we observe a full transition from a hard to a soft spectral state that occurs within less than 2.5 hours - shorter than previously reported for any other similar Cyg X-1 transition. We describe the spectra with a phenomenological model dominated by a broken power law, and we fit the X-ray variability power spectra with a combination of a cutoff power law and Lorentzian components. The spectral and timing properties are correlated: the power spectrum Lorentzian components have an energy-dependent amplitude, and their peak frequencies increase with photon spectral index. Averaged over 3.2-10 Hz, the time lag between the variability in the 4.5-5.7 keV and 9.5-15 keV bands increases with decreasing hardness when the variability is dominated by the Lorentzian components during the hard state. The lag is small when there is a large power law noise contribution, shortly after the transition to the soft state. Interestingly, the soft state not only shows the shortest lags, but also the longest lags when the spectrum is at its softest and faintest. We discuss our results in terms of emission models for black hole binaries.

We present a timing and spectral analysis of the variable X-ray source 3A 1954+319. Our analysis is mainly based on an outburst serendipitously observed during INTEGRAL Key Program observations of the Cygnus region in 2008 fall and on the Swift/BAT longterm light curve. Previous observations, though sparse, have identified the source to be one of only nine known symbiotic X-ray binaries, i.e., systems composed of an accreting neutron star orbiting in a highly inhomogeneous medium around an M-giant companion. The spectrum of 3A 1954+319 above 20 keV can be best described by a broken power law model. The extremely long pulse period of ~5.3 hours is clearly visible in the INTEGRAL/ISGRI light curve and confirmed through an epoch folding period search. Furthermore, the light curve allows us to determine a very strong spin up of -2x10^-4 h/h during the outburst. This spin up is confirmed by the pulse period evolution calculated from Swift/BAT data. The Swift/BAT data also show a long spin-down trend prior to the 2008 outburst, which is confirmed in archival INTEGRAL/ISGRI data. We discuss possible accretion models and geometries allowing for the transfer of such large amounts of angular momentum and investigate the harder spectrum of this outburst compared to previously published results.

We report on a Spitzer/IRS (mid-infrared), RXTE/PCA+HEXTE (X-ray), and Ryle (radio) simultaneous multiwavelength study of the microquasar Cygnus X-1, which aimed at an investigation of the origin of its mid-infrared emission. Compact jets were present in two out of three observations, and we show that they strongly contribute to the mid-infrared continuum. During the first observation, we detect the spectral break - where the transition from the optically thick to the optically thin regime takes place - at about 2.9e13 Hz. We then show that the jet's optically thin synchrotron emission accounts for the Cygnus X-1's emission beyond 400 keV, although it cannot alone explain its 3-200 keV continuum. A compact jet was also present during the second observation, but we do not detect the break, since it has likely shifted to higher frequencies. In contrast, the compact jet was absent during the last observation, and we show that the 5-30 micron mid-infrared continuum of Cygnus X-1 stems from the blue supergiant companion star HD 226868. Indeed, the emission can then be understood as the combination of the photospheric Raleigh-Jeans tail and the bremsstrahlung from the expanding stellar wind. Moreover, the stellar wind is found to be clumpy, with a filling factor f {\infty}~0.09-0.10. Its bremsstrahlung emission is likely anti-correlated to the soft X-ray emission, suggesting an anti-correlation between the mass-loss and mass-accretion rates. Nevertheless, we do not detect any mid-infared spectroscopic evidence of interaction between the jets and the Cygnus X-1's environment and/or companion star's stellar wind.

Because of their inherently high flux allowing the detection of clear signals, black hole X-ray binaries are interesting candidates for polarization studies, even if no polarization signals have been observed from them before. Such measurements would provide further detailed insight into these sources' emission mechanisms. We measured the polarization of the gamma-ray emission from the black hole binary system Cygnus X-1 with the INTEGRAL/IBIS telescope. Spectral modeling of the data reveals two emission mechanisms: The 250-400 keV data are consistent with emission dominated by Compton scattering on thermal electrons and are weakly polarized. The second spectral component seen in the 400keV-2MeV band is by contrast strongly polarized, revealing that the MeV emission is probably related to the jet first detected in the radio band.

We present broad band analysis of the Be/X-ray transient 1A 1118-61 by Suzaku at the peak of its 3rd observed outburst in January 2009 and 2 weeks later when the source flux had decayed by an order of magnitude. The continuum was modeled with a \texttt{cutoffpl} model as well as a compTT model, with both cases requiring an additional black body component at lower energies. We confirm the detection of a cyclotron line at ~5 keV and discuss the possibility of a first harmonic at ~110 keV. Pulse profile comparisons show a change in the profile structure at lower energies, an indication for possible changes in the accretion geometry. Phase resolved spectroscopy in the outburst data show a change in the continuum throughout the pulse period. The decrease in the CRSF centroid energy also indicates that the viewing angle on the accretion column is changing throughout the pulse period.

In 2005 February we observed Cygnus X-1 over a period of 10 days quasi-continuously with the Rossi X-ray Timing Explorer and the Ryle telescope. We present the results of the spectral and timing analysis on a timescale of 90 min and show that the behavior of Cyg X-1 is similar to that found during our years long monitoring campaign. As a highlight we present evidence for a full transition from the hard to the soft state that happened during less than three hours. The observation provided a more complete picture of a state transition than before, especially concerning the evolution of the time lags, due to unique transition coverage and analysis with high time resolution.

The flight calibration of the spectral response of CCD instruments below 1.5 keV is difficult in general because of the lack of strong lines in the on-board calibration sources typically available. We have been using 1E 0102.2-7219, the brightest supernova remnant in the Small Magellanic Cloud, to evaluate the response models of the ACIS CCDs on the Chandra X-ray Observatory (CXO), the EPIC CCDs on the XMM-Newton Observatory, the XIS CCDs on the Suzaku Observatory, and the XRT CCD on the Swift Observatory. E0102 has strong lines of O, Ne, and Mg below 1.5 keV and little or no Fe emission to complicate the spectrum. The spectrum of E0102 has been well characterized using high-resolution grating instruments, namely the XMM-Newton RGS and the CXO HETG, through which a consistent spectral model has been developed that can then be used to fit the lower-resolution CCD spectra. We have also used the measured intensities of the lines to investigate the consistency of the effective area models for the various instruments around the bright O (~570 eV and 654 eV) and Ne (~910 eV and 1022 eV) lines. We find that the measured fluxes of the O VII triplet, the O VIII Ly-alpha line, the Ne IX triplet, and the Ne X Ly-alpha line generally agree to within +/-10 % for all instruments, with 28 of our 32 fitted normalizations within +/-10% of the RGS-determined value. The maximum discrepancies, computed as the percentage difference between the lowest and highest normalization for any instrument pair, are 23% for the O VII triplet, 24% for the O VIII Ly-alpha line, 13% for the Ne IX triplet, and 19% for the Ne X Ly-alpha line. If only the CXO and XMM are compared, the maximum discrepancies are 22% for the O VII triplet, 16% for the O VIII Ly-alpha line, 4% for the Ne IX triplet, and 12% for the Ne X Ly-alpha line.

We present a detailed analysis of observations of the high mass X-ray binary Cen X-3 spanning two consecutive binary orbits performed with the RXTE satellite in early March 1997. The PCA and HEXTE light curves both show a clear reduction in count rate after mid-orbit for both binary revolutions. We therefore analyze two broad band spectra for each orbit, before and after mid-orbit. Consistent with earlier observations these four joint PCA and HEXTE spectra can be well described using a phenomenological pulsar continuum model, including an iron emission line and a cyclotron resonance scattering feature. While no strong spectral variations were detected, the second half of orbit 2 shows a tendency toward being softer and more strongly absorbed. In order to follow the orbital phase-dependent evolution of the spectrum in greater detail, we model spectra for shorter exposures, confirming that most spectral parameters show either a gradual or sudden change for the second half of the second orbit. A comparison with a simple wind model indicates the existence of an accretion wake in this system. We also present and discuss high resolution pulse profiles for several different energy bands, as well as their hardness ratios. PCA and HEXTE spectra were created for 24 phase bins and fitted using the same model as in the phase averaged case. Systematic pulse phase-dependent variations of several continuum and cyclotron line parameters were detected, most notably a significant increase of the cyclotron line energy during the early rise of the main peak, followed by a gradual decrease. We show that applying a simple dipole model for the magnetic field is not sufficient to describe our data.

We present an evaluation of a simulated cosmic ray shower, based on {\sc geant4} and {\sc top-c}, which tracks all the particles in the shower. {\sc top-c} (Task Oriented Parallel C) provides a framework for parallel algorithm development which makes tractable the problem of following each particle. This method is compared with a simulation program which employs the Hillas thinning algorithm.

{\sc top-c} (Task Oriented Parallel C) is a freely available package for parallel computing. It is designed to be easy to learn and to have good tolerance for the high latencies that are common in commodity networks of computers. It has been successfully used in a wide range of examples, providing linear speedup with the number of computers. A brief overview of {\sc top-c} is provided, along with recent experience with cosmic ray physics simulations.

We describe a scalable parallelization of Geant4 using commodity hardware in a collaborative effort between the College of Computer Science and the Department of Physics at Northeastern University. The system consists of a Beowulf cluster of 32 Pentium II processors with 128 MBytes of memory each, connected via ATM and fast Ethernet. The bulk of the parallelization is done using TOP-C (Task Oriented Parallel C), software widely used in the computational algebra community. TOP-C provides a flexible and powerful framework for parallel algorithm development, is easy to learn, and is available at no cost. Its task oriented nature allows one to parallelize legacy code while hiding the details of interprocess communications. Applications include fast interactive simulation of computationally intensive processes such as electromagnetic showers. General results motivate wider applications of TOP-C to other simulation problems as well as to pattern recognition in high energy physics.