Search

Gamma-ray bursts (GRBs) are traditionally classified as either short GRBs with durations $\lesssim 2$ s that are powered by compact object mergers, or long GRBs with durations $\gtrsim 2$ s that powered by the deaths of massive stars. Recent results, however, have challenged this dichotomy and suggest that there exists a population of merger-driven long bursts. One such example, GRB 191019A, has a $t_{90} \approx 64$ s but many of its other properties -- including its host galaxy, afterglow luminosity and lack of associated supernova -- are more consistent with a short GRB. Here we propose an alternative interpretation: that GRB 191019A (which is located in the nucleus of its host) is an atypical jetted tidal disruption event (TDE). In particular, we suggest the short timescale and rapid decline, not expected for standard TDEs, are the result of an "ultra-deep" encounter, in which the star came well within the tidal radius of the black hole and promptly self-intersected, circularised, accreted, and launched a relativistic outflow. This model reproduces the timescale and luminosity through a prompt super-Eddington accretion phase and accounts for the lack of late optical emission. This would make GRB 191019A only the fifth jetted TDE and the first discovered ultra-deep TDE. The ultra-deep TDE model can be distinguished from merger-driven long GRBs via the soft X-ray flash that results from prompt self-intersection of the debris stream; the detection of this flash will be possible with wide-field and soft-X-ray satellites such as $\textit{Einstein Probe}$ or $\textit{SVOM}$., Comment: 10 pages, 2 figures and 1 table. Accepted to ApJL

We present the first release of the Gravitational Wave AfterglowPy Analysis (GWAPA) webtool. GWAPA is designed to provide the community with an interactive tool for rapid analysis of gravitational wave afterglow counterparts and can be extended to the general case of gamma-ray burst afterglows seen at different angles. It is based on the afterglowpy package and allows users to upload observational data and vary afterglow parameters to infer the properties of the explosion. Multiple jet structures, including top hat, Gaussian and power laws, in addition to a spherical outflow model are implemented. A Python script for MCMC fitting is also available to download, with initial guesses taken from GWAPA., Comment: Published in RNAAS. GWAPA will be available at https://gwapa.web.roma2.infn.it/ once essential server repairs are complete

In recent years, searches of archival X-ray data have revealed galaxies exhibiting nuclear quasi-periodic eruptions with periods of several hours. These are reminiscent of the tidal disruption of a star by a supermassive black hole, and the repeated, partial stripping of a white dwarf in an eccentric orbit around a ~10^5 solar mass black hole provides an attractive model. A separate class of periodic nuclear transients, with significantly longer timescales, have recently been discovered optically, and may arise from the partial stripping of a main-sequence star by a ~10^7 solar mass black hole. No clear connection between these classes has been made. We present the discovery of an X-ray nuclear transient which shows quasi-periodic outbursts with a period of weeks. We discuss possible origins for the emission, and propose that this system bridges the two existing classes outlined above. This discovery was made possible by the rapid identification, dissemination and follow up of an X-ray transient found by the new live \swift-XRT transient detector, demonstrating the importance of low-latency, sensitive searches for X-ray transients., Comment: To be published in Nature Astronomy at 1600 BST on September 7th. This version for arXiv includes the main article, Methods and Supplementary Information combined into a single file

Nebular-phase observations of peculiar Type Ia supernovae (SNe Ia) provide important constraints on progenitor scenarios and explosion dynamics for both these rare SNe and the more common, cosmologically useful SNe Ia. We present observations from an extensive ground-based and space-based follow-up campaign to characterize SN 2022pul, a "super-Chandrasekhar" mass SN Ia (alternatively "03fg-like" SN), from before peak brightness to well into the nebular phase across optical to mid-infrared (MIR) wavelengths. The early rise of the light curve is atypical, exhibiting two distinct components, consistent with SN Ia ejecta interacting with dense carbon-oxygen rich circumstellar material (CSM). In the optical, SN 2022pul is most similar to SN 2012dn, having a low estimated peak luminosity ($M_{B}=-18.9$ mag) and high photospheric velocity relative to other 03fg-like SNe. In the nebular phase, SN 2022pul adds to the increasing diversity of the 03fg-like subclass. From 168 to 336 days after peak $B$-band brightness, SN 2022pul exhibits asymmetric and narrow emission from [O I] $\lambda\lambda 6300,\ 6364$ (${\rm FWHM} \approx 2{,}000$ km s$^{-1}$), strong, broad emission from [Ca II] $\lambda\lambda 7291,\ 7323$ (${\rm FWHM} \approx 7{,}300$ km s$^{-1}$), and a rapid Fe III to Fe II ionization change. Finally, we present the first-ever optical-to-mid-infrared (MIR) nebular spectrum of an 03fg-like SN Ia using data from JWST. In the MIR, strong lines of neon and argon, weak emission from stable nickel, and strong thermal dust emission (with $T \approx 500$ K), combined with prominent [O I] in the optical, suggest that SN 2022pul was produced by a white dwarf merger within carbon/oxygen-rich CSM., Comment: 23 pages, 11 figures, submitted to ApJ

