We report on INTEGRAL observations of the soft $\gamma$-ray repeater SGR 1935+2154 performed between 2020 April 28 and May 3. Several short bursts with fluence of $\sim10^{-7}-10^{-6}$ erg cm$^{-2}$ were detected by the IBIS instrument in the 20-200 keV range. The burst with the hardest spectrum, discovered and localized in real time by the INTEGRAL Burst Alert System, was spatially and temporally coincident with a short and very bright radio burst detected by the CHIME and STARE2 radio telescopes at 400-800 MHz and 1.4 GHz, respectively. Its lightcurve shows three narrow peaks separated by $\sim$29 ms time intervals, superimposed on a broad pulse lasting $\sim$0.6 s. The brightest peak had a delay of 6.5$\pm$1.0 ms with respect to the 1.4 GHz radio pulse (that coincides with the second and brightest component seen at lower frequencies). The burst spectrum, an exponentially cut-off power law with photon index $\Gamma=0.7_{-0.2}^{+0.4}$ and peak energy $E_p=65\pm5$ keV, is harder than those of the bursts usually observed from this and other magnetars. By the analysis of an expanding dust scattering ring seen in X-rays with the {\it Neil Gehrels Swift Observatory} XRT instrument, we derived a distance of 4.4$_{-1.3}^{+2.8}$ kpc for SGR 1935+2154, independent of its possible association with the supernova remnant G57.2+0.8. At this distance, the burst 20-200 keV fluence of $(6.1\pm 0.3)\times10^{-7}$ erg cm$^{-2}$ corresponds to an isotropic emitted energy of $\sim1.4\times10^{39}$ erg. This is the first burst with a radio counterpart observed from a soft $\gamma$-ray repeater and it strongly supports models based on magnetars that have been proposed for extragalactic fast radio bursts.

Type I X-ray bursts in GS 1826-24, and in several other systems, may induce cooling of the hot inner accretion flow that surrounds the bursting neutron star. Given that GS 1826-24 remained persistently in the hard state over the period 2003-2008 and presented regular bursting properties, we stacked the spectra of the X-ray bursts detected by INTEGRAL (JEM-X and ISGRI) and XMM-Newton (RGS) during that period to study the effect of the burst photons on the properties of the Comptonizing medium. The extended energy range provided by these instruments allows the simultaneous observation of the burst and persistent emission spectra. We detect an overall change in the shape of the persistent emission spectrum in response to the burst photon shower. For the first time, we observe simultaneously a drop in the hard X-ray emission, together with a soft X-ray excess with respect to the burst blackbody emission. The hard X-ray drop can be explained by burst-induced coronal cooling, while the bulk of the soft X-ray excess can be described by fitting the burst emission with an atmosphere model, instead of a simple blackbody model. Traditionally, the persistent emission was assumed to be invariant during X-ray bursts, and more recently to change only in normalization but not in spectral shape; the observed change in the persistent emission level during X-ray bursts may thus trigger the revision of existing neutron star mass-radius constraints, as the derived values rely on the assumption that the persistent emission does not change during X-ray bursts. The traditional burst fitting technique leads to up to a 10% overestimation of the bolometric burst flux in GS 1826-24, which significantly hampers the comparisons of the KEPLER and MESA model against this 'textbook burster'.

During the June 2015 outburst of the black hole binary V404 Cyg, rapid changes in the X-ray brightness and spectra were common. The INTEGRAL monitoring campaign detected spectacular Eddington-limited X-ray flares, but also rapid variations at much lower flux levels. On 2015 June 21 at 20 h 50 min, the 3-10 keV JEM-X data as well as simultaneous optical data started to display a gradual brightening from one of these low-flux states. This was followed 15 min later by an order-of-magnitude increase of flux in the 20-40 keV IBIS/ISGRI light curve in just 15 s. The best-fitting model for both the pre- and post-transition spectra required a Compton-thick partially covering absorber. The absorber parameters remained constant, but the spectral slope varied significantly during the event, with the photon index decreasing from $\Gamma \approx 3.7$ to $\Gamma \approx 2.3$. We propose that the rapid 20-40 keV flux increase was either caused by a spectral state transition that was hidden from our direct view, or that there was a sudden reduction in the amount of Compton down-scattering of the primary X-ray emission in the disk outflow.

