Search

We present optical, radio and X-ray observations of a rapidly-evolving transient AT2019wxt (PS19hgw), discovered during the search for an electromagnetic (EM) counterpart to the gravitational-wave (GW) trigger S191213g (LIGO Scientific Collaboration & Virgo Collaboration 2019a). Although S191213g was not confirmed as a significant GW event in the off-line analysis of LIGO-Virgo data, AT2019wxt remained an interesting transient due its peculiar nature. The optical/NIR light curve of AT2019wxt displayed a double-peaked structure evolving rapidly in a manner analogous to currently know ultra-stripped supernovae (USSNe) candidates. This double-peaked structure suggests presence of an extended envelope around the progenitor, best modelled with two-components: i) early-time shock-cooling emission and ii) late-time radioactive $^{56}$Ni decay. We constrain the ejecta mass of AT2019wxt at $M_{ej} \approx{0.20 M_{\odot}}$ which indicates a significantly stripped progenitor that was possibly in a binary system. We also followed-up AT2019wxt with long-term Chandra and Jansky Very Large Array observations spanning $\sim$260 days. We detected no definitive counterparts at the location of AT2019wxt in these long-term X-ray and radio observational campaigns. We establish the X-ray upper limit at $9.93\times10^{-17}$ erg cm$^{-2}$ s$^{-1}$ and detect an excess radio emission from the region of AT2019wxt. However, there is little evidence for SN1993J- or GW170817-like variability of the radio flux over the course of our observations. A substantial host galaxy contribution to the measured radio flux is likely. The discovery and early-time peak capture of AT2019wxt in optical/NIR observation during EMGW follow-up observations highlights the need of dedicated early, multi-band photometric observations to identify USSNe., Comment: 20 pages, 14 figures, Submitted to ApJ

Several new radio facilities have a field of view and sensitivity well suited for transient searches. This makes it more important than ever to accurately determine transient rates in radio surveys. The work presented here seeks to do this task by using Monte-Carlo simulations. In particular, the user inputs either a real or simulated observational setup, and the simulations code calculates transient rate as a function of transient duration and peak flux. These simulations allow for simulating a wide variety of scenarios including observations with varying sensitivities and durations, multiple overlapping telescope pointings, and a wide variety of light curve shapes with the user having the ability to easily add more. While the current scientific focus is on the radio regime, with examples given here from the MeerKAT telescope in South Africa, the simulations code can be easily adapted to other wavelength regimes., Comment: Minor updates to plots

