This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, demonstrate the EVN's potential to generate new and exciting results that will transform our view of the cosmos. Together with e-MERLIN, the EVN provides a range of baseline lengths that permit unique studies of faint radio sources to be made over a wide range of spatial scales. The science cases are reviewed in six chapters that cover the following broad areas: cosmology, galaxy formation and evolution, innermost regions of active galactic nuclei, explosive phenomena and transients, stars and stellar masers in the Milky Way, celestial reference frames and space applications. The document concludes with identifying the synergies with other radio, as well as multi-band/multi-messenger instruments, and provide the recommendations for future improvements. The appendices briefly describe other radio VLBI arrays, the technological framework for EVN developments, and a selection of spectral lines of astrophysical interest below 100 GHz. The document includes a glossary for non-specialists, and a list of acronyms at the end.

We report the first plausible optical electromagnetic (EM) counterpart to a (candidate) binary black hole (BBH) merger. Detected by the Zwicky Transient Facility (ZTF), the EM flare is consistent with expectations for a kicked BBH merger in the accretion disk of an active galactic nucleus (AGN), and is unlikely ($<O(0.01\%$)) due to intrinsic variability of this source. The lack of color evolution implies that it is not a supernovae and instead is strongly suggestive of a constant temperature shock. Other false-positive events, such as microlensing or a tidal disruption event, are ruled out or constrained to be $<O(0.1\%$). If the flare is associated with S190521g, we find plausible values of: total mass $ M_{\rm BBH} \sim 100 M_{\odot}$, kick velocity $v_k \sim 200\, {\rm km}\, {\rm s}^{-1}$ at $\theta \sim 60^{\circ}$ in a disk with aspect ratio $H/a \sim 0.01$ (i.e., disk height $H$ at radius $a$) and gas density $\rho \sim 10^{-10}\, {\rm g}\, {\rm cm}^{-3}$. The merger could have occurred at a disk migration trap ($a \sim 700\, r_{g}$; $r_g \equiv G M_{\rm SMBH} / c^2$, where $M_{\rm SMBH}$ is the mass of the AGN supermassive black hole). The combination of parameters implies a significant spin for at least one of the black holes in S190521g. The timing of our spectroscopy prevents useful constraints on broad-line asymmetry due to an off-center flare. We predict a repeat flare in this source due to a re-encountering with the disk in $\sim 1.6\, {\rm yr}\, (M_{\rm SMBH}/10^{8}M_{\odot})\, (a/10^{3}r_{g})^{3/2}$.

Radio observations of tidal disruption events (TDEs) probe material ejected by the disruption of stars by supermassive black holes (SMBHs), uniquely tracing the formation and evolution of jets and outflows, revealing details of the disruption hydrodynamics, and illuminating the environments around previously-dormant SMBHs. To date, observations reveal a surprisingly diverse population. A small fraction of TDEs (at most a few percent) have been observed to produce radio-luminous mildly relativistic jets. The remainder of the population are radio quiet, producing less luminous jets, non-relativistic outflows or, possibly, no radio emission at all. Here, we review the radio observations that have been made of TDEs to date and discuss possible explanations for their properties, focusing on detected sources and, in particular, on the two best-studied events: Sw J1644+57 and ASASSN-14li. We also discuss what we have learned about the host galaxies of TDEs from radio observations and review constraints on the rates of bright and faint radio outflows in TDEs. Upcoming X-ray, optical, near-IR, and radio surveys will greatly expand the sample of TDEs, and technological advances open the exciting possibility of discovering a sample of TDEs in the radio band unbiased by host galaxy extinction.

Cosmic neutrinos provide a unique window into the otherwise-hidden mechanism of particle acceleration in astrophysical objects. A flux of high-energy neutrinos was discovered in 2013, and the IceCube Collaboration recently associated one high-energy neutrino with a flare from the relativistic jet of an active galaxy pointed towards the Earth. However a combined analysis of many similar active galaxies revealed no excess from the broader population, leaving the vast majority of the cosmic neutrino flux unexplained. Here we present the association of a radio-emitting tidal disruption event (AT2019dsg) with another high-energy neutrino, identified as part of our systematic search for optical counterparts to high-energy neutrinos with the Zwicky Transient Facility (ZTF). The probability of finding any radio-emitting tidal disruption event by chance is 0.5%, while the probability of finding one as bright in bolometric energy flux as AT2019dsg is 0.2%. Our electromagnetic observations can be explained through a multi-zone model, with radio analysis revealing a central engine, embedded in a UV photosphere, that powers an extended synchrotron-emitting outflow. This provides an ideal site for PeV neutrino production. The association suggests that tidal disruption events contribute to the cosmic neutrino flux. Unlike previous work which considered the rare subset of tidal disruption events with relativistic jets, our observations of AT2019dsg suggest an empirical model with a mildly-relativistic outflow.

As part of an on-going study of radio transients in Epoch 1 (2017-2019) of the Very Large Array Sky Survey (VLASS), we have discovered a sample of 0.2<z<3.2 active galactic nuclei (AGN) selected in the optical/infrared that have recently brightened dramatically in the radio. These sources would have previously been classified as radio-quiet based on upper limits from the Faint Images of the Radio Sky at Twenty-centimeters (FIRST; 1993-2011) survey; however, they are now consistent with radio-loud quasars. We present a quasi-simultaneous, multi-band (1-18 GHz) VLA follow-up campaign of our sample of AGN with extreme radio variability. We conclude that the radio properties are most consistent with AGN that have recently launched jets within the past few decades, potentially making them among the youngest radio AGN known.

We present nearly 500 days of observations of the tidal disruption event ASASSN-18pg, spanning from 54 days before peak light to 441 days after peak light. Our dataset includes X-ray, UV, and optical photometry, optical spectroscopy, radio observations, and the first published spectropolarimetric observations of a TDE. ASASSN-18pg was discovered on 2018 July 11 by the All-Sky Automated Survey for Supernovae (ASAS-SN) at a distance of $d=78.6$ Mpc, and with a peak UV magnitude of $m\simeq14$ it is both one of the nearest and brightest TDEs discovered to-date. The photometric data allow us to track both the rise to peak and the long-term evolution of the TDE. ASASSN-18pg peaked at a luminosity of $L\simeq2.2\times10^{44}$ erg s$^{-1}$, and its late-time evolution is shallower than a flux $\propto t^{-5/3}$ power-law model, similar to what has been seen in other TDEs. ASASSN-18pg exhibited Balmer lines and spectroscopic features consistent with Bowen fluorescence prior to peak which remained detectable for roughly 225 days after peak. Analysis of the two-component H$\alpha$ profile indicates that, if they are the result of reprocessing of emission from the accretion disk, the different spectroscopic lines may be coming from regions between $\sim10$ and $\sim60$ light-days from the black hole. No X-ray emission is detected from the TDE and there is no evidence of a jet or strong outflow detected in the radio. Our spectropolarimetric observations give no strong evidence for significant asphericity in the emission region, with the emission region having an axis ratio of at least $\sim0.65$.

We explore the galaxy structure of four tidal disruption event (TDE) host galaxies on 30 pc to kpc scales using HST WFC3 multi-band imaging. The star formation histories of these hosts are diverse, including one post-starburst galaxy (ASASSN-14li), two hosts with recent weak starbursts (ASASSN-14ae and iPTF15af), and one early type (PTF09ge). Compared to early type galaxies of similar stellar masses, the TDE hosts have higher central surface brightnesses and stellar mass surface densities on 30-100 pc scales. The TDE hosts do not show the large, kpc-scale tidal disruptions seen in some post-starburst galaxies; the hosts have low morphological asymmetries similar to those of early type galaxies. The lack of strong asymmetries are inconsistent with a recent major (~1:1 mass) merger, although minor ($\lesssim$1:3) mergers are possible. Given the time elapsed since the end of the starbursts in the three post-burst TDE hosts and the constraints on the merger mass ratios, it is unlikely that a bound supermassive black hole binary (SMBHB) has had time to coalesce. The TDE hosts have low central (<140 pc) ellipticities compared to early type galaxies. The low central ellipticities disfavor a strong radial anisotropy as the cause for the enhanced TDE rate, although we cannot rule out eccentric disks at the scale of the black hole gravitational radius of influence (~1 pc). These observations suggest that the high central stellar densities are a more important driver than SMBHBs or radial anisotropies in increasing the TDE rate in galaxies with recent starbursts.