We present a multiwavelength analysis of the simultaneous optical and X-ray light curves of the microquasar V404 Cyg during the June 2015 outburst. We have performed a comprehensive analysis of all the INTEGRAL/IBIS, JEM-X, and OMC observations during the brightest epoch of the outburst, along with complementary NuSTAR, AAVSO, and VSNET data, to examine the timing relationship between the simultaneous optical and X-ray light curves, in order to understand the emission mechanisms and physical locations. We have identified all optical flares which have simultaneous X-ray observations, and performed cross-correlation analysis to estimate the time delays between the optical and soft and hard X-ray emission. We have also compared the evolution of the optical and X-ray emission with the hardness-ratios. We have identified several types of behaviour during the outburst. On many occasions, the optical flares occur simultaneously with X-ray flares, but at other times positive and negative time delays between the optical and X-ray emission are measured. We conclude that the observed optical variability is driven by different physical mechanisms, including reprocessing of X-rays in the accretion disc and/or the companion star, interaction of the jet ejections with surrounding material or with previously ejected blobs, and synchrotron emission from the jet.

SAX J1748.9-2021 is an accreting X-ray millisecond pulsar observed in outburst five times since its discovery in 1998. In early October 2017, the source started its sixth outburst, which lasted only ~13 days, significantly shorter than the typical 30 days duration of the previous outbursts. It reached a 0.3-70 keV unabsorbed peak luminosity of $\sim3\times10^{36}$ erg/s. This is the weakest outburst ever reported for this source to date. We analyzed almost simultaneous XMM-Newton, NuSTAR and INTEGRAL observations taken during the decaying phase of its 2017 outburst. We found that the spectral properties of SAX J1748.9-2021 are consistent with an absorbed Comptonization plus a blackbody component. The former, characterized by an electron temperature of ~20 keV, a photon index of ~1.6-1.7 keV and seed photon temperature of 0.44 keV, can be associated to a hot corona or the accretion column, while the latter is more likely originating from the neutron star surface (kT$_{bb}\sim0.6$ keV, R$_{bb}\sim2.5$ km). These findings suggest that SAX J1748.9-2021 was observed in a $hard$ spectral state, as it is typically the case for accreting millisecond pulsars in outburst.

In December 2015 the black hole binary V404 Cyg underwent a secondary outburst after the main June 2015 event. We monitored this re-brightening with the INTEGRAL and Swift satellites, and in this paper we report the results of the time-resolved spectral analysis of these data. The December outburst shared similar characteristics with the June one. The well sampled INTEGRAL light curve shows up to 10 Crab flares, which are separated by relatively weak non-flaring emission phases when compared to the June outburst. The spectra are well described by absorbed Comptonization models, with hard photon indices, $\Gamma \lesssim 2$, and significant detections of a high-energy cut-off only during the bright flares. This is in contrast to the June outburst, where the Comptonization models gave electron temperatures mostly in the 30--50 keV range, while some spectra were soft ($\Gamma \sim 2.5$) without signs of any spectral cut-off. Similarly to the June outburst, we see clear sings of a variable local absorber in the soft energy band covered by Swift/XRT and INTEGRAL/JEM-X, which causes rapid spectral variations observed during the flares. During one flare, both Swift and INTEGRAL captured V404 Cyg in a state where the absorber was nearly Compton thick, $N_\textrm{H} \approx 10^{24}\,\textrm{cm}^{-2}$, and the broad band spectrum was similar to obscured AGN spectra, as seen during the "plateaus" in the June outburst. We conclude that the spectral behaviour of V404 Cyg during the December outburst was analogous with the first few days of the June outburst, both having hard X-ray flares that were intermittently influenced by obscuration due to nearly Compton-thick outflows launched from the accretion disc.

