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.

PSR B1259-63 is a gamma-ray binary system hosting a radio pulsar orbiting around a O9.5Ve star, LS 2883, with a period of ~3.4 years. The interaction of the pulsar wind with the LS 2883 outflow leads to unpulsed broad band emission in the radio, X-rays, GeV and TeV domains. While the radio, X-ray and TeV light curves show rather similar behaviour, the GeV light curve appears very different with a huge outburst about a month after a periastron. The energy release during this outburst seems to significantly exceed the spin down luminosity of the pulsar and the GeV light curve and energy release varies from one orbit to the next. In this paper we present for the first time the results of optical observations of the system in 2017, and also reanalyze the available X-ray and GeV data. We present a new model in which the GeV data are explained as a combination of the bremsstrahlung and inverse Compton emission from the unshocked and weakly shocked electrons of the pulsar wind. The X-ray and TeV emission is produced by synchrotron and inverse Compton emission of energetic electrons accelerated on a strong shock arising due to stellar/pulsar winds collision. The brightness of the GeV flare is explained in our model as a beaming effect of the energy released in a cone oriented, during the time of flare, in the direction of the observer.

The Laser Interferometer Gravitational Wave Observatory (LIGO) and Virgo Collaboration's Observing Run 3 has demanded the development of widely-applicable tools for gravitational wave follow-up. These tools must address the main challenges of the multi-messenger era, namely covering large localisation regions and quickly identifying decaying transients. To address these challenges, we present a public web interface to assist astronomers in conducting galaxy-targeted follow-up of gravitational wave events by offering a fast and public list of targets post-gravitational wave trigger. After a gravitational wave trigger, the back-end galaxy retrieval algorithm identifies and scores galaxies based on the LIGO and Virgo computed probabilities and properties of the galaxies taken from the Galaxy List for the Advanced Detector Era (GLADE) V2 galaxy catalogue. Within minutes, the user can retrieve, download, and limit ranked galaxy lists from the web application. The algorithm and website have been tested on past gravitational wave events, and execution times have been analysed. The algorithm is being triggered automatically during Observing Run 3 and its features will be extended if needed. The web application was developed using the Python based Flask web framework. The web application is freely available and publicly accessible at gwtool.watchertelescope.ie.

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes.

We present the results from Nordic Optical Telescope and X-shooter follow-up campaigns of the tidal disruption event (TDE) iPTF16fnl, covering the first $\sim$100 days after the transient discovery. We followed the source photometrically until the TDE emission was no longer detected above the host galaxy light. The bolometric luminosity evolution of the TDE is consistent with an exponential decay with e-folding constant t$_{0}$=17.6$\pm$0.2 days. The early time spectra of the transient are dominated by broad He II $\lambda$4686, H$\beta$, H$\alpha$ and N III $\lambda$4100 emission lines. The latter is known to be produced together with the N III $\lambda$4640 in the Bowen fluorescence mechanism. Thanks to the medium resolution X-shooter spectra we have been able to separate the Bowen blend contribution from the broad He II emission line. The detection of the Bowen fluorescence lines in iPTF16fnl place this transient among the N-rich TDE subset. In the late-time X-shooter spectra, narrow emission lines of [O III] and [N II] originating from the host galaxy are detected, suggesting that the host galaxy harbours a weak AGN in its core. The properties of all broad emission lines evolve with time. The equivalent widths follow an exponential decay compatible with the bolometric luminosity evolution. The full-width a half maximum of the broad lines decline with time and the line profiles develop a narrow core at later epochs. Overall, the optical emission of iPTF16fnl can be explained by being produced in an optically thick region in which high densities favour the Bowen fluorescence mechanism and where multiple electron scatterings are responsible for the line broadening.

The discovery of gravitational waves from merging compact objects has opened up a new window to the Universe. Planned third-generation gravitational-wave detectors such as Einstein Telescope will potentially deliver hundreds of such events at redshifts below $z\sim0.1$. Finding electromagnetic counterparts to these events will be a major observational challenge. We demonstrate how Einstein Telescope will provide advance warning of such events on a timescale of hours, based on the low frequency emission from the pre-merger system. In addition, we suggest how this early warning enables prompt identification of any electromagnetic counterpart using the Large Synoptic Survey Telescope.