We present an analysis of ground-based and JWST observations of SN~2022pul, a peculiar "03fg-like" (or "super-Chandrasekhar") Type Ia supernova (SN Ia), in the nebular phase at 338d post explosion. Our combined spectrum continuously covers 0.4--14 $\mu$m and includes the first mid-infrared spectrum of an 03fg-like SN Ia. Compared to normal SN Ia 2021aefx, SN 2022pul exhibits a lower mean ionization state, asymmetric emission-line profiles, stronger emission from the intermediate-mass elements (IMEs) argon and calcium, weaker emission from iron-group elements (IGEs), and the first unambiguous detection of neon in a SN Ia. Strong, broad, centrally peaked [Ne II] line at 12.81 $\mu$m was previously predicted as a hallmark of "violent merger'' SN Ia models, where dynamical interaction between two sub-$M_{ch}$ white dwarfs (WDs) causes disruption of the lower mass WD and detonation of the other. The violent merger scenario was already a leading hypothesis for 03fg-like SNe Ia; in SN 2022pul it can explain the large-scale ejecta asymmetries seen between the IMEs and IGEs and the central location of narrow oxygen and broad neon. We modify extant models to add clumping of the ejecta to better reproduce the optical iron emission, and add mass in the innermost region ($< 2000$ km s$^{-1}$) to account for the observed narrow [O I]~$\lambda\lambda6300$, 6364 emission. A violent WD-WD merger explains many of the observations of SN 2022pul, and our results favor this model interpretation for the subclass of 03fg-like SN Ia., Comment: 20 pages, 10 figures, published in ApJ

We report the discovery of Swift J221951-484240 (hereafter: J221951), a luminous slow-evolving blue transient that was detected by the Neil Gehrels Swift Observatory Ultra-violet/Optical Telescope (Swift/UVOT) during the follow-up of Gravitational Wave alert S190930t, to which it is unrelated. Swift/UVOT photometry shows the UV spectral energy distribution of the transient to be well modelled by a slowly shrinking black body with an approximately constant temperature of T~2.5x10^4 K. At a redshift z=0.5205, J221951 had a peak absolute magnitude of M_u,AB = -23 mag, peak bolometric luminosity L_max=1.1x10^45 erg s^-1 and a total radiated energy of E>2.6x10^52 erg. The archival WISE IR photometry shows a slow rise prior to a peak near the discovery date. Spectroscopic UV observations display broad absorption lines in N V and O VI, pointing toward an outflow at coronal temperatures. The lack of emission in the higher H~Lyman lines, N I and other neutral lines is consistent with a viewing angle close to the plane of the accretion or debris disc. The origin of J221951 can not be determined with certainty but has properties consistent with a tidal disruption event and the turn-on of an active galactic nucleus., Comment: 37 pages (25 main + 12 supplementary), submitted to MNRAS

Type Ia supernovae remain poorly understood despite decades of investigation. Massive computationally intensive hydrodynamic simulations have been developed and run to model an ever-growing number of proposed progenitor channels. Further complicating the matter, a large number of sub-types of Type Ia supernovae have been identified in recent decades. Due to the massive computational load required, inference of the internal structure of Type Ia supernovae ejecta directly from observations using simulations has previously been computationally intractable. However, deep-learning emulators for radiation transport simulations have alleviated such barriers. We perform abundance tomography on 40 Type Ia supernovae from optical spectra using the radiative transfer code TARDIS accelerated by the probabilistic DALEK deep-learning emulator. We apply a parametric model of potential ejecta structures to comparatively investigate abundance distributions and internal ionization fractions of intermediate-mass elements between normal and 1991T-like Type Ia supernovae. Our inference shows that 1991T-like Type Ia supernovae are under-abundant in the typical intermediate mass elements that heavily contribute to the spectral line formation seen in normal Type Ia supernovae at early times. Additionally, we find that the intermediate-mass elements present in 1991T-like Type Ia supernovae are highly ionized compared to those in the normal Type Ia population. Finally, we conclude that the transition between normal and 1991T-like Type Ia supernovae appears to be continuous observationally and that the observed differences come out of a combination of both abundance and ionization fractions in these supernovae populations., Comment: 15 Pages, 6 Figures, Article, Submitted to ApJ

