Black-hole binary (BHB) systems comprise a stellar-mass black hole and a closely orbiting companion star. Matter is transferred from the companion to the black hole, forming an accretion disk, corona and jet structures. The resulting release of gravitational energy leads to emission of X-rays. The radiation is affected by special/general relativistic effects, and can serve as a probe of the properties of the black hole and surrounding environment, if the accretion geometry is properly identified. Two competing models describe the disk-corona geometry for the hard spectral state of BHBs, based on spectral and timing measurements. Measuring the polarization of hard X-rays reflected from the disk allows the geometry to be determined. The extent of the corona differs between the two models, affecting the strength of relativistic effects (e.g., enhancement of polarization fraction and rotation of polarization angle). Here, we report observational results on linear polarization of hard X-ray (19-181 keV) emission from a BHB, Cygnus X-1, in the hard state. The low polarization fraction, <8.6% (upper limit at 90% confidence level), and the alignment of the polarization angle with the jet axis show that the dominant emission is not influenced by strong gravity. When considered together with existing spectral and timing data, our result reveals that the accretion corona is either an extended structure, or is located far from the black hole in the hard state of Cygnus X-1.

We report a polarimetric constraint on the hard X-ray synchrotron jet emission from the Cygnus X-1 black-hole binary system. The observational data were obtained using the PoGO+ hard X-ray polarimeter in July 2016, when Cygnus X-1 was in the hard state. We have previously reported that emission from an extended corona with a low polarization fraction is dominating, and that the polarization angle is perpendicular to the disk surface. In the soft gamma-ray regime, a highly-polarized synchrotron jet is reported with INTEGRAL observations. To constrain the polarization fraction and flux of such a jet component in the hard X-ray regime, we now extend analyses through vector calculations in the Stokes QU plane, where the dominant corona emission and the jet component are considered simultaneously. The presence of another emission component with different polarization angle could partly cancel out the net polarization. The 90% upper limit of the polarization fraction for the additional synchrotron jet component is estimated as <10%, <5%, and <5% for polarization angle perpendicular to the disk surface, parallel to the surface, and aligned with the emission reported by INTEGRAL data, respectively. From the 20-180 keV total flux of 2.6 x 10^-8 erg s^-1 cm^-2, the upper limit of the polarized flux is estimated as <3 x 10^-9 erg s^-1 cm^-2.

The linear polarisation fraction and angle of the hard X-ray emission from the Crab provide unique insight into high energy radiation mechanisms, complementing the usual imaging, timing and spectroscopic approaches. Results have recently been presented by two missions operating in partially overlapping energy bands, PoGO+ (18-160 keV) and AstroSat CZTI (100-380 keV). We previously reported PoGO+ results on the polarisation parameters integrated across the light-curve and for the entire nebula-dominated off-pulse region. We now introduce finer phase binning, in light of the AstroSat CZTI claim that the polarisation fraction varies across the off-pulse region. Since both missions are operating in a regime where errors on the reconstructed polarisation parameters are non-Gaussian, we adopt a Bayesian approach to compare results from each mission. We find no statistically significant variation in off-pulse polarisation parameters, neither when considering the mission data separately nor when they are combined. This supports expectations from standard high-energy emission models.

Strong magnetic fields, synchrotron emission, and Compton scattering are omnipresent in compact celestial X-ray sources. Emissions in the X-ray energy band are consequently expected to be linearly polarized. X-ray polarimetry provides a unique diagnostic to study the location and fundamental mechanisms behind emission processes. The polarization of emissions from a bright celestial X-ray source, the Crab, is reported here for the first time in the hard X-ray band (~20-160 keV). The Crab is a complex system consisting of a central pulsar, a diffuse pulsar wind nebula, as well as structures in the inner nebula including a jet and torus. Measurements are made by a purpose-built and calibrated polarimeter, PoGO+. The polarization vector is found to be aligned with the spin axis of the pulsar for a polarization fraction, PF = (20.9 $\pm$ 5.0)%. This is higher than that of the optical diffuse nebula, implying a more compact emission site, though not as compact as, e.g., the synchrotron knot. Contrary to measurements at higher energies, no significant temporal evolution of phase-integrated polarisation parameters is observed. The polarization parameters for the pulsar itself are measured for the first time in the X-ray energy band and are consistent with observations at optical wavelengths.

Polarimetric observations of celestial sources in the hard X-ray band stand to provide new information on emission mechanisms and source geometries. PoGO+ is a Compton scattering polarimeter (20-150 keV) optimised for the observation of the Crab (pulsar and wind nebula) and Cygnus X-1 (black hole binary), from a stratospheric balloon-borne platform launched from the Esrange Space Centre in summer 2016. Prior to flight, the response of the polarimeter has been studied with polarised and unpolarised X-rays allowing a Geant4-based simulation model to be validated. The expected modulation factor for Crab observations is found to be $M_{\mathrm{Crab}}=(41.75\pm0.85)\%$, resulting in an expected Minimum Detectable Polarisation (MDP) of $7.3\%$ for a 7 day flight. This will allow a measurement of the Crab polarisation parameters with at least $5\sigma$ statistical significance assuming a polarisation fraction $\sim20\%$ $-$ a significant improvement over the PoGOLite Pathfinder mission which flew in 2013 and from which the PoGO+ design is developed.