While tidal disruption events (TDEs) have long been heralded as laboratories for the study of quiescent black holes, the small number of known TDEs and uncertainties in their emission mechanism have hindered progress towards this promise. Here present 17 new TDEs that have been detected recently by the Zwicky Transient Facility along with Swift UV and X-ray follow-up observations. Our homogeneous analysis of the optical/UV light curves, including 22 previously known TDEs from the literature, reveals a clean separation of light curve properties with spectroscopic class. The TDEs with Bowen fluorescence features in their optical spectra have smaller blackbody radii, as well as longer rise times and higher disruption rates compared to the rest of the sample. The Bowen fluorescence mechanism requires a high density which can be reached at smaller radii, which in turn yields longer diffusion timescales. Thus, the difference in rise times suggests the pre-peak TDE light curves are governed not by the fallback timescale, but instead by the diffusion of photons through the tidal debris. The small subset of TDEs that show only helium emission lines in their spectra have the longest rise times, the highest luminosities and the lowest rates. We also report, for the first time, the detection of soft X-ray flares from a TDE on day timescales. Based on the fact the flares peak at a luminosity similar to the optical/UV blackbody luminosity, we attribute them to brief glimpses through a reprocessing layer that otherwise obscures the inner accretion flow.

We present Chandra X-ray observations of four optically-selected tidal disruption events (TDEs) obtained 4-9 years after discovery. Three sources were detected with luminosities between 9X10^40 and 3X10^42 erg/s. The spectrum of PTF09axc is consistent with a power law of index 2.5+-0.1, whereas the spectrum of PTF09ge is very soft. The power law spectrum of PTF09axc and prior literature findings, provide evidence that TDEs transition from an early-time soft state to a late-time hard state many years after disruption. We propose that the time to peak luminosity for optical and X-ray emission may differ substantially in TDEs, with X-rays being produced or becoming observable later. This delay helps explain the differences in observed properties such as L_opt/L_ X of optically and X-ray selected TDEs. We update TDE rate predictions for the eROSITA instrument: it ranges from 3 per yr to 990 per yr, depending sensitively on the distribution of black hole spins and the time delay between disruption and peak X-ray brightness. We further predict an asymmetry in the number of retrograde and prograde disks in samples of optically and X-ray selected TDEs. The details of the observational biases can contribute to observed differences between optically and X-ray selected TDEs (with optically selected TDEs being fainter in X-rays for retrograde TDE disks).

We report the discovery of six active galactic nuclei (AGN) caught "turning on" during the first nine months of the Zwicky Transient Facility (ZTF) survey. The host galaxies were classified as LINERs by weak narrow forbidden line emission in their archival SDSS spectra, and detected by ZTF as nuclear transients. In five of the cases, we found via follow-up spectroscopy that they had transformed into broad-line AGN, reminiscent of the changing-look LINER iPTF 16bco. In one case, ZTF18aajupnt/AT2018dyk, follow-up HST UV and ground-based optical spectra revealed the transformation into a narrow-line Seyfert 1 (NLS1) with strong [Fe VII, X, XIV] and He II 4686 coronal lines. Swift monitoring observations of this source reveal bright UV emission that tracks the optical flare, accompanied by a luminous soft X-ray flare that peaks ~60 days later. Spitzer follow-up observations also detect a luminous mid-infrared flare implying a large covering fraction of dust. Archival light curves of the entire sample from CRTS, ATLAS, and ASAS-SN constrain the onset of the optical nuclear flaring from a prolonged quiescent state. Here we present the systematic selection and follow-up of this new class of changing-look LINERs, compare their properties to previously reported changing-look Seyfert galaxies, and conclude that they are a unique class of transients well-suited to test the uncertain physical processes associated with the LINER accretion state.

The origin of thermal optical and UV emission from stellar tidal disruption flares (TDFs) remains an open question. We present Hubble Space Telescope far-UV (FUV) observations of eight optical/UV selected TDFs 5-10 years post-peak. Six sources are cleanly detected, showing point-like FUV emission from the centers of their host galaxies. We discover that the light curves of TDFs from low-mass black holes ($<10^{6.5} M_\odot$) show significant late-time flattening. Conversely, FUV light curves from high-mass black hole TDFs are generally consistent with an extrapolation from the early-time light curve. The observed late-time emission cannot be explained by existing models for early-time TDF light curves (i.e. reprocessing or circularization shocks), but is instead consistent with a viscously spreading, unobscured accretion disk. These disk models can only reproduce the observed FUV luminosities, however, if they are assumed to be thermally and viscously stable, in contrast to the simplest predictions of alpha-disk theory. For one TDF in our sample, we measure an upper limit to the UV luminosity that is significantly lower than expectations from theoretical modeling and an extrapolation of the early-time light curve. This dearth of late-time emission could be due to a disk instability/state change absent in the rest of the sample. The disk models that explain the late-time UV detections solve the TDF "missing energy problem" by radiating a rest-mass energy of ~0.1 solar mass over a period of decades, primarily in extreme UV wavelengths.

We present new medium resolution, optical long-slit spectra of a sample of 6 UV/optical and 15 X-ray selected tidal disruption event candidate host galaxies. We measure emission line ratios from the optical spectra, finding that the large majority of hosts are quiescent galaxies, while those displaying emission lines are generally consistent with star-formation dominated environments; only 3 sources show clear evidence of nuclear activity. We measure bulge velocity dispersions using absorption lines and infer host black hole (BH) masses using the M-$\sigma$ relation. While the optical and X-ray host BH masses are statistically consistent with coming from the same parent distribution, the optical host distribution has a visible peak near $M_{\rm BH} \sim 10^6 M_\odot$, whereas the X-ray host distribution appears flat in $M_{\rm BH}$. We find a subset of X-ray selected candidates that are hosted in galaxies significantly less luminous (M$_{\rm g}$$\sim$-16) and less massive (stellar mass$\sim$10$^{8.5-9}$M$_{\odot}$) than those of optical events. Using statistical tests we find suggestive evidence that, in terms of black hole mass, stellar mass and absolute magnitude, the hard X-ray hosts differ from the UV/optical and soft X-ray samples. Similar to individual studies, we find that the emission region size for the X-ray sample is much smaller than the optical emission region, consistent with a compact accretion disk. We find a typical Eddington ratio of the X-ray emission is $\sim$0.01, as opposed to the optical events which have L$_{\rm BB}$$\sim$L$_{\rm Edd}$.

We present optical spectroscopic and Swift UVOT/XRT observations of the X-ray and UV/optical bright tidal disruption event (TDE) AT 2018fyk/ASASSN-18ul discovered by ASAS-SN. The Swift lightcurve is atypical for a TDE, entering a plateau after $\sim$40 days of decline from peak. After 80 days the UV/optical lightcurve breaks again to decline further, while the X-ray emission becomes brighter and harder. In addition to broad H, He and potentially O/Fe lines, narrow emission lines emerge in the optical spectra during the plateau phase. We identify both high ionisation (O III) and low ionisation (Fe II) lines, which are visible for $\sim$45 days. We similarly identify Fe II lines in optical spectra of ASASSN-15oi 330 d after discovery, indicating that a class of Fe-rich TDEs exists. The spectral similarity between AT 2018fyk, narrow-line Seyfert 1 galaxies and some extreme coronal line emitters suggests that TDEs are capable of creating similar physical conditions in the nuclei of galaxies. The Fe II lines can be associated with the formation of a compact accretion disk, as the emergence of low ionisation emission lines requires optically thick, high density gas. Taken together with the plateau in X-ray and UV/optical luminosity this indicates that emission from the central source is efficiently reprocessed into UV/optical wavelengths. Such a two-component lightcurve is very similar to that seen in the TDE candidate ASASSN-15lh, and is a natural consequence of a highly relativistic orbital pericenter.

