Showing papers in "arXiv: High Energy Astrophysical Phenomena in 2018"

Binary Black Hole Population Properties Inferred from the First and Second Observing Runs of Advanced LIGO and Advanced Virgo

Abstract: We present results on the mass, spin, and redshift distributions with phenomenological population models using the ten binary black hole mergers detected in the first and second observing runs completed by Advanced LIGO and Advanced Virgo. We constrain properties of the binary black hole (BBH) mass spectrum using models with a range of parameterizations of the BBH mass and spin distributions. We find that the mass distribution of the more massive black hole in such binaries is well approximated by models with no more than 1% of black holes more massive than $45\,M_\odot$, and a power law index of $\alpha = {1.3}^{+1.4}_{-1.7}$ (90% credibility). We also show that BBHs are unlikely to be composed of black holes with large spins aligned to the orbital angular momentum. Modelling the evolution of the BBH merger rate with redshift, we show that it is at or increasing with redshift with 93% probability. Marginalizing over uncertainties in the BBH population, we find robust estimates of the BBH merger rate density of $R = {53.2}^{+55.8}_{-28.2}$ Gpc$^{-3}$ yr$^{-1}$ (90% credibility). As the BBH catalog grows in future observing runs, we expect that uncertainties in the population model parameters will shrink, potentially providing insights into the formation of black holes via supernovae, binary interactions of massive stars, stellar cluster dynamics, and the formation history of black holes across cosmic time.

Multi-messenger observations of a flaring blazar coincident with high-energy neutrino IceCube-170922A

Abstract: Individual astrophysical sources previously detected in neutrinos are limited to the Sun and the supernova 1987A, whereas the origins of the diffuse flux of high-energy cosmic neutrinos remain unidentified. On 22 September 2017 we detected a high-energy neutrino, IceCube-170922A, with an energy of approximately 290 TeV. Its arrival direction was consistent with the location of a known gamma-ray blazar TXS 0506+056, observed to be in a flaring state. An extensive multi-wavelength campaign followed, ranging from radio frequencies to gamma-rays. These observations characterize the variability and energetics of the blazar and include the first detection of TXS 0506+056 in very-high-energy gamma-rays. This observation of a neutrino in spatial coincidence with a gamma-ray emitting blazar during an active phase suggests that blazars may be a source of high-energy neutrinos.

GWTC-1: A Gravitational-Wave Transient Catalog of Compact Binary Mergers Observed by LIGO and Virgo during the First and Second Observing Runs

Abstract: We present the results from three gravitational-wave searches for coalescing compact binaries with component masses above 1$\mathrm{M}_\odot$ during the first and second observing runs of the Advanced gravitational-wave detector network. During the first observing run (O1), from September $12^\mathrm{th}$, 2015 to January $19^\mathrm{th}$, 2016, gravitational waves from three binary black hole mergers were detected. The second observing run (O2), which ran from November $30^\mathrm{th}$, 2016 to August $25^\mathrm{th}$, 2017, saw the first detection of gravitational waves from a binary neutron star inspiral, in addition to the observation of gravitational waves from a total of seven binary black hole mergers, four of which we report here for the first time: GW170729, GW170809, GW170818 and GW170823. For all significant gravitational-wave events, we provide estimates of the source properties. The detected binary black holes have total masses between $18.6_{-0.7}^{+3.2}\mathrm{M}_\odot$, and $84.4_{-11.1}^{+15.8} \mathrm{M}_\odot$, and range in distance between $320_{-110}^{+120}$ Mpc and $2840_{-1360}^{+1400}$ Mpc. No neutron star - black hole mergers were detected. In addition to highly significant gravitational-wave events, we also provide a list of marginal event candidates with an estimated false alarm rate less than 1 per 30 days. From these results over the first two observing runs, which include approximately one gravitational-wave detection per 15 days of data searched, we infer merger rates at the 90% confidence intervals of $110\, -\, 3840$ $\mathrm{Gpc}^{-3}\,\mathrm{y}^{-1}$ for binary neutron stars and $9.7\, -\, 101$ $\mathrm{Gpc}^{-3}\,\mathrm{y}^{-1}$ for binary black holes assuming fixed population distributions, and determine a neutron star - black hole merger rate 90% upper limit of $610$ $\mathrm{Gpc}^{-3}\,\mathrm{y}^{-1}$.

The binary neutron star event ligo/virgo gw170817 a hundred and sixty days after merger: Synchrotron emission across the electromagnetic spectrum

