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.

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 a systematic search for periodically varying quasar and supermassive black hole binary (SMBHB) candidates in the Pan-STARRS1 Medium Deep Survey. From $\sim9,000$ color-selected quasars in a $\sim50$ deg$^{2}$ sky area, we initially identify $26$ candidates with more than $1.5$ cycles of variation. We extend the baseline of observations via our imaging campaign with the Discovery Channel Telescope and the Las Cumbres Observatory network and reevaluate the candidates using a more rigorous, maximum likelihood method. Using a range of statistical criteria and assuming the Damped Random Walk model for normal quasar variability, we identify one statistically significant periodic candidate. We also investigate the capabilities of detecting SMBHBs by the Large Synoptic Survey Telescope using our study with MDS as a benchmark and explore any complementary, multiwavelength evidence for SMBHBs in our sample.

Time domain and multi-messenger astrophysics are growing and important modes of observational astronomy that will help define astrophysics in the 2020s. Significant effort is being put into developing the components of a follow-up system for dynamically turning survey alerts into data. This system consists of: 1) brokers that will aggregate, classify, and filter alerts; 2) Target Observation Managers (TOMs) for prioritizing targets and managing observations and data; and 3) observatory interfaces, schedulers, and facilities along with data reduction software and science archives. These efforts need continued community support and funding in order to complete and maintain them. Many of the efforts can be community open-source software projects but they will benefit from the leadership of professional software developers. The coordination should be done by institutions that are involved in the follow-up system such as the national observatories (e.g. LSST/Gemini/NOAO Mid-scale/Community Science and Data Center) or a new MMA institute. These tools will help the community to produce the most science from new facilities and will provide new capabilities for all users of the facilities that adopt them.

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 analyze extensive spectroscopic and photometric data of the hypervariable quasar SDSS J131424+530527 (RMID 017) at z=0.456, an optical "changing look" quasar from the Sloan Digital Sky Survey Reverberation Mapping project that increased in optical luminosity by a factor of 10 between 2014 and 2017. The observed broad emission lines all respond in luminosity and width to the changing optical continuum, as expected for photoionization in a stratified, virialized broad emission line region. The luminosity changes therefore result from intrinsic changes in accretion power rather than variable obscuration. The variability is continuous and apparently stochastic, disfavoring an origin as a discrete event such as a tidal disruption flare or microlensing event. It is coordinated on day timescales with blue leading red, consistent with reprocessing powering the entire optical SED. We show that this process cannot work in a standard thin disk geometry on energetic grounds, and would instead require a large covering factor reprocessor. Disk instability models could potentially also explain the data, provided that the instability sets in near the inner radius of a geometrically thick accretion disk.

Time-domain Astrophysics (TDA), a foundation of Astronomy, has become a major part of current and projected (2020s) astrophysics. While much has been derived from temporal measures of flux and color, the real physics comes from spectroscopy. With LSST coming on line in 2022, with TDA as one of its original drivers, the deluge of Transients and new types of variables will be truly astronomical. With multi-wavelength targeted EM surveys and multi-messenger (e.g. LIGO-international and advanced ICECUBE), and the possibility of full-sky and full-time X-ray imagers, the discovery of new Transients and Variables will flood telescopes on the ground and in space, and this just for multi-band imaging without spectroscopy. In this White Paper we briefly summarize several long-standing major science objectives that can be realized with TDA imaging and spectroscopy (near UV to mid IR) from space. We provide a brief description of how these can be achieved with the Time-domain Spectroscopic Observatory (TSO), a Probe-Class mission concept that ELTs on the ground and Flagship missions in space can not achieve on their own.

Long gamma-ray bursts (GRBs) are the most luminous known electromagnetic radiation sources in the Universe for the 3 to 300 sec of their prompt flashes (isotropic X/ gamma-ray luminosities up to 10^53 ergs/sec). Their afterglows have first day rest-frame UV/optical absolute magnitudes AB = -30 to -23. This luminous continuum nUV-nIR back-light provides the ultimate probe of the SFR(z) back to the first Pop III to II.5 stars, expected to be massive and GRB progenitors. GRB afterglow spectra in the first 1 to 3 hours will directly measure their host galaxy ionization fraction x_i vs. z in the Epoch of Reionization (EOR), tracing the growth of structure. Only 28% of Swift GRBs have measured redshifts due to limited followup at R, J >21. Some ~25% of GRBs are optically dark due to dust absorption in their host galaxies, but those with low NH in their X-ray spectra are likely at z >7. Current 8-10m telescopes and coming ELTs cannot pursue optically dark GRBs promptly, nor can JWST or WFIRST slew within 0.5 to 1 days of a GRB. The Time-domain Spectroscopic Observatory (TSO) is a proposed Probe-class 1.3m telescope at L2, with imaging and spectroscopy (R = 200, 1800) in 4 bands (0.3 to 5um) and rapid slew capability to 90% of sky. TSO would finally utilize z > 6 to 12 GRBs as the most direct probe of the SFR(z), EOR(z), and possibly the first direct detection of the core collapse of the very first (Pop III) stars.

With the epochal first detection of gravitational waves from a binary neutron star (NS) merger with the GW170817 event, and its direct confirmation that NS-NS mergers are significant sources of the of the r-process nucleosynthesis of heavy elements, an immense new arena for prompt EM (X-rays through IR and radio) studies of fundamental physics has been opened. Over the next decade, GW observatories will expand in scale and sensitivity so the need for facilities that can provide prompt, high sensitivity, broad-band EM followup becomes more urgent. NS-NS or NS-black hole (BH) mergers will be instantly recognized (and announced) by the LIGO-international collaboration. LSST will be a prime resource for rapid tiling of what will usually be large (~10-100 degree squared) error boxes. X-ray through IR Telescopes in space with (nearly) full-sky access that can rapidly image and tile are crucial for providing the earliest imaging and spectroscopic studies of the kilonova emission immediately following NS-NS mergers. The Time-domain Spectroscopic Observatory (TSO) is a proposed Probe-class 1.3 m telescope at L2, with imaging and spectroscopy (R = 200, 1800) in 4 bands (0.3 - 5 micron) and rapid slew capability to 90% of sky. TSO nUV-mid-IR spectra will enable new constraints on NS structure and nucleosynthesis.

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.

The ubiquitous variability of quasars across a wide range of wavelengths and timescales encodes critical information about the structure and dynamics of the circumnuclear emitting regions that are too small to be directly resolved, as well as detailed underlying physics of accretion and feedback processes in these active supermassive black holes. We emphasize the importance of studying quasar variability with time-domain spectroscopy, focusing on two science cases: (1) reverberation mapping (RM) to measure the broad-line region sizes and black hole masses in distant quasars; (2) spectroscopic follow-up of extreme variability quasars that dramatically change their continuum and broad-line flux within several years. We highlight the need for dedicated optical-infrared spectroscopic survey facilities in the coming decades to accompany wide-area time-domain imaging surveys, including: (1) the next phase of the Sloan Digital Sky Survey (SDSS-V; ~2020-2025), an all-sky, time-domain multi-object spectroscopic survey with 2.5m-class telescopes; (2) the planned Maunakea Spectroscopic Explorer, a dedicated 10m-class spectroscopic survey telescope with a 1.5 sq. deg field-of-view and multiplex of thousands of fibers in both optical and near-IR (J+H) to begin operations in 2029; (3) the Time-domain Spectroscopic Observatory (TSO), a proposed Probe-class ~1.3m telescope at L2, with imaging and spectroscopy (R=200, 1800) in 4 bands (0.3 - 5 micron) and rapid slew capability to 90% of sky, which will extend the coverage of Hbeta to z=8.

The tidal disruption of stars by (super-)massive black holes in galactic nuclei has been discussed in theoretical terms for about 30 years but only in the past decade have we been able to detect such events in substantial numbers. Thus, we are now starting to carry out observational tests of models for the disruption. We are also formulating expectations for the inspiral and disruption of white dwarfs by intermediate-mass black holes with masses $< 10^5\;{\rm M}_\odot$. Such events are very rich with information and open a new window to intermediate-mass black holes, thought to live in dwarf galaxies and star clusters. They can inform us of the demographics of intermediate-mass black holes, stellar populations and dynamics in their immediate vicinity, and the physics of accretion of hydrogen-deficient material. The combination of upcoming transient surveys using ground-based, electromagnetic observatories and low-frequency gravitational wave observations is ideal for exploiting tidal disruptions of white dwarfs. The detection rate of gravitational wave signals, optimistically, may reach a few dozen per year in a volume up to $z\approx 0.1$. Gravitational wave observations are particularly useful because they yield the masses of the objects involved and allow determination of the spin of the black hole, affording tests of physical models for black hole formation and growth. They also give us advance warning of the electromagnetic flares by weeks or more. The right computing infrastructure for modern models for the disruption process and event rates will allow us to make the most of the upcoming observing facilities.

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

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.

Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ \Delta ra, \Delta dec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.

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.