Both multi-messenger astronomy and new high-throughput wide-field surveys require flexible tools for the selection and analysis of astrophysical transients. We here introduce the Alert Management, Photometry and Evaluation of Lightcurves (AMPEL) system, an analysis framework designed for high-throughput surveys and suited for streamed data. AMPEL combines the functionality of an alert broker with a generic framework capable of hosting user-contributed code, that encourages provenance and keeps track of the varying information states that a transient displays. The latter concept includes information gathered over time and data policies such as access or calibration levels. We describe a novel ongoing real-time multi-messenger analysis using AMPEL to combine IceCube neutrino data with the alert streams of the Zwicky Transient Facility (ZTF). We also reprocess the first four months of ZTF public alerts, and compare the yields of more than 200 different transient selection functions to quantify efficiencies for selecting Type Ia supernovae that were reported to the Transient Name Server (TNS). We highlight three channels suitable for (1) the collection of a complete sample of extragalactic transients, (2) immediate follow-up of nearby transients and (3) follow-up campaigns targeting young, extragalactic transients. We confirm ZTF completeness in that all TNS supernovae positioned on active CCD regions were detected. AMPEL can assist in filtering transients in real time, running alert reaction simulations, the reprocessing of full datasets as well as in the final scientific analysis of transient data. This text introduces how users can design their own channels for inclusion in the AMPEL live instance that parses the ZTF stream and the real-time submission of high quality extragalactic supernova candidates to the TNS.

We report the discovery of non-stellar hydrogen Balmer and metastable helium absorption lines accompanying a transient, high-velocity (0.05$c$) broad absorption line (BAL) system in the optical spectra of the tidal disruption event (TDE) AT2018zr ($z=0.071$). In the HST UV spectra, absorption of high- and low-ionization lines are also present at this velocity, making AT2018zr resemble a low-ionization broad absorption line (LoBAL) QSO. We conclude that these transient absorption features are more likely to arise in fast outflows produced by the TDE than absorbed by the unbound debris. In accordance with the outflow picture, we are able to reproduce the flat-topped H$\alpha$ emission in a spherically expanding medium, without invoking the typical prescription of an elliptical disk. We also report the appearance of narrow ($\sim$1000~km~s$^{-1}$) NIII$\lambda$4640, HeII$\lambda4686$, H$\alpha$, and H$\beta$, emission in the late-time optical spectra of AT2018zr, which may be a result of UV continuum hardening at late time as observed by Swift. Including AT2018zr, we find a high association rate (3 out of 4) of BALs in the UV spectra of TDEs. This suggests that outflows may be ubiquitous among TDEs and may be less sensitive to viewing angle effects compared to QSO outflows.

We present results of the Sky2Night project: a systematic, unbiased search for fast optical transients with the Palomar Transient Factory. We have observed 407 deg$^2$ in $R$-band for 8 nights at a cadence of 2 hours. During the entire duration of the project, the 4.2m William Herschel Telescope on La Palma was dedicated to obtaining identification spectra for the detected transients. During the search, we found 12 supernovae, 10 outbursting cataclysmic variables, 9 flaring M-stars, 3 flaring active Galactic nuclei and no extragalactic fast optical transients. Using this systematic survey for transients, we have calculated robust observed rates for the detected types of transients, and upper limits of the rate of extragalactic fast optical transients of $\mathcal{R}<37\times 10^{-4}$deg$^{-2}$d$^{-1}$ and $\mathcal{R}<9.3\times 10^{-4}$deg$^{-2}$d$^{-1}$ for timescales of 4h and 1d and a limiting magnitude of $R\approx19.7$. We use the results of this project to determine what kind of and how many astrophysical false positives we can expect when following up gravitational wave detections in search for kilonovae.

We present Zwicky Transient Facility (ZTF) observations of the tidal disruption flare AT2018zr/PS18kh reported by Holoien et al. and detected during ZTF commissioning. The ZTF light curve of the tidal disruption event (TDE) samples the rise-to-peak exceptionally well, with 50 days of g- and r-band detections before the time of maximum light. We also present our multi-wavelength follow-up observations, including the detection of a thermal (kT~100 eV) X-ray source that is two orders of magnitude fainter than the contemporaneous optical/UV blackbody luminosity, and a stringent upper limit to the radio emission. We use observations of 128 known active galactic nuclei (AGN) to assess the quality of the ZTF astrometry, finding a median host-flare distance of 0.2" for genuine nuclear flares. Using ZTF observations of variability from known AGN and supernovae we show how these sources can be separated from TDEs. A combination of light-curve shape, color, and location in the host galaxy can be used to select a clean TDE sample from multi-band optical surveys such as ZTF or LSST.

The Zwicky Transient Facility (ZTF) is a new optical time-domain survey that uses the Palomar 48-inch Schmidt telescope. A custom-built wide-field camera provides a 47 deg$^2$ field of view and 8 second readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory (PTF). We describe the design and implementation of the camera and observing system. The ZTF data system at the Infrared Processing and Analysis Center provides near-real-time reduction to identify moving and varying objects. We outline the analysis pipelines, data products, and associated archive. Finally, we present on-sky performance analysis and first scientific results from commissioning and the early survey. ZTF's public alert stream will serve as a useful precursor for that of the Large Synoptic Survey Telescope.

The Zwicky Transient Facility is a large optical survey in multiple filters producing hundreds of thousands of transient alerts per night. We describe here various machine learning (ML) implementations and plans to make the maximal use of the large data set by taking advantage of the temporal nature of the data, and further combining it with other data sets. We start with the initial steps of separating bogus candidates from real ones, separating stars and galaxies, and go on to the classification of real objects into various classes. Besides the usual methods (e.g., based on features extracted from light curves) we also describe early plans for alternate methods including the use of domain adaptation, and deep learning. In a similar fashion we describe efforts to detect fast moving asteroids. We also describe the use of the Zooniverse platform for helping with classifications through the creation of training samples, and active learning. Finally we mention the synergistic aspects of ZTF and LSST from the ML perspective.

The Zwicky Transient Facility (ZTF), a public-private enterprise, is a new time domain survey employing a dedicated camera on the Palomar 48-inch Schmidt telescope with a 47 deg$^2$ field of view and 8 second readout time. It is well positioned in the development of time domain astronomy, offering operations at 10% of the scale and style of the Large Synoptic Survey Telescope (LSST) with a single 1-m class survey telescope. The public surveys will cover the observable northern sky every three nights in g and r filters and the visible Galactic plane every night in g and r. Alerts generated by these surveys are sent in real time to brokers. A consortium of universities which provided funding ("partnership") are undertaking several boutique surveys. The combination of these surveys producing one million alerts per night allows for exploration of transient and variable astrophysical phenomena brighter than r $\sim$ 20.5 on timescales of minutes to years. We describe the primary science objectives driving ZTF including the physics of supernovae and relativistic explosions, multi-messenger astrophysics, supernova cosmology, active galactic nuclei and tidal disruption events, stellar variability, and Solar System objects.