GRB 201015A is a peculiarly low luminosity, spectrally soft gamma-ray burst (GRB), with $T_{\rm 90} = 9.8 \pm 3.5$ s (time interval of detection of 90\% of photons from the GRB), and an associated supernova (likely to be type Ic or Ic-BL). GRB 201015A has an isotropic energy $E_{\gamma,\rm iso} = 1.75 ^{+0.60} _{-0.53} \times 10^{50}$ erg, and photon index $\Gamma = 3.00 ^{+0.50} _{-0.42}$ (15-150 keV). It follows the Amati relation, a correlation between $E_{\gamma,\rm iso}$ and spectral peak energy $E_{\rm p}$ followed by long GRBs. It appears exceptionally soft based on $\Gamma$, the hardness ratio of HR = $0.47 \pm 0.24$, and low-$E_{\rm p}$, so we have compared it to other GRBs sharing these properties. These events can be explained by shock breakout, poorly collimated jets, and off-axis viewing. Follow-up observations of the afterglow taken in the X-ray, optical, and radio, reveal a surprisingly late flattening in the X-ray from $t = (2.61 \pm 1.27)\times 10^4$ s to $t = 1.67 ^{+1.14} _{-0.65} \times 10^6$ s. We fit the data to closure relations describing the synchrotron emission, finding the electron spectral index to be $p = 2.42 ^{+0.44} _{-0.30}$, and evidence of late-time energy injection with coefficient $q = 0.24 ^{+0.24} _{-0.18}$. The jet half opening angle lower limit ($\theta_{j} \ge 16^{\circ}$) is inferred from the non-detection of a jet break. The launch of SVOM and Einstein Probe in 2023, should enable detection of more low luminosity events like this, providing a fuller picture of the variety of GRBs., Comment: 15 pages, 4 figures

Type Iax supernovae (SN Iax) are the largest known class of peculiar white dwarf supernovae, distinct from normal Type Ia supernovae (SN Ia). The unique properties of SN Iax, especially their strong photospheric lines out to extremely late times, allow us to model their optical spectra and derive physical parameters for the long-lasting photosphere. We present an extensive spectral timeseries, including 21 new spectra, of SN Iax 2014dt from +11 to +562 days after maximum light. We are able to reproduce the entire timeseries with a self-consistent, nearly unaltered deflagration explosion model from Fink et al. (2014) using TARDIS, an open-source radiative transfer code (Kerzendorf & Sim 2014; Kerzendorf et al. 2023). We find that the photospheric velocity of SN 2014dt slows its evolution between +64 and +148 days, which closely overlaps the phase when we see SN 2014dt diverge from the normal spectral evolution of SN Ia (+90 to +150 days). The photospheric velocity at these epochs, ~400$-$1000 km s$^{-1}$, may demarcate a boundary within the ejecta below which the physics of SN Iax and normal SN Ia differ. Our results suggest that SN 2014dt is consistent with a weak deflagration explosion model that leaves behind a bound remnant and drives an optically thick, quasi-steady-state wind creating the photospheric lines at late times. The data also suggest that this wind may weaken at epochs past +450 days, perhaps indicating a radioactive power source that has decayed away., Comment: Accepted to ApJ, 22 pages, 8 figures, 3 tables

We present JWST near- and mid-infrared spectroscopic observations of the nearby normal Type Ia supernova SN 2021aefx in the nebular phase at $+255$ days past maximum light. Our Near Infrared Spectrograph (NIRSpec) and Mid Infrared Instrument (MIRI) observations, combined with ground-based optical data from the South African Large Telescope (SALT), constitute the first complete optical $+$ NIR $+$ MIR nebular SN Ia spectrum covering 0.3$-$14 $\mu$m. This spectrum unveils the previously unobserved 2.5$-$5 $\mu$m region, revealing strong nebular iron and stable nickel emission, indicative of high-density burning that can constrain the progenitor mass. The data show a significant improvement in sensitivity and resolution compared to previous Spitzer MIR data. We identify numerous NIR and MIR nebular emission lines from iron-group elements and as well as lines from the intermediate-mass element argon. The argon lines extend to higher velocities than the iron-group elements, suggesting stratified ejecta that are a hallmark of delayed-detonation or double-detonation SN Ia models. We present fits to simple geometric line profiles to features beyond 1.2 $\mu$m and find that most lines are consistent with Gaussian or spherical emission distributions, while the [Ar III] 8.99 $\mu$m line has a distinctively flat-topped profile indicating a thick spherical shell of emission. Using our line profile fits, we investigate the emissivity structure of SN 2021aefx and measure kinematic properties. Continued observations of SN 2021aefx and other SNe Ia with JWST will be transformative to the study of SN Ia composition, ionization structure, density, and temperature, and will provide important constraints on SN Ia progenitor and explosion models., Comment: published in ApJ Letters, 17 pages, 12 figures