Abstract: We report deep Chandra, HST and VLA observations of the binary neutron star event GW170817 at $t<160$ d after merger. These observations show that GW170817 has been steadily brightening with time and might have now reached its peak, and constrain the emission process as non-thermal synchrotron emission where the cooling frequency $ u_c$ is above the X-ray band and the synchrotron frequency $ u_m$ is below the radio band. The very simple power-law spectrum extending for eight orders of magnitude in frequency enables the most precise measurement of the index $p$ of the distribution of non-thermal relativistic electrons $N(\gamma)\propto \gamma^{-p}$ accelerated by a shock launched by a NS-NS merger to date. We find $p=2.17\pm0.01$, which indicates that radiation from ejecta with $\Gamma\sim3-10$ dominates the observed emission. While constraining the nature of the emission process, these observations do \emph{not} constrain the nature of the relativistic ejecta. We employ simulations of explosive outflows launched in NS ejecta clouds to show that the spectral and temporal evolution of the non-thermal emission from GW170817 is consistent with both emission from radially stratified quasi-spherical ejecta traveling at mildly relativistic speeds, \emph{and} emission from off-axis collimated ejecta characterized by a narrow cone of ultra-relativistic material with slower wings extending to larger angles. In the latter scenario, GW170817 harbored a normal SGRB directed away from our line of sight. Observations at $t\le 200$ days are unlikely to settle the debate as in both scenarios the observed emission is effectively dominated by radiation from mildly relativistic material.

The NANOGrav 11-year Data Set: Pulsar-timing Constraints On The Stochastic Gravitational-wave Background

Abstract: We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released $11$-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no significant evidence for a GWB, we place constraints on a GWB from a population of supermassive black-hole binaries, cosmic strings, and a primordial GWB. For the first time, we find that the GWB upper limits and detection statistics are sensitive to the Solar System ephemeris (SSE) model used, and that SSE errors can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first PTA constraints that are robust against SSE uncertainties. We thus place a $95\%$ upper limit on the GW strain amplitude of $A_\mathrm{GWB}<1.45\times 10^{-15}$ at a frequency of $f=1$ yr$^{-1}$ for a fiducial $f^{-2/3}$ power-law spectrum, and with inter-pulsar correlations modeled. This is a factor of $\sim 2$ improvement over the NANOGrav $9$-year limit, calculated using the same procedure. Previous PTA upper limits on the GWB will need revision in light of SSE systematic uncertainties. We use our constraints to characterize the combined influence on the GWB of the stellar mass-density in galactic cores, the eccentricity of SMBH binaries, and SMBH--galactic-bulge scaling relationships. We constrain cosmic-string tension using recent simulations, yielding an SSE-marginalized $95\%$ upper limit on the cosmic string tension of $G\mu < 5.3\times 10^{-11}$---a factor of $\sim 2$ better than the published NANOGrav $9$-year constraints. Our SSE-marginalized $95\%$ upper limit on the energy density of a primordial GWB (for a radiation-dominated post-inflation Universe) is $\Omega_\mathrm{GWB}(f)h^2<3.4\times10^{-10}$.

The construction and use of LISA sensitivity curves

Abstract: The Laser Interferometer Space Antenna (LISA) will open the mHz band of the gravitational wave spectrum for exploration. Sensitivity curves are a useful tool for surveying the types of sources that can be detected by the LISA mission. Here we describe how the sensitivity curve is constructed, and how it can be used to compute the signal-to-noise ratio for a wide range of binary systems. We adopt the 2018 LISA Phase-0 reference design parameters. We consider both sky-averaged sensitivities, and the sensitivity to sources at particular sky locations. The calculations are included in a publicly available {\em Python} notebook.

Relativistic Jets in Active Galactic Nuclei

Abstract: The nuclei of most normal galaxies contain supermassive black holes, which can accrete gas through a disk and become active. These Active Galactic Nuclei, AGN, can form jets which are observed on scales from AU to Mpc and from meter wavelengths to TeV gamma ray energies. High resolution radio imaging and multi-wavelength/messenger campaigns are elucidating the conditions under which this happens. Evidence is presented that: AGN jets are formed when the black hole spins and the accretion disk is strongly magnetized, perhaps on account of gas accretingat high latitude beyond the black hole sphere of influence; AGN jets are collimated close to the black hole by magnetic stress associated with a disk wind; higher power jets can emerge from their galactic nuclei in a relativistic, supersonic and proton-dominated state and they terminate in strong, hot spot shocks; lower power jets are degraded to buoyant plumes and bubbles; jets may accelerate protons to EeV energies which contribute to the cosmic ray spectrum and which may initiate pair cascades that can efficiently radiate synchrotron gamma rays; jets were far more common when the universe was a few billion years old and black holes and massive galaxies were growing rapidly; jets can have a major influence on their environments, stimulating and limiting the growth of galaxies. The observational prospects for securing our understanding of AGN jets are bright.

Highest-frequency detection of FRB 121102 at 4-8 GHz using the Breakthrough Listen Digital Backend at the Green Bank Telescope