Tidal disruption events (TDEs) are rare, 10^(-7)/yr/Mpc^3 (Hung et al. 2018), yet the large survey volume of LSST implies a very large detection rate of 200/yr/(1000 deg^2) (van Velzen et al. 2011), a factor of 250 increase in the detection capability of the current generation of optical synoptic surveys, e.g. ZTF, ASAS-SN, Pan-STARRS, and ATLAS. The goal of this LSST cadence white paper is to determine which survey strategy will ensure the efficient selection and characterization of TDEs in the LSST Wide-Fast-Deep Survey transient alert stream. We conclude that the baseline cadence design fails to 1) measure the u-r color and color evolution of transients, a critical parameter for distinguishing TDEs from SNe, and to 2) catch the pre-peak light curves of transients, an essential measurement for probing their rise times, which are in turn a probe of black hole mass in TDEs. If we do not harvest the fruits of the LSST transient alert stream with photometric classification and early detections, both TDE and SN science will be greatly limited. Hence, we propose for a "smart" and "colorful" rolling cadence in the Wide-Fast Deep (WFD) Survey, that allows for efficient photometric transient classification from well sampled multi-band light curves, with the 20,000 deg^2 survey divided into eight 2500 deg^2 strips each observed for one year in Years 2-9, with the full WFD area observed in Years 1 & 10. This will yield a legacy sample of 200 TDEs per year with early detections in u, g, and r bands for efficient classification and full light curve characterization.

The tidal disruption and subsequent accretion of a star by a supermassive black hole can be used as a laboratory to study the physics of relativistic jets. The ngVLA is the only planned instrument that can both discover and characterize a large number of these short-lived radio sources. In particular the high-frequency capabilities of the ngVLA enable this important leap forward. Multi-frequency radio follow-up observations (3 - 100 GHz) of tidal disruption events found in optical or X-ray surveys will provide a measurement of the jet efficiency as a function of black hole spin, thus enabling a direct test of the prediction that relativistic jets require high spin. Hundreds of tidal disruption jets will be discovered in a blind ngVLA survey for radio transients. By including VLBI observations with the ngVLA Long Baseline Array, we can resolve some of these sources, obtaining a robust measurement of the jet launch date and the magnetic field strength. From the thermal emission of the tidal disruption flare we can measure the accretion rate at this launch date, thus providing another unique opportunity to identify the conditions that lead to jet production.

We present results from a systematic selection of tidal disruption events (TDEs) in a wide-area (4800~deg$^2$), $g+R$ band, Intermediate Palomar Transient Factory (iPTF) experiment. Our selection targets typical optically-selected TDEs: bright ($>$60\% flux increase) and blue transients residing in the center of red galaxies. Using photometric selection criteria to down-select from a total of 493 nuclear transients to a sample of 26 sources, we then use follow-up UV imaging with the Neil Gehrels Swift Telescope, ground-based optical spectroscopy, and light curve fitting to classify them as 14 Type Ia supernovae (SNe Ia), 9 highly variable active galactic nuclei (AGNs), 2 confirmed TDEs, and 1 potential core-collapse supernova. We find it possible to filter AGNs by employing a more stringent transient color cut ($g-r <$ $-$0.2 mag); further, UV imaging is the best discriminator for filtering SNe, since SNe Ia can appear as blue, optically, as TDEs in their early phases. However, when UV-optical color is unavailable, higher precision astrometry can also effectively reduce SNe contamination in the optical. Our most stringent optical photometric selection criteria yields a 4.5:1 contamination rate, allowing for a manageable number of TDE candidates for complete spectroscopic follow-up and real-time classification in the ZTF era. We measure a TDE per galaxy rate of 1.7$^{+2.9}_{-1.3}$ $\times$10$^{-4}$ gal$^{-1}$ yr$^{-1}$ (90\% CL in Poisson statistics). This does not account for TDEs outside our selection criteria, thus may not reflect the total TDE population, which is yet to be fully mapped.

The tidal disruption of a star by a supermassive black hole can result in transient radio emission. The electrons producing these synchrotron radio flares could either be accelerated inside a relativistic jet or externally by shocks resulting from an outflow interacting with the circumnuclear medium. Until now, evidence for the internal emission mechanism has been lacking; nearly all tidal disruption flare studies have adopted the external shock model to explain the observed properties of radio flares. Here we report a result that presents a challenge to external emission models: we discovered a cross-correlation between the soft X-ray (0.3-1 keV) and 16 GHz radio flux of the tidal disruption flare ASASSN-14li. Variability features in the X-ray light curve appear again in the radio light curve, but after a time lag of about 13 days. This demonstrates that soft X-ray emitting accretion disk regulates the radio emission. This coupling appears to be inconsistent with all previous external emission models for this source but is naturally explained if the radio emission originates from a freely expanding jet. We show that emission internal to an adiabatically expanding jet can also reproduce the observed evolution of the radio spectral energy distribution. Furthermore, both the correlation between X-ray and radio luminosity as well as our radio spectral modeling imply an approximately linear coupling between the accretion rate and jet power.

A tidal disruption event (TDE) ensues when a star passes too close to the supermassive black hole (SMBH) in a galactic center and is ripped apart by the tidal field of the SMBH. The gaseous debris produced in a TDE can power a bright electromagnetic flare as it is accreted by the SMBH; so far, several dozen TDE candidates have been observed. For SMBHs with masses above $\sim 10^7 M_\odot$, the tidal disruption of solar-type stars occurs within ten gravitational radii of the SMBH, implying that general relativity (GR) is needed to describe gravity. Three promising signatures of GR in TDEs are: (1) a super-exponential cutoff in the volumetric TDE rate for SMBH masses above $\sim 10^8 M_\odot$ due to direct capture of tidal debris by the event horizon, (2) delays in accretion disk formation (and a consequent alteration of the early-time light curve) caused by the effects of relativistic precession on stream circularization, and (3) quasi-periodic modulation of X-ray emission due to global precession of misaligned accretion disks and the jets they launch. We review theoretical models and simulations of TDEs in Newtonian gravity, then describe how relativistic modifications give rise to these proposed observational signatures, as well as more speculative effects of GR. We conclude with a brief summary of TDE observations and the extent to which they show indications of these predicted relativistic signatures.

The tidal disruption of a star by a massive black hole is expected to yield a luminous flare of thermal emission. About two dozen of these stellar tidal disruption flares (TDFs) may have been detected in optical transient surveys. However, explaining the observed properties of these events within the tidal disruption paradigm is not yet possible. This theoretical ambiguity has led some authors to suggest that optical TDFs are due to a different process, such as a nuclear supernova or accretion disk instabilities. Here we present a test of a fundamental prediction of the tidal disruption event scenario: a suppression of the flare rate due to the direct capture of stars by the black hole. Using a recently compiled sample of candidate TDFs with black hole mass measurements, plus a careful treatment of selection effects in this flux-limited sample, we confirm that the dearth of observed TDFs from high-mass black holes is statistically significant. All the TDF impostor models we consider fail to explain the observed mass function; the only scenario that fits the data is a suppression of the rate due to direct captures. We find that this suppression can explain the low volumetric rate of the luminous TDF candidate ASASSN-15lh, thus supporting the hypothesis that this flare belongs to the TDF family. Our work is the first to present the optical TDF luminosity function. A steep power law is required to explain the observed rest-frame g-band luminosity, $dN/dL_{g} \propto L_{g}^{-2.5}$. The mean event rate of the flares in our sample is about $1 \times10^{-4}$ per galaxy per year, consistent with the theoretically expected tidal disruption rate.

We report on late time radio and X-ray observations of the tidal disruption event candidate ASASSN-14li, covering the first 1000 days of the decay phase. For the first $\sim200$ days the radio and X-ray emission fade in concert. This phase is better fit by an exponential decay at X-ray wavelengths, while the radio emission is well described by either an exponential or the canonical $t^{-5/3}$ decay assumed for tidal disruption events. The correlation between radio and X-ray emission during this period can be fit as $L_{R}\propto L_{X}^{1.9\pm0.2}$. After 400 days the radio emission at $15.5\,\textrm{GHz}$ has reached a plateau level of $244\pm8\,\mu\textrm{Jy}$ which it maintains for at least the next 600 days, while the X-ray emission continues to fade exponentially. This steady level of radio emission is likely due to relic radio lobes from the weak AGN-like activity implied by historical radio observations. We note that while most existing models are based upon the evolution of ejecta which are decoupled from the central black hole, the radio : X-ray correlation during the declining phase is also consistent with core jet emission coupled to a radiatively efficient accretion flow.

