• ### Stellar Tidal Disruption Events in General Relativity(1801.10180)

Jan. 30, 2018 gr-qc, astro-ph.HE
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.
• ### On the mass and luminosity functions of tidal disruption flares: rate suppression due to black hole event horizons(1707.03458)

Jan. 10, 2018 astro-ph.GA, astro-ph.HE
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.
• ### Black hole masses of tidal disruption event host galaxies(1706.08965)

June 27, 2017 astro-ph.GA, astro-ph.HE
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}$$\leq22 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\times10^5 M_{\odot}$$\leq$M$_{\rm BH}$$\leq2\times10^7 M_{\odot}. The TDE host galaxy sample is dominated by low mass black holes (\sim10^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}$$\geq10^{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.
• ### Optical/UV-to-X-Ray Echoes from the Tidal Disruption Flare ASASSN-14li(1703.07024)

March 21, 2017 astro-ph.HE
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.
• ### An enhanced rate of tidal disruptions in the centrally overdense E+A galaxy NGC 3156(1604.02056)

Dec. 29, 2016 astro-ph.GA, astro-ph.HE
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.
• ### Discovery of transient infrared emission from dust heated by stellar tidal disruption flares(1605.04304)

Sept. 13, 2016 astro-ph.GA, astro-ph.HE
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.
• ### Sensitive Radio Survey of Obscured Quasar Candidates(1603.07325)

March 23, 2016 astro-ph.GA
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.
• ### A radio jet from the optical and x-ray bright stellar tidal disruption flare ASASSN-14li(1511.08803)

Jan. 15, 2016 astro-ph.GA, astro-ph.HE
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.
• ### Gaia transient detection efficiency: hunting for nuclear transients(1510.07767)

Oct. 27, 2015 astro-ph.GA, astro-ph.IM
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.
• ### Nature and evolution of powerful radio galaxies and their link with the quasar luminosity function(1411.2968)

Nov. 11, 2014 astro-ph.GA, astro-ph.HE
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.
• ### Measurement of the rate of stellar tidal disruption flares(1407.6425)

Aug. 26, 2014 astro-ph.GA, astro-ph.HE
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 contribution of spin to jet-disk coupling in black holes(1308.1437)

Aug. 6, 2013 astro-ph.CO, astro-ph.HE
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
• ### Constraints on off-axis jets from stellar tidal disruption flares(1210.0022)

March 4, 2013 astro-ph.CO, astro-ph.HE
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.
• ### The rate of stellar tidal disruption flares from SDSS data(1210.1246)

Oct. 3, 2012 astro-ph.CO
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).
• ### Radio galaxies of the local universe: all-sky catalog, luminosity functions, and clustering(1206.0031)

July 20, 2012 astro-ph.CO