Multiwavelength monitoring of Be/X-ray binaries is crucial to understand the mechanisms producing their outbursts. H 1145-619 is one of these systems, which has recently displayed X-ray activity. We investigate the correlation between the optical emission and the X-ray activity to predict the occurrence of new X-ray outbursts from the inferred state of the circumstellar disc. We have performed a multiwavelength study of H 1145-619 from 1973 to 2017 and present here a global analysis of its variability over the last 40 years. We have used optical spectra from the SAAO, SMARTS and SALT telescopes and optical photometry from INTEGRAL/OMC and ASAS. We also used X-ray observations from INTEGRAL/JEM-X, and IBIS to generate the light curves and combined them with Swift/XRT to extract the X-ray spectra. In addition, we have compiled archival observations and measurements from the literature to complement these data. Comparing the evolution of the optical continuum emission with the Halpha line variability, we have identified three different patterns of optical variability: Global increases and decreases of the optical brightness (observed from 1982 to 1994 and from 2009 to 2017) that can be explained by the dissipation and replenishment of the circumstellar disc, superorbital variations with a period of Psuperorb~590 days (observed in 2002-2009) which seems to be related with the circumstellar disc, and optical outbursts (observed in 1998-1999 and 2002-2005) that we interpret as mass ejections. We have discovered the presence of a retrograde one-armed density wave, which appeared in 2016 and is still present in the circumstellar disc. We have carried out the most complete long-term optical study of the Be/X-ray binary H 1145-619 in correlation with its X-ray activity. We found for the fist time the presence of a retrograde density perturbation in the circumstellar disc of a Be/X-ray binary.

When the upper layer of an accreting neutron star experiences a thermonuclear runaway of helium and hydrogen, it exhibits an X-ray burst of a few keV with a cool-down phase of typically 1~minute. When there is a surplus of hydrogen, hydrogen fusion is expected to simmer during that same minute due to the rp process, which consists of rapid proton captures and slow beta-decays of proton-rich isotopes. We have analyzed the high-quality light curves of 1254 X-ray bursts, obtained with the Proportional Counter Array on the Rossi X-ray Timing Explorer between 1996 and 2012, to systematically study the cooling and rp process. This is a follow-up of a study on a selection of 37 bursts from systems that lack hydrogen and show only cooling during the bursts. We find that the bolometric light curves are well described by the combination of a power law and a one-sided Gaussian. The power-law decay index is between 1.3 and 2.1 and similar to that for the 37-bursts sample. There are individual bursters with a narrower range. The Gaussian is detected in half of all bursts, with a typical standard deviation of 50~s and a fluence ranging up to 60% of the total fluence. The Gaussian appears consistent with being due to the rp process. The Gaussian fluence fraction suggests that the layer where the rp process is active is underabundant in H by a factor of at least five with respect to cosmic abundances. Ninety-four percent of all bursts from ultracompact X-ray binaries lack the Gaussian component, and the remaining 6% are marginal detections. This is consistent with a hydrogen deficiency in these binaries. We find no clear correlation between the power law and Gaussian light-curve components.