The progenitors of astronomical transients are linked to a specific stellar population and galactic environment, and observing their host galaxies hence constrains the physical nature of the transient itself. Here, we use imaging from the Hubble Space Telescope, and spatially-resolved, medium resolution spectroscopy from the Very Large Telescope obtained with X-Shooter and MUSE to study the host of the very luminous transient ASASSN-15lh. The dominant stellar population at the transient site is old (around 1 to 2 Gyr), without signs of recent star-formation. We also detect emission from ionized gas, originating from three different, time-invariable, narrow components of collisionally-excited metal and Balmer lines. The ratios of emission lines in the Baldwin-Phillips-Terlevich diagnostic diagram indicate that the ionization source is a weak Active Galactic Nucleus with a black hole mass of $M_\bullet = 5_{-3}^{+8}\cdot10^{8} M_\odot$, derived through the $M_\bullet$-$\sigma$ relation. The narrow line components show spatial and velocity offsets on scales of 1 kpc and 500 km/s, respectively; these offsets are best explained by gas kinematics in the narrow-line region. The location of the central component, which we argue is also the position of the supermassive black hole, aligns with that of the transient within an uncertainty of 170 pc. Using this positional coincidence as well as other similarities with the hosts of Tidal Disruption Events, we strengthen the argument that the transient emission observed as ASASSN-15lh is related to the disruption of a star around a supermassive black hole, most probably spinning with a Kerr parameter $a_\bullet\gtrsim0.5$.

The mass of the central black hole in a galaxy that hosted a tidal disruption event (TDE) is an important parameter in understanding its energetics and dynamics. We present the first homogeneously measured black hole masses of a complete sample of 12 optically/UV selected TDE host galaxies (down to $g_{host}$$\leq$22 mag and $z$=0.37) in the Northern sky. The mass estimates are based on velocity dispersion measurements, performed on late time optical spectroscopic observations. We find black hole masses in the range 3$\times$10$^5$ M$_{\odot}$$\leq$M$_{\rm BH}$$\leq$2$\times$10$^7$ M$_{\odot}$. The TDE host galaxy sample is dominated by low mass black holes ($\sim$10$^6$ M$_{\odot}$), as expected from theoretical predictions. The blackbody peak luminosity of TDEs with M$_{\rm BH}$$\leq$10$^{7.1}$ M$_{\odot}$ is consistent with the Eddington limit of the SMBH, whereas the two TDEs with M$_{\rm BH}$$\geq$10$^{7.1}$ M$_{\odot}$ have peak luminosities below their SMBH Eddington luminosity, in line with the theoretical expectation that the fallback rate for M$_{\rm BH}$$\geq$10$^{7.1}$ M$_{\odot}$ is sub-Eddington. In addition, our observations suggest that TDEs around lower mass black holes evolve faster. These findings corroborate the standard TDE picture in 10$^6$ M$_{\odot}$ black holes. Our results imply an increased tension between observational and theoretical TDE rates. By comparing the blackbody emission radius with theoretical predictions, we conclude that the optical/UV emission is produced in a region consistent with the stream self-intersection radius of shallow encounters, ruling out a compact accretion disk as the direct origin of the blackbody radiation at peak brightness.

We carried out the first multi-wavelength (optical/UV and X-ray) photometric reverberation mapping of a tidal disruption flare (TDF) ASASSN-14li. We find that its X-ray variations are correlated with and lag the optical/UV fluctuations by 32$\pm$4 days. Based on the direction and the magnitude of the X-ray time lag, we rule out X-ray reprocessing and direct emission from a standard circular thin disk as the dominant source of its optical/UV emission. The lag magnitude also rules out an AGN disk-driven instability as the origin of ASASSN-14li and thus strongly supports the tidal disruption picture for this event and similar objects. We suggest that the majority of the optical/UV emission likely originates from debris stream self-interactions. Perturbations at the self-interaction sites produce optical/UV variability and travel down to the black hole where they modulate the X-rays. The time lag between the optical/UV and the X-rays variations thus correspond to the time taken by these fluctuations to travel from the self-interaction site to close to the black hole. We further discuss these time lags within the context of the three variants of the self-interaction model. High-cadence monitoring observations of future TDFs will be sensitive enough to detect these echoes and would allow us to establish the origin of optical/UV emission in TDFs in general.

We report on the discovery in the LOFAR Multifrequency Snapshot Sky Survey (MSSS) of a giant radio galaxy (GRG) with a projected size of $2.56 \pm 0.07$ Mpc projected on the sky. It is associated with the galaxy triplet UGC 9555, within which one is identified as a broad-line galaxy in the Sloan Digital Sky Survey (SDSS) at a redshift of $0.05453 \pm 1 \times 10^{-5} $, and with a velocity dispersion of $215.86 \pm 6.34$ km/s. From archival radio observations we see that this galaxy hosts a compact flat-spectrum radio source, and we conclude that it is the active galactic nucleus (AGN) responsible for generating the radio lobes. The radio luminosity distribution of the jets, and the broad-line classification of the host AGN, indicate this GRG is orientated well out of the plane of the sky, making its physical size one of the largest known for any GRG. Analysis of the infrared data suggests that the host is a lenticular type galaxy with a large stellar mass ($\log~\mathrm{M}/\mathrm{M}_\odot = 11.56 \pm 0.12$), and a moderate star formation rate ($1.2 \pm 0.3~\mathrm{M}_\odot/\mathrm{year}$). Spatially smoothing the SDSS images shows the system around UGC 9555 to be significantly disturbed, with a prominent extension to the south-east. Overall, the evidence suggests this host galaxy has undergone one or more recent moderate merger events and is also experiencing tidal interactions with surrounding galaxies, which have caused the star formation and provided the supply of gas to trigger and fuel the Mpc-scale radio lobes.

High-redshift quasars are important to study galaxy and active galactic nuclei (AGN) evolution, test cosmological models, and study supermassive black hole growth. Optical searches for high-redshift sources have been very successful, but radio searches are not hampered by dust obscuration and should be more effective at finding sources at even higher redshifts. Identifying high-redshift sources based on radio data is, however, not trivial. Here we report on new multi-frequency Giant Metrewave Radio Telescope (GMRT) observations of eight z>4.5 sources previously studied at high angular resolution with very long baseline interferometry (VLBI). Combining these observations with those from the literature, we construct broad-band radio spectra of all 30 z>4.5 sources that have been observed with VLBI. In the sample we found flat, steep and peaked spectra in approximately equal proportions. Despite several selection effects, we conclude that the z>4.5 VLBI (and likely also non-VLBI) sources have diverse spectra and that only about a quarter of the sources in the sample have flat spectra. Previously, the majority of high-redshift radio sources were identified based on their ultra-steep spectra (USS). Recently a new method has been proposed to identify these objects based on their megahertz-peaked spectra (MPS). Neither method would have identified more than 18% of the high-redshift sources in this sample. More effective methods are necessary to reliably identify complete samples of high-redshift sources based on radio data.

Time domain optical surveys have discovered roughly a dozen candidate stellar tidal disruption flares in the last five years, and future surveys like the {\it Large Synoptic Survey Telescope} will likely find hundreds to thousands more. These tidal disruption events (TDEs) present an interesting puzzle: a majority of the current TDE sample is hosted by rare post-starburst galaxies, and tens of percent are hosted in even rarer E+A galaxies, which make up $\sim 0.1\%$ of all galaxies in the local universe. E+As are therefore overrepresented among TDE hosts by 1-2 orders of magnitude, a discrepancy unlikely to be accounted for by selection effects. We analyze {\it Hubble Space Telescope} photometry of one of the nearest E+A galaxies, NGC~3156, to estimate the rate of stellar tidal disruption produced as two-body relaxation diffuses stars onto orbits in the loss cone of the central supermassive black hole. The rate of TDEs produced by two-body relaxation in NGC~3156 is large when compared to other galaxies with similar black hole mass: $\dot{N}_{\rm TDE}\sim 1\times 10^{-3}~{\rm yr}^{-1}$. This suggests that the preference of TDEs for E+A hosts may be due to central stellar overdensities produced in recent starbursts.