We present simulated optical light curves of super-Eddington tidal disruption events (TDEs) using the zero-Bernoulli accretion (ZEBRA) flow model, which proposes that during the super-Eddington phase, the disc is quasi-spherical, radiation-pressure dominated, and accompanied by the production of strong jets. We construct light curves for both on- and off-axis (with respect to the jet) observers to account for the anisotropic nature of the jetted emission. We find that at optical wavelengths, emission from the accretion flow is orders of magnitude brighter than that produced by the jet, even with boosting from synchrotron self-Compton. Comparing to the observed jetted TDE Swift J2058.4+0516, we find that the ZEBRA model accurately captures the timescale for which accretion remains super-Eddington and reproduces the luminosity of the transient. However, we find the shape of the light curves deviate at early times and the radius and temperature of our modelled ZEBRA are $\sim2.7 - 4.1$ times smaller and $\sim1.4 - 2.3$ times larger, respectively, than observed. We suggest that this indicates the ZEBRA inflates more, and more rapidly, than currently predicted by the model, and we discuss possible extensions to the model to account for this. Such refinements, coupled with valuable new data from upcoming large scale surveys, could help to resolve the nature of super-Eddington TDEs and how they are powered., Comment: Accepted to MNRAS. 9 pages and 5 figures

The Second Swift-XRT Point Source catalogue offers a combination of sky coverage and sensitivity and presents an invaluable opportunity for transient discovery. We search the catalogue at the positions of inactive and active galaxies, and identify transient candidates by comparison with XMM-Newton and ROSAT. We recover 167 previously known transients and find 19 sources consistent with being new sources, estimating a completeness of $\sim65\%$. These 19 new sources are split approximately equally between inactive and active hosts and their peak X-ray luminosities span $\sim 10^{42} - 10^{47}$ erg s$^{-1}$. We find eight are best fit with non-thermal spectral models and one with a blackbody. We also discuss our methodology and its application to the forthcoming Living Swift-XRT Point Source catalogue for the potential near real time serendipitous discovery of $\sim$ few new X-ray transients per year., Comment: 12 pages, 7 figures, 8 tables plus appendix with 2 additional figures. Accepted by MNRAS

The progenitor systems and explosion mechanism of Type Ia supernovae are still unknown. Currently favoured progenitors include double-degenerate systems consisting of two carbon-oxygen white dwarfs with thin helium shells. In the double-detonation scenario, violent accretion leads to a helium detonation on the more massive primary white dwarf that turns into a carbon detonation in its core and explodes it. We investigate the fate of the secondary white dwarf, focusing on changes of the ejecta and observables of the explosion if the secondary explodes as well rather than survives. We simulate a binary system of a $1.05\,\mathrm{M_\odot}$ and a $0.7\,\mathrm{M_\odot}$ carbon-oxygen white dwarf with $0.03\,\mathrm{M_\odot}$ helium shells each. We follow the system self-consistently from inspiral to ignition, through the explosion, to synthetic observables. We confirm that the primary white dwarf explodes self-consistently. The helium detonation around the secondary white dwarf, however, fails to ignite a carbon detonation. We restart the simulation igniting the carbon detonation in the secondary white dwarf by hand and compare the ejecta and observables of both explosions. We find that the outer ejecta at $v~>~15\,000$\,km\,s$^{-1}$ are indistinguishable. Light curves and spectra are very similar until $\sim~40\,\mathrm{d}$ after explosion and the ejecta are much more spherical than violent merger models. The inner ejecta differ significantly slowing down the decline rate of the bolometric light curve after maximum of the model with a secondary explosion by $\sim20$ per cent. We expect future synthetic 3D nebular spectra to confirm or rule out either model., Comment: 12 pages, 7 figures, accepted to MNRAS, comments welcome

Understanding the evolution of massive binary stars requires accurate estimates of their masses. This understanding is critically important because massive star evolution can potentially lead to gravitational wave sources such as binary black holes or neutron stars. For Wolf-Rayet stars with optically thick stellar winds, their masses can only be determined with accurate inclination angle estimates from binary systems which have spectroscopic $M \sin i$ measurements. Orbitally-phased polarization signals can encode the inclination angle of binary systems, where the Wolf-Rayet winds act as scattering regions. We investigated four Wolf-Rayet + O star binary systems, WR 42, WR 79, WR 127, and WR 153, with publicly available phased polarization data to estimate their masses. To avoid the biases present in analytic models of polarization while retaining computational expediency, we used a Monte Carlo radiative transfer model accurately emulated by a neural network. We used the emulated model to investigate the posterior distribution of parameters of our four systems. Our mass estimates calculated from the estimated inclination angles put strong constraints on existing mass estimates for three of the systems, and disagrees with the existing mass estimates for WR 153. We recommend a concerted effort to obtain polarization observations that can be used to estimate the masses of Wolf-Rayet binary systems and increase our understanding of their evolutionary paths., Comment: 24 pages, 13 figures, accepted to ApJ