Wide-field optical surveys have begun to uncover large samples of fast (t_rise < 5d), luminous (M_peak < -18), blue transients. While commonly attributed to the breakout of a supernova shock into a dense wind, the great distances to the transients of this class found so far have hampered detailed investigation of their properties. We present photometry and spectroscopy from a comprehensive worldwide campaign to observe AT2018cow (ATLAS18qqn), the first fast-luminous optical transient to be found in real time at low redshift. Our first spectra (<2 days after discovery) are entirely featureless. A very broad absorption feature suggestive of near-relativistic velocities develops between 3-8 days, then disappears. Broad emission features of H and He develop after >10 days. The spectrum remains extremely hot throughout its evolution, and the photospheric radius contracts with time (receding below R<10^14 cm after 1 month). This behaviour does not match that of any known supernova, although a relativistic jet within a fallback supernova could explain some of the observed features. Alternatively, the transient could originate from the disruption of a star by an intermediate-mass black hole, although this would require long-lasting emission of highly super-Eddington thermal radiation. In either case, AT2018cow suggests that the population of fast luminous transients represents a new class of astrophysical event. Intensive follow-up of this event in its late phases, and of any future events found at comparable distance, will be essential to better constrain their origins.

This white paper specifies the footprints, cadence requirements, and total-depth requirements needed to allow the most-successful AGN studies in the four currently selected LSST Deep-Drilling Fields (DDFs): ELAIS-S1, XMM-LSS, CDF-S, and COSMOS. The information provided on cadence and total-depth requirements will also likely be applicable to enabling effective AGN science in any additional DDFs that are chosen.

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 past year has witnessed discovery of the first identified counterparts to a gravitational wave transient (GW 170817A) and a very high-energy neutrino (IceCube-170922A). These source identifications, and ensuing detailed studies, have realized longstanding dreams of astronomers and physicists to routinely carry out observations of cosmic sources by other than electromagnetic means, and inaugurated the era of "multi-messenger" astronomy. While this new era promises extraordinary physical insights into the universe, it brings with it new challenges, including: highly heterogeneous, high-volume, high-velocity datasets; globe-spanning cross-disciplinary teams of researchers, regularly brought together into transient collaborations; an extraordinary breadth and depth of domain-specific knowledge and computing resources required to anticipate, model, and interpret observations; and the routine need for adaptive, distributed, rapid-response observing campaigns to fully exploit the scientific potential of each source. We argue, therefore, that the time is ripe for the community to conceive and propose an Institute for Multi-Messenger Astrophysics that would coordinate its resources in a sustained and strategic fashion to efficiently address these challenges, while simultaneously serving as a center for education and key supporting activities. In this fashion, we can prepare now to realize the bright future that we see, beyond, through these newly opened windows onto the universe.

Graham et al. (2015a) reported a periodically varying quasar and supermassive black hole binary candidate, PG1302-102 (hereafter PG1302), which was discovered in the Catalina Real-Time Transient Survey (CRTS). Its combined Lincoln Near-Earth Asteroid Research (LINEAR) and CRTS optical light curve is well fitted to a sinusoid of an observed period of $\approx 1,884$ days and well modeled by the relativistic Doppler boosting of the secondary mini-disk (D'Orazio et al. 2015). However, the LINEAR+CRTS light curve from MJD $\approx 52700$ to MJD $\approx 56400$ covers only $\sim 2$ cycles of periodic variation, which is a short baseline that can be highly susceptible to normal, stochastic quasar variability (Vaughan et al. 2016). In this Letter, we present a re-analysis of PG1302, using the latest light curve from the All-Sky Automated Survey for Supernovae (ASAS-SN), which extends the observational baseline to the present day (MJD $\approx 58200$), and adopting a maximum likelihood method which searches for a periodic component in addition to stochastic quasar variability. When the ASAS-SN data are combined with the previous LINEAR+CRTS data, the evidence for periodicity decreases. For genuine periodicity one would expect that additional data would strengthen the evidence, so the decrease in significance may be an indication that the binary model is disfavored.

We performed a wide-area (2000 deg$^{2}$) g and I band experiment as part of a two month extension to the Intermediate Palomar Transient Factory. We discovered 36 extragalactic transients including iPTF17lf, a highly reddened local SN Ia, iPTF17bkj, a new member of the rare class of transitional Ibn/IIn supernovae, and iPTF17be, a candidate luminous blue variable outburst. We do not detect any luminous red novae and place an upper limit on their rate. We show that adding a slow-cadence I band component to upcoming surveys such as the Zwicky Transient Facility will improve the photometric selection of cool and dusty transients.

We present late-time observations by Swift and XMM-Newton of the tidal disruption event (TDE) ASASSN-15oi that reveal that the source brightened in the X-rays by a factor of $\sim10$ one year after its discovery, while it faded in the UV/optical by a factor of $\sim 100$. The XMM-Newton observations measure a soft X-ray blackbody component with $kT_{\rm bb} \sim 45$ eV, corresponding to radiation from several gravitational radii of a central $\sim 10^6 M_\odot$ black hole. The last Swift epoch taken almost 600 days after discovery shows that the X-ray source has faded back to its levels during the UV/optical peak. The timescale of the X-ray brightening suggests that the X-ray emission could be coming from delayed accretion through a newly forming debris disk, and that the prompt UV/optical emission is from the prior circularization of the disk through stream-stream collisions. The lack of spectral evolution during the X-ray brightening disfavors ionization breakout of a TDE "veiled" by obscuring material. This is the first time a TDE has been shown to have a delayed peak in soft X-rays relative to the UV/optical peak, which may be the first clear signature of the real-time assembly of a nascent accretion disk, and provides strong evidence for the origin of the UV/optical emission from circularization, as opposed to reprocessed emission of accretion radiation.

We present an optical to X-ray spectral analysis of four 'hypervariable' AGN (HVAs) discovered by comparing Pan-STARRS data to that from SDSS over a 10 year baseline (Lawrence et al 2016). There is some evidence that these objects are X-ray loud for their corresponding UV luminosities, but given that we measured them in a historic high state, it is not clear whether to take the high-state or low-state as typical of the properties of these HVAs. We estimate black hole masses based on MgII and H$\alpha$ emission line profiles, and either the high and low state luminosities, finding mass ranges $\log(M_{\rm BH}/M_{\odot}) = 8.2-8.8$ and $\log(M_{\rm BH}/M_{\odot}) = 7.9-8.3$ respectively. We then fit energy conserving models to the SEDs, obtaining strong constraints on the bolometric luminosity and $\alpha_{\rm OX}$. We compare the SED properties with a larger, X-ray selected AGN sample for both of these scenarios, and observe distinct groupings in spectral shape versus luminosity parameter space. In general, the SED properties are closer to normal if we assume that the low-state is representative. This supports the idea that the large slow outbursts may be due to extrinsic effects (for example microlensing) as opposed to accretion rate changes, but a larger sample of HVAs is needed to be confident of this conclusion.

We present analysis of the first Chandra observation of PSO J334.2028+01.4075 (PSO J334), targeted as a binary-AGN candidate based on periodic variations of the optical flux. With no prior targeted X-ray coverage for PSO J334, our new 40 ksec Chandra observation allows for the opportunity to differentiate between a single or binary-AGN system, and if a binary, can characterize the mode of accretion. Simulations show that the two expected accretion disk morphologies for binary-AGN systems are (i) a "cavity", where the inner region of the accretion disk is mostly empty and emission is truncated blueward of the wavelength associated with the temperature of the innermost ring, or (ii) "minidisks", where there is substantial accretion from the cirumbinary disk onto one or both of the members of the binary, each with their own shock-heated thin-disk accretion system. We find the X-ray emission to be well-fit with an absorbed power-law, incompatible with the simple cavity scenario. Further, we construct an SED of PSO J334 by combining radio through X-ray observations and find that the SED agrees well with that of a normal AGN, most likely incompatible with the minidisk scenario. Other analyses, such as locating the quasar on IR color-color diagrams and analyzing the quasar mass predicted by the fundamental plane of black hole activity, further highlight the similarity of PSO J334 with respect to normal AGN. On the multi-wavelength fronts we investigated, we find no evidence supporting PSO J334 as a binary-AGN system, though our analysis remains insensitive to some binary configurations.

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 present ground-based and \textit{Swift} observations of iPTF16fnl, a likely tidal disruption event (TDE) discovered by the intermediate Palomar Transient Factory (iPTF) survey at 66.6 Mpc. The lightcurve of the object peaked at absolute $M_g=-17.2$ mag. The maximum bolometric luminosity (from optical and UV) was $L_p~\simeq~(1.0\,\pm\,0.15) \times 10^{43}$ erg/s, an order of magnitude fainter than any other optical TDE discovered so far. The luminosity in the first 60 days is consistent with an exponential decay, with $L \propto e^{-(t-t_0)/\tau}$, where $t_0$=~57631.0 (MJD) and $\tau\simeq 15$ days. The X-ray shows a marginal detection at $L_X=2.4^{1.9}_{-1.1}\times 10^{39}$ erg/s (\textit{Swift} X-ray Telescope). No radio counterpart was detected down to 3$\sigma$, providing upper limits for monochromatic radio luminosity of $\nu L_{\nu} < 2.3\times10^{36}$ erg/s and $\nu L_{\nu}<1.7\times 10^{37}$ erg/s (VLA, 6.1 and 22 GHz). The blackbody temperature, obtained from combined \textit{Swift} UV and optical photometry, shows a constant value of 19,000 K. The transient spectrum at peak is characterized by broad He II and H$\alpha$ emission lines, with an FWHM of about 14,000 km/s and 10,000 km/s respectively. He I lines are also detected at $\lambda\lambda$ 5875 and 6678. The spectrum of the host is dominated by strong Balmer absorption lines, which are consistent with a post-starburst (E+A) galaxy with an age of $\sim$650 Myr and solar metallicity. The characteristics of iPTF16fnl make it an outlier on both luminosity and decay timescales, as compared to other optically selected TDEs. The discovery of such a faint optical event suggests a higher rate of tidal disruptions, as low luminosity events may have gone unnoticed in previous searches.

We report the discovery by the intermediate Palomar Transient Factory (iPTF) of a candidate tidal disruption event (TDE) iPTF16axa at $z=0.108$, and present its broadband photometric and spectroscopic evolution from 3 months of follow-up observations with ground-based telescopes and Swift. The light curve is well fitted with a $t^{-5/3}$ decay, and we constrain the rise-time to peak to be $<$49 rest-frame days after disruption, which is roughly consistent with the fallback timescale expected for the $\sim 5\times$10$^{6}$ $M_\odot$ black hole inferred from the stellar velocity dispersion of the host galaxy. The UV and optical spectral energy distribution (SED) is well described by a constant blackbody temperature of T$\sim$ 3$\times$10$^4$ K over the monitoring period, with an observed peak luminosity of 1.1$\times$10$^{44}$ erg s$^{-1}$. The optical spectra are characterized by a strong blue continuum and broad HeII and H$\alpha$ lines characteristic of TDEs. We compare the photometric and spectroscopic signatures of iPTF16axa with 11 TDE candidates in the literature with well-sampled optical light curves. Based on a single-temperature fit to the optical and near-UV photometry, most of these TDE candidates have peak luminosities confined between log(L [erg s$^{-1}$]) = 43.4-44.4, with constant temperatures of a few $\times 10^{4}$ K during their power-law declines, implying blackbody radii on the order of ten times the tidal disruption radius, that decrease monotonically with time. For TDE candidates with hydrogen and helium emission, the high helium-to-hydrogen ratios suggest that the emission arises from high-density gas, where nebular arguments break down. We find no correlation between the peak luminosity and the black hole mass, contrary to the expectations for TDEs to have $\dot{M} \propto M_{\rm BH}^{-1/2}$.

We analyze four extreme active galactic nuclei (AGN) transients to explore the possibility that they are caused by rare, high-amplitude microlensing events. These previously unknown type-I AGN are located in the redshift range 0.6-1.1 and show changes of > 1.5 mag in the g-band on a time-scale of ~years. Multi-epoch optical spectroscopy, from the William Herschel Telescope, shows clear differential variability in the broad line fluxes with respect to the continuum changes and also evolution in the line profiles. In two cases, a simple point-source, point-lens microlensing model provides an excellent match to the long-term variability seen in these objects. For both models, the parameter constraints are consistent with the microlensing being due to an intervening stellar mass object but as yet there is no confirmation of the presence of an intervening galaxy. The models predict a peak amplification of 10.3/13.5 and an Einstein time-scale of 7.5/10.8 yr, respectively. In one case, the data also allow constraints on the size of the CIII] emitting region, with some simplifying assumptions, to to be ~1.0-6.5 light-days and a lower limit on the size of the MgII emitting region to be > 9 light-days (half-light radii). This CIII] radius is perhaps surprisingly small. In the remaining two objects, there is spectroscopic evidence for an intervening absorber but the extra structure seen in the light curves requires a more complex lensing scenario to adequately explain.