Abstract: We report the first detections of the repeating fast radio burst source FRB 121102 above 5.2 GHz. Observations were performed using the 4$-$8 GHz receiver of the Robert C. Byrd Green Bank Telescope with the Breakthrough Listen digital backend. We present the spectral, temporal and polarization properties of 21 bursts detected within the first 60 minutes of a total 6-hour observations. These observations comprise the highest burst density yet reported in the literature, with 18 bursts being detected in the first 30 minutes. A few bursts clearly show temporal sub-structures with distinct spectral properties. These sub-structures superimpose to provide enhanced peak signal-to-noise ratio at higher trial dispersion measures. Broad features occur in $\sim 1$ GHz wide subbands that typically differ in peak frequency between bursts within the band. Finer-scale structures ($\sim 10-50$ MHz) within these bursts are consistent with that expected from Galactic diffractive interstellar scintillation. The bursts exhibit nearly 100% linear polarization, and a large average rotation measure of 9.359$\pm$0.012 $\times$ 10$^{\rm 4}$ rad m$^{\rm -2}$ (in the observer's frame). No circular polarization was found for any burst. We measure an approximately constant polarization position angle in the 13 brightest bursts. The peak flux densities of the reported bursts have average values (0.2$\pm$0.1 Jy), similar to those seen at lower frequencies ($<3$ GHz), while the average burst widths (0.64$\pm$0.46 ms) are relatively narrower.

The Structure and Dynamics of the Sub-parsec Scale Jet in M87 Based on 50 VLBA Observations Over 17 Years at 43 GHz

Abstract: The central radio source in M87 provides the best opportunity to study jet formation because it has a large angular size for the gravitational radius of the black hole and has a bright jet that is well resolved by VLBI observations. We present intensive monitoring observations from 2007 and 2008, plus roughly annual observations that span 17 years, all made with the the Very Long Baseline Array at 43 GHz with a resolution of about 30 by 60 Rs. Our high-dynamic-range images clearly show the wide-opening-angle structure and the counter-jet. The jet and counter-jet are nearly symmetric in the inner 1.5 milli-arcseconds (mas; 0.12 pc in projection) with both being edge brightened. Both show deviations from parabolic shape in the form of an initial rapid expansion and subsequent contraction followed by further rapid expansion and, beyond the visible counter-jet, subsequent collimation. Proper motions and counter-jet/jet intensity ratios both indicate acceleration from apparent speeds of $\lesssim 0.5c$ to $\gtrsim 2c$ in the inner about 2 mas and suggest a helical flow. The jet displays a sideways shift with an approximately 8 to 10 year quasi-periodicity. The shift propagates outwards non-ballistically and significantly more slowly than the flow speed revealed by the fastest moving components. Polarization data show a systematic structure with magnetic field vectors that suggest a toroidal field close to the core.

The Giant Radio Array for Neutrino Detection (GRAND): Science and Design

Abstract: The Giant Radio Array for Neutrino Detection (GRAND) is a planned large-scale observatory of ultra-high-energy (UHE) cosmic particles, with energies exceeding 10^8 GeV. Its goal is to solve the long-standing mystery of the origin of UHE cosmic rays. To do this, GRAND will detect an unprecedented number of UHE cosmic rays and search for the undiscovered UHE neutrinos and gamma rays associated to them with unmatched sensitivity. GRAND will use large arrays of antennas to detect the radio emission coming from extensive air showers initiated by UHE particles in the atmosphere. Its design is modular: 20 separate, independent sub-arrays, each of 10 000 radio antennas deployed over 10 000 km^2. A staged construction plan will validate key detection techniques while achieving important science goals early. Here we present the science goals, detection strategy, preliminary design, performance goals, and construction plans for GRAND.

Discovery of a Candidate Black Hole - Giant Star Binary System in the Galactic Field

Abstract: We report the discovery of the first likely black hole in a non-interacting binary system with a field red giant. By combining radial velocity measurements from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) with photometric variability data from the All-Sky Automated Survey for Supernovae (ASAS-SN), we identified the bright rapidly-rotating giant 2MASS J05215658+4359220 as a binary system with a massive unseen companion. Subsequent radial velocity measurements reveal a system with an orbital period of 83 days and near-zero eccentricity. The photometric variability period of the giant is consistent with the orbital period, indicative of star spots and tidal synchronization. Constraints on the giant's mass and radius from its luminosity, surface gravity, and temperature imply an unseen companion with mass of $3.3^{+2.8}_{-0.7}$ M$_\odot$, indicating a low-mass black hole or an exceedingly massive neutron star. Measurement of the astrometric binary motion by {\it Gaia} will further characterize the system. This discovery demonstrates the potential of massive spectroscopic surveys like APOGEE and all-sky, high-cadence photometric surveys like ASAS-SN to revolutionize our understanding of the compact object mass function, and to test theories of binary star evolution and the supernova mechanism.

An embedded X-ray source shines through the aspherical AT2018cow: revealing the inner workings of the most luminous fast-evolving optical transients