The black-hole binary V404 Cyg entered the outburst phase in June 2015 after 26 years of X-ray quiescence, and with its behaviour broke the outburst evolution pattern typical of most black-hole binaries. We observed the entire outburst with the Swift satellite and performed time-resolved spectroscopy of its most active phase, obtaining over a thousand spectra with exposures from tens to hundreds of seconds. All the spectra can be fitted with an absorbed power law model, which most of the time required the presence of a partial covering. A blue-shifted iron-Kalpha line appears in 10% of the spectra together with the signature of high column densities, and about 20% of the spectra seem to show signatures of reflection. None of the spectra showed the unambiguous presence of soft disk-blackbody emission, while the observed bolometric flux exceeded the Eddington value in 3% of the spectra. Our results can be explained assuming that the inner part of the accretion flow is inflated into a slim disk that both hides the innermost (and brightest) regions of the flow, and produces a cold, clumpy, high-density outflow that introduces the high-absorption and fast spectral variability observed. We argue that the black hole in V404 Cyg might have been accreting erratically or even continuously at Eddington/Super-Eddington rates - thus sustaining a surrounding slim disk - while being partly or completely obscured by the inflated disk and its outflow. Hence, the largest flares produced by the source might not be accretion-driven events, but instead the effects of the unveiling of the extremely bright source hidden within the system.

The cooling of type-I X-ray bursts can be used to probe the nuclear burning conditions in neutron star envelopes. The flux decay of the bursts has been traditionally modelled with an exponential, even if theoretical considerations predict power-law-like decays. We have analysed a total of 540 type-I X-ray bursts from five low-mass X-ray binaries observed with the Rossi X-ray Timing Explorer. We grouped the bursts according to the source spectral state during which they were observed (hard or soft), flagging those bursts that showed signs of photospheric radius expansion (PRE). The decay phase of all the bursts were then fitted with a dynamic power-law index method. This method provides a new way of probing the chemical composition of the accreted material. Our results show that in the hydrogen-rich sources the power-law decay index is variable during the burst tails and that simple cooling models qualitatively describe the cooling of presumably helium-rich sources 4U 1728-34 and 3A 1820-303. The cooling in the hydrogen-rich sources 4U 1608-52, 4U 1636-536, and GS 1826-24, instead, is clearly different and depends on the spectral states and whether PRE occurred or not. Especially the hard state bursts behave differently than the models predict, exhibiting a peculiar rise in the cooling index at low burst fluxes, which suggests that the cooling in the tail is much faster than expected. Our results indicate that the drivers of the bursting behaviour are not only the accretion rate and chemical composition of the accreted material, but also the cooling that is somehow linked to the spectral states. The latter suggests that the properties of the burning layers deep in the neutron star envelope might be impacted differently depending on the spectral state.

Typical black hole binaries in outburst show spectral states and transitions, characterized by a clear connection between the inflow onto the black hole and outflow from its vicinity. The transient stellar mass black hole binary V404 Cyg apparently does not fit in this picture. Its outbursts are characterized by intense flares and intermittent low-flux states, with a dynamical range of several orders of magnitude on timescales of hours. During the 2015 June-July X-ray outburst a joint Swift and INTEGRAL observing campaign captured V404 Cyg in one of these low-flux states. The simultaneous Swift/XRT and INTEGRAL/JEM-X/ISGRI spectrum is reminiscent of that of obscured/absorbed AGN. It can be modeled as a Comptonization spectrum, heavily absorbed by a partial covering, high-column density material ($N_\textrm{H} \approx 1.4\times10^{24}\,\textrm{cm}^{-2}$), and a dominant reflection component, including a narrow Iron-K$\alpha$ line. Such spectral distribution can be produced by a geometrically thick accretion flow able to launch a clumpy mass outflow, likely responsible for both the high intrinsic absorption and the intense reflection emission observed. Similarly to what happens in certain obscured AGN, the low-flux states might not be solely related to a decrease in the intrinsic luminosity, but could instead be caused by an almost complete obscuration of the inner accretion flow.

In June 2015, the black hole X-ray binary (BHXRB) V404 Cygni went into outburst for the first time since 1989. Here, we present a comprehensive search for quasi-periodic oscillations (QPOs) of V404 Cygni during its recent outburst, utilizing data from six instruments on board five different X-ray missions: Swift/XRT, Fermi/GBM, Chandra/ACIS, INTEGRAL's IBIS/ISGRI and JEM-X, and NuSTAR. We report the detection of a QPO at 18 mHz simultaneously with both Fermi/GBM and Swift/XRT, another example of a rare but slowly growing new class of mHz-QPOs in BHXRBs linked to sources with a high orbital inclination. Additionally, we find a duo of QPOs in a Chandra/ACIS observation at 73 mHz and 1.03 Hz, as well as a QPO at 136 mHz in a single Swift/XRT observation that can be interpreted as standard Type-C QPOs. Aside from the detected QPOs, there is significant structure in the broadband power, with a strong feature observable in the Chandra observations between 0.1 and 1 Hz. We discuss our results in the context of current models for QPO formation.