We present a radio-quiet quasar at z=0.237 discovered "turning on" by the intermediate Palomar Transient Factory (iPTF). The transient, iPTF 16bco, was detected by iPTF in the nucleus of a galaxy with an archival SDSS spectrum with weak narrow-line emission characteristic of a low-ionization emission line region (LINER). Our follow-up spectra show the dramatic appearance of broad Balmer lines and a power-law continuum characteristic of a luminous (L_bol~10^45 erg/s) type 1 quasar 12 years later. Our photometric monitoring with PTF from 2009-2012, and serendipitous X-ray observations from the XMM-Newton Slew Survey in 2011 and 2015, constrain the change of state to have occurred less than 500 days before the iPTF detection. An enhanced broad Halpha to [OIII]5007 line ratio in the type 1 state relative to other changing-look quasars also is suggestive of the most rapid change of state yet observed in a quasar. We argue that the >10 increase in Eddington ratio inferred from the brightening in UV and X-ray continuum flux is more likely due to an intrinsic change in the accretion rate of a pre-existing accretion disk, than an external mechanism such as variable obscuration, microlensing, or the tidal disruption of a star. However, further monitoring will be helpful in better constraining the mechanism driving this change of state. The rapid "turn on" of the quasar is much shorter than the viscous infall timescale of an accretion disk, and requires a disk instability that can develop around a ~10^8 M_sun black hole on timescales less than a year.

We present a systematic search for periodically varying quasars and supermassive black hole binary (SMBHB) candidates in the Pan-STARRS1 (PS1) Medium Deep Survey's MD09 field. From a color-selected sample of 670 quasars extracted from a multi-band deep-stack catalog of point sources, we locally select variable quasars and look for coherent periods with the Lomb-Scargle periodogram. Three candidates from our sample demonstrate strong variability for more than ~3 cycles, and their PS1 light curves are well fitted to sinusoidal functions. We test the persistence of the candidates' apparent periodic variations detected during the 4.2 years of the PS1 survey with archival photometric data from the SDSS Stripe 82 survey or new monitoring with the Large Monolithic Imager at the Discovery Channel Telescope. None of the three periodic candidates (including PSO J334.2028+1.4075) remain persistent over the extended baseline of 7 - 14 years, corresponding to a detection rate of < 1 in 670 quasars in a search area of 5 deg^2. Even though SMBHBs should be a common product of the hierarchal growth of galaxies, and periodic variability in SMBHBs has been theoretically predicted, a systematic search for such signatures in a large optical survey is strongly limited by its temporal baseline and the "red noise" associated with normal quasar variability. We show that follow-up long-term monitoring (>5 cycles) is crucial to our search for these systems.

Super-luminous supernovae that radiate more than 10^44 ergs per second at their peak luminosity have recently been discovered in faint galaxies at redshifts of 0.1-4. Some evolve slowly, resembling models of 'pair-instability' supernovae. Such models involve stars with original masses 140-260 times that of the Sun that now have carbon-oxygen cores of 65-30 solar masses. In these stars, the photons that prevent gravitational collapse are converted to electron-positron pairs, causing rapid contraction and thermonuclear explosions. Many solar masses of 56Ni are synthesized; this isotope decays to 56Fe via 56Co, powering bright light curves. Such massive progenitors are expected to have formed from metal-poor gas in the early Universe. Recently, supernova 2007bi in a galaxy at redshift 0.127 (about 12 billion years after the Big Bang) with a metallicity one-third that of the Sun was observed to look like a fading pair-instability supernova. Here we report observations of two slow-to-fade super-luminous supernovae that show relatively fast rise times and blue colours, which are incompatible with pair-instability models. Their late-time light-curve and spectral similarities to supernova 2007bi call the nature of that event into question. Our early spectra closely resemble typical fast-declining super-luminous supernovae, which are not powered by radioactivity. Modelling our observations with 10-16 solar masses of magnetar-energized ejecta demonstrates the possibility of a common explosion mechanism. The lack of unambiguous nearby pair-instability events suggests that their local rate of occurrence is less than 6x10^-6 times that of the core-collapse rate.

We analyze the wavelength-dependent variability of a sample of spectroscopically confirmed active galactic nuclei (AGN) selected from near-UV ($NUV$) variable sources in the GALEX Time Domain Survey that have a large amplitude of optical variability (difference-flux S/N $>$ 3) in the Pan-STARRS1 Medium Deep Survey (PS1 MDS). By matching GALEX and PS1 epochs in 5 bands ($NUV$, $g_{P1}$, $r_{P1}$, $i_{P1}$, $z_{P1}$) in time, and taking their flux difference, we create co-temporal difference-flux spectral energy distributions ($\Delta f$SEDs) using two chosen epochs for each of the 23 objects in our sample on timescales of about a year. We confirm the "bluer-when-brighter" trend reported in previous studies, and measure a median spectral index of the $\Delta f$SEDs of $\alpha_{\lambda}$ = 2.1 that is consistent with an accretion disk spectrum. We further fit the $\Delta f$SEDs of each source with a standard accretion disk model in which the accretion rate changes from one epoch to the other. In our sample, 17 out of 23 ($\sim$74 %) sources are well described by this variable accretion-rate disk model, with a median average characteristic disk temperature $\bar{T}^*$ of $1.2\times 10^5$~K that is consistent with the temperatures expected given the distribution of accretion rates and black hole masses inferred for the sample. Our analysis also shows that the variable accretion rate model is a better fit to the $\Delta f$SEDs than a simple power law.

We discuss 76 large amplitude transients (Delta-m>1.5) occurring in the nuclei of galaxies, nearly all with no previously known Active Galactic Nucleus (AGN). They have been discovered as part of the Pan-STARRS1 (PS1) 3pi survey, by comparison with SDSS photometry a decade earlier, and then monitored with the Liverpool Telescope, and studied spectroscopically with the William Herschel Telescope (WHT). Based on colours, light curve shape, and spectra, these transients fall into four groups. A few are misclassified stars or objects of unknown type. Some are red/fast transients and are known or likely nuclear supernovae. A few are either radio sources or erratic variables and so likely blazars. However the majority (~66%) are blue and evolve slowly, on a timescale of years. Spectroscopy shows them to be AGN at z~ 0.3 - 1.4, which must have brightened since the SDSS photometry by around an order of magnitude. It is likely that these objects were in fact AGN a decade ago, but too weak to be recognised by SDSS; they could then be classed as "hypervariable" AGN. By searching the SDSS Stripe 82 quasar database, we find 15 similar objects. We discuss several possible explanations for these slow blue hypervariables - (i) unusually luminous tidal disruption events; (ii) extinction events; (iii) changes in accretion state; and (iv) large amplitude microlensing by stars in foreground galaxies. A mixture of explanations (iii) and (iv) seems most likely. Both hold promise of considerable new insight into the AGN phenomenon.