Abstract: We present the first extensive radio to gamma-ray observations of a fast-rising blue optical transient (FBOT), AT2018cow, over its first ~100 days. AT2018cow rose over a few days to a peak luminosity $L_{pk}\sim4\times 10^{44}$ erg/s exceeding those of superluminous supernovae (SNe), before declining as $\propto t^{-2}$. Initial spectra at $\lesssim 15$ days were mostly featureless and indicated large expansion velocities v~0.1c and temperatures reaching 30000 K. Later spectra revealed a persistent optically-thick photosphere and the emergence of H and He emission features with v~sim 4000 km/s with no evidence for ejecta cooling. Our broad-band monitoring revealed a hard X-ray spectral component at $E\ge 10$ keV, in addition to luminous and highly variable soft X-rays, with properties unprecedented among astronomical transients. An abrupt change in the X-ray decay rate and variability appears to accompany the change in optical spectral properties. AT2018cow showed bright radio emission consistent with the interaction of a blastwave with $v_{sh}$~0.1c with a dense environment ($\dot M\sim10^{-3}-10^{-4}\,M_{\odot}yr^{-1}$ for $v_w=1000$ km\s). While these properties exclude Ni-powered transients, our multi-wavelength analysis instead indicates that AT2018cow harbored a "central engine", either a compact object (magnetar or black hole) or an embedded internal shock produced by interaction with a compact, dense circumstellar medium. The engine released $\sim10^{50}-10^{51.5}$ erg over $\sim10^3-10^5$ s and resides within low-mass fast-moving material with equatorial-polar density asymmetry ($M_{ej,fast}\lesssim0.3\,\rm{M_{\odot}}$). Successful SNe from low-mass H-rich stars (like electron-capture SNe) or failed explosions from blue supergiants satisfy these constraints. Intermediate-mass black-holes are disfavored by the large environmental density probed by the radio observations.

Superluminous supernovae

Abstract: Superluminous supernovae are a new class of supernovae that were recognized about a decade ago. Both observational and theoretical progress has been significant in the last decade. In this review, we first briefly summarize the observational properties of superluminous supernovae. We then introduce the three major suggested luminosity sources to explain the huge luminosities of superluminous supernovae, i.e., the nuclear decay of 56Ni, the interaction between supernova ejecta and dense circumstellar media, and the spin down of magnetars. We compare these models and discuss their strengths and weaknesses.

Peering into the dark side: magnesium lines establish a massive neutron star in PSR J2215+5135

Abstract: New millisecond pulsars (MSPs) in compact binaries provide a good opportunity to search for the most massive neutron stars. Their main-sequence companion stars are often strongly irradiated by the pulsar, displacing the effective center of light from their barycenter and making mass measurements uncertain. We present a series of optical spectroscopic and photometric observations of PSR J2215+5135, a "redback" binary MSP in a 4.14 hr orbit, and measure a drastic temperature contrast between the dark/cold (T$_\mathrm{N}$=5660$^{+260}_{-380}$ K) and bright/hot (T$_\mathrm{D}$=8080$^{+470}_{-280}$ K) sides of the companion star. We find that the radial velocities depend systematically on the atmospheric absorption lines used to measure them. Namely, the semi-amplitude of the radial velocity curve of J2215 measured with magnesium triplet lines is systematically higher than that measured with hydrogen Balmer lines, by 10%. We interpret this as a consequence of strong irradiation, whereby metallic lines dominate the dark side of the companion (which moves faster) and Balmer lines trace its bright (slower) side. Further, using a physical model of an irradiated star to fit simultaneously the two-species radial velocity curves and the three-band light curves, we find a center-of-mass velocity of K$_2$=412.3$\pm$5.0 km s$^{-1}$ and an orbital inclination i=63.9$^\circ$$^{+2.4}_{-2.7}$. Our model is able to reproduce the observed fluxes and velocities without invoking irradiation by an extended source. We measure masses of M$_1$=2.27$^{+0.17}_{-0.15}$ M$_\odot$ and M$_2$=0.33$^{+0.03}_{-0.02}$ M$_\odot$ for the neutron star and the companion star, respectively. If confirmed, such a massive pulsar would rule out some of the proposed equations of state for the neutron star interior.

Indication of anisotropy in arrival directions of ultra-high-energy cosmic rays through comparison to the flux pattern of extragalactic gamma-ray sources

Abstract: A new analysis of the dataset from the Pierre Auger Observatory provides evidence for anisotropy in the arrival directions of ultra-high-energy cosmic rays on an intermediate angular scale, which is indicative of excess arrivals from strong, nearby sources. The data consist of 5514 events above 20 EeV with zenith angles up to 80 deg recorded before 2017 April 30. Sky models have been created for two distinct populations of extragalactic gamma-ray emitters: active galactic nuclei from the second catalog of hard Fermi-LAT sources (2FHL) and starburst galaxies from a sample that was examined with Fermi-LAT. Flux-limited samples, which include all types of galaxies from the Swift-BAT and 2MASS surveys, have been investigated for comparison. The sky model of cosmic-ray density constructed using each catalog has two free parameters, the fraction of events correlating with astrophysical objects and an angular scale characterizing the clustering of cosmic rays around extragalactic sources. A maximum-likelihood ratio test is used to evaluate the best values of these parameters and to quantify the strength of each model by contrast with isotropy. It is found that the starburst model fits the data better than the hypothesis of isotropy with a statistical significance of 4.0 sigma, the highest value of the test statistic being for energies above 39 EeV. The three alternative models are favored against isotropy with 2.7-3.2 sigma significance. The origin of the indicated deviation from isotropy is examined and prospects for more sensitive future studies are discussed.