In this paper, we report on the observational performance of the Swift Ultra-violet/Optical Telescope (UVOT) in response to the Gravitational Wave alerts announced by the Advanced Laser Interferometer Gravitational Wave Observatory and the Advanced Virgo detector during the O3 period. We provide the observational strategy for follow-up of GW alerts and provide an overview of the processing and analysis of candidate optical/UV sources. For the O3 period, we also provide a statistical overview and report on serendipitous sources discovered by Swift/UVOT. Swift followed 18 gravitational-wave candidate alerts, with UVOT observing a total of 424 deg^2. We found 27 sources that changed in magnitude at the 3 sigma level compared with archival u or g-band catalogued values. Swift/UVOT also followed up a further 13 sources reported by other facilities during the O3 period. Using catalogue information, we divided these 40 sources into five initial classifications: 11 candidate active galactic nuclei (AGN)/quasars, 3 Cataclysmic Variables (CVs), 9 supernovae, 11 unidentified sources that had archival photometry and 6 uncatalogued sources for which no archival photometry was available. We have no strong evidence to identify any of these transients as counterparts to the GW events. The 17 unclassified sources are likely a mix of AGN and a class of fast-evolving transient, and one source may be a CV., Comment: 25 pages, 6 figures and 5 tables. Submitted to MNRAS. Supplementary contains 23 pages with 8 figures and 1 table

Manual fits to spectral times series of Type Ia supernovae have provided a method of reconstructing the explosion from a parametric model but due to lack of information about model uncertainties or parameter degeneracies direct comparison between theory and observation is difficult. In order to mitigate this important problem we present a new way to probabilistically reconstruct the outer ejecta of the normal Type Ia supernova SN 2002bo. A single epoch spectrum, taken 10 days before maximum light, is fit by a 13-parameter model describing the elemental composition of the ejecta and the explosion physics (density, temperature, velocity, and explosion epoch). Model evaluation is performed through the application of a novel rapid spectral synthesis technique in which the radiative transfer code, TARDIS, is accelerated by a machine-learning framework. Analysis of the posterior distribution reveals a complex and degenerate parameter space and allows direct comparison to various hydrodynamic models. Our analysis favors detonation over deflagration scenarios and we find that our technique offers a novel way to compare simulation to observation., Comment: 12 Pages, 4 Figures, submitted to AAS Journals, comments and constructive criticism welcome

THESEUS is a medium size space mission of the European Space Agency, currently under evaluation for a possible launch in 2032. Its main objectives are to investigate the early Universe through the observation of gamma-ray bursts and to study the gravitational waves electromagnetic counterparts and neutrino events. On the other hand, its instruments, which include a wide field of view X-ray (0.3-5 keV) telescope based on lobster-eye focusing optics and a gamma-ray spectrometer with imaging capabilities in the 2-150 keV range, are also ideal for carrying out unprecedented studies in time domain astrophysics. In addition, the presence onboard of a 70 cm near infrared telescope will allow simultaneous multi-wavelegth studies. Here we present the THESEUS capabilities for studying the time variability of different classes of sources in parallel to, and without affecting, the gamma-ray bursts hunt., Comment: Submitted to Experimental Astronomy

Multi-messenger astrophysics is becoming a major avenue to explore the Universe, with the potential to span a vast range of redshifts. The growing synergies between different probes is opening new frontiers, which promise profound insights into several aspects of fundamental physics and cosmology. In this context, THESEUS will play a central role during the 2030s in detecting and localizing the electromagnetic counterparts of gravitational wave and neutrino sources that the unprecedented sensitivity of next generation detectors will discover at much higher rates than the present. Here, we review the most important target signals from multi-messenger sources that THESEUS will be able to detect and characterize, discussing detection rate expectations and scientific impact., Comment: Submitted to Experimental Astronomy

