We report results from a multi-mission observational campaign of the transient X-ray pulsar 2S 1417-624 during its 2025 outburst, using data from NICER, IXPE, and NuSTAR. Phase-averaged and phase-resolved spectroscopy with NICER and NuSTAR reveal that a typical broken power-law model with a high-energy cut-off well describes the broadband spectra. Several spectral parameters, however, show clear and systematic modulations with pulse phase, indicating variations in the physical conditions of the emitting plasma over the neutron star's rotation. IXPE provides the first polarimetric measurements of this source, yielding a phase-averaged polarization degree (PD) of 4.8 ± 1.2
We present a first look at the Galactic population of heretofore known HMXBs as observed by SRG/eROSITA during its first four surveys. eROSITA's sensitivity of ∼10^-13 erg s^-1 cm^-2, translating to 10^32-10^34,erg s^-1 in luminosity for most known HMXBs in the Milky Way, has thus far never been reached by any wide-area survey instrument. We present the extended log N-log L distribution of known HMXBs reaching down to 10^32 erg s^-1 using eROSITA, and show the large scatter that can be induced by source intrinsic variability. We present sub-type resolved luminosity distributions, showing that the Supergiant X-ray binaries (SgXBs) and Be X-ray binaries (BeXRBs) occupy different parts of the overall distribution, and reanalyse RXTE/ASM data and MAXI for comparison to eROSITA. The luminosity regime uncovered by eROSITA allows a systematic study of the "transient" BeXRBs, which are typically below the detection threshold of monitors outside of outburst, and whose low luminosity behavior has been a longstanding question. Signatures of stable accretion at low luminosities have been observed with pointed instruments for a fraction of the overall sample, so far. With the eROSITA results, we posit that accretion outside of outburst is likely the norm, since a vast majority (> 80
Context. Cyg X-1 is one of the brightest X-ray binaries and has been observed multiple times with the Imaging X-ray Polarimetry Explorer (IXPE). Recent studies report tentative evidence for a polarisation modulation with the orbital period P, but a half-period (P/2) signal, expected from reflection off the companion star and its stellar wind, has not been reported. Aims. We aim to quantify the reflection-induced variations of the polarisation degree PD and polarisation angle PA in Cyg X-1 as a function of orbital phase and energy, and interpret these in terms of binary geometry and wind structure. Methods. We set up a radiative transfer model combining a general relativistic description of the polarised source emission (kerrC) with a focussed stellar wind model for the binary medium. Using the 3D X-ray radiative transfer code SKIRT, we simulate broadband Stokes I, Q, and U fluxes, surface brightness maps, and linear polarisation maps over one binary orbit. Results. We find a prominent double-peaked (P/2) polarisation modulation, with a peak-to-peak PD amplitude of 0.25, 0.81, and 1.24 percentage points in the 2-4, 4-6, and 6-8 keV bands, respectively, with a strong energy dependence. The PA modulation is more modest, with |ΔPA| < 4.6°. Crucially, X-ray reprocessing reduces the overall PD relative to the source polarisation. Conclusions. The modulation is driven by reflection off the companion star and the focussed wind, which induces a polarisation signal that alternately reinforces and counteracts the source polarisation throughout the orbit. The diffuse scattering halo surrounding the source systematically reduces the PD, an effect that should be accounted for in all wind-fed XRBs. The PD amplitude increases with energy as absorption disproportionately attenuates the distant-reflection signal; as the extinction drops, the reflection signal becomes increasingly important.
Recent X-ray and soft gamma-ray spectro-polarimetric observations of black hole X-ray binaries, such as Cygnus X-1, reveal discrepancies that challenge current models of accretion and high-energy emission. We present a general framework for broadband spectro-polarimetric modeling that simultaneously fits spectral and polarization data from independent instruments across different energy bands. The method robustly combines measurements with different energy binning and statistical significance while accounting for energy-dependent polarization mixing between model components. Applied to Cygnus X-1 over the 2 keV-2 MeV range, it provides an improved characterization of the soft gamma-ray hard-tail and constrains the optically thin synchrotron cutoff energy to (3.9^+0.6_-0.5) × 10^2 keV. We discuss the implications for non-thermal electron acceleration in the context of Bohm diffusion and synchrotron cooling, and propose a new explanation for the observed polarization misalignment based on magnetic field helicity and Doppler boosting. The framework is broadly applicable to multi-instrument spectro-polarimetric analyses.
We present a comprehensive analysis of the X-ray spectrum of the highly obscured high-mass X-ray binary IGR J16318-4848 based on similar to 250 ks Chandra/HETGS observations covering different orbital phases. We detected a strong narrow Fe K alpha emission line and interpret our finding as the iron line originating in the stellar wind close to the compact object. We tentatively identify a Compton shoulder associated with the Fe K alpha line. The ratio of the Compton-shoulder flux to that of the Fe K alpha core exhibits orbital-phase variability, reaching a maximum of similar to 10% around orbital phase phi orb similar or equal to 0.33. In addition, we do not detect ionized lines, implying that the plasma in the system is predominantly cold. We preliminarily report that the Fe K alpha line centroid exhibits an orbital modulation with an amplitude of 167 +/- 32 km s-1, adopting the best-known circular orbit in the source. This study presents the first detailed view of the high-resolution X-ray spectra of IGR J16318-4848 across different orbital phases, providing timely context for the newly taken XRISM observations of the source.
Cygnus X-1 is a persistent, high-mass black hole X-ray binary (BHXRB) which in the hard state shows many similar properties to transient BHXRBs, along with intriguing differences, such as the lack of quasi-periodic oscillations. Here, we compare for the first time the detailed spectral-timing properties of Cyg X-1 with a transient BHXRB, MAXI J1820+070, combining data from XMM-Newton and NICER (Neutron Star Interior Composition Explorer) with contemporaneous INTEGRAL (International Gamma-Ray Astrophysics Laboratory) data to study the power spectra, rms spectra, and time lags over a broad 0.5-200 keV range. We select bright hard state MAXI J1820+070 data with similar power-spectral shapes to the Cyg X-1 data, to compare the source behaviours while accounting for the evolution of spectral-timing properties, notably the lags, through the hard state. Cyg X-1 shows no evidence for soft lags in the 1-10 Hz frequency range where they are clearly detected for MAXI J1820+070. Furthermore, the low-frequency hard lags and rms spectra evolve much more strongly during the hard state of Cyg X-1 than for MAXI J1820+070. We argue that these differences cannot be explained by the different black hole masses of these systems, but may be related to their different accretion rates and corresponding locations on the hardness-intensity diagram. We conjecture that there is a significant luminosity-dependence of coronal geometry in the hard state of BHXRBs, rather than an intrinsic difference between Cyg X-1 and transient BHXRBs. This possibility has also been suggested to explain a common time-lag feature that appears in the hard intermediate states of Cyg X-1 and transient BHXRBs.
We present the first X-ray spectropolarimetric results for Cygnus X-1 in its soft state from a campaign of five IXPE observations conducted during 2023 May-June. Companion multiwavelength data during the campaign are likewise shown. The 2-8 keV X-rays exhibit a net polarization degree PD=1.99%+/-0.13% (68% confidence). The polarization signal is found to increase with energy across IXPE's 2-8 keV bandpass. The polarized X-rays exhibit an energy-independent polarization angle of PA=-25.7+/-1.8 deg. East of North (68% confidence). This is consistent with being aligned to Cyg X-1's AU-scale compact radio jet and its pc-scale radio lobes. In comparison to earlier hard-state observations, the soft state exhibits a factor of 2 lower polarization degree, but a similar trend with energy and a similar (also energy-independent) position angle. When scaling by the natural unit of the disk temperature, we find the appearance of a consistent trendline in the polarization degree between soft and hard states. Our favored polarimetric model indicates Cyg X-1's spin is likely high (a* above ~0.96). The substantial X-ray polarization in Cyg X-1's soft state is most readily explained as resulting from a large portion of X-rays emitted from the disk returning and reflecting off the disk surface, generating a high polarization degree and a polarization direction parallel to the black hole spin axis and radio jet. In IXPE's bandpass, the polarization signal is dominated by the returning reflection emission. This constitutes polarimetric evidence for strong gravitational lensing of X-rays close to the black hole.
Aims. Soft γ-ray emission (100 keV–10 MeV) has previously been detected in the hard state of several microquasars. In some sources, this emission was found to be highly polarized and was suggested to be emitted at the base of the jet. Until now, no γ-ray polarization had been found in any other state. Methods. Using INTEGRAL/IBIS, we studied the soft γ-ray spectral and polarization properties of Swift J1727.8−1613 throughout its outburst. Results. We detect a highly polarized spectral component in both the hard intermediate state and the early stages of the soft intermediate state above 210 keV. In the hard intermediate state, the polarization angle significantly deviates from the compact jet angle projected onto the sky, whereas in the soft intermediate they are closely aligned. This constitutes the first detection of jet-aligned polarization in the soft γ-ray for a microquasar. We attribute this polarized spectral component to synchrotron emission from the jet, which indicates that some of the jet might persist into the softer states.
The Neutron Star Interior Composition Explorer (NICER) monitoring campaign of Cyg X-1 allows us to study its spectral-timing behavior at energies <1 keV across all states. The hard state power spectrum can be decomposed into two main broad Lorentzians with a transition at around 1 Hz. The lower-frequency Lorentzian is the dominant component at low energies. The higher-frequency Lorentzian begins to contribute significantly to the variability above 1.5 keV and dominates at high energies. We show that the low- and high-frequency Lorentzians likely represent individual physical processes. The lower-frequency Lorentzian can be associated with a (possibly Comptonized) disk component, while the higher-frequency Lorentzian is clearly associated with the Comptonizing plasma. At the transition of these components, we discover a low-energy timing phenomenon characterized by an abrupt lag change of hard (≳2 keV) with respect to soft (≲1.5 keV) photons, accompanied by a drop in coherence, and a reduction in amplitude of the second broad Lorentzian. The frequency of the phenomenon increases with the frequencies of the Lorentzians as the source softens and cannot be seen when the power spectrum is single-humped. A comparison to transient low-mass X-ray binaries shows that this feature does not only appear in Cyg X-1, but that it is a general property of accreting black hole binaries. In Cyg X-1, we find that the variability at low and high energies is overall highly coherent in the hard and intermediate states. The high coherence shows that there is a process at work which links the variability, suggesting a physical connection between the accretion disk and Comptonizing plasma. This process fundamentally changes in the soft state, where strong red noise at high energies is incoherent to the variability at low energies.
Strong stellar winds are an important feature in wind-accreting high-mass X-ray binary (HMXB) systems, providing insights into stellar evolution and their impact on surrounding environments. However, the long-term evolution and temporal variability of these winds are not fully understood. This work probes the archetypal wind-accreting HMXB Vela X-1 using MAXI observations over 14 years, focusing on orbit-to-orbit absorption variability in the 2-10 keV band. Additionally, the relation between hardness ratio trends in binary orbits and neutron star spin states is investigated. We calculate hardness ratios to track absorption variability, comparing flux changes across energy bands, as the effect of absorption on the flux is energy-dependent. Variability is analyzed by comparing hardness ratio trends across binary orbits to the MAXI long-term averaged evolution. The long-term averaged hardness ratio evolution displays a stable pattern. Yet, individual binary orbits reveal different hardness ratio evolutions between consecutive orbits with no evident periodicity. Less than half of the binary orbits align with the long-term evolution. Moreover, neutron star spin-up episodes exhibit harder-than-average hardness trends compared to spin-down episodes, although their distributions overlap considerably. The long-term averaged hardness ratio dispersion is consistent with absorption column densities reported in literature from shorter observations, suggesting that heterogeneous wind structures, including accretion wakes and wind clumps, drive observed variations. The orbit-to-orbit variability indicates that pointed X-ray observations provide limited insight into wind structure. The link between neutron star spin states and hardness trends underscores the influence of accretion on absorption, with variability tied to stellar wind density fluctuations.
Travel to academic conferences -- where international flights are the norm -- is responsible for a sizeable fraction of the greenhouse gas (GHG) emissions associated with academic work. In order to provide a benchmark for comparison with other fields, as well as for future reduction strategies and assessments, we estimate the CO2-equivalent emissions for conference travel in the field of astronomy for the prepandemic year 2019. The GHG emission of the international astronomical community's 362 conferences and schools in 2019 amounted to 42,500 tCO2e, assuming a radiative-forcing index factor of 1.95 for air travel. This equates to an average of 1.0 $\pm$ 0.6 tCO2e per participant per meeting. The total travel distance adds up to roughly 1.5 Astronomical Units, that is, 1.5 times the distance between the Earth and the Sun. We present scenarios for the reduction of this value, for instance with virtual conferencing or hub models, while still prioritizing the benefits conferences bring to the scientific community.
Context. Cygnus X-1 (Cyg X-1) is a high-mass X-ray binary where accretion onto the black hole (BH) is mediated by the stellar wind from the blue supergiant companion star HDE 226868. Due to its inclination, the system is a perfect laboratory to study the not yet well-understood stellar wind structure. In fact, depending on the position of the BH along the orbit, X-ray observations can probe different layers of the stellar wind. Deeper wind layers can be investigated at superior conjunction (i.e. null orbital phases). Aims. We aim to characterise the stellar wind in the Cyg X-1/HDE 226868 system, analysing one passage at superior conjunction covered by XMM-Newton during the 'Cyg X-1 Hard state Observations of a Complete Binary Orbit in X-rays' (CHOCBOX) campaign. Methods. To analyse the properties of the stellar wind, we computed colour-colour diagrams. Since X-ray absorption is energy-dependent, colour indices provide information on the parameters of the stellar wind, such as the column density, N-H,N- w, and the covering factor, f(c). We fitted colour-colour diagrams with models that include both a continuum and a stellar wind component. We used the kernel density estimation method to infer the unknown probability distribution of the data points in the colour-colour diagram, and selected the model corresponding to the highest likelihood. In order to study the temporal evolution of the wind around superior conjunction, we extracted and fitted time-resolved colour-colour diagrams. Results. We found that the model that best describes the shape of the colour-colour diagram of Cyg X-1 at superior conjunction requires the wind to be partially ionised. The shape of the colour-colour diagram strongly varies during the analysed observation, due to concurrent changes of the mean N-H,N- w and the f(c) of the wind. Our results suggest the existence of a linear scaling between the rapid variability amplitude of N-H,N- w (on timescales between 10 s and 11 ks) and its long-term variations (on timescales > 11 ks). Using the inferred best-fit values, we estimated the stellar mass loss rate to be similar to 7 x 10(-6) M-circle dot yr(-1) and the clumps to have a characteristic mass of similar to 10(17) g.
Stellar winds of massive stars are known to be driven by line absorption of UV photons, a mechanism which is prone to instabilities, causing the wind to be clumpy. The clumpy structure hampers wind mass-loss estimates, limiting our understanding of massive star evolution. The wind structure also impacts accretion in high-mass X-ray binary (HMXB) systems. We analyse the wavelength-dependent variability of X-ray absorption in the wind to study its structure. Such an approach is possible in HMXBs, where the compact object serves as an X-ray backlight. We probe different parts of the wind by analysing data taken at superior and inferior conjunction. We apply excess variance spectroscopy to study the wavelength-dependent soft X-ray variability of the HMXB Cygnus X-1 in the low/hard spectral state. Excess variance spectroscopy quantifies the variability of an object above the statistical noise as a function of wavelength, which allows us to study the variability of individual spectral lines. As one of the first studies, we apply this technique to high-resolution gratings spectra provided by Chandra, accounting for various systematic effects. The frequency dependence is investigated by changing the time binning. The strong orbital phase dependence we observe in the excess variance is consistent with column density variations predicted by a simple model for a clumpy wind. We identify spikes of increased variability with spectral features found by previous spectroscopic analyses of the same data set, most notably from silicon in over-dense clumps in the wind. In the silicon line region, the variability power is redistributed towards lower frequencies, hinting at increased line variability in large clumps. In prospect of the microcalorimetry missions that are scheduled to launch within the next decade, excess variance spectra present a promising approach to constrain the wind structure.
In black hole X-ray binaries, a nonthermal high-energy component is sometimes detected at energies above 200 keV. The origin of this high-energy component is debated and distinct spectral modelizations can lead to different interpretations. High-energy polarimetry measurements with the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) enable new diagnostics on the physics responsible for the MeV spectral component in black hole X-ray binaries. In this work, we aim to investigate the high-energy behavior of three bright sources discovered by the Monitor of All-sky X-ray Image: MAXI J1535−571, MAXI J1820+070, and MAXI J1348−630. We took advantage of their brightness to investigate their soft γ-ray (0.1–2 MeV) properties with INTEGRAL. We used both spectral and polarimetric approaches to probe their high-energy emission with the aim of bringing new constraints on the MeV emission in black hole X-ray binaries. We first studied the spectral characteristics of the sources in the 3–2000 keV range using JEM-X, IBIS, and SPI, with a semi-phenomenological description of the data. We then used IBIS as a Compton telescope in order to evaluate the polarization properties of the sources above 300 keV. A high-energy component was detected during the hard-intermediate state and soft-intermediate state of MAXI J1535−571, the low-hard state of MAXI J1820+070, and the low-hard state of MAXI J1348−630. The components detected in MAXI J1820+070 and MAXI J1348 − 630 were polarized with a polarization fraction of 26 ± 9° and > 56% in the 300–1000 keV range, respectively. With no polarization information for MAXI J1535−571, the component detected could either come from the jets or the corona. In the case of MAXI J1820+070, the extrapolation of the synchrotron spectrum measured in the infrared indicates that the component is likely due to a nonthermal distribution of electrons from a hybrid corona. For MAXI J1348−630, the high fraction of polarization points toward a jets origin; however, we cannot formally conclude this without any infrared data giving information on the optically thin part of the synchrotron spectrum.
Victoria Grinberg, who is the first “liaison scientist” at the European Space Agency in the Netherlands, talks to Hamish Johnston about her interest in binary stars, why she wants to make astronomy more sustainable and her love of scientific illustrations.
ABSTRACT Neutron stars accreting from OB supergiants are often divided between persistently and transiently accreting systems, called supergiant X-ray binaries (SgXBs) and supergiant fast X-ray transients (SFXTs). This dichotomy in accretion behaviour is typically attributed to systematic differences in the massive stellar wind, binary orbit, or magnetic field configuration, but direct observational evidence for these hypotheses remains sparse. To investigate their stellar winds, we present the results of pilot 100-GHz observations of one SFXT and one SgXB with the Northern Extended Millimetre Array. The SFXT, IGR J18410-0535, is detected as a point source at 63.4 ± 9.6 μJy, while the SgXB, IGR J18410-0535 remains undetected. Radio observations of IGR J18410-0535 imply a flat or inverted low-frequency spectrum, arguing for wind emission and against non-thermal flaring. Due to the uncertain SFXT distance, however, the observations do not necessarily imply a difference between the wind properties of the SFXT and SgXB. We compare the mm constraints with other HMXBs and isolated OB supergiants, before considering how future mm campaigns can constrain HMXB wind properties by including X-ray measurements. Specifically, we discuss caveats and future steps to successfully measure wind mass-loss rates and velocities in HMXBs with coordinated mm, radio, and X-ray campaigns.
Black hole X-ray binaries are ideal environments to test the accretion phenomena in the presence of strong gravitational potentials. KERRBB held an important place in the X-ray spectral continuum method for measuring the black hole spin modeling the emission from the innermost regions of the accretion disk. In this work, we present the results of X-ray spectral analysis using publicly available RXTE data of GRO J1655-40 obtained during the 2005 outburst with the two relativistic accretion disk models, KERRBB and KYNBB. Our analysis showed that both models provide identical results with black hole spin measurements, disk temperature, and disk luminosity when the inner edge of the accretion disk is set at the innermost stable circular orbit (ISCO) for the same accretion rates. We could not obtain reasonable fits for similar to 89% of the observations with a fixed black hole spin value at a(*) = 0.7 using both models. Allowing the spin parameter to vary improved the fit statistic significantly with reduced chi(2) values being reduced from -10 to 100 to below 2. Both models revealed black hole spin values varying between 0.52 < a(*) < 0.94, which can be interpreted as a variable inner edge of the disk throughout different accretion states.
4U 1957+11 is a black hole candidate system that has been in a soft X-ray spectral state since its discovery. We present analyses of recent joint NICER and NuSTAR spectra, which are extremely well described by a highly inclined disk accreting into a near maximally spinning black hole. Owing to the broad X-ray coverage of NuSTAR, the fitted spin and inclination are strongly constrained for our hypothesized disk models. The faintest spectra are observed out to 20 keV, even though their hard tail components are almost absent when described with a simple corona. The hard tail increases with luminosity, but shows clear two-track behavior with one track having appreciably stronger tails. The disk spectrum color-correction factor is anticorrelated with the strength of the hard tail (e.g., as measured by the Compton y parameter). Although the spin and inclination parameters are strongly constrained for our chosen model, the mass and distance are degenerate parameters. We use our spectral fits, along with a theoretical prior on color-correction, an observational prior on likely fractional Eddington luminosity, and an observational prior on distance obtained from Gaia studies, to present mass and distance contours for this system. The most likely parameters, given our presumed disk model, suggest a 4.6 M _⊙ black hole at 7.8 kpc observed at luminosities ranging from ≈1.7% to 9% of Eddington. This would place 4U 1957+11 as one of the few actively accreting sources within the mass gap of ≈2–5 M _⊙ where there are few known massive neutron stars or low-mass black holes. Higher mass and distance, however, remain viable.
Winds of massive stars are suspected to be inhomogeneous (or clumpy), which biases the measures of their mass-loss rates. In high-mass X-ray binaries (HMXBs), the compact object can be used as an orbiting X-ray point source to probe the wind and constrain its clumpiness. We perform a spectrotiming analysis of the HMXB OAO 1657–415 with nonsimultaneous NuSTAR and NICER observations. We compute the hardness ratio from the energy-resolved light curves, and, using an adaptive rebinning technique, we thus select appropriate time segments to search for rapid spectral variations on timescales of a few hundred to thousands of seconds. The column density and intensity of the iron Kα line were strongly correlated, and the recorded spectral variations were consistent with accretion from a clumpy wind. We also illustrate a novel framework to measure clump sizes and masses in HMXBs more accurately based on the absorption measurements and orbital parameters of the source. We then discuss the limitations posed by current X-ray spacecraft in such measurements and present prospects with future X-ray missions. We find that the source pulse profiles show a moderate dependence on energy. We identify a previously undetected dip in the pulse profile visible throughout the NuSTAR observation near spin phase 0.15 possibly caused by intrinsic changes in accretion geometry close to the neutron star. We do not find any evidence for the debated cyclotron line at ∼36 keV in the time-averaged or phase-resolved spectra with NuSTAR.
ABSTRACT Black hole X-ray binaries are ideal environments to study the accretion phenomena in strong gravitational potentials. These systems undergo dramatic accretion state transitions and analysis of the X-ray spectra is used to probe the properties of the accretion disc and its evolution. In this work, we present a systematic investigation of ∼1800 spectra obtained by Rossi X-Ray Timing Explorer Proportional Counter Array observations of GRO J1655−40 and LMC X-3 to explore the nature of the accretion disc via non-relativistic and relativistic disc models describing the thermal emission in black hole X-ray binaries. We demonstrate that the non-relativistic modelling throughout an outburst with the phenomenological multicolour disc model DISKBB yields significantly lower and often unphysical inner disc radii and correspondingly higher (∼50–60 per cent) disc temperatures compared to its relativistic counterparts KYNBB and KERRBB. We obtained the dimensionless spin parameters of a* = 0.774 ± 0.069 and a* = 0.752 ± 0.061 for GRO J1655−40 with KERRBB and KYNBB, respectively. We report a spin value of a* = 0.098 ± 0.063 for LMC X-3 using the updated black hole mass of 6.98 M⊙. Both measurements are consistent with the previous studies. Using our results, we highlight the importance of self-consistent modelling of the thermal emission, especially when estimating the spin with the continuum-fitting method which assumes the disc terminates at the innermost stable circular orbit at all times.