Search for Spatial Extension in High-Latitude Sources Detected by the Fermi Large Area Telescope

Abstract: We present a search for spatial extension in high-latitude ($|b|>5^\circ$) sources in recent Fermi point source catalogs. The result is the Fermi High-Latitude Extended Sources Catalog, which provides source extensions (or upper limits thereof) and likelihood profiles for a suite of tested source morphologies. We find 24 extended sources, 19 of which were not previously characterized as extended. These include sources that are potentially associated with supernova remnants and star forming regions. We also found extended $\gamma$-ray emission in the vicinity of the Cen A radio lobes and - at GeV energies for the first time - spatially coincident with the radio emission of the SNR CTA 1, as well as from the Crab Nebula. We also searched for halos around active galactic nuclei, which are predicted from electromagnetic cascades that are induced by the $e^+e^-$ pairs that are deflected in intergalactic magnetic fields. These are produced when $\gamma$-rays interact with background radiation fields. We do not find evidence for extension in individual sources or in stacked source samples. This enables us to place limits on the flux of the extended source components, which are then used to constrain the intergalactic magnetic field a coherence length $\lambda \gtrsim 10\,$kpc, even when conservative assumptions on the source duty cycle are made. This improves previous limits by several orders of magnitude.

A magnetar origin for the kilonova ejecta in GW170817

Abstract: The neutron star (NS) merger GW170817 was followed over several days by optical-wavelength ("blue") kilonova (KN) emission likely powered by the radioactive decay of light r-process nuclei synthesized by ejecta with a low neutron abundance (electron fraction Ye ~ 0.25-0.35). While the composition and high velocities of the blue KN ejecta are consistent with shock-heated dynamical material, the large quantity is in tension with the results of numerical simulations. We propose an alternative ejecta source: the neutrino-heated, magnetically-accelerated wind from the strongly-magnetized hypermassive NS (HMNS) remnant. A rapidly-spinning HMNS with an ordered surface magnetic field of strength B ~ 1-3e14 G and lifetime t_rem ~ 0.1-1 s can simultaneously explain the velocity, total mass, and electron fraction of the blue KN ejecta. The inferred HMNS lifetime is close to its Alfven crossing time, suggesting global magnetic torques could be responsible for bringing the HMNS into solid body rotation and instigating its gravitational collapse. Different origins for the KN ejecta may be distinguished by their predictions for the emission in the first hours after the merger, when the luminosity is enhanced by heating from internal shocks; the latter are likely generic to any temporally-extended ejecta source (e.g. magnetar or accretion disk wind) and are not unique to the emergence of a relativistic jet. The same shocks could mix and homogenizes the composition to a low but non-zero lanthanide mass fraction, X_La ~ 1e-3, as advocated by some authors, but only if the mixing occurs after neutrons are consumed in the r-process on a timescale >~ 1 s.

A Decline in the X-ray through Radio Emission from GW170817 Continues to Support an Off-Axis Structured Jet.

Abstract: We present new observations of the binary neutron star merger GW170817 at $\Delta t\approx 220-290$ days post-merger, at radio (Karl G. Jansky Very Large Array; VLA), X-ray (Chandra X-ray Observatory) and optical (Hubble Space Telescope; HST) wavelengths. These observations provide the first evidence for a turnover in the X-ray light curve, mirroring a decline in the radio emission at $\gtrsim5\sigma$ significance. The radio-to-X-ray spectral energy distribution exhibits no evolution into the declining phase. Our full multi-wavelength dataset is consistent with the predicted behavior of our previously published models of a successful structured jet expanding into a low-density circumbinary medium, but pure cocoon models with a choked jet cannot be ruled out. If future observations continue to track our predictions, we expect that the radio and X-ray emission will remain detectable until $\sim 1000$ days post-merger.

The Cow: discovery of a luminous, hot and rapidly evolving transient