THESEUS, one of the two space mission concepts being studied by ESA as candidates for next M5 mission within its Comsic Vision programme, aims at fully exploiting Gamma-Ray Bursts (GRB) to solve key questions about the early Universe, as well as becoming a cornerstone of multi-messenger and time-domain astrophysics. By investigating the first billion years of the Universe through high-redshift GRBs, THESEUS will shed light on the main open issues in modern cosmology, such as the population of primordial low mass and luminosity galaxies, sources and evolution of cosmic re-ionization, SFR and metallicity evolution up to the "cosmic dawn" and across Pop-III stars. At the same time, the mission will provide a substantial advancement of multi-messenger and time-domain astrophysics by enabling the identification, accurate localisation and study of electromagnetic counterparts to sources of gravitational waves and neutrinos, which will be routinely detected in the late '20s and early '30s by the second and third generation Gravitational Wave (GW) interferometers and future neutrino detectors, as well as of all kinds of GRBs and most classes of other X/gamma-ray transient sources. In all these cases, THESEUS will provide great synergies with future large observing facilities in the multi-messenger domain. A Guest Observer programme, comprising Target of Opportunity (ToO) observations, will expand the science return of the mission, to include, e.g., solar system minor bodies, exoplanets, and AGN., Comment: Submitted to Experimental Astronomy. arXiv admin note: text overlap with arXiv:2102.08702

The proposed THESEUS mission will vastly expand the capabilities to monitor the high-energy sky, and will exploit large samples of gamma-ray bursts to probe the early Universe back to the first generation of stars, and to advance multi-messenger astrophysics by detecting and localizing the counterparts of gravitational waves and cosmic neutrino sources. The combination and coordination of these activities with multi-wavelength, multi-messenger facilities expected to be operating in the thirties will open new avenues of exploration in many areas of astrophysics, cosmology and fundamental physics, thus adding considerable strength to the overall scientific impact of THESEUS and these facilities. We discuss here a number of these powerful synergies., Comment: Revised version after submission to Experimental Astronomy

At peak, long-duration gamma-ray bursts are the most luminous sources of electromagnetic radiation known. Since their progenitors are massive stars, they provide a tracer of star formation and star-forming galaxies over the whole of cosmic history. Their bright power-law afterglows provide ideal backlights for absorption studies of the interstellar and intergalactic medium back to the reionization era. The proposed THESEUS mission is designed to detect large samples of GRBs at $z>6$ in the 2030s, at a time when supporting observations with major next generation facilities will be possible, thus enabling a range of transformative science. THESEUS will allow us to explore the faint end of the luminosity function of galaxies and the star formation rate density to high redshifts; constrain the progress of re-ionisation beyond $z\gtrsim6$; study in detail early chemical enrichment from stellar explosions, including signatures of Population III stars; and potentially characterize the dark energy equation of state at the highest redshifts., Comment: Submitted to Experimental Astronomy

THESEUS is a space mission concept, currently under Phase A study by ESA as candidate M5 mission, aiming at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. In addition to fully exploiting high-redshift GRBs for cosmology (pop-III stars, cosmic re-ionization, SFR and metallicity evolution up to the "cosmic dawn"), THESEUS will allow the identification and study of the electromagnetic counterparts to sources of gravitational waves which will be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/aVirgo, LISA, KAGRA, and Einstein Telescope (ET), as well as of most classes of X/gamma-ray transient sources, thus providing an ideal synergy with the large e.m. facilities of the near future like, e.g., LSST, ELT, TMT, SKA, CTA, ATHENA. These breakthrough scientific objectives will be achieved by an unprecedented combination of X/gamma-ray monitors, providing the capabilities of detecting and accurately localize and kind of GRBs and may classes of transient in an energy band as large as 0.1 keV - 10 MeV, with an on-board NIR telescope providing detection, localization (arcsec) and redshift measurement of the NIR counterpart. A Guest Observer programme, further improving the scientific return and community involvement is also envisaged. We summarize the main scientific requirements of the mission and provide an overview of the updated concept, design (instruments and spacecraft) and mission profile., Comment: Proceedings of the SPIE 2020, paper 11444-302

On 2020 February 24, during their third observing run ("O3"), the Laser Interferometer Gravitational-wave Observatory and Virgo Collaboration (LVC) detected S200224ca: a candidate gravitational wave (GW) event produced by a binary black hole (BBH) merger. This event was one of the best-localized compact binary coalescences detected in O3 (with 50%/90% error regions of 13/72 deg$^2$), and so the Neil Gehrels Swift Observatory performed rapid near-UV/X-ray follow-up observations. Swift-XRT and UVOT covered approximately 79.2% and 62.4% (respectively) of the GW error region, making S200224ca the BBH event most thoroughly followed-up in near-UV (u-band) and X-ray to date. No likely EM counterparts to the GW event were found by the Swift BAT, XRT, or UVOT, nor by other observatories. Here we report on the results of our searches for an EM counterpart, both in the BAT data near the time of the merger, and in follow-up UVOT/XRT observations. We also discuss the upper limits we can place on EM radiation from S200224ca, and the implications these limits have on the physics of BBH mergers. Namely, we place a shallow upper limit on the dimensionless BH charge, $\hat{q} < 1.4 \times10^{-4}$, and an upper limit on the isotropic-equivalent energy of a blast wave $E < 4.1\times10^{51}$ erg (assuming typical GRB parameters)., Comment: 14 pages, 6 figures, accepted for publication in ApJ