Long gamma-ray bursts are produced during the dramatic deaths of massive stars with very short lifetimes, meaning that they explode close to the birth place of their progenitors. During a short period they become the most luminous objects observable in the Universe, being perfect beacons to study high-redshift star-forming regions. To use the afterglow of GRB 161023A at a redshift $z=2.710$ as a background source to study the environment of the explosion and the intervening systems along its line-of-sight. r the first time, we complement UV/Optical/NIR spectroscopy with millimetre spectroscopy using ALMA, which allows us to probe the molecular content of the host galaxy. The X-shooter spectrum shows a plethora of absorption features including fine-structure and metastable transitions of Fe, Ni, Si, C and O. We present photometry ranging from 43 s to over 500 days after the burst. We infer a host-galaxy metallicity of [Zn/H] $=-1.11\pm0.07$, which corrected for dust depletion results in [X/H] $=-0.94\pm0.08$. We do not detect molecular features in the ALMA data, but we derive limits on the molecular content of $log(N_{CO}/cm^{-2})<15.7$ and $log(N_{HCO+}/cm^{-2})<13.2$, which are consistent with those that we obtain from the optical spectra, $log(N_{H_2}/cm^{-2})<15.2$ and $log(N_{CO}/cm^{-2})<14.5$. Within the host galaxy we detect three velocity systems through UV/Optical/NIR absorption spectroscopy, all with levels that were excited by the GRB afterglow. We determine the distance from these systems to the GRB to be in the range between 0.7 and 1.0 kpc. The sight-line to GRB 161023A shows 9 independent intervening systems, most of them with multiple components. (Abridged)

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA).

Thanks to its high orbit and a set of complementary detectors providing continuous coverage of the whole sky, the INTEGRAL satellite has unique capabilities for the identification and study of the electromagnetic radiation associated to gravitational waves signals and, more generally, for multi-messenger astrophysics. Here we briefly review the results obtained during the first two observing runs of the advanced LIGO/Virgo interferometers.

Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa <= 0.5 cm2/gm at a velocity of 0.2 +/- 0.1c. The power source is constrained to have a power law slope of beta = -1.2 +/- 0.3, consistent with radioactive powering from r-process nuclides. We identify line features in the spectra that are consistent with light r-process elements (90 < A < 140). As it fades, the transient rapidly becomes red, and emission may have contribution by a higher opacity, lanthanide-rich ejecta component. This indicates that neutron star mergers produce gravitational waves, radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements.

We report the e INTernational Gamma-ray Astrophysics Laboratory (INTEGRAL) detection of the short gamma-ray burst GRB 170817A (discovered by Fermi-GBM) with a signal-to-noise ratio of 4.6, and, for the first time, its association with the gravitational waves (GWs) from binary neutron star (BNS) merging event GW170817 detected by the LIGO and Virgo observatories. The significance of association between the gamma-ray burst observed by INTEGRAL and GW170817 is 3.2 $\sigma$, while the association between the Fermi-GBM and INTEGRAL detections is 4.2 $\sigma$. GRB 170817A was detected by the SPI-ACS instrument about 2 s after the end of the gravitational wave event. We measure a fluence of $(1.4 \pm 0.4 \pm 0.6) \times$10$^{-7}$ erg cm$^{-2})$ (75--2000 keV), where, respectively, the statistical error is given at the 1 $\sigma$ confidence level, and the systematic error corresponds to the uncertainty in the spectral model and instrument response. We also report on the pointed follow-up observations carried out by INTEGRAL, starting 19.5 h after the event, and lasting for 5.4 days. We provide a stringent upper limit on any electromagnetic signal in a very broad energy range, from 3 keV to 8 MeV, constraining the soft gamma-ray afterglow flux to $<7.1\times$10$^{-11}$ erg cm$^{-2}$ s$^{-1}$ (80--300 keV). Exploiting the unique capabilities of INTEGRAL, we constrained the gamma-ray line emission from radioactive decays that are expected to be the principal source of the energy behind a kilonova event following a BNS coalescence. Finally, we put a stringent upper limit on any delayed bursting activity, for example from a newly formed magnetar.

We present a comprehensive temporal and spectral analysis of the long Swift GRB 120327A afterglow data to investigate the possible causes of the observed early time colour variations. We collected data from various instruments/telescopes in different bands (X-rays, ultra- violet, optical and near-infrared) and determined the shapes of the afterglow early-time light curves. We studied the overall temporal behaviour and the spectral energy distributions from early to late times. The ultra-violet, optical, and near-infrared light curves can be modelled with a single power-law component between 200 and 2e4 s after the burst event. The X-ray light curve shows a canonical steep-shallow-steep behaviour, typical of long gamma-ray bursts. At early times a colour variation is observed in the ultra-violet/optical bands, while at very late times a hint of a re-brightening is visible. The observed early time colour change can be explained as a variation in the intrinsic optical spectral index, rather than an evolution of the optical extinction.