Abstract: We present the ATLAS discovery and initial analysis of the first 18 days of the unusual transient event, ATLAS18qqn/AT2018cow. It is characterized by a high peak luminosity ($\sim$1.7 $\times$ 10$^{44}$ erg s$^{-1}$), rapidly evolving light curves ($>$5 mag rise to peak in $\sim$3.5 days), and hot blackbody spectra, peaking at $\sim$27000 K that are relatively featureless and unchanging over the first two weeks. The bolometric light curve cannot be powered by radioactive decay under realistic assumptions. The detection of high-energy emission may suggest a central engine as the powering source. Using a magnetar model, we estimated an ejected mass of $0.1-0.4$ \msol, which lies between that of low-energy core-collapse events and the kilonova, AT2017gfo. The spectra cooled rapidly from 27000 to 15000 K in just over 2 weeks but remained smooth and featureless. Broad and shallow emission lines appear after about 20 days, and we tentatively identify them as He I although they would be redshifted from their rest wavelengths. We rule out that there are any features in the spectra due to intermediate mass elements up to and including the Fe-group. The presence of r-process elements cannot be ruled out. If these lines are due to He, then we suggest a low-mass star with residual He as a potential progenitor. Alternatively, models of magnetars formed in neutron-star mergers give plausible matches to the data.

A Strong Jet Signature in the Late-Time Lightcurve of GW170817

Abstract: We present new 0.6-10 GHz observations of the binary neutron star merger GW170817 covering the period up to 300 days post-merger, taken with the Karl G. Jansky Very Large Array, the Australia Telescope Compact Array, the Giant Metrewave Radio Telescope and the MeerKAT telescope. We use these data to precisely characterize the decay phase of the late-time radio light curve. We find that the temporal decay is consistent with a power-law slope of t^-2.2, and that the transition between the power-law rise and decay is relatively sharp. Such a slope cannot be produced by a quasi-isotropic (cocoon-dominated) outflow, but is instead the classic signature of a relativistic jet. This provides strong observational evidence that GW170817 produced a successful jet, and directly demonstrates the link between binary neutron star mergers and short-hard GRBs. Using simple analytical arguments, we derive constraints on the geometry and the jet opening angle of GW170817. These results are consistent with those from our companion Very Long Baseline Interferometry (VLBI) paper, reporting superluminal motion in GW170817.

A Living Theory Catalogue for Fast Radio Bursts

Abstract: At present, we have almost as many theories to explain Fast Radio Bursts as we have Fast Radio Bursts observed. This landscape will be changing rapidly with CHIME/FRB, recently commissioned in Canada, and HIRAX, under construction in South Africa. This is an opportune time to review existing theories and their observational consequences, allowing us to efficiently curtail viable astrophysical models as more data becomes available. In this article we provide a currently up to date catalogue of the numerous and varied theories proposed for Fast Radio Bursts so far. We also launch an online evolving repository for the use and benefit of the community to dynamically update our theoretical knowledge and discuss constraints and uses of Fast Radio Bursts.

The evolution of the X-ray afterglow emission of GW 170817 / GRB 170817A in XMM-Newton observations

Abstract: We report our observation of the short GRB 170817A, associated to the binary neutron star merger event GW 170817, perfomed in the X-ray band with XMM-Newton 135 d after the event (on the 29th December 2017). We find evidence for a flattening of the X-ray light curve with respect to the previously observed brightening. This is also supported by a nearly simultaneous optical Hubble Space Telescope and successive X-ray Chandra and low-frequency radio observations recently reported in the literature. Since the optical-to-X-ray spectral slope did not change with respect to previous observations, we exclude that the change in the temporal evolution of the light curve is due to the passage of the cooling frequency: its origin must be geometric or dynamical. We interpret all the existing afterglow data with two models: i) a structured jet and ii) a jet-less isotropic fireball with some stratification in its radial velocity structure. Both models fit the data and predict that the radio flux must decrease simultaneously with the optical and the X-ray one, making hard to distinguish between them at the present stage. Polarimetric measures and the rate of short GRB-GW association in future LIGO/Virgo runs will be key to disentangle these two geometrically different scenarios.

Numerical simulations of the jet dynamics and synchrotron radiation of binary neutron star merger event GW170817/GRB170817A

Abstract: We present numerical simulations of energetic flows propagating through the debris cloud of a binary neutron star (BNS) merger. Starting from the scale of the central engine, we use a moving-mesh hydrodynamics code to simulate the complete dynamical evolution of the produced relativistic jets. We compute synchrotron emission directly from the simulations and present multi-band light curves of the early (sub-day) through late (weeks to years) afterglow stages. Our work systematically compares two distinct models for the central engine, referred to as the narrow and wide engine scenario, which is associated with a successful structured jet and a quasi-isotropic explosion respectively. Both engine models naturally evolve angular and radial structure through hydrodynamical interaction with the merger debris cloud. They both also result in a relativistic blast wave capable of producing the observed multi-band afterglow data. However, we find that the narrow and wide engine scenario might be differentiated by a new emission component that we refer to as a merger flash. This component is a consequence of applying the synchrotron radiation model to the shocked optically thin merger cloud. Such modeling is appropriate if injection of non-thermal electrons is sustained in the breakout relativistic shell, for example by internal shocks or magnetic reconnection. The rapidly declining signature may be detectable for future BNS mergers during the first minutes to day following the GW chirp. Furthermore, its non-detection for the GRB170817A event may disfavor the wide, quasi-isotropic explosion model.