The Neil Gehrels Swift Observatory followed up 18 gravitational wave (GW) triggers from the LIGO/Virgo collaboration during the O3 observing run in 2019/2020, performing approximately 6500 pointings in total. Of these events, four were finally classified (if real) as binary black hole (BH) triggers, six as binary neutron star (NS) events, two each of NSBH and Mass Gap triggers, one an unmodelled (Burst) trigger, and the remaining three were subsequently retracted. Thus far, four of these O3 triggers have been formally confirmed as real gravitational wave events. While no likely electromagnetic counterparts to any of these GW events have been identified in the X-ray data (to an average upper limit of 3.60 x 10^{-12} erg cm^{-2} s^{-1} over 0.3-10 keV), or at other wavelengths, we present a summary of all the Swift-XRT observations performed during O3, together with typical upper limits for each trigger observed. The majority of X-ray sources detected during O3 were previously uncatalogued; while some of these will be new (transient) sources, others are simply too faint to have been detected by earlier survey missions such as ROSAT. The all-sky survey currently being performed by eROSITA will be a very useful comparison for future observing runs, reducing the number of apparent candidate X-ray counterparts by up to 95 per cent., Comment: 23 pages (including 4 pages of references, and a 4 page table in the appendix), 5 figures (4 in colour), accepted for publication in MNRAS. (Replaced due to annoying spelling typo in the abstract.)

We report the results of optical follow-up observations of 29 gravitational-wave triggers during the first half of the LIGO-Virgo Collaboration (LVC) O3 run with the Gravitational-wave Optical Transient Observer (GOTO) in its prototype 4-telescope configuration (GOTO-4). While no viable electromagnetic counterpart candidate was identified, we estimate our 3D (volumetric) coverage using test light curves of on- and off-axis gamma-ray bursts and kilonovae. In cases where the source region was observable immediately, GOTO-4 was able to respond to a GW alert in less than a minute. The average time of first observation was $8.79$ hours after receiving an alert ($9.90$ hours after trigger). A mean of $732.3$ square degrees were tiled per event, representing on average $45.3$ per cent of the LVC probability map, or $70.3$ per cent of the observable probability. This coverage will further improve as the facility scales up alongside the localisation performance of the evolving gravitational-wave detector network. Even in its 4-telescope prototype configuration, GOTO is capable of detecting AT2017gfo-like kilonovae beyond 200~Mpc in favourable observing conditions. We cannot currently place meaningful electromagnetic limits on the population of distant ($\hat{D}_L = 1.3$~Gpc) binary black hole mergers because our test models are too faint to recover at this distance. However, as GOTO is upgraded towards its full 32-telescope, 2 node (La Palma \& Australia) configuration, it is expected to be sufficiently sensitive to cover the predicted O4 binary neutron star merger volume, and will be able to respond to both northern and southern triggers., Comment: 15 pages, 7 figures, 3 tables. Accepted for publication in MNRAS. Author's final submitted version

On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. Preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical/near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS and VINROUGE projects also carried out a search on this event. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN) possibly generated by this NS-BH merger, and for the strategy of future searches. Altogether, our observations allow us to exclude a KN with large ejecta mass $M\gtrsim 0.1\,\mathrm{M_\odot}$ to a high ($>90\%$) confidence, and we can exclude much smaller masses in a subsample of our observations. This disfavours the tidal disruption of the neutron star during the merger. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundreds Mpc will be detected only by large facilities with both high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event., Comment: 52 pages, revised version now accepted for publication in A&A. Abstract abridged to meet arXiv requirements

The Neil Gehrels Swift Observatory carried out prompt searches for gravitational wave (GW) events detected by the LIGO/Virgo Collaboration (LVC) during the second observing run ("O2"). Swift performed extensive tiling of eight LVC triggers, two of which had very low false-alarm rates (GW 170814 and the epochal GW 170817), indicating a high confidence of being astrophysical in origin; the latter was the first GW event to have an electromagnetic counterpart detected. In this paper we describe the follow-up performed during O2 and the results of our searches. No GW electromagnetic counterparts were detected; this result is expected, as GW 170817 remained the only astrophysical event containing at least one neutron star after LVC's later retraction of some events. A number of X-ray sources were detected, with the majority of identified sources being active galactic nuclei. We discuss the detection rate of transient X-ray sources and their implications in the O2 tiling searches. Finally, we describe the lessons learned during O2, and how these are being used to improve the \swift\ follow-up of GW events. In particular, we simulate a population of GRB afterglows to evaluate our source ranking system's ability to differentiate them from unrelated and uncatalogued X-ray sources. We find that $\approx$60-70% of afterglows whose jets are oriented towards Earth will be given high rank (i.e., "interesting" designation) by the completion of our second follow-up phase (assuming their location in the sky was observed), but that this fraction can be increased to nearly 100% by performing a third follow-up observation of sources exhibiting fading behavior., Comment: 18 pages, 4 figures, 2 tables, accepted for publication in ApJS

