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.

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.

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.

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.

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.

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.

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.

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.

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 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.

The spin of supermassive black holes could power jets from active galactic nuclei (AGN), although direct observational evidence for this conjecture is sparse. The accretion disk luminosity and jet power, on the other hand, have long been observed to follow a linear correlation. If jet power is coupled to black hole spin, deviations from the jet-disk correlation for a sample of AGN can be used to probe the dispersion of the spin parameter within this sample. To exploit this idea, we matched double-lobed radio sources from the FIRST survey to spectroscopically confirmed quasars from SDSS. We obtain 763 FR-II quasars with a median redshift of z=1.2 A tight correlation between the optical luminosity of the accretion disk and the lobe radio luminosity is observed. We estimate that 5-20% of the bolometric disk luminosity is converted to jet power. Most of the scatter to the optical-radio correlation is due to environment; deviations from jet-disk coupling due to internal factors (e.g., spin) contribute at most 0.2 dex. Under the assumption that the Blandford-Znajek mechanism operates in AGN, we obtain an upper limit of 0.1 dex to the dispersion of the product of the spin and the magnetic flux threading the horizon. If black hole spin determines the AGN jet efficiency, then our observations imply that all FR-II quasars have very similar spin. In addition, the quasar spin distribution needs to have a wide gap to explain the radio-quiet population. The alternative, and perhaps more likely, interpretation of the tight jet-disk correlation is that black hole spin is not relevant for powering AGN jets.

Contributions of the Pierre Auger Collaboration to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013

Many decades of observations of active galactic nuclei and X-ray binaries have shown that relativistic jets are ubiquitous when compact objects accrete. One could therefore anticipate the launch of a jet after a star is disrupted and accreted by a massive black hole. This birth of a relativistic jet may have been observed recently in two stellar tidal disruption flares (TDFs), which were discovered in gamma-rays by Swift. Yet no transient radio emission has been detected from the tens of TDF candidates that were discovered at optical to soft X-ray frequencies. Because the sample that was followed-up at radio frequencies is small, the non-detections can be explained by Doppler boosting, which reduces the jet flux for off-axis observers. And since the existing follow-up observation are mostly within ~10 months of the discovery, the non-detections can also be due to a delay of the radio emission with respect to the time of disruption. To test the conjecture that all TDFs launch jets, we obtained 5 GHz follow-up observations with the Jansky VLA of seven known TDFs. To avoid missing delayed jet emission, our observations probe 1-8 years since the estimated time of disruption. None of the sources are detected, with very deep upper limits at the 10 micro Jansky level. These observations rule out the hypothesis that these TDFs launched jets similar to radio-loud quasars. We also constrain the possibility that the flares hosted a jet identical to Sw 1644+57, the first and best-sampled relativistic TDF. We thus obtain evidence for a dichotomy in the stellar tidal disruption population, implying that the jet launching mechanism is sensitive to the parameters of the disruption.

We report an observational estimate of the rate of stellar tidal disruption flares (TDFs), based on our (successful) search for these events in archival Sloan Digital Sky Survey (SDSS) multi-epoch imaging data. Our pipeline took advantage of the excellent astrometry of SDSS to separate nuclear flares from supernovae. The 10 year baseline and the high cadence of the observations facilitate a clear-cut identification of variable active galactic nuclei. We found 186 nuclear flares, of which two are excellent TDF candidates. To compute the rate of (optical) stellar tidal disruption events, we simulate our entire pipeline to obtain the efficiency of detection for a given light curve. Using a simple model to extrapolate the observed light curves forward and backward in time, we find our best-estimate of the TDF rate: 3x10-5 per galaxy per year. In addition, we give a model-independent upper limit to the TDF rate: < 3x10-4 per galaxy per year (90% CL).

To understand the feedback of black holes on their environment or the acceleration of ultra-high energy cosmic rays in the present cosmic epoch, a systematic, all-sky inventory of radio galaxies in the local universe is needed. Here we present the first catalog of radio-emitting galaxies that meets this requirement. Our catalog allows the selection of volume-limited subsamples containing all low-power radio galaxies, similar to the prototypical low-power radio galaxies Cen A or M87, within some hundred Mpc. It is constructed by matching radio emission from the NVSS and SUMSS surveys to galaxies of the 2MASS Redshift Survey (2MRS) using an image-level algorithm that properly treats the extended structure of radio sources. The sample contains 575 radio-emitting galaxies with a flux greater than 213 mJy at 1.4 GHz. Over 30% of the galaxies in our catalog are not contained in existing large-area extra-galactic radio samples. We compute the optical and radio luminosity functions and the fraction of radio galaxies as a function of galaxy luminosity. We find that the galaxy density around radio galaxies is significantly higher than around non-radio galaxies of the same luminosity and morphology. This enhanced clustering suggests a causal relation between external galaxy properties, such as environment or merger history, and the formation of powerful jets in the present universe. Since the enhancement is observed with respect to galaxies of the same luminosity and Hubble type, it is not primarily driven by black hole mass. Our automated matching procedure is found to select radio-emitting galaxies with high efficiency (99%) and purity (91%), which is key for future processing of deeper, larger samples.

Using archival SDSS multi-epoch imaging data (Stripe 82), we have searched for the tidal disruption of stars by super-massive black holes in non-active galaxies. Two candidate tidal disruption events (TDEs) are identified. They have optical black-body temperatures 2 10^4 K and observed peak luminosities M_g=-18.3 and -20.4; their cooling rates are very low, qualitatively consistent with expectations for tidal disruption flares. Their properties are examined using i) SDSS imaging to compare them to other flares observed in the search, ii) UV emission measured by GALEX and iii) spectra of the hosts and of one of the flares. Our pipeline excludes optically identifiable AGN hosts, and our variability monitoring over 9 years provides strong evidence that these are not flares in hidden AGNs. The spectra and color evolution of the flares are unlike any SN observed to date, their strong late-time UV emission is particularly distinctive, and they are nuclear at high resolution, arguing against their being first cases of a previously-unobserved class of SNe or more extreme examples of known SN types. Taken together, the observed properties are difficult to reconcile with a SN or AGN-flare explanation, although an entirely new process specific to the inner few-hundred parsecs of non-active galaxies cannot be excluded. Our observed rate and method show the feasibility of obtaining a candidate TDE sample of hundreds of events and O(1) purity, using geometric resolution and host and flare color alone. A by-product of this work is quantification of the power-spectrum of extreme flares in AGNs.

A star that passes too close to a massive black hole will be torn apart by tidal forces. The flare of photons emitted during the accretion of the stellar debris is predicted to be observable and candidates of such events have been observed at optical to X-ray frequencies. If a fraction of the accreted material is fed into a jet, tidal flares should be detectable at radio frequencies too, thus comprising a new class of rare radio transients. Using the well-established scaling between accretion power and jet luminosity and basic synchrotron theory, we construct an empirically-rooted model to predict the jet luminosity for a time-dependent accretion rate. We apply this model to stellar tidal disruptions and predict the snapshot rate of these events. For a small angle between the observer and the jet, our model reproduces the observed radio flux of the tidal flare candidate GRB 110328A. We find that future radio surveys will be able to test whether the majority of tidal disruptions are accompanied by a jet.