When a star passes within the tidal radius of a supermassive black hole, it will be torn apart. For a star with the mass of the Sun ($M_\odot$) and a non-spinning black hole with a mass $<10^8 M_\odot$, the tidal radius lies outside the black hole event horizon and the disruption results in a luminous flare. Here we report observations over a period of 10 months of a transient, hitherto interpreted as a superluminous supernova. Our data show that the transient rebrightened substantially in the ultraviolet and that the spectrum went through three different spectroscopic phases without ever becoming nebular. Our observations are more consistent with a tidal disruption event than a superluminous supernova because of the temperature evolution, the presence of highly ionised CNO gas in the line of sight and our improved localisation of the transient in the nucleus of a passive galaxy, where the presence of massive stars is highly unlikely. While the supermassive black hole has a mass $> 10^8 M_\odot$, a star with the same mass as the Sun could be disrupted outside the event horizon if the black hole were spinning rapidly. The rapid spin and high black hole mass can explain the high luminosity of this event.

Stars that pass within the Roche radius of a supermassive black hole will be tidally disrupted, yielding a sudden injection of gas close to the black hole horizon which produces an electromagnetic flare. A few dozen of these flares have been discovered in recent years, but current observations provide poor constraints on the bolometric luminosity and total accreted mass of these events. Using images from the Wide-field Infrared Survey Explorer (WISE), we have discovered transient 3.4 micron emission from several previously known tidal disruption flares. The observations can be explained by dust heated to its sublimation temperature due to the intense radiation of the tidal flare. From the break in the infrared light curve we infer that this hot dust is located ~0.1 pc from the supermassive black hole. Since the dust has been heated by absorbing UV and (potentially) soft X-ray photons of the flare, the reprocessing light curve yields an estimate of the bolometric flare luminosity. For the flare PTF-09ge, we infer that the most likely value of the luminosity integrated over frequencies at which dust can absorb photons is $8\times 10^{44}$ erg/s, with a factor of 3 uncertainty due to the unknown temperature of the dust. This bolometric luminosity is a factor ~10 larger than the observed black body luminosity. Our work is the first to probe dust in the nuclei of non-active galaxies on sub-parsec scales. The observed infrared luminosity implies a covering factor ~1% for the nuclear dust in the host galaxies.

The Lockman Hole is a well-studied extragalactic field with extensive multi-band ancillary data covering a wide range in frequency, essential for characterising the physical and evolutionary properties of the various source populations detected in deep radio fields (mainly star-forming galaxies and AGNs). In this paper we present new 150-MHz observations carried out with the LOw Frequency ARray (LOFAR), allowing us to explore a new spectral window for the faint radio source population. This 150-MHz image covers an area of 34.7 square degrees with a resolution of 18.6$\times$14.7 arcsec and reaches an rms of 160 $\mu$Jy beam$^{-1}$ at the centre of the field. As expected for a low-frequency selected sample, the vast majority of sources exhibit steep spectra, with a median spectral index of $\alpha_{150}^{1400}=-0.78\pm0.015$. The median spectral index becomes slightly flatter (increasing from $\alpha_{150}^{1400}=-0.84$ to $\alpha_{150}^{1400}=-0.75$) with decreasing flux density down to $S_{150} \sim$10 mJy before flattening out and remaining constant below this flux level. For a bright subset of the 150-MHz selected sample we can trace the spectral properties down to lower frequencies using 60-MHz LOFAR observations, finding tentative evidence for sources to become flatter in spectrum between 60 and 150 MHz. Using the deep, multi-frequency data available in the Lockman Hole, we identify a sample of 100 Ultra-steep spectrum (USS) sources and 13 peaked spectrum sources. We estimate that up to 21 percent of these could have $z>4$ and are candidate high-$z$ radio galaxies, but further follow-up observations are required to confirm the physical nature of these objects.

We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.

The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy -- corrected for geometrical effects -- is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.

The Pierre Auger Observatory has begun a major Upgrade of its already impressive capabilities, with an emphasis on improved mass composition determination using the surface detectors of the Observatory. Known as AugerPrime, the upgrade will include new 4 m$^2$ plastic scintillator detectors on top of all 1660 water-Cherenkov detectors, updated and more flexible surface detector electronics, a large array of buried muon detectors, and an extended duty cycle for operations of the fluorescence detectors. This Preliminary Design Report was produced by the Collaboration in April 2015 as an internal document and information for funding agencies. It outlines the scientific and technical case for AugerPrime. We now release it to the public via the arXiv server. We invite you to review the large number of fundamental results already achieved by the Observatory and our plans for the future.

Megahertz peaked-spectrum (MPS) sources have spectra that peak at frequencies below 1 GHz in the observer's frame and are believed to be radio-loud active galactic nuclei (AGN). We recently presented a new method to search for high-redshift AGN by identifying unusually compact MPS sources. In this paper, we present European VLBI Network (EVN) observations of 11 MPS sources which we use to determine their sizes and investigate the nature of the sources with ~10 mas resolution. Of the 11 sources, we detect nine with the EVN. Combining the EVN observations with spectral and redshift information, we show that the detected sources are all AGN with linear sizes smaller than 1.1 kpc and are likely young. This shows that low-frequency colour-colour diagrams are an easy and efficient way of selecting small AGN and explains our high detection fraction (82%) in comparison to comparable surveys. Finally we argue that the detected sources are all likely compact symmetric objects and that none of the sources are blazars.

We study the radio properties of moderately obscured quasars over a range of redshifts to understand the role of radio activity in accretion using the Jansky Very Large Array (JVLA) at 6.0GHz and 1.4GHz. Our z~2.5 sample consists of optically-selected obscured quasar candidates, all of which are radio-quiet, with typical radio luminosities of $\nu L_{\nu}$[1.4 GHz] < $10^{40}$ erg s$^{-1}$. Only a single source is individually detected in our deep (rms~10 $\mu$Jy) exposures. This population would not be identified by radio-based selection methods used for distinguishing dusty star-forming galaxies and obscured active nuclei. In our pilot A-array study of z~0.5 radio-quiet quasars, we spatially resolve four of five objects on scales ~ 5 kpc and find they have steep spectral indices. Therefore, radio emission in these sources could be due to jet-driven or radiatively driven bubbles interacting with interstellar material on the scale of the host galaxy. Finally, we also study the population of ~ 200 faint (~40 $\mu$Jy - 40 mJy) radio sources observed over ~ 120 arcmin$^2$ of our data. 60% of these detections are matched in the SDSS and/or WISE and are, in roughly equal shares, active nuclei at a broad range of redshifts, passive galaxies with no other signs of nuclear activity and IR-bright but optically faint sources. Spectroscopically or photometrically confirmed star-forming galaxies constitute only a small minority of the matches. Such sensitive radio surveys allow us to address important questions of AGN evolution and evaluate the AGN contribution to the radio-quiet sky.

The tidal disruption of a star by a supermassive black hole leads to a short-lived thermal flare. Despite extensive searches, radio follow-up observations of known thermal stellar tidal disruption flares (TDFs) have not yet produced a conclusive detection. We present a detection of variable radio emission from a thermal TDF, which we interpret as originating from a newly-launched jet. The multi-wavelength properties of the source present a natural analogy with accretion state changes of stellar mass black holes, suggesting all TDFs could be accompanied by a jet. In the rest frame of the TDF, our radio observations are an order of magnitude more sensitive than nearly all previous upper limits, explaining how these jets, if common, could thus far have escaped detection.