We study the properties of 975 active galactic nuclei (AGN) selected by variability in the Pan-STARRS1 Medium-Deep Survey. Using complementary multi wavelength data from the ultraviolet to the far-infrared, we use SED fitting to determine the AGN and host properties at $z<1$, and compare to a well-matched control sample. We confirm the trend previously observed that the variability amplitude decreases with AGN luminosity, but on the other hand, we observe that the slope of this relation steepens with wavelength resulting in a "redder when brighter" trend at low luminosities. Our results show that AGN are hosted by more massive hosts than control sample galaxies, while the restframe, dust-corrected $NUV-r$ color distribution of AGN hosts is similar to control galaxies. We find a positive correlation between the AGN luminosity and star formation rate (SFR), independent of redshift. AGN hosts populate the whole range of SFRs within and outside the Main Sequence of star forming galaxies. Comparing the distribution of AGN hosts and control galaxies, we show that AGN hosts are less likely to be hosted by quiescent galaxies, but more likely to be hosted by Main Sequence or starburst galaxies.

We apply a combination of a Genetic Algorithms (GA) and Support Vector Machines (SVM) machine learning algorithm to solve two important problems faced by the astronomical community: star/galaxy separation, and photometric redshift estimation of galaxies in survey catalogs. We use the GA to select the relevant features in the first step, followed by optimization of SVM parameters in the second step to obtain an optimal set of parameters to classify or regress, in process of which we avoid over-fitting. We apply our method to star/galaxy separation in Pan-STARRS1 data. We show that our method correctly classifies 98% of objects down to i_P1= 24.5, with a completeness (or true positive rate) of 99% for galaxies, and 88% for stars. By combining colors with morphology, our star/classification method yields better results than the new SExtractor classifier spread_model in particular at the faint end (i_P1>22). We also use our method to derive photometric redshifts for galaxies in the COSMOS bright multi-wavelength dataset down to an error in (1+z) of sigma=0.013, which compares well with estimates from SED fitting on the same data (sigma=0.007) while making a significantly smaller number of assumptions.

We present a Hubble Space Telescope STIS spectrum of ASASSN-14li, the first rest-frame UV spectrum of a tidal disruption flare (TDF). The underlying continuum is well fit by a blackbody with $T_{\mathrm{UV}} = 3.5 \times 10^{4}$ K, an order of magnitude smaller than the temperature inferred from X-ray spectra (and significantly more precise than previous efforts based on optical and near-UV photometry). Super-imposed on this blue continuum, we detect three classes of features: narrow absorption from the Milky Way (probably a high-velocity cloud), and narrow absorption and broad (FWHM $\approx 2000$-8000 km s$^{-1}$) emission lines at/near the systemic host velocity. The absorption lines are blueshifted with respect to the emission lines by $\Delta v = -(250$-400) km s$^{-1}$. Due both to this velocity offset and the lack of common low-ionization features (Mg II, Fe II), we argue these arise from the same absorbing material responsible for the low-velocity outflow discovered at X-ray wavelengths. The broad nuclear emission lines display a remarkable abundance pattern: N III], N IV], He II are quite prominent, while the common quasar emission lines of C III] and Mg II are weak or entirely absent. Detailed modeling of this spectrum will help elucidate fundamental questions regarding the nature of the emission process(es) at work in TDFs, while future UV spectroscopy of ASASSN-14li would help to confirm (or refute) the previously proposed connection between TDFs and "N-rich" quasars.

Extreme coronal-line emitter (ECLE) SDSSJ095209.56+214313.3, known by its strong, fading, high ionization lines, has been a long standing candidate for a tidal disruption event, however a supernova origin has not yet been ruled out. Here we add several new pieces of information to the puzzle of the nature of the transient that powered its variable coronal lines: 1) an optical light curve from the Lincoln Near Earth Asteroid Research (LINEAR) survey that serendipitously catches the optical flare, and 2) late-time observations of the host galaxy with the Swift Ultraviolet and Optical Telescope (UVOT) and X-ray telescope (XRT) and the ground-based Mercator telescope. The well-sampled, $\sim10$-year long, unfiltered LINEAR light curve constrains the onset of the flare to a precision of $\pm5$ days and enables us to place a lower limit on the peak optical magnitude. Difference imaging allows us to estimate the location of the flare in proximity of the host galaxy core. Comparison of the \textsl{GALEX} data (early 2006) with the recently acquired Swift UVOT (June 2015) and Mercator observations (April 2015) demonstrate a decrease in the UV flux over a $\sim 10$ year period, confirming that the flare was UV-bright. The long-lived UV-bright emission, detected 1.8 rest-frame years after the start of the flare, strongly disfavors a SN origin. These new data allow us to conclude that the flare was indeed powered by the tidal disruption of a star by a supermassive black hole and that TDEs are in fact capable of powering the enigmatic class of ECLEs.

We present late-time observations of the tidal disruption event candidate PS1-10jh. UV and optical imaging with HST/WFC3 localize the transient to be coincident with the host galaxy nucleus to an accuracy of 0.023 arcsec, corresponding to 66 pc. The UV flux in the F225W filter, measured 3.35 rest-frame years after the peak of the nuclear flare, is consistent with a decline that continues to follow a $t^{-5/3}$ power-law with no spectral evolution. Late epochs of optical spectroscopy obtained with MMT ~ 2 and 4 years after the peak, enable a clean subtraction of the host galaxy from the early spectra, revealing broad helium emission lines on top of a hot continuum, and placing stringent upper limits on the presence of hydrogen line emission. We do not measure Balmer H\delta absorption in the host galaxy strong enough to be indicative of a rare, post-starburst "E+A" galaxy as reported by Arcavi et al. (2014). The light curve of PS1-10jh over a baseline of 3.5 yr is best modeled by fallback accretion of a tidally disrupted star. Its strong broad helium emission relative to hydrogen (He II \lambda 4686/H\alpha > 5) could be indicative of either the hydrogen-poor chemical composition of the disrupted star, or certain conditions in the tidal debris of a solar-composition star in the presence of an optically-thick, extended reprocessing envelope.

Tidal forces close to massive black holes can violently disrupt stars that make a close approach. These extreme events are discovered via bright X-ray and optical/UV flares in galactic centers. Prior studies based on modeling decaying flux trends have been able to estimate broad properties, such as the mass accretion rate. Here we report the detection of flows of highly ionized X-ray gas in high-resolution X-ray spectra of a nearby tidal disruption event. Variability within the absorption-dominated spectra indicates that the gas is relatively close to the black hole. Narrow line widths indicate that the gas does not stretch over a large range of radii, giving a low volume filling factor. Modest outflow speeds of a few hundred kilometers per second are observed, significantly below the escape speed from the radius set by variability. The gas flow is consistent with a rotating wind from the inner, super-Eddington region of a nascent accretion disk, or with a filament of disrupted stellar gas near to the apocenter of an elliptical orbit. Flows of this sort are predicted by fundamental analytical theory and more recent numerical simulations.

For the first time in history, humans have reached the point where it is possible to construct a revolutionary space-based observatory that has the capability to find dozens of Earth-like worlds, and possibly some with signs of life. This same telescope, designed as a long-lived facility, would also produce transformational scientific advances in every area of astronomy and astrophysics from black hole physics to galaxy formation, from star and planet formation to the origins of the Solar System. The Association of Universities for Research in Astronomy (AURA) commissioned a study on a next-generation UVOIR space observatory with the highest possible scientific impact in the era following JWST. This community-based study focuses on the future space-based options for UV and optical astronomy that significantly advance our understanding of the origin and evolution of the cosmos and the life within it. The committee concludes that a space telescope equipped with a 12-meter class primary mirror can find and characterize dozens of Earth-like planets and make fundamental advances across nearly all fields of astrophysics. The concept is called the High Definition Space Telescope (HDST). The telescope would be located at the Sun-Earth L2 point and would cover a spectral range that, at a minimum, runs from 0.1 to 2 microns. Unlike JWST, HDST will not need to operate at cryogenic temperatures. HDST can be made to be serviceable on orbit but does not require servicing to complete its primary scientific objectives. We present the scientific and technical requirements for HDST and show that it could allow us to determine whether or not life is common outside the Solar System. We do not propose a specific design for such a telescope, but show that designing, building and funding such a facility is feasible beginning in the next decade - if the necessary strategic investments in technology begin now.