We present the full panchromatic afterglow light curve data of GW170817, including new radio data as well as archival optical and X-ray data, between 0.5 and 940 days post-merger. By compiling all archival data, and reprocessing a subset of it, we have evaluated the impact of differences in data processing or flux determination methods used by different groups, and attempted to mitigate these differences to provide a more uniform dataset. Simple power-law fits to the uniform afterglow light curve indicate a $t^{0.86\pm0.04}$ rise, a $t^{-1.92\pm0.12}$ decline, and a peak occurring at $155\pm4$ days. The afterglow is optically thin throughout its evolution, consistent with a single spectral index ($-0.584\pm0.002$) across all epochs. This gives a precise and updated estimate of the electron power-law index, $p=2.168\pm0.004$. By studying the diffuse X-ray emission from the host galaxy, we place a conservative upper limit on the hot ionized ISM density, $<$0.01 cm$^{-3}$, consistent with previous afterglow studies. Using the late-time afterglow data we rule out any long-lived neutron star remnant having magnetic field strength between 10$^{10.4}$ G and 10$^{16}$ G. Our fits to the afterglow data using an analytical model that includes VLBI proper motion from Mooley et al. (2018), and a structured jet model that ignores the proper motion, indicates that the proper motion measurement needs to be considered while seeking an accurate estimate of the viewing angle., Comment: Accepted by ApJ. Revised X-ray and up-to-date afterglow measurements are available in ASCII format on GitHub (https://github.com/kmooley/GW170817/) and mirrored on tauceti.caltech.edu (http://www.tauceti.caltech.edu/kunal/gw170817/). The online dataset will be continuously updated as new measurements become available

Several on-going or planned synoptic optical surveys are offering or will soon be offering an unprecedented opportunity for discovering larger samples of the rarest types of stripped-envelope core-collapse supernovae (SNe), such as those associated with relativistic jets, mildly-relativistic ejecta, or strong interaction with the circumstellar medium (CSM). Observations at radio wavelengths are a useful tool to probe the fastest moving ejecta, as well as denser circumstellar environments, and can thus help us identify the rarest type of core-collapse explosions. Here, we discuss how to set up an efficient radio follow-up program to detect and correctly identify radio-emitting stripped-envelope core-collapse explosions. We use a method similar to the one described in \citealt{Carbone2018}, and determine the optimal timing of GHz radio observations assuming a sensitivity comparable to that of the Karl G. Jansky Very Large Array. The optimization is done so as to ensure that the collected radio observations can identify the type of explosion powering the radio counterpart by using the smallest possible amount of telescope time. We also present a previously unpublished upper-limit on the late-time radio emission from supernova iPTF17cw. Finally, we conclude by discussing implications for follow-up in the X-rays., Comment: 9 pages, 2 figures, 7 tables, submitted to ApJ

We performed simulations of a large number of so-called very faint X-ray transient sources from surveys obtained using the X-ray telescope aboard the Neil Gehrels \emph{Swift} Observatory on two Galactic globular clusters, and the Galactic Center. We calculated the ratio between the duty cycle we input in our simulations and the one we measure after the simulations. We found that fluctuations in outburst duration and recurrence times affect our estimation of the duty cycle more than non detected outbursts. This biases our measures to overestimate the simulated duty cycle of sources. Moreover, we determined that compact surveys are necessary to detect outbursts with short duration because they could fall in gaps between observations, if such gaps are longer than their duration. On the other hand, long surveys are necessary to detect sources with low duty cycle because the smallest duty cycle a survey can observe is given by the ratio between the shortest outburst duration and the total length of the survey. If one has a limited amount of observing time, these two effects are competing, and a compromise is required which is set by the goals of the proposed survey. We have also performed simulations with several artificial survey strategies in order to evaluate the optimal observing campaign aimed at detecting transients as well as at having the most accurate estimates of the duty cycle. As expected, the best campaign would be a regular and dense monitoring that extends for a very long period. The closest real example of such a dataset is the monitoring of the Galactic Centre., Comment: 13 pages, 6 figures. Accepted for publication in MNRAS

GHz radio astronomy has played a fundamental role in the recent dazzling discovery of GW170817, a neutron star (NS)-NS merger observed in both gravitational waves (GWs) and light at all wavelengths. Here we show how the expected progress in sensitivity of ground-based GW detectors over the next decade calls for U.S.-based GHz radio arrays to be improved beyond current levels. We discuss specifically how several new scientific opportunities would emerge in multi-messenger time-domain astrophysics if a next generation GHz radio facility with sensitivity and resolution $10\times$ better than the current Jansky Very Large Array (VLA) were to work in tandem with ground-based GW detectors. These opportunities include probing the properties, structure, and size of relativistic jets and wide-angle ejecta from NS-NS mergers, as well as unraveling the physics of their progenitors via host galaxy studies., Comment: 8 pages, 5 figures, Submitted to the Astro2020 Decadal Survey

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network., Comment: White paper submitted to the Astro2020 Decadal Survey

In light of the recent dazzling discovery of GW170817, we discuss several new scientific opportunities that would emerge in multi-messenger time-domain astrophysics if a facility like the next generation Very Large Array (ngVLA) were to work in tandem with ground-based gravitational wave (GW) detectors. These opportunities include probing wide-angle ejecta and off-axis afterglows of neutron star (NS)-NS mergers; enabling direct size measurements of radio ejecta from NS-NS mergers; and unraveling the physics behind the progenitors of compact binary mergers via host galaxy studies at radio wavelengths. Our results show that, thanks to its unprecedented sensitivity and resolution, the ngVLA will enable transformational results in the multi-messenger exploration of the transient radio sky., Comment: To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA). arXiv admin note: substantial text overlap with arXiv:1711.06873