We present the results of a four-month campaign searching for low-frequency radio transients near the North Celestial Pole with the Low-Frequency Array (LOFAR), as part of the Multifrequency Snapshot Sky Survey (MSSS). The data were recorded between 2011 December and 2012 April and comprised 2149 11-minute snapshots, each covering 175 deg^2. We have found one convincing candidate astrophysical transient, with a duration of a few minutes and a flux density at 60 MHz of 15-25 Jy. The transient does not repeat and has no obvious optical or high-energy counterpart, as a result of which its nature is unclear. The detection of this event implies a transient rate at 60 MHz of 3.9 (+14.7, -3.7) x 10^-4 day^-1 deg^-2, and a transient surface density of 1.5 x 10^-5 deg^-2, at a 7.9-Jy limiting flux density and ~10-minute time-scale. The campaign data were also searched for transients at a range of other time-scales, from 0.5 to 297 min, which allowed us to place a range of limits on transient rates at 60 MHz as a function of observation duration.

Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time-structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for "Earth-skimming" neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an $E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times 10^{19}$ eV is $E_\nu^2 dN_\nu/dE_\nu < 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~ s^{-1}~ sr^{-1}}$.

A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.

Joint contributions of the Pierre Auger Collaboration and the Telescope Array Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.

We have conducted three searches for correlations between ultra-high energy cosmic rays detected by the Telescope Array and the Pierre Auger Observatory, and high-energy neutrino candidate events from IceCube. Two cross-correlation analyses with UHECRs are done: one with 39 cascades from the IceCube `high-energy starting events' sample and the other with 16 high-energy `track events'. The angular separation between the arrival directions of neutrinos and UHECRs is scanned over. The same events are also used in a separate search using a maximum likelihood approach, after the neutrino arrival directions are stacked. To estimate the significance we assume UHECR magnetic deflections to be inversely proportional to their energy, with values $3^\circ$, $6^\circ$ and $9^\circ$ at 100 EeV to allow for the uncertainties on the magnetic field strength and UHECR charge. A similar analysis is performed on stacked UHECR arrival directions and the IceCube sample of through-going muon track events which were optimized for neutrino point-source searches.

We present a study of the detectability of transient events associated with galaxies for the Gaia European Space Agency astrometric mission. We simulated the on-board detections, and on-ground processing for a mock galaxy catalogue to establish the properties required for the discovery of transient events by Gaia, specifically tidal disruption events (TDEs) and supernovae (SNe). Transients may either be discovered by the on-board detection of a new source or by the brightening of a previously known source. We show that Gaia transients can be identified as new detections on-board for offsets from the host galaxy nucleus of 0.1--0.5,arcsec, depending on magnitude and scanning angle. The Gaia detection system shows no significant loss of SNe at close radial distances to the nucleus. We used the detection efficiencies to predict the number of transients events discovered by Gaia. For a limiting magnitude of 19, we expect around 1300 SNe per year: 65% SN Ia, 28% SN II and 7% SN Ibc, and ~20 TDEs per year.

The existence of double-double radio galaxies (DDRGs) is evidence for recurrent jet activity in AGN, as expected from standard accretion models. A detailed study of these rare sources provides new perspectives for investigating the AGN duty cycle, AGN-galaxy feedback, and accretion mechanisms. Large catalogues of radio sources provide statistical information about the evolution of the radio-loud AGN population out to high redshifts. Using wide-field imaging with the LOFAR telescope, we study both a well-known DDRG as well as a large number of radio sources in the field of view. We present a high resolution image of the DDRG B1834+620 obtained at 144 MHz using LOFAR commissioning data. Our image covers about 100 square degrees and contains over 1000 sources. The four components of the DDRG B1834+620 have been resolved for the first time at 144 MHz. Inner lobes were found to point towards the direction of the outer lobes, unlike standard FR~II sources. Polarized emission was detected in the northern outer lobe. The high spatial resolution allows the identification of a large number of small double-lobed radio sources; roughly 10% of all sources in the field are doubles with a separation smaller than 1 arcmin. The spectral fit of the four components is consistent with a scenario in which the outer lobes are still active or the jets recently switched off, while emission of the inner lobes is the result of a mix-up of new and old jet activity. From the presence of the newly extended features in the inner lobes of the DDRG, we can infer that the mechanism responsible for their formation is the bow shock that is driven by the newly launched jet. We find that the density of the small doubles exceeds the density of FR-II sources with similar properties at 1.4 GHz, but this difference becomes smaller for low flux densities.

Contributions of the Pierre Auger Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.

We report the discovery and characterisation of a deeply eclipsing AM CVn-system, Gaia14aae (= ASSASN-14cn). Gaia14aae was identified independently by the All-Sky Automated Survey for Supernovae (ASAS-SN; Shappee et al. 2014) and by the Gaia Science Alerts project, during two separate outbursts. A third outburst is seen in archival Pan-STARRS-1 (PS1; Schlafly et al. 2012; Tonry et al. 2012; Magnier et al. 2013) and ASAS-SN data. Spectroscopy reveals a hot, hydrogen-deficient spectrum with clear double-peaked emission lines, consistent with an accreting double degenerate classification. We use follow-up photometry to constrain the orbital parameters of the system. We find an orbital period of 49.71 min, which places Gaia14aae at the long period extremum of the outbursting AM CVn period distribution. Gaia14aae is dominated by the light from its accreting white dwarf. Assuming an orbital inclination of 90 degrees for the binary system, the contact phases of the white dwarf lead to lower limits of 0.78 M solar and 0.015 M solar on the masses of the accretor and donor respectively and a lower limit on the mass ratio of 0.019. Gaia14aae is only the third eclipsing AM CVn star known, and the first in which the WD is totally eclipsed. Using a helium WD model, we estimate the accretor's effective temperature to be 12900+-200 K. The three out-burst events occurred within 4 months of each other, while no other outburst activity is seen in the previous 8 years of Catalina Real-time Transient Survey (CRTS; Drake et al. 2009), Pan-STARRS-1 and ASAS-SN data. This suggests that these events might be rebrightenings of the first outburst rather than individual events.

We use the integrated polarized radio emission at 1.4 GHz ($\Pi_{\rm 1.4\,GHz}$) from a large sample of AGN (796 sources at redshifts $z<0.7$) to study the large-scale magnetic field properties of radio galaxies in relation to the host galaxy accretion state. We find a fundamental difference in $\Pi_{\rm 1.4\,GHz}$ between radiative-mode AGN (i.e. high-excitation radio galaxies, HERGs, and radio-loud QSOs) and jet-mode AGN (i.e. low-excitation radio galaxies, LERGs). While LERGs can achieve a wide range of $\Pi_{\rm 1.4\,GHz}$ (up to $\sim$$30\%$), the HERGs and radio-loud QSOs are limited to $\Pi_{\rm 1.4\,GHz} \lesssim 15\%$. A difference in $\Pi_{\rm 1.4\,GHz}$ is also seen when the sample is divided at 0.5% of the total Eddington-scaled accretion rate, where the weakly accreting sources can attain higher values of $\Pi_{\rm 1.4\,GHz}$. We do not find any clear evidence that this is driven by intrinsic magnetic field differences of the different radio morphological classes. Instead, we attribute the differences in $\Pi_{\rm 1.4\,GHz}$ to the local environments of the radio sources, in terms of both the ambient gas density and the magnetoionic properties of this gas. Thus, not only are different large-scale gaseous environments potentially responsible for the different accretion states of HERGs and LERGs, we argue that the large-scale magnetised environments may also be important for the formation of powerful AGN jets. Upcoming high angular resolution and broadband radio polarization surveys will provide the high precision Faraday rotation measure and depolarization data required to robustly test this claim.