Supermassive black hole binaries (SMBHBs) should be an inevitable consequence of the hierarchical growth of massive galaxies through mergers, and the strongest sirens of gravitational waves (GWs) in the cosmos. And yet, their direct detection has remained elusive due to the compact (sub-parsec) orbital separations of gravitationally bound SMBHBs. Here we exploit a theoretically predicted signature of a SMBHB in the time domain: periodic variability caused by a mass accretion rate that is modulated by the binary's orbital motion. We report our first significant periodically varying quasar detection from the systematic search in the Pan-STARRS1 (PS1) Medium Deep Survey. Our SMBHB candidate, PSO J334.2028+01.4075, is a luminous radio-loud quasar at $z=2.060$, with extended baseline photometry from the Catalina Real-Time Transient Survey, as well as archival spectroscopy from the FIRST Bright Quasar Survey. The observed period ($542 \pm 15$ days) and estimated black hole mass ($\log (M_{\rm BH}/M_\odot) = 9.97 \pm 0.50$), correspond to an orbital separation of $7^{+8}_{-4}$ Schwarzschild radii ($\sim 0.006^{+0.007}_{-0.003}$ pc), assuming the rest-frame period of the quasar variability traces the orbital period of the binary. This SMBHB candidate, discovered at the peak redshift for SMBH mergers, is in a physically stable configuration for a circumbinary accretion disk, and within the regime of GW-driven orbital decay. Our search with PS1 is a benchmark study for the exciting capabilities of LSST, which will have orders of magnitude larger survey power, and will potentially pinpoint the locations of thousands of SMBHBs in the variable night sky.

The radius and surface composition of an exploding massive star,as well as the explosion energy per unit mass, can be measured using early UV observations of core collapse supernovae (SNe). We present the first results from a simultaneous GALEX/PTF search for early UV emission from SNe. Six Type II SNe and one Type II superluminous SN (SLSN-II) are clearly detected in the GALEX NUV data. We compare our detection rate with theoretical estimates based on early, shock-cooling UV light curves calculated from models that fit existing Swift and GALEX observations well, combined with volumetric SN rates. We find that our observations are in good agreement with calculated rates assuming that red supergiants (RSGs) explode with fiducial radii of 500 solar, explosion energies of 10^51 erg, and ejecta masses of 10 solar masses. Exploding blue supergiants and Wolf-Rayet stars are poorly constrained. We describe how such observations can be used to derive the progenitor radius, surface composition and explosion energy per unit mass of such SN events, and we demonstrate why UV observations are critical for such measurements. We use the fiducial RSG parameters to estimate the detection rate of SNe during the shock-cooling phase (<1d after explosion) for several ground-based surveys (PTF, ZTF, and LSST). We show that the proposed wide-field UV explorer ULTRASAT mission, is expected to find >100 SNe per year (~0.5 SN per deg^2), independent of host galaxy extinction, down to an NUV detection limit of 21.5 mag AB. Our pilot GALEX/PTF project thus convincingly demonstrates that a dedicated, systematic SN survey at the NUV band is a compelling method to study how massive stars end their life.

We present the GALEX detection of a UV burst at the time of explosion of an optically normal Type II-P supernova (PS1-13arp) from the Pan-STARRS1 survey at z=0.1665. The temperature and luminosity of the UV burst match the theoretical predictions for shock breakout in a red supergiant, but with a duration a factor of ~50 longer than expected. We compare the $NUV$ light curve of PS1-13arp to previous GALEX detections of Type IIP SNe, and find clear distinctions that indicate that the UV emission is powered by shock breakout, and not by the subsequent cooling envelope emission previously detected in these systems. We interpret the ~ 1 d duration of the UV signal with a shock breakout in the wind of a red supergiant with a pre-explosion mass-loss rate of ~ 10^-3 Msun yr^-1. This mass-loss rate is enough to prolong the duration of the shock breakout signal, but not enough to produce an excess in the optical plateau light curve or narrow emission lines powered by circumstellar interaction. This detection of non-standard, potentially episodic high mass-loss in a RSG SN progenitor has favorable consequences for the prospects of future wide-field UV surveys to detect shock breakout directly in these systems, and provide a sensitive probe of the pre-explosion conditions of SN progenitors.

In recent years, wide-field sky surveys providing deep multi-band imaging have presented a new path for indirectly characterizing the progenitor populations of core-collapse supernovae (SN): systematic light curve studies. We assemble a set of 76 grizy-band Type IIP SN light curves from Pan-STARRS1, obtained over a constant survey program of 4 years and classified using both spectroscopy and machine learning-based photometric techniques. We develop and apply a new Bayesian model for the full multi-band evolution of each light curve in the sample. We find no evidence of a sub-population of fast-declining explosions (historically referred to as "Type IIL" SNe). However, we identify a highly significant relation between the plateau phase decay rate and peak luminosity among our SNe IIP. These results argue in favor of a single parameter, likely determined by initial stellar mass, predominantly controlling the explosions of red supergiants. This relation could also be applied for supernova cosmology, offering a standardizable candle good to an intrinsic scatter of 0.2 mag. We compare each light curve to physical models from hydrodynamic simulations to estimate progenitor initial masses and other properties of the Pan-STARRS1 Type IIP SN sample. We show that correction of systematic discrepancies between modeled and observed SN IIP light curve properties and an expanded grid of progenitor properties, are needed to enable robust progenitor inferences from multi-band light curve samples of this kind. This work will serve as a pathfinder for photometric studies of core-collapse SNe to be conducted through future wide field transient searches.

We present a novel method for the light-curve characterization of Pan-STARRS1 Medium Deep Survey (PS1 MDS) extragalactic sources into stochastic variables (SV) and burst-like (BL) transients, using multi-band image-differencing time-series data. We select detections in difference images associated with galaxy hosts using a star/galaxy catalog extracted from the deep PS1 MDS stacked images, and adopt a maximum a posteriori formulation to model their difference-flux time-series in four Pan-STARRS1 photometric bands g,r,i, and z. We use three deterministic light-curve models to fit burst-like transients and one stochastic light curve model, the Ornstein-Uhlenbeck process, in order to fit variability that is characteristic of active galactic nuclei (AGN). We assess the quality of fit of the models band-wise source-wise, using their estimated leave-out-one cross-validation likelihoods and corrected Akaike information criteria. We then apply a K-means clustering algorithm on these statistics, to determine the source classification in each band. The final source classification is derived as a combination of the individual filter classifications. We use our clustering method to characterize 4361 extragalactic image difference detected sources in the first 2.5 years of the PS1 MDS, into 1529 BL, and 2262 SV, with a purity of 95.00% for AGN, and 90.97% for SN based on our verification sets. We combine our light-curve classifications with their nuclear or off-nuclear host galaxy offsets, to define a robust photometric sample of 1233 active galactic nuclei and 812 supernovae. We use these samples to identify simple photometric priors that would enable their real-time identification in future wide-field synoptic surveys.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

In the past decade, several rapidly-evolving transients have been discovered whose timescales and luminosities are not easily explained by traditional supernovae (SN) models. The sample size of these objects has remained small due, at least in part, to the challenge of detecting short timescale transients with traditional survey cadences. Here we present the results from a search within the Pan-STARRS1 Medium Deep Survey (PS1-MDS) for rapidly-evolving and luminous transients. We identify 10 new transients with a time above half-maximum of less than 12 days and -16.5 > M > -20 mag. This increases the number of known events in this region of SN phase space by roughly a factor of three. The median redshift of the PS1-MDS sample is z=0.275 and they all exploded in star forming galaxies. In general, the transients possess faster rise than decline timescale and blue colors at maximum light (g - r < -0.2). Best fit blackbodies reveal photospheric temperatures/radii that expand/cool with time and explosion spectra taken near maximum light are dominated by a blue continuum, consistent with a hot, optically thick, ejecta. We find it difficult to reconcile the short timescale, high peak luminosity (L > 10^43 erg/s), and lack of UV line blanketing observed in many of these transients with an explosion powered mainly by the radioactive decay of Ni-56. Rather, we find that many are consistent with either (1) cooling envelope emission from the explosion of a star with a low-mass extended envelope which ejected very little (<0.03 M_sun) radioactive material, or (2) a shock breakout within a dense, optically thick, wind surrounding the progenitor star. After calculating the detection efficiency for objects with rapid timescales in the PS1-MDS we find a volumetric rate of 4800 - 8000 events/yr/Gpc^3 (4-7% of the core-collapse SN rate at z=0.2).

We present the Pan-STARRS1 discovery of the long-lived and blue transient PS1-11af, which was also detected by GALEX with coordinated observations in the near-ultraviolet (NUV) band. PS1-11af is associated with the nucleus of an early-type galaxy at redshift z=0.4046 that exhibits no evidence for star formation or AGN activity. Four epochs of spectroscopy reveal a pair of transient broad absorption features in the UV on otherwise featureless spectra. Despite the superficial similarity of these features to P-Cygni absorptions of supernovae (SNe), we conclude that PS1-11af is not consistent with the properties of known types of SNe. Blackbody fits to the spectral energy distribution are inconsistent with the cooling, expanding ejecta of a SN, and the velocities of the absorption features are too high to represent material in homologous expansion near a SN photosphere. However, the constant blue colors and slow evolution of the luminosity are similar to previous optically-selected tidal disruption events (TDEs). The shape of the optical light curve is consistent with models for TDEs, but the minimum accreted mass necessary to power the observed luminosity is only ~0.002M_sun, which points to a partial disruption model. A full disruption model predicts higher bolometric luminosities, which would require most of the radiation to be emitted in a separate component at high energies where we lack observations. In addition, the observed temperature is lower than that predicted by pure accretion disk models for TDEs and requires reprocessing to a constant, lower temperature. Three deep non-detections in the radio with the VLA over the first two years after the event set strict limits on the production of any relativistic outflow comparable to Swift J1644+57, even if off-axis.