Aims. INTEGRAL has been monitoring the Galactic center region for more than a decade. Over this time INTEGRAL has detected hundreds of type-I X-ray bursts from the neutron star low-mass X-ray binary 4U 1728-34, a.k.a. "the slow burster". Our aim is to study the connection between the persistent X-ray spectra and the X-ray burst spectra in a broad spectral range. Methods. We performed spectral modeling of the persistent emission and the X-ray burst emission of 4U 1728-34 using data from the INTEGRAL JEM-X and IBIS/ISGRI instruments. Results. We constructed a hardness intensity diagram to track spectral state variations. In the soft state the energy spectra are characterized by two thermal components - likely from the accretion disc and the boundary/spreading layer - together with a weak hard X-ray tail that we detect in 4U 1728-34 for the first time in the 40 to 80 keV range. In the hard state the source is detected up to 200 keV and the spectrum can be described by a thermal Comptonization model plus an additional component: either a powerlaw tail or reflection. By stacking 123 X-ray bursts in the hard state, we detect emission up to 80 keV during the X-ray bursts. We find that during the bursts the emission above 40 keV decreases by a factor of about three with respect to the persistent emission level. Conclusions. Our results suggest that the enhanced X-ray burst emission changes the spectral properties of the accretion disc in the hard state. The likely cause is an X-ray burst induced cooling of the electrons in the inner hot flow near the neutron star.

GS 0836-429 is a neutron star X-ray transient that displays Type-I X-ray bursts. In 2003 and 2004 it experienced two outbursts in X-rays. We present here an analysis of the system bursting properties during these outbursts. We studied the evolution of the 2003-2004 outbursts in soft X-rays using RXTE (2.5-12 keV; ASM), and in hard X-rays with INTEGRAL (17-80 keV, IBIS/ISGRI). Using data from the JEM-X monitor onboard INTEGRAL we detected 61 Type-I X-ray bursts, and confirm that the source displayed a quasi-periodic burst recurrence time of about 2.3 hours. We improve the characterization of the fuel composition, as well as the description of the typical burst durations and fluences. We estimate the average value of $\alpha$ to be $49\pm\,3$. This value together with the observed burst profiles indicate a regime of a mixed He/H runaway triggered by unstable helium ignition. In addition, we report the detection of four series of double bursts, with burst recurrence times of $\leq\,20$ minutes. The measured recurrence time in double bursts is too short to allow the accretion of enough fresh material, necessary to trigger a Type-I X-ray burst. This suggests the presence of left-over, unburned material from the preceding burst which gets ignited in a time scale of minutes. The energies and time scales of the secondary bursts suggest a lower fraction of hydrogen compared to that estimated for the primary bursts. The persistent emission was roughly constant during the period when the Type I X-ray bursts were detected. We derive an average accretion rate during our observations of $\dot{m}\sim\,8\,\%\,\dot{m}_{Edd}$. The spectrum of the persistent emission can be fit with a non-thermal component, indicative for the source to be in a hard state when the INTEGRAL observations were performed.

In this work we present a system for the automatic classification of the light curves of eclipsing binaries. This system is based on a classification scheme that aims to separate eclipsing binary sistems according to their geometrical configuration in a modified version of the traditional classification scheme. The classification is performed by a Bayesian ensemble of neural networks trained with {\em Hipparcos} data of seven different categories including eccentric binary systems and two types of pulsating light curve morphologies.