We present a 324.5MHz image of the NOAO Bo\"otes field that was made using Very Large Array (VLA) P-band observations. The image has a resolution of 5.6x5.1arcsec, a radius of $2.05^\circ$ and a central noise of ~0.2mJy\beam. Both the resolution and noise of the image are an order of magnitude better than what was previously available at this frequency and will serve as a valuable addition to the already extensive multiwavelength data that are available for this field. The final source catalogue contains 1370 sources and has a median 325 to 1400MHz spectral index of -0.72. Using a radio colour-colour diagram of the unresolved sources in our catalogue, we identify 33 megahertz peaked-spectrum (MPS) sources. Based on the turnover frequency linear size relation for the gigahertz peaked-spectrum (GPS) and compact steep-spectrum (CSS) sources, we expect that the MPS sources that are compact on scales of tens of milliarcseconds should be young radio loud active galactic nuclei at high (z>2) redshifts. Of the 33 MPS sources, we were able to determine redshifts for 24, with an average redshift of 1.3. Given that five of the sources are at z>2, that the four faint sources for which we could not find redshifts are likely at even higher redshifts and that we could only select sources that are compact on a scale of ~5arcsec, there is encouraging evidence that the MPS method can be used to search for high-redshift sources.

Observational consequences of the tidal disruption of stars by supermassive black holes (SMBHs) can enable us to discover quiescent SMBHs and constrain their mass function. Moreover, observing jetted TDEs (from previously non-active galaxies) provides us with a new means of studying the early phases of jet formation and evolution in an otherwise "pristine" environment. Although several (tens) TDEs have been discovered since 1999, only two jetted TDEs have been recently discovered in hard X-rays, and only one, Swift J1644+57, has a precise localization which further supports the TDE interpretation. These events alone are not sufficient to address those science issues, which require a substantial increase of the current sample. Despite the way they were discovered, the highest discovery potential for {\em jetted} TDEs is not held by current and up-coming X-ray instruments, which will yield only a few to a few tens events per year. In fact, the best strategy is to use the Square Kilometer Array to detect TDEs and trigger multi-wavelength follow-ups, yielding hundreds candidates per year, up to $z \sim 2$. Radio and X-ray synergy, however, can in principle constrain important quantities such as the absolute rate of jetted TDEs, their jet power, bulk Lorentz factor, the black hole mass function, and perhaps discover massive black holes (MBH) with $<10^{5} M_{\odot}$. Finally, when comparing SKA results with information from optical surveys like LSST, one can more directly constrain the efficiency of jet production.

Current wide-area radio surveys are dominated by active galactic nuclei, yet many of these sources have no identified optical counterparts. Here we investigate whether one can constrain the nature and properties of these sources, using Fanaroff-Riley type II (FRII) radio galaxies as probes. These sources are easy to identify since the angular separation of their lobes remains almost constant at some tens of arcseconds for z>1. Using a simple algorithm applied to the FIRST survey, we obtain the largest FRII sample to date, containing over ten thousand double-lobed sources. A subset of 459 sources is matched to SDSS quasars. This sample yields a statistically meaningful description of the fraction of quasars with lobes as a function of redshift and luminosity. This relation is combined with the bolometric quasar luminosity function, as derived from surveys at IR to hard X-ray frequencies, and a disc-lobe correlation to obtain a robust prediction for the density of FRIIs on the radio sky. We find that the observed density can be explained by the population of known quasars, implying that the majority of powerful jets originate from a radiatively efficient accretion flow with a linear jet-disc coupling. Finally, we show that high-redshift jets are more often quenched within 100 kpc, suggesting a higher efficiency of jet-induced feedback into their host galaxies.

We present LOFAR Low Band observations of the Bootes and 3C295 fields. Our images made at 34, 46, and 62 MHz reach noise levels of 12, 8, and 5 mJy beam$^{-1}$, making them the deepest images ever obtained in this frequency range. In total, we detect between 300 and 400 sources in each of these images, covering an area of 17 to 52 deg$^{2}$. From the observations we derive Euclidean-normalized differential source counts. The 62 MHz source counts agree with previous GMRT 153 MHz and VLA 74 MHz differential source counts, scaling with a spectral index of $-0.7$. We find that a spectral index scaling of $-0.5$ is required to match up the LOFAR 34 MHz source counts. This result is also in agreement with source counts from the 38 MHz 8C survey, indicating that the average spectral index of radio sources flattens towards lower frequencies. We also find evidence for spectral flattening using the individual flux measurements of sources between 34 and 1400 MHz and by calculating the spectral index averaged over the source population. To select ultra-steep spectrum ($\alpha < -1.1$) radio sources, that could be associated with massive high redshift radio galaxies, we compute spectral indices between 62 MHz, 153 MHz and 1.4 GHz for sources in the Bo\"otes field. We cross-correlate these radio sources with optical and infrared catalogues and fit the spectral energy distribution to obtain photometric redshifts. We find that most of these ultra-steep spectrum sources are located in the $ 0.7 \lesssim z \lesssim 2.5$ range.

Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage. This can be achieved by combining data from observatories located in both the northern and southern hemispheres. To this end, a joint analysis using data recorded at the Telescope Array and the Pierre Auger Observatory above $10^{19}$ eV is presented in this work. The resulting multipolar expansion of the flux of cosmic rays allows us to perform a series of anisotropy searches, and in particular to report on the angular power spectrum of cosmic rays above $10^{19}$ eV. No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole.

We report an observational estimate of the rate of stellar tidal disruption flares (TDFs) in inactive galaxies, based on a successful search for these events among transients in galaxies using archival SDSS multi-epoch imaging data (Stripe 82). This search yielded 186 nuclear flares in galaxies, of which two are excellent TDF candidates. Because of the systematic nature of the search, the very large number of galaxies, the long time of observation, and the fact that non-TDFs were excluded without resorting to assumptions about TDF characteristics, this study provides an unparalleled opportunity to measure the TDF rate. To compute the rate of optical stellar tidal disruption events, we simulate our entire pipeline to obtain the efficiency of detection. The rate depends on the light curves of TDFs, which are presently still poorly constrained. Using only the observed part of the SDSS light curves gives a model-independent upper limit to the optical TDF rate: < 2 10^-4 per year per galaxy (90% CL). We develop three empirical models of the light curves, based on the two SDSS light curves and two more recent and better-sampled Pan-STARRS TDF light curves, leading to our best-estimate of the rate: (1.5 - 2.0)_{-1.3}^{+2.7} 10^-5 per year per galaxy. We explore the modeling uncertainties by considering two theoretically motivated light curve models, as well as two different relationships between black hole mass and galaxy luminosity, and two different treatments of the cutoff in the visibility of TDFs at large black hole mass. From this we conclude that these sources of uncertainty are not significantly larger than the statistical ones. Our results are applicable for galaxies hosting black holes with mass in the range of few million to 10^8 solar masses, and translates to a volumetric TDF rate of (4 - 8) 10^-8 per year per cubic Mpc.

We describe the method devised to reconstruct inclined cosmic-ray air showers with zenith angles greater than $60^\circ$ detected with the surface array of the Pierre Auger Observatory. The measured signals at the ground level are fitted to muon density distributions predicted with atmospheric cascade models to obtain the relative shower size as an overall normalization parameter. The method is evaluated using simulated showers to test its performance. The energy of the cosmic rays is calibrated using a sub-sample of events reconstructed with both the fluorescence and surface array techniques. The reconstruction method described here provides the basis of complementary analyses including an independent measurement of the energy spectrum of ultra-high energy cosmic rays using very inclined events collected by the Pierre Auger Observatory.

A flux of neutrons from an astrophysical source in the Galaxy can be detected in the Pierre Auger Observatory as an excess of cosmic-ray air showers arriving from the direction of the source. To avoid the statistical penalty for making many trials, classes of objects are tested in combinations as nine "target sets", in addition to the search for a neutron flux from the Galactic Center or from the Galactic Plane. Within a target set, each candidate source is weighted in proportion to its electromagnetic flux, its exposure to the Auger Observatory, and its flux attenuation factor due to neutron decay. These searches do not find evidence for a neutron flux from any class of candidate sources. Tabulated results give the combined p-value for each class, with and without the weights, and also the flux upper limit for the most significant candidate source within each class. These limits on fluxes of neutrons significantly constrain models of EeV proton emission from non-transient discrete sources in the Galaxy.