During the first observing run of LIGO, two gravitational wave events and one lower-significance trigger (LVT151012) were reported by the LIGO/Virgo collaboration. At the time of LVT151012, the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) was pointing at a region of the sky coincident with the high localization probability area of the event and thus permitted us to search for its electromagnetic counterpart (both prompt and afterglow emission). The imaging instruments on-board INTEGRAL (IBIS/ISGRI, IBIS/PICsIT, SPI, and the two JEM-X modules) have been exploited to attempt the detection of any electromagnetic emission associated with LVT151012 over 3 decades in energy (from 3 keV to 8 MeV). The omni-directional instruments on-board the satellite, i.e. the SPI-ACS and IBIS monitored the entire LVT151012 localization region at energies above 75 keV. We did not find any significant transient source that was spatially and/or temporally coincident with LVT151012, obtaining tight upper limits on the associated hard X-ray and $\gamma$-ray radiation. For typical spectral models, the upper limits on the fluence of the emission from any 1 s long-lasting counterpart of LVT151012 ranges from $F_{\gamma}=$3.5$\times$10$^{-8}$ erg cm$^{-2}$ (20 - 200 keV) to $F_{\gamma}$=7.1$\times$10$^{-7}$ erg cm$^{-2}$ (75 - 2000 keV), constraining the ratio of the isotropic equivalent energy released in the electromagnetic emission to the total energy of the gravitational waves: $E_{75-2000~keV}/E_{GW}<$4.4$\times$10$^{-5}$. Finally, we provide an exhaustive summary of the capabilities of all instruments on-board INTEGRAL to hunt for $\gamma$-ray counterparts of gravitational wave events, exploiting both serendipitous and pointed follow-up observations. This will serve as a reference for all future searches.

A long and intense gamma-ray burst (GRB) was detected by INTEGRAL on July 11 2012 with a duration of ~115s and fluence of 2.8x10^-4 erg cm^-2 in the 20 keV-8 MeV energy range. GRB 120711A was at z~1.405 and produced soft gamma-ray emission (>20 keV) for at least ~10 ks after the trigger. The GRB was observed by several ground-based telescopes that detected a powerful optical flash peaking at an R-band brightness of ~11.5 mag at ~126 s after the trigger. We present a comprehensive temporal and spectral analysis of the long-lasting soft gamma-ray emission detected in the 20-200 keV band with INTEGRAL, the Fermi/LAT post-GRB detection above 100 MeV, the soft X-ray afterglow from XMM-Newton, Chandra, and Swift and the optical/NIR detections from Watcher, Skynet, GROND, and REM. We modelled the long-lasting soft gamma-ray emission using the standard afterglow scenario, which indicates a forward shock origin. The combination of data extending from the NIR to GeV energies suggest that the emission is produced by a broken power-law spectrum consistent with synchrotron radiation. The afterglow is well modelled using a stratified wind-like environment with a density profile k~1.2, suggesting a massive star progenitor (i.e. Wolf-Rayet). The analysis of the reverse and forward shock emission reveals an initial Lorentz factor of ~120-340, a jet half-opening angle of ~2deg-5deg, and a baryon load of ~10^-5-10^-6 Msun consistent with the expectations of the fireball model when the emission is highly relativistic. Long-lasting soft gamma-ray emission from other INTEGRAL GRBs with high peak fluxes, such as GRB 041219A, was not detected, suggesting that a combination of high Lorentz factor, emission above 100 MeV, and possibly a powerful reverse shock are required. Similar long-lasting soft gamma-ray emission has recently been observed from the nearby and extremely bright Fermi/LAT burst GRB 130427A.