PoGOLite is a balloon-borne hard X-ray polarimeter dedicated to the study of point sources. Compton scattered events are registered using an array of plastic scintillator units to determine the polarisation of incident X-rays in the energy range 20 - 240 keV. In 2013, a near circumpolar balloon flight of 14 days duration was completed after launch from Esrange, Sweden, resulting in a measurement of the linear polarisation of the Crab emission. Building on the experience gained from this Pathfinder flight, the polarimeter is being modified to improve performance for a second flight in 2016. Such optimisations, based on Geant4 Monte Carlo simulations, take into account the source characteristics, the instrument response and the background environment which is dominated by atmospheric neutrons. This paper describes the optimisation of the polarimeter and details the associated increase in performance. The resulting design, PoGO+, is expected to improve the Minimum Detectable Polarisation (MDP) for the Crab from 19.8% to 11.1% for a 5 day flight. Assuming the same Crab polarisation fraction as measured during the 2013 flight, this improvement in MDP will allow a 5{\sigma} constrained result. It will also allow the study of the nebula emission only (Crab off-pulse) and Cygnus X-1 if in the hard state.

We have measured the linear polarisation of hard X-ray emission from the Crab in a previously unexplored energy interval, 20-120 keV. The introduction of two new observational parameters, the polarisation fraction and angle stands to disentangle geometrical and physical effects, thereby providing information on the pulsar wind geometry and magnetic field environment. Measurements are conducted using the PoGOLite Pathfinder - a balloon-borne polarimeter. Polarisation is determined by measuring the azimuthal Compton scattering angle of incident X-rays in an array of plastic scintillators housed in an anticoincidence well. The polarimetric response has been characterised prior to flight using both polarised and unpolarised calibration sources. We address possible systematic effects through observations of a background field. The measured polarisation fraction for the integrated Crab light-curve is ($18.4^{+9.8}_{-10.6}$)%, corresponding to an upper limit (99% credibility) of 42.4%, for a polarisation angle of ($149.2\pm16.0)^\circ$.

In the 50 years since the advent of X-ray astronomy there have been many scientific advances due to the development of new experimental techniques for detecting and characterising X-rays. Observations of X-ray polarisation have, however, not undergone a similar development. This is a shortcoming since a plethora of open questions related to the nature of X-ray sources could be resolved through measurements of the linear polarisation of emitted X-rays. The PoGOLite Pathfinder is a balloon-borne hard X-ray polarimeter operating in the 25 - 240 keV energy band from a stabilised observation platform. Polarisation is determined using coincident energy deposits in a segmented array of plastic scintillators surrounded by a BGO anticoincidence system and a polyethylene neutron shield. The PoGOLite Pathfinder was launched from the SSC Esrange Space Centre in July 2013. A near-circumpolar flight was achieved with a duration of approximately two weeks. The flight performance of the Pathfinder design is discussed for the three Crab observations conducted. The signal-to-background ratio for the observations is shown to be 0.25$\pm$0.03 and the Minimum Detectable Polarisation (99% C.L.) is (28.4$\pm$2.2)%. A strategy for the continuation of the PoGOLite programme is outlined based on experience gained during the 2013 maiden flight.

Polarimetric studies of astrophysical sources can make important contributions to resolve the geometry of the emitting region and determine the photon emission mechanism. PoGOLite is a balloon-borne polarimeter operating in the hard X-ray band (25-240 keV), with a Pathfinder mission focussing on Crab observations. Within the polarimeter, the distribution of Compton scattering angles is used to determine the polarisation fraction and angle of incident photons. To assure an unbiased measurement of the polarisation during a balloon flight it is crucial to characterise the performance of the instrument before the launch. This paper presents the results of the PoGOLite calibration tests and simulations performed before the 2013 balloon flight. The tests performed confirm that the polarimeter does not have any intrinsic asymmetries and therefore does not induce bias into the measurements. Generally, good agreement is found between results from test data and simulations which allows the polarimeter performance to be estimated for Crab observations.

PoGOLino is a balloon-borne scintillator-based experiment developed to study the largely unexplored high altitude neutron environment at high geomagnetic latitudes. The instrument comprises two detectors that make use of LiCAF, a novel neutron sensitive scintillator, sandwiched by BGO crystals for background reduction. The experiment was launched on March 20th 2013 from the Esrange Space Centre, Northern Sweden (geomagnetic latitude of $65^\circ$), for a three hour flight during which the instrument took data up to an altitude of 30.9 km. The detector design and ground calibration results are presented together with the measurement results from the balloon flight.