With the arrival of a number of wide-field snapshot image-plane radio transient surveys, there will be a huge influx of images in the coming years making it impossible to manually analyse the datasets. Automated pipelines to process the information stored in the images are being developed, such as the LOFAR Transients Pipeline, outputting light curves and various transient parameters. These pipelines have a number of tuneable parameters that require training to meet the survey requirements. This paper utilises both observed and simulated datasets to demonstrate different machine learning strategies that can be used to train these parameters. The datasets used are from LOFAR observations and we process the data using the LOFAR Transients Pipeline; however the strategies developed are applicable to any light curve datasets at different frequencies and can be adapted to different automated pipelines. These machine learning strategies are publicly available as Python tools that can be downloaded and adapted to different datasets (https://github.com/AntoniaR/TraP_ML_tools)., Comment: Astronomy & Computing Accepted, 25 pages, 20 figures

We present late-time radio observations of GW170817, the first binary neutron star merger discovered through gravitational waves by the advanced LIGO and Virgo detectors. Our observations, carried out with the Karl G. Jansky Very Large Array, were optimized to detect polarized radio emission, and thus to constrain the linear polarization fraction of GW170817. At an epoch of ~244 days after the merger, we rule out linearly polarized emission above a fraction of ~12% at a frequency of 2.8 GHz (99% confidence). Within the structured jet scenario (a.k.a. successful jet plus cocoon system) for GW170817, the derived upper-limit on the radio continuum linear polarization fraction strongly constrains the magnetic field configuration in the shocked ejecta. We show that our results for GW170817 are compatible with the low level of linear polarization found in afterglows of cosmological long gamma-ray bursts. Finally, we discuss our findings in the context of future expectations for the study of radio counterparts of binary neutron star mergers identified by ground-based gravitational-wave detectors., Comment: 5 pages, 2 figures, 1 table

Motivated by the recent discovery of the binary neutron-star (BNS) merger GW170817, we determine the optimal observational setup for detecting and characterizing radio counterparts of nearby ($d_L\sim40$\,Mpc) BNS mergers. We simulate GW170817-like radio transients, and radio afterglows generated by fast jets with isotropic energy $E_{\rm iso}\sim 10^{50}$\,erg, expanding in a low-density interstellar medium (ISM; $n_{\rm ISM}=10^{-4}-10^{-2}$\,cm$^{-3}$), observed from different viewing angles (from slightly off-axis to largely off-axis). We then determine the optimal timing of GHz radio observations following the precise localization of the BNS radio counterpart candidate, assuming a sensitivity comparable to that of the Karl G. Jansky Very Large Array. The optimization is done so as to ensure that properties such as viewing angle and circumstellar density can be correctly reconstructed with the minimum number of observations. We show that radio is the optimal band to explore the fastest ejecta from BNSs in low-density ISM, since the optical emission is likelyto be dominated by the so-called `kilonova' component, while X-rays from the jet are detectable only for a small subset of the BNS models considered here. Finally, we discuss how future radio arrays like the next generation VLA (ngVLA) would improve the detectability of BNS mergers with physical parameters similar to the ones here explored., Comment: Accepted for Publication in ApJ; 8 pages, 2 Figures, 8 Tables

We have developed a method to determine the transient surface density and transient rate for any given survey, using Monte-Carlo simulations. This method allows us to determine the transient rate as a function of both the flux and the duration of the transients in the whole flux-duration plane rather than one or a few points as currently available methods do. It is applicable to every survey strategy that is monitoring the same part of the sky, regardless the instrument or wavelength of the survey, or the target sources. We have simulated both top-hat and Fast Rise Exponential Decay light curves, highlighting how the shape of the light curve might affect the detectability of transients. Another application for this method is to estimate the number of transients of a given kind that are expected to be detected by a survey, provided that their rate is known., Comment: 12 pages, 6 figures, accepted for publication in MNRAS

Current and future astronomical survey facilities provide a remarkably rich opportunity for transient astronomy, combining unprecedented fields of view with high sensitivity and the ability to access previously unexplored wavelength regimes. This is particularly true of LOFAR, a recently-commissioned, low-frequency radio interferometer, based in the Netherlands and with stations across Europe. The identification of and response to transients is one of LOFAR's key science goals. However, the large data volumes which LOFAR produces, combined with the scientific requirement for rapid response, make automation essential. To support this, we have developed the LOFAR Transients Pipeline, or TraP. The TraP ingests multi-frequency image data from LOFAR or other instruments and searches it for transients and variables, providing automatic alerts of significant detections and populating a lightcurve database for further analysis by astronomers. Here, we discuss the scientific goals of the TraP and how it has been designed to meet them. We describe its implementation, including both the algorithms adopted to maximize performance as well as the development methodology used to ensure it is robust and reliable, particularly in the presence of artefacts typical of radio astronomy imaging. Finally, we report on a series of tests of the pipeline carried out using simulated LOFAR observations with a known population of transients., Comment: 30 pages, 11 figures; Accepted for publication in Astronomy & Computing; Code at https://github.com/transientskp/tkp