Testing the universality of free fall by tracking a pulsar in a stellar triple system

Abstract: Einstein's theory of gravity, general relativity, has passed stringent tests in laboratories, elsewhere in the Solar Sytem, and in pulsar binaries. Nevertheless it is known to be incompatible with quantum mechanics and must differ from the true behaviour of matter in strong fields and at small spatial scales. A key aspect of general relativity to test is the strong equivalence principle (SEP), which states that all freely falling objects, regardless of how strong their gravity, experience the same acceleration in the same gravitational field. Essentially all alternatives to general relativity violate this principle at some level. Previous direct tests of the SEP are limited by the weak gravity of the bodies in the Earth-Moon-Sun system or by the weak gravitational pull of the Galaxy on pulsar-white dwarf binaries. PSR~J0337+1715 is a hierarchical stellar triple system, where the inner binary consists of a millisecond radio pulsar in a $1.6$-day orbit with a white dwarf. This inner binary is in a $327$-day orbit with another white dwarf. In this system, the pulsar and the inner companion fall toward the outer companion with an acceleration about $10^8$ times greater than that produced by falling in the Galactic potential, and the pulsar's gravitational binding energy is roughly $10\%$ of its mass. Here we report that in spite of the pulsar's strong gravity, the accelerations experienced by it and the inner white dwarf differ by a fraction of no more than $2.6\times 10^{-6}$ ($95\%$ confidence level). We can roughly compare this to other SEP tests by using the strong-field Nordtvedt parameter $\hat\eta_N$. Our limit on $\hat\eta_N$ is a factor of ten smaller than that obtained from (weak-field) Solar-System SEP tests and a factor of almost a thousand smaller than that obtained from other strong-field SEP tests.

Fast Radio Burst 121102 Pulse Detection and Periodicity: A Machine Learning Approach

Abstract: We report the detection of 72 new pulses from the repeating fast radio burst FRB 121102 in Breakthrough Listen C-band (4-8 GHz) observations at the Green Bank Telescope. The new pulses were found with a convolutional neural network in data taken on August 26, 2017, where 21 bursts have been previously detected. Our technique combines neural network detection with dedispersion verification. For the current application we demonstrate its advantage over a traditional brute-force dedis- persion algorithm in terms of higher sensitivity, lower false positive rates, and faster computational speed. Together with the 21 previously reported pulses, this observa- tion marks the highest number of FRB 121102 pulses from a single observation, total- ing 93 pulses in five hours, including 45 pulses within the first 30 minutes. The number of data points reveal trends in pulse fluence, pulse detection rate, and pulse frequency structure. We introduce a new periodicity search technique, based on the Rayleigh test, to analyze the time of arrivals, with which we exclude with 99% confidence pe- riodicity in time of arrivals with periods larger than 5.1 times the model-dependent time-stamp uncertainty. In particular, we rule out constant periods >10 ms in the barycentric arrival times, though intrinsic periodicity in the time of emission remains plausible.

Hydrodynamics of circumbinary accretion: Angular momentum transfer and binary orbital evolution

Abstract: We carry out 2D viscous hydrodynamical simulations of circumbinary accretion using the AREPO code. We self-consistently compute the accretion flow over a wide range of spatial scales, from the circumbinary disk (CBD) far from the central binary, through accretion streamers, to the disks around individual binary components, resolving the flow down to 2% of the binary separation. We focus on equal-mass binaries with arbitrary eccentricities. We evolve the flow over long (viscous) timescales until a quasi-steady is reached, in which the mass supply rate at large distances $\dot{M}_0$ (assumed constant) equals the time-averaged mass transfer rate across the disk and the total mass accretion rate onto the binary components. This quasi-steady state allows us to compute the secular angular momentum transfer rate onto the binary, $\langle\dot{J}_b\rangle$, and the resulting orbital evolution. Through direct computation of the gravitational and accretion torques on the binary, we find that $\langle\dot{J}_b\rangle$ is consistently positive (i.e., the binary gains angular momentum), with $l_0\equiv\langle\dot{J}_b\rangle/\dot M_0$ in the range of $(0.4-0.8)a_b^2\Omega_b$, depending on the binary eccentricity (where $a_b,~\Omega_b$ are the binary semi-major axis and angular frequency); we also find that this $\langle\dot{J}_b\rangle$ is equal to the net angular momentum current across the CBD, indicating that global angular momentum balance is achieved in our simulations. We compute the time-averaged rate of change of the binary orbital energy for eccentric binaries, and thus obtain the secular rates $\langle\dot a_b\rangle$ and $\langle \dot{e}_b\rangle$. In all cases, $\langle\dot{a}_b\rangle$ is positive, i.e., the binary expands while accreting. We discuss the implications of our results for the merger of supermassive binary black holes and for the formation of close stellar binaries.