We present X-ray and optical observations of the short duration gamma-ray burst GRB 071227 and its host at $z=0.381$, obtained using \textit{Swift}, Gemini South and the Very Large Telescope. We identify a short-lived and moderately bright optical transient, with flux significantly in excess of that expected from a simple extrapolation of the X-ray spectrum at 0.2-0.3 days after burst. We fit the SED with afterglow models allowing for high extinction and thermal emission models that approximate a kilonova to assess the excess' origins. While some kilonova contribution is plausible, it is not favoured due to the low temperature and high luminosity required, implying superluminal expansion and a large ejecta mass of $\sim 0.1$ M$_{\odot}$. We find, instead, that the transient is broadly consistent with power-law spectra with additional dust extinction of $E(B-V)\sim0.4$ mag, although a possibly thermal excess remains in the \textit{z}-band. We investigate the host, a spiral galaxy with an edge-on orientation, resolving its spectrum along its major axis to construct the galaxy rotation curve and analyse the star formation and chemical properties. The integrated host emission shows evidence for high extinction, consistent with the afterglow findings. The metallicity and extinction are consistent with previous studies of this host and indicate the galaxy is a typical, but dusty, late-type SGRB host., Comment: 16 pages, 13 figures, accepted to MNRAS

We report our identification of the optical afterglow and host galaxy of the short-duration gamma-ray burst GRB 160821B. The spectroscopic redshift of the host is $z=0.162$, making it one of the lowest redshift sGRBs identified by Swift. Our intensive follow-up campaign using a range of ground-based facilities as well as HST, XMM and Swift, shows evidence for a late-time excess of optical and near-infrared emission in addition to a complex afterglow. The afterglow light-curve at X-ray frequencies reveals a narrow jet, $\theta_j\sim1.9^{+0.10}_{-0.03}$ deg, that is refreshed at $>1$ day post-burst by a slower outflow with significantly more energy than the initial outflow that produced the main GRB. Observations of the 5 GHz radio afterglow shows a reverse shock into a mildly magnetised shell. The optical and near-infrared excess is fainter than AT2017gfo associated with GW170817, and is well explained by a kilonova with dynamic ejecta mass $M_{\rm dyn}=(1.0\pm0.6)\times10^{-3}$ M$_{\odot}$ and a secular (postmerger) ejecta mass with $M_{\rm pm}=(1.0\pm0.6)\times10^{-2}$ M$_\odot$, consistent with a binary neutron star merger resulting in a short-lived massive neutron star. This optical and near-infrared dataset provides the best-sampled kilonova light-curve without a gravitational wave trigger to date., Comment: 17 pages, 6 figures, Version accepted by ApJ

We report the X-ray data analysis of two transient ultraluminous X-ray sources (ULXs, hereafter X1 and X2) located in the nearby galaxy NGC 7090. While they were not detected in the 2004 XMM-Newton and 2005 Chandra observations, their 0.3-10 keV X-ray luminosities reached $>3\times10^{39}\,\mathrm{erg\,s^{-1}}$ in later XMM-Newton or Swift observations, showing increases in flux by a factor of $>80$ and $>300$ for X1 and X2, respectively. X1 showed indications of spectral variability: at the highest luminosity, its X-ray spectra can be fitted with a powerlaw ($\Gamma=1.55\pm0.15$), or a multicolour disc model with $T_{\mathrm{in}}=2.07^{+0.30}_{-0.23}$ keV; the X-ray spectrum became softer ($\Gamma=2.67^{+0.69}_{-0.64}$), or cooler ($T_\mathrm{in}=0.64^{+0.28}_{-0.17}$ keV) at lower luminosity. No strong evidence for spectral variability was found for X2. Its X-ray spectra can be fitted with a simple powerlaw model ($\Gamma=1.61^{+0.55}_{-0.50}$), or a multicolour disc model ($1.69^{+1.17}_{-0.48}$ keV). A possible optical counterpart for X1 is revealed in HST imaging. No optical variability is found, indicating that the optical radiation may be dominated by the companion star. Future X-ray and optical observations are necessary to determine the true nature of the compact object., Comment: Accepted for publication in MNRAS, 7 pages, 3 figures