Nowadays, analyzing and reducing the ever larger astronomical datasets is becoming a crucial challenge, especially for long cumulated observation times. The INTEGRAL/SPI X-gamma-ray spectrometer is an instrument for which it is essential to process many exposures at the same time in order to increase the low signal-to-noise ratio of the weakest sources. In this context, the conventional methods for data reduction are inefficient and sometimes not feasible at all. Processing several years of data simultaneously requires computing not only the solution of a large system of equations, but also the associated uncertainties. We aim at reducing the computation time and the memory usage. Since the SPI transfer function is sparse, we have used some popular methods for the solution of large sparse linear systems; we briefly review these methods. We use the Multifrontal Massively Parallel Solver (MUMPS) to compute the solution of the system of equations. We also need to compute the variance of the solution, which amounts to computing selected entries of the inverse of the sparse matrix corresponding to our linear system. This can be achieved through one of the latest features of the MUMPS software that has been partly motivated by this work. In this paper we provide a brief presentation of this feature and evaluate its effectiveness on astrophysical problems requiring the processing of large datasets simultaneously, such as the study of the entire emission of the Galaxy. We used these algorithms to solve the large sparse systems arising from SPI data processing and to obtain both their solutions and the associated variances. In conclusion, thanks to these newly developed tools, processing large datasets arising from SPI is now feasible with both a reasonable execution time and a low memory usage.

The INTEGRAL/SPI, X-gamma-ray spectrometer (20 keV - 8 MeV) is an instrument for which recovering source intensity variations is not straightforward and can constitute a difficulty for data analysis. In most cases, determining the source intensity changes between exposures is largely based on a priori information. We propose techniques that help to overcome the difficulty related to source intensity variations, which make this step more rational. In addition, the constructed "synthetic" light curves should permit us to obtain a sky model that describes the data better and optimizes the source signal-to-noise ratios. For this purpose, the time intensity variation of each source was modeled as a combination of piecewise segments of time during which a given source exhibits a constant intensity. To optimize the signal-to-noise ratios, the number of segments was minimized. We present a first method that takes advantage of previous time series that can be obtained from another instrument on-board the INTEGRAL observatory. A data segmentation algorithm was then used to synthesize the time series into segments. The second method no longer needs external light curves, but solely SPI raw data. For this, we developed a specific algorithm that involves the SPI transfer function. The time segmentation algorithms that were developed solve a difficulty inherent to the SPI instrument, which is the intensity variations of sources between exposures, and it allows us to obtain more information about the sources' behavior.

We present recent improvements in polarization analysis with the INTEGRAL SPI data. The SPI detector plane consists of 19 independent Ge crystals and can operate as a polarimeter. The anisotropy characteristics of Compton diffusions can provide information on the polarization parameters of the incident flux. By including the physics of the polarized Compton process in the instrument simulation, we are able to determine the instrument response for a linearly polarized emission at any position angle. We compare the observed data with the simulation sets by a minimum \chi^2 technique to determine the polarization parameters of the source (angle and fraction). We have tested our analysis procedure with Crab nebula observations and find a position angle similar to those previously reported in the literature, with a comfortable significance. Since the instrument response depends on the incident angle, each exposure in the SPI data requires its own set of simulations, calculated for 18 polarization angles (from 0{\deg} to 170{\deg} in steps of 10{\deg}) and unpolarized emission. The analysis of a large amount of observations for a given source, required to obtain statistically significant results, represents a large amount of computing time, but it is the only way to access this complementary information in the hard X-ray regime. Indeed, major scientific advances are expected from such studies since the observational results will help to discriminate between the different models proposed for the high energy emission of compact objects like X-ray binaries and active galactic nuclei or gamma-ray bursts.

Operational since 2002 on-board the INTEGRAL observatory, the SPI spectrometer can be used to perform polarization measurements in the hard X-ray/soft gamma-ray domain (~ 130 keV - 8 MeV). However, this phenomenon is complex to measure at high energy and requires high fluxes. Cyg X-1 appears as the best candidate amongst the X-ray binaries since it is one of the brightest persistent sources in this energy domain. Furthermore, a polarized component has recently been reported above 400 keV from IBIS data. We have therefore dedicated our efforts to develop the required tools to study the polarization in the INTEGRAL SPI data and have first applied them to 2.6 Ms of Cyg X-1 observations, covering 6.5 years of the INTEGRAL mission. We have found that the high energy emission of Cyg X-1 is indeed polarized, with a mean polarization fraction of 76 % +/- 15 % at a position angle estimated to 42 +/- 3 degrees, for energies above 230 keV. The polarization fraction clearly increases with energy. In the 130-230 keV band, the polarization fraction is lower than 20 %, but exceeds 75 % between 370 and 850 keV, with the (total) emission vanishing above this energy. This result strongly suggests that the emission originates from the jet structure known to emit in the radio domain. The same synchrotron process could be responsible for the emission from radio to MeV, implying the presence of high energy electrons. This illustrates why the polarization of the high energy emission in compact objects is an increasingly important observational objective.