We present the selection and classification of over a thousand ultraviolet (UV) variable sources discovered in ~ 40 deg^2 of GALEX Time Domain Survey (TDS) NUV images observed with a cadence of 2 days and a baseline of observations of ~ 3 years. The GALEX TDS fields were designed to be in spatial and temporal coordination with the Pan-STARRS1 Medium Deep Survey, which provides deep optical imaging and simultaneous optical transient detections via image differencing. We characterize the GALEX photometric errors empirically as a function of mean magnitude, and select sources that vary at the 5\sigma level in at least one epoch. We measure the statistical properties of the UV variability, including the structure function on timescales of days and years. We report classifications for the GALEX TDS sample using a combination of optical host colors and morphology, UV light curve characteristics, and matches to archival X-ray, and spectroscopy catalogs. We classify 62% of the sources as active galaxies (358 quasars and 305 active galactic nuclei), and 10% as variable stars (including 37 RR Lyrae, 53 M dwarf flare stars, and 2 cataclysmic variables). The remaining unclassified sources include UV-bright extragalactic transients, two of which have been spectroscopically confirmed to be a young core-collapse supernova and a flare from the tidal disruption of a star by dormant supermassive black hole. We calculate a surface density for variable sources in the UV with NUV < 23 mag and |\Delta m| > 0.2 mag of ~ 8.0, 7.7, and 1.8 deg^-2 for quasars, AGNs, and RR Lyrae stars, respectively. We also calculate a surface density rate in the UV for transient sources, using the effective survey time at the cadence appropriate to each class, of ~ 15 and 52 deg^-2 yr^-1 for M dwarfs and extragalactic transients, respectively. (Abridged.)

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