We present analysis of 100 ks contiguous XMM-Newton data of the prototypical wind accretor Vela X-1. The observation covered eclipse egress between orbital phases 0.134 and 0.265, during which a giant flare took place, enabling us to study the spectral properties both outside and during the flare. This giant flare with a peak luminosity of $3.92^{+0.42}_{-0.09} \times 10^{37}$ erg s$^{-1}$ allows estimates of the physical parameters of the accreted structure with a mass of $\sim$ $10^{21}$ g. We have been able to model several contributions to the observed spectrum with a phenomenological model formed by three absorbed power laws plus three emission lines. After analysing the variations with orbital phase of the column density of each component, as well as those in the Fe and Ni fluorescence lines, we provide a physical interpretation for each spectral component. Meanwhile, the first two components are two aspects of the principal accretion component from the surface of the neutron star, and the third component seems to be the \textit{X-ray light echo} formed in the stellar wind of the companion.

INTEGRAL's JEM-X instrument offers a very rare opportunity to observe the full prompt X-ray emission from GRBs and the transition to the afterglow phase. A study of prompt X-ray flares in some bursts from the INTEGRAL GRB sample and the early X-ray post-GRB emission from 3-35 keV is presented here. Significant post-GRB emission above 10 keV is observed for GRB 041219A and GRB 081003A.

Reliable timing calibration is essential for the accurate comparison of XMM-Newton light curves with those from other observatories, to ultimately use them to derive precise physical quantities. The XMM-Newton timing calibration is based on pulsar analysis. However, as pulsars show both timing noise and glitches, it is essential to monitor these calibration sources regularly. To this end, the XMM-Newton observatory performs observations twice a year of the Crab pulsar to monitor the absolute timing accuracy of the EPIC-pn camera in the fast Timing and Burst modes. We present the results of this monitoring campaign, comparing XMM-Newton data from the Crab pulsar (PSR B0531+21) with radio measurements. In addition, we use five pulsars (PSR J0537-69, PSR B0540-69, PSR B0833-45, PSR B1509-58 and PSR B1055-52) with periods ranging from 16 ms to 197 ms to verify the relative timing accuracy. We analysed 38 XMM-Newton observations (0.2-12.0 keV) of the Crab taken over the first ten years of the mission and 13 observations from the five complementary pulsars. All the data were processed with the SAS, the XMM-Newton Scientific Analysis Software, version 9.0. Epoch folding techniques coupled with \chi^{2} tests were used to derive relative timing accuracies. The absolute timing accuracy was determined using the Crab data and comparing the time shift between the main X-ray and radio peaks in the phase folded light curves. The relative timing accuracy of XMM-Newton is found to be better than 10^{-8}. The strongest X-ray pulse peak precedes the corresponding radio peak by 306\pm9 \mus, which is in agreement with other high energy observatories such as Chandra, INTEGRAL and RXTE. The derived absolute timing accuracy from our analysis is \pm48 \mus.

(Abridged) GRB 061122 is one of the brightest GRBs detected within INTEGRAL's field of view to date. The two gamma-ray detectors on INTEGRAL were used to investigate the spectral characteristics of GRB 061122. A search for linear polarisation in the prompt emission was carried out using the SPI multiple event data in the energy range 100 keV-1 MeV. The prompt spectrum was best fit by a combination of a blackbody and a power--law model (the quasithermal model), with evidence for high energy emission continuing above 8 MeV. A pseudo-redshift value of pz = 0.95 +/- 0.18 was determined using the spectral fit parameters. The isotropic energy at this pseudo-redshift is 8.5 x 10^{52} erg. The jet opening angle was estimated to be smaller than 2.8 deg or larger than 11.9 deg from the X-ray lightcurve. An upper limit of 60% polarisation was determined for the prompt emission of GRB 061122, using the multiple event data. The high energy emission observed in the spectrum may be due to the reverse shock interacting with the GRB ejecta when it is decelerated by the circumburst medium. This behaviour has been observed in a small fraction of GRBs to date, but is expected to be more commonly observed by the Fermi Gamma-ray Space Telescope. The conditions for polarisation are met if the jet opening angle is less than 2.8 deg, but further constraints on the level of polarisation are not possible.

We present results of a timing analysis of various isolated pulsars using ESA's \emph{XMM-Newton} observatory. Isolated pulsars are useful for calibration purposes because of their stable emission. We have analyzed six pulsars with different pulse profiles in a range of periods between 15 and 200 ms. All observations were made using the \emph{EPIC-pn camera} in its faster modes (Small window, Timing and Burst modes). We investigate the relative timing accuracy of the camera by comparing the pulse periods determined from the \emph{EPIC-pn camera} observations with those from radio observations. As a result of our analysis we conclude that the relative timing accuracy of the \emph{EPIC-pn camera} is of the order of $1\times 10^{-8}$.