A Multimessenger Picture of the Flaring Blazar TXS 0506+056: implications for High-Energy Neutrino Emission and Cosmic Ray Acceleration

Abstract: Detection of the IceCube-170922A neutrino coincident with the flaring blazar TXS 0506+056, the first and only 3-sigma high-energy neutrino source association to date, offers a potential breakthrough in our understanding of high-energy cosmic particles and blazar physics. We present a comprehensive analysis of TXS 0506+056 during its flaring state, using newly collected Swift, NuSTAR, and X-shooter data with Fermi observations and numerical models to constrain the blazar's particle acceleration processes and multimessenger (electromagnetic and high-energy neutrino) emissions. Accounting properly for electromagnetic cascades in the emission region, we find a physically-consistent picture only within a hybrid leptonic scenario, with gamma-rays produced by external inverse-Compton processes and high-energy neutrinos via a radiatively-subdominant hadronic component. We derive robust constraints on the blazar's neutrino and cosmic-ray emissions and demonstrate that, because of cascade effects, the 0.1-100keV emissions of TXS 0506+056 serve as a better probe of its hadronic acceleration and high-energy neutrino production processes than its GeV-TeV emissions. If the IceCube neutrino association holds, physical conditions in the TXS 0506+056 jet must be close to optimal for high-energy neutrino production, and are not favorable for ultra-high-energy cosmic-ray acceleration. Alternatively, the challenges we identify in generating a significant rate of IceCube neutrino detections from TXS 0506+056 may disfavor single-zone models. In concert with continued operations of the high-energy neutrino observatories, we advocate regular X-ray monitoring of TXS 0506+056 and other blazars in order to test single-zone blazar emission models, clarify the nature and extent of their hadronic acceleration processes, and carry out the most sensitive possible search for additional multimessenger sources.

Depth of Ultra High Energy Cosmic Ray Induced Air Shower Maxima Measured by the Telescope Array Black Rock and Long Ridge FADC Fluorescence Detectors and Surface Array in Hybrid Mode

Abstract: The Telescope Array observatory utilizes fluorescence detectors and surface detectors to observe air showers produced by ultra high energy cosmic rays in the Earth's atmosphere Cosmic ray events observed in this way are termed hybrid data The depth of air shower maximum is related to the mass of the primary particle that generates the shower This paper reports on shower maxima data collected over 85 years using the Black Rock Mesa and Long Ridge fluorescence detectors in conjunction with the array of surface detectors We compare the means and standard deviations of the observed $X_{\mathrm{max}}$ distributions with Monte Carlo $X_{\mathrm{max}}$ distributions of unmixed protons, helium, nitrogen, and iron, all generated using the QGSJet~II-04 hadronic model We also perform an unbinned maximum likelihood test of the observed data, which is subjected to variable systematic shifting of the data $X_{\mathrm{max}}$ distributions to allow us to test the full distributions, and compare them to the Monte Carlo to see which elements are not compatible with the observed data For all energy bins, QGSJet~II-04 protons are found to be compatible with Telescope Array hybrid data at the 95% confidence level after some systematic $X_{\mathrm{max}}$ shifting of the data Three other QGSJet~II-04 elements are found to be compatible using the same test procedure in an energy range limited to the highest energies where data statistics are sparse

Large-scale cosmic-ray anisotropies above 4 EeV measured by the Pierre Auger Observatory

Abstract: We present a detailed study of the large-scale anisotropies of cosmic rays with energies above 4 EeV measured using the Pierre Auger Observatory. For the energy bins [4,8] EeV and $E\geq 8$ EeV, the most significant signal is a dipolar modulation in right ascension at energies above 8 EeV, as previously reported. In this paper we further scrutinize the highest-energy bin by splitting it into three energy ranges. We find that the amplitude of the dipole increases with energy above 4 EeV. The growth can be fitted with a power law with index $\beta=0.79\pm 0.19$. The directions of the dipoles are consistent with an extragalactic origin of these anisotropies at all the energies considered. Additionally we have estimated the quadrupolar components of the anisotropy: they are not statistically significant. We discuss the results in the context of the predictions from different models for the distribution of ultrahigh-energy sources and cosmic magnetic fields.

A turnover in the radio light curve of GW170817

Abstract: We present 2-9 GHz radio observations of GW170817 covering the period 125-200 days post-merger, taken with the Australia Telescope Compact Array and the Karl G. Jansky Very Large Array. Our observations demonstrate that the radio afterglow peaked at $149\pm2$ days post-merger and is now declining in flux density. We see no evidence for evolution in the radio-only spectral index, which remains consistent with optically-thin synchrotron emission connecting the radio, optical, and X-ray regimes. The peak implies a total energy in the synchrotron-emitting component of a ${\rm few}\times 10^{50}\,$erg. The temporal decay rate is most consistent with mildly- or non-relativistic material and we do not see evidence for a very energetic off-axis jet, but we cannot distinguish between a lower-energy jet and more isotropic emission.