The flare of radiation from the tidal disruption and accretion of a star can be used as a marker for supermassive black holes that otherwise lie dormant and undetected in the centres of distant galaxies. Previous candidate flares have had declining light curves in good agreement with expectations, but with poor constraints on the time of disruption and the type of star disrupted, because the rising emission was not observed. Recently, two `relativistic' candidate tidal disruption events were discovered, each of whose extreme X-ray luminosity and synchrotron radio emission were interpreted as the onset of emission from a relativistic jet. Here we report the discovery of a luminous ultraviolet-optical flare from the nuclear region of an inactive galaxy at a redshift of 0.1696. The observed continuum is cooler than expected for a simple accreting debris disk, but the well-sampled rise and decline of its light curve follows the predicted mass accretion rate, and can be modelled to determine the time of disruption to an accuracy of two days. The black hole has a mass of about 2 million solar masses, modulo a factor dependent on the mass and radius of the star disrupted. On the basis of the spectroscopic signature of ionized helium from the unbound debris, we determine that the disrupted star was a helium-rich stellar core.

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.

We present the discovery of two ultra-luminous supernovae (SNe) at z ~ 0.9 with the Pan-STARRS1 Medium-Deep Survey. These SNe, PS1-10ky and PS1-10awh, are amongst the most luminous SNe ever discovered, comparable to the unusual transients SN 2005ap and SCP 06F6. Like SN 2005ap and SCP 06F6, they show characteristic high luminosities (M_bol ~ -22.5 mag), blue spectra with a few broad absorption lines, and no evidence for H or He. We have constructed a full multi-color light curve sensitive to the peak of the spectral energy distribution in the rest-frame ultraviolet, and we have obtained time-series spectroscopy for these SNe. Given the similarities between the SNe, we combine their light curves to estimate a total radiated energy over the course of explosion of (0.9-1.4) x 10^51 erg. We find photospheric velocities of 12,000-19,000 km/s with no evidence for deceleration measured across ~3 rest-frame weeks around light-curve peak, consistent with the expansion of an optically-thick massive shell of material. We show that, consistent with findings for other ultra-luminous SNe in this class, radioactive decay is not sufficient to power PS1-10ky, and we discuss two plausible origins for these events: the initial spin-down of a newborn magnetar in a core-collapse SN, or SN shock breakout from the dense circumstellar wind surrounding a Wolf-Rayet star.

We report on the discovery and observations of the extremely luminous optical transient CSS100217:102913+404220 (CSS100217 hereafter). Spectroscopic observations show this transient was coincident with a galaxy at redshift z=0.147, and reached an apparent magnitude of V ~ 16.3. After correcting for foreground Galactic extinction we determine the absolute magnitude to be M_V =-22.7 approximately 45 days after maximum light. Based on our unfiltered optical photometry the peak optical emission was L = 1.3 x 10^45 erg s^-1, and over a period of 287 rest-frame days had an integrated bolometric luminosity of 1.2 x 10^52 erg. Analysis of the pre-outburst SDSS spectrum of the source shows features consistent with a Narrow-line Seyfert1 (NLS1) galaxy. High-resolution HST and Keck followup observations show the event occurred within 150pc of nucleus of the galaxy, suggesting a possible link to the active nuclear region. However, the rapid outburst along with photometric and spectroscopic evolution are much more consistent with a luminous supernova. Line diagnostics suggest that the host galaxy is undergoing significant star formation. We use extensive follow-up of the event along with archival CSS and SDSS data to investigate the three most likely sources of such an event; 1) an extremely luminous supernova; 2) the tidal disruption of a star by the massive nuclear black hole; 3) variability of the central AGN. We find that CSS100217 was likely an extremely luminous type IIn supernova that occurred within range of the narrow-line region of an AGN. We discuss how similar events may have been missed in past supernova surveys because of confusion with AGN activity.

We report the discovery of a supernova (SN) with the highest apparent energy output to date and conclude that it represents an extreme example of the Type IIn subclass. The SN, which was discovered behind the Large Magellanic Cloud at z = 0.289 by the SuperMACHO microlensing survey, peaked at M_R = -21.5 mag and only declined by 2.9 mag over 4.7 years after the peak. Over this period, SN 2003ma had an integrated bolometric luminosity of 4 x 10^51 ergs, more than any other SN to date. The radiated energy is close to the limit allowed by conventional core-collapse explosions. Optical spectra reveal that SN 2003ma has persistent single-peaked intermediate-width hydrogen lines, a signature of interaction between the SN and a dense circumstellar medium. The light curves show further evidence for circumstellar interaction, including a long plateau with a shape very similar to the classic SN IIn 1988Z -- however, SN 2003ma is ten times more luminous at all epochs. The fast velocity measured for the intermediate-width H_alpha component (~6000 km/s) points towards an extremely energetic explosion (> 10^52 ergs), which imparts a faster blast-wave speed to the post-shock material and a higher luminosity from the interaction than is observed in typical SNe IIn. Mid-infrared observations of SN 2003ma suggest an infrared light echo is produced by normal interstellar dust at a distance ~0.5 pc from the SN.

Recent searches by unbiased, wide-field surveys have uncovered a group of extremely luminous optical transients. The initial discoveries of SN 2005ap by the Texas Supernova Search and SCP-06F6 in a deep Hubble pencil beam survey were followed by the Palomar Transient Factory confirmation of host redshifts for other similar transients. The transients share the common properties of high optical luminosities (peak magnitudes ~ -21 to -23), blue colors, and a lack of H or He spectral features. The physical mechanism that produces the luminosity is uncertain, with suggestions ranging from jet-driven explosion to pulsational pair-instability. Here we report the most detailed photometric and spectral coverage of an ultra-bright transient (SN 2010gx) detected in the Pan-STARRS 1 sky survey. In common with other transients in this family, early-time spectra show a blue continuum, and prominent broad absorption lines of O II. However, about 25d after discovery, the spectra developed type Ic supernova features, showing the characteristic broad Fe II and Si II absorption lines. Detailed, post-maximum follow-up may show that all SN 2005ap and SCP-06F6 type transients are linked to supernovae Ic. This poses problems in understanding the physics of the explosions: there is no indication from late-time photometry that the luminosity is powered by 56Ni, the broad lightcurves suggest very large ejected masses, and the slow spectral evolution is quite different from typical Ic timescales. The nature of the progenitor stars and the origin of the luminosity are intriguing and open questions.

We present the early UV and optical light curve of Type IIP supernova (SN) 2010aq at z=0.0862, and compare it to analytical models for thermal emission following SN shock breakout in a red supergiant star. SN 2010aq was discovered in joint monitoring between the Galaxy Evolution Explorer (GALEX) Time Domain Survey (TDS) in the NUV and the Pan-STARRS1 Medium Deep Survey (PS1 MDS) in the g, r, i, and z bands. The GALEX and Pan-STARRS1 observations detect the SN less than 1 day after shock breakout, measure a diluted blackbody temperature of 31,000 +/- 6,000 K 1 day later, and follow the rise in the UV/optical light curve over the next 2 days caused by the expansion and cooling of the SN ejecta. The high signal-to-noise ratio of the simultaneous UV and optical photometry allows us to fit for a progenitor star radius of 700 +/- 200 R_sun, the size of a red supergiant star. An excess in UV emission two weeks after shock breakout compared to SNe well fitted by model atmosphere-code synthetic spectra with solar metallicity, is best explained by suppressed line blanketing due to a lower metallicity progenitor star in SN 2010aq. Continued monitoring of PS1 MDS fields by the GALEX TDS will increase the sample of early UV detections of Type II SNe by an order of magnitude, and probe the diversity of SN progenitor star properties.

We present photometric and spectroscopic observations of a luminous type IIP Supernova 2009kf discovered by the Pan-STARRS 1 (PS1) survey and detected also by GALEX. The SN shows a plateau in its optical and bolometric light curves, lasting approximately 70 days in the rest frame, with absolute magnitude of M_V = -18.4 mag. The P-Cygni profiles of hydrogen indicate expansion velocities of 9000km/s at 61 days after discovery which is extremely high for a type IIP SN. SN 2009kf is also remarkably bright in the near-ultraviolet (NUV) and shows a slow evolution 10-20 days after optical discovery. The NUV and optical luminosity at these epochs can be modelled with a black-body with a hot effective temperature (T ~16,000 K) and a large radius (R ~1x10^{15} cm). The bright bolometric and NUV luminosity, the lightcurve peak and plateau duration, the high velocities and temperatures suggest that 2009kf is a type IIP SN powered by a larger than normal explosion energy. Recently discovered high-z SNe (0.7 < z < 2.3) have been assumed to be IIn SNe, with the bright UV luminosities due to the interaction of SN ejecta with a dense circumstellar medium (CSM). UV bright SNe similar to SN 2009kf could also account for these high-z events, and its absolute magnitude M_NUV = -21.5 +/- 0.5 mag suggests such SNe could be discovered out to z ~2.5 in the PS1 survey.

A dormant supermassive black hole lurking in the center of a galaxy will be revealed when a star passes close enough to be torn apart by tidal forces, and a flare of electromagnetic radiation is emitted when the bound fraction of the stellar debris falls back onto the black hole and is accreted. Here we present the third candidate tidal disruption event discovered in the GALEX Deep Imaging Survey: a 1.6x10^{43} erg s^{-1} UV/optical flare from a star-forming galaxy at z=0.1855. The UV/optical SED during the peak of the flare measured by GALEX and Palomar LFC imaging can be modeled as a single temperature blackbody with T_{bb}=1.7x10^{5} K and a bolometric luminosity of 3x10^{45} erg s^{-1}, assuming an internal extinction with E(B-V)_{gas}=0.3. The Chandra upper limit on the X-ray luminosity during the peak of the flare, L_{X}(2-10 keV)< 10^{41} erg s^{-1}, is 2 orders of magnitude fainter than expected from the ratios of UV to X-ray flux density observed in active galaxies. We compare the light curves and broadband properties of all three tidal disruption candidates discovered by GALEX, and find that (1) the light curves are well fitted by the power-law decline expected for the fallback of debris from a tidally disrupted solar-type star, and (2) the UV/optical SEDs can be attributed to thermal emission from an envelope of debris located at roughly 10 times the tidal disruption radius of a ~10^{7} M_sun central black hole. We use the observed peak absolute optical magnitudes of the flares (-17.5 > M_{g} > -18.9) to predict the detection capabilities of upcoming optical synoptic surveys. (Abridged)

We discuss the central role played by X-ray studies to reconstruct the past history of formation and evolution of supermassive Black Holes (BHs), and the role they played in shaping the properties of their host galaxies. We shortly review the progress in this field contributed by the current X-ray and multiwavelength surveys. Then, we focus on the outstanding scientific questions that have been opened by observations carried out in the last years and that represent the legacy of Chandra and XMM, as for X-ray observations, and the legacy of the SDSS, as for wide area surveys: 1) When and how did the first supermassive black holes form? 2) How does cosmic environment regulate nuclear activity (and star formation) across cosmic time? 3) What is the history of nuclear activity in a galaxy lifetime? We show that the most efficient observational strategy to address these questions is to carry out a large-area X-ray survey, reaching a sensitivity comparable to that of deep Chandra and XMM pointings, but extending over several thousands of square degrees. Such a survey can only be carried out with a Wide-Field X-ray Telescope (WFXT) with a high survey speed, due to the combination of large field of view and large effective area, i.e., grasp, and sharp PSF. We emphasize the important synergies that WFXT will have with a number of future groundbased and space telescopes, covering from the radio to the X-ray bands and discuss the immense legacy value that such a mission will have for extragalactic astronomy at large.

Tidal disruption events provide a unique probe of quiescent black holes in the nuclei of distant galaxies. The next generation of synoptic surveys will yield a large sample of flares from the tidal disruption of stars by massive black holes that will give insights to four key science questions: 1) What is the assembly history of massive black holes in the universe? 2) Is there a population of intermediate mass black holes that are the primordial seeds of supermassive black holes? 3) How can we increase our understanding of the physics of accretion onto black holes? 4) Can we localize sources of gravitational waves from the detection of tidal disruption events around massive black holes and recoiling binary black hole mergers?

Over the next decade, we can expect time domain astronomy to flourish at optical and radio wavelengths. In parallel with these efforts, a dedicated transient "machine" operating at higher energies (X-ray band through soft gamma-rays) is required to reveal the unique subset of events with variable emission predominantly visible above 100 eV. Here we focus on the transient phase space never yet sampled due to the lack of a sensitive, wide-field and triggering facility dedicated exclusively to catching high energy transients and enabling rapid coordinated multi-wavelength follow-up. We first describe the advancements in our understanding of known X-ray transients that can only be enabled through such a facility and then focus on the classes of transients theoretically predicted to be out of reach of current detection capabilities. Finally there is the exciting opportunity of revealing new classes of X-ray transients and unveiling their nature through coordinated follow-up observations at longer wavelengths.

SDSSJ092712.65+294344.0 was identified by the SDSS as a quasar, but has the unusual property of having two emission line systems offset by 2650 km/s. One of these contains the usual combination of broad and narrow lines, the other only narrow lines. In the first paper commenting on this system (Komossa et al. 2008), it was interpreted as a galaxy in which a pair of black holes had merged, imparting a several thousand km/s recoil to the new, larger black hole. In two other papers (Bogdanovic, Eracleous & Sigurdsson 2008; Dotti et al. 2008), it was interpreted as a small-separation binary black hole. We propose a new interpretation: that this system is a more distant analog of NGC1275, a large and a small galaxy interacting near the center of a rich cluster.

We report the discovery by the Robotic Optical Transient Experiment (ROTSE-IIIb) telescope of SN 2008es, an overluminous supernova (SN) at z=0.205 with a peak visual magnitude of -22.2. We present multiwavelength follow-up observations with the Swift satellite and several ground-based optical telescopes. The ROTSE-IIIb observations constrain the time of explosion to be 23+/-1 rest-frame days before maximum. The linear decay of the optical light curve, and the combination of a symmetric, broad H\alpha emission line profile with broad P Cygni H\beta and Na I \lambda5892 profiles, are properties reminiscent of the bright Type II-L SNe 1979C and 1980K, although SN 2008es is greater than 10 times more luminous. The host galaxy is undetected in pre-supernova Sloan Digital Sky Survey images, and similar to Type II-L SN 2005ap (the most luminous SN ever observed), the host is most likely a dwarf galaxy with M_r > -17. Swift Ultraviolet/Optical Telescope observations in combination with Palomar photometry measure the SED of the SN from 200 to 800 nm to be a blackbody that cools from a temperature of 14,000 K at the time of the optical peak to 6400 K 65 days later. The inferred blackbody radius is in good agreement with the radius expected for the expansion speed measured from the broad lines (10,000 km/s). The bolometric luminosity at the optical peak is 2.8 x 10^44 erg/s, with a total energy radiated over the next 65 days of 5.6 x 10^50 erg. We favor a model in which the exceptional peak luminosity is a consequence of the core-collapse explosion of a progenitor star with a low-mass extended hydrogen envelope and a stellar wind with a density close to the upper limit on the mass-loss rate measured from the lack of an X-ray detection by the Swift X-Ray Telescope. (Abridged).

We present a multi-epoch quantitative spectroscopic analysis of the Type IIn SN 1994W, an event interpreted by Chugai et al. as stemming from the interaction between the ejecta of a SN and a 0.4Msun circumstellar shell ejected 1.5yr before core collapse. During the brightening phase, our models suggest that the source of optical radiation is not unique, perhaps associated with an inner optically-thick Cold Dense Shell (CDS) and outer optically-thin shocked material. During the fading phase, our models support a single source of radiation, an hydrogen-rich optically-thick layer with a near-constant temperature of ~7000K that recedes from a radius of 4.3x10^15 at peak to 2.3x10^15cm 40 days later. We reproduce the hybrid narrow-core broad-wing line profile shapes of SN 1994W at all times, invoking an optically-thick photosphere exclusively (i.e., without any external optically-thick shell). In SN 1994W, slow expansion makes scattering with thermal electrons a key escape mechanism for photons trapped in optically-thick line cores, and allows the resulting broad incoherent electron-scattering wings to be seen around narrow line cores. In SNe with larger expansion velocities, the thermal broadening due to incoherent scattering is masked by the broad profile and the dominant frequency redshift occasioned by bulk motions. Given the absence of broad lines at all times and the very low 56Ni yields, we speculate whether SN 1994W could have resulted from an interaction between two ejected shells without core collapse. The high conversion efficiency of kinetic to thermal energy may not require a SN-like energy budget for SN1994W.

We report the serendipitous detection by GALEX of fast (<1 day) rising (>1 mag) UV emission from two Type II plateau (II-P) supernovae (SNe) at z=0.185 and 0.324 discovered by the Supernova Legacy Survey. Optical photometry and VLT spectroscopy 2 weeks after the GALEX detections link the onset of UV emission to the time of shock breakout. Using radiation hydrodynamics and non-LTE radiative transfer simulations, and starting from a standard red supergiant (RSG; Type II-P SN progenitor) star evolved self-consistently from the main sequence to iron core collapse, we model the shock breakout phase and the 55 hr that follow. The small scale height of our RSG atmosphere model suggests that the breakout signature is a thermal soft X-ray burst (lambda_peak ~ 90\AA) with a duration of <~ 2000 s. Longer durations are possible but require either an extended and tenuous non-standard envelope, or an unusually dense RSG wind with \dot{M} ~ 10^(-3) Msun yr^(-1). The GALEX observations miss the peak of the luminous (M_FUV ~ -20) UV burst but unambiguously capture the rise of the emission and a subsequent 2 day long plateau. The postbreakout, UV-bright plateau is a prediction of our model in which the shift of the peak of the spectral energy distribution (SED) from ~100 to ~1000\AA and the ejecta expansion both counteract the decrease in bolometric luminosity from ~10^11 to ~10^9 L_sun over that period. Based on the observed detection efficiency of our study we make predictions for the breakout detection rate of the GALEX Time Domain Survey.

We present two luminous UV/optical flares from the nuclei of apparently inactive early-type galaxies at z=0.37 and 0.33 that have the radiative properties of a flare from the tidal disruption of a star. In this paper we report the second candidate tidal disruption event discovery in the UV by the GALEX Deep Imaging Survey, and present simultaneous optical light curves from the CFHTLS Deep Imaging Survey for both UV flares. The first few months of the UV/optical light curves are well fitted with the canonical t^(-5/3) power-law decay predicted for emission from the fallback of debris from a tidally disrupted star. Chandra ACIS X-ray observations during the flares detect soft X-ray sources with T_bb= (2-5) x 10^5 K or Gamma > 3 and place limits on hard X-ray emission from an underlying AGN down to L_X (2-10 keV) <~ 10^41 ergs/s. Blackbody fits to the UV/optical spectral energy distributions of the flares indicate peak flare luminosities of > 10^44-10^45 ergs/s. The temperature, luminosity, and light curves of both flares are in excellent agreement with emission from a tidally disrupted main sequence star onto a central black hole of several times 10^7 msun. The observed detection rate of our search over ~ 2.9 deg^2 of GALEX Deep Imaging Survey data spanning from 2003 to 2007 is consistent with tidal disruption rates calculated from dynamical models, and we use these models to make predictions for the detection rates of the next generation of optical synoptic surveys.

We present up to two decades of spectroscopic monitoring of the double-peaked broad Halpha emission lines of seven broad-line radio galaxies. These extremely broad, double-peaked lines are successfully modeled by emission from gas in the outer accretion disk, and their profiles vary on timescales of months to years. This variability is a valuable probe for understanding the structure and dynamics of the accretion disk. We characterize the long-term broad-line profile variability in detail, and compare the evolution of the line profiles with predictions from a few simple physical models for dynamical processes in the accretion disk. We find no evidence for persistent, periodic variability that would be indicative of a precessing elliptical disk, or a circular disk with a long-lived, single-armed spiral or warp. We do find transient, periodic variability on the orbital timescale, possibly resulting from shocks induced by tidal perturbations, and variability in the wings of the profile which suggests changes in the emissivity of the inner accretion flow. Dramatic but localized profile changes are observed during flares and high-states of broad-line flux. In 3C 332, these changes can be explained by a slow, smooth, secular change in disk illumination. In Arp 102B and 3C 390.3, however, a simple disk model cannot explain the asymmetric profile changes, suggesting a need for more sophisticated models. We also observe sharp, transient features that appear in the blue peak of the objects, which require a well-organized velocity field.

A supermassive black hole in the nucleus of a galaxy will be revealed when a star passes close enough to be torn apart by tidal forces and a flare of radiation is emitted by the stream of stellar debris that plunges into the black hole. Since common active galactic nuclei have accreting black holes that can also produce flares, a convincing demonstration that a stellar tidal disruption has occurred generally begins with a ``normal'' galaxy that has no evidence of prior nuclear activity. Here we report a luminous UV flare from an elliptical galaxy at z = 0.37 in the Groth field of the GALEX Deep Imaging Survey that has no evidence of a Seyfert nucleus from optical spectroscopy and X-ray imaging obtained during the flare. Multiwavelength data collected at the time of the event, and for 2 years following, allow us to constrain, for the first time, the spectral energy distribution of a candidate tidal disruption flare from optical through X-rays. The luminosity and temperature of the radiation and the decay curve of the flare are in excellent agreement with theoretical predictions for the tidal disruption of a star, and provide the strongest empirical evidence for a stellar disruption event to date.

Large-amplitude, high-luminosity soft X-ray flares were detected by the ROSAT All-Sky Survey in several galaxies with no evidence of Seyfert activity in their ground-based optical spectra. These flares had the properties predicted for a tidal disruption of a star by a central supermassive black hole. We report Chandra observations of three of these galaxies taken a decade after their flares that reveal weak nuclear X-ray sources that are from 240 to 6000 times fainter than their luminosities at peak, supporting the theory that these were special events and not ongoing active galactic nucleus (AGN) variability. The decline of RX J1624.9+7554 by a factor of 6000 is consistent with the (t-t_D)^(-5/3) decay predicted for the fall-back phase of a tidal disruption event, but only if ROSAT was lucky enough to catch the event exactly at its peak in 1990 October. RX J1242.6-1119A has declined by a factor of 240, also consistent with (t-t_D)^(-5/3). In the H II galaxy NGC 5905 we find only resolved, soft X-ray emission that is undoubtedly associated with starburst activity. When accounting for the starburst component, the ROSAT observations of NGC 5905, as well as the Chandra upper limit on its nuclear flux, are consistent with a (t-t_D)^(-5/3) decay by at least a factor of 1000. Although we found weak Seyfert~2 emission lines in Hubble Space Telescope spectra of NGC 5905, indicating that a low-luminosity AGN was present prior to the X-ray flare, we favor a tidal disruption explanation for the flare itself.

Large amplitude, high luminosity soft X-ray flares were detected by the ROSAT All-Sky Survey in several galaxies with no evidence for Seyfert activity in their ground-based optical spectra. These flares had the properties predicted for a tidal disruption event by a central supermassive black hole: soft X-ray spectrum, time-scale of months, and large X-ray luminosity (10^42 - 10^44 ergs s^-1). In order to evaluate the alternative hypothesis that these flares could have been some form of extreme AGN variability, we obtained follow-up optical spectroscopy of three of the flaring galaxies a decade later with the Space Telescope Imaging Spectrograph (STIS) and a narrow slit to search for or place stringent limits on the presence of any persistent Seyfert-like emission in their nuclei. Two of the galaxies, RX J1624.9+7554 and RX J1242.6-1119, show no evidence for emission lines or a non-stellar continuum in their Hubble Space Telescope (HST) nuclear spectra, consistent with their ground-based classification as inactive galaxies. They remain the most convincing as hosts of tidal disruption events. NGC 5905, previously known as a starburst H II galaxy due to its strong emission lines, has in its inner 0."1 a nucleus with narrow emission-line ratios requiring a Seyfert 2 classification. This weak Seyfert 2 nucleus in NGC 5905, which was masked by the many surrounding H II regions in ground-based spectra, requires a low level of prior non-stellar photoionization, thus raising some uncertainty about the nature of its X-ray flare, which may or may not have been a tidal disruption event. The absence of both broad emission lines and nuclear X-ray absorption in NGC 5905 also characterizes it as a ``true'' Seyfert 2 galaxy, yet one that has varied by more than a factor of 100 in X-rays.

We use unpublished and published VLBI results to investigate the geometry and the statistical properties of the velocity field traced by H2O masers in five galactic regions of star formation -- Sgr B2(M), W49N, W51(MAIN), W51N, and W3(OH). In all sources the angular distribution of the H2O hot spots demonstrates approximate self-similarity (fractality) over almost four orders of magnitude in scale, with the calculated fractal dimension d between (approximately) 0.2 and 1.0. In all sources, the lower order structure functions for the line-of-sight component of the velocity field are satisfactorily approximated by power laws, with the exponents near their classic Kolmogorov values for the high-Reynolds-number incompressible turbulence. These two facts, as well as the observed significant excess of large deviations of the two-point velocity increments from their mean values, strongly suggest that the H2O masers in regions of star formation trace turbulence. We propose a new conceptual model of these masers in which maser hot spots originate at the sites of ultimate dissipation of highly supersonic turbulence produced in the ambient gas by the intensive gas outflow from a newly-born star. Due to the high brightness and small angular sizes of masing hot spots and the possibility of measuring their positions and velocities with high precision, they become a unique probe of supersonic turbulence.

We present K-band $\lambda/\Delta\lambda$ ~ 2600 spectroscopy of five stars (K ~ 14 - 16 mag) within 0.''5 of Sgr A*, the radio source associated with the compact massive object suspected to be a 2.6 x 10$^{6}$ \msun black hole at the center of our Galaxy. High spatial resolution of ~ 0.''09, and good strehl ratios of ~ 0.2 achieved with adaptive optics on the 10-meter Keck telescope make it possible to measure moderate-resolution spectra of these stars individually for the first time. Two stars (S0-17 and S0-18) are identified as late-type stars by the detection of CO bandhead absorption in their spectra. Their absolute K magnitudes and CO bandhead absorption strengths are consistent with early K giants. Three stars (S0-1, S0-2, and S0-16), with r$_{proj}$ $<$ 0.0075 pc (~ 0.''2) from Sgr A*, lack CO bandhead absorption, confirming the results of earlier lower spectral and lower spatial resolution observations that the majority of the stars in the Sgr A* Cluster are early-type stars. The absolute K magnitudes of the early-type stars suggest that they are late O - early B main sequence stars of ages $<$ 20 Myr. The presence of young stars in the Sgr A* Cluster, so close to the central supermassive black hole, poses the intriguing problem of how these stars could have formed, or could have been brought, within its strong tidal field.