A possible close supermassive black-hole binary in a quasar with optical periodicity

Quasars have long been known to be variable sources at all wavelengths. Their optical variability is stochastic, can be due to a variety of physical mechanisms, and is well-described statistically in terms of a damped random walk model. The recent availability of large collections of astronomical time series of flux measurements (light curves) offers new data sets for a systematic exploration of quasar variability. Here we report on the detection of a strong, smooth periodic signal in the optical variability of the quasar PG 1302-102 with a mean observed period of 1,884 $\pm$ 88 days. It was identified in a search for periodic variability in a data set of light curves for 247,000 known, spectroscopically confirmed quasars with a temporal baseline of $\sim9$ years. While the interpretation of this phenomenon is still uncertain, the most plausible mechanisms involve a binary system of two supermassive black holes with a subparsec separation. Such systems are an expected consequence of galaxy mergers and can provide important constraints on models of galaxy formation and evolution.Comment: 19 pages, 6 figures. Published online by Nature on 7 January 201

Optical Discovery of Probable Stellar Tidal Disruption Flares

Using archival Sloan Digital Sky Survey (SDSS) multi-epoch imaging data (Stripe 82), we have searched for the tidal disruption of stars by supermassive black holes in non-active galaxies. Two candidate tidal disruption events (TDEs) are identified. The TDE flares have optical blackbody temperatures of 2 × 10^4 K and observed peak luminosities of M_g = –18.3 and –20.4 (νL_ν = 5 × 10^(42), 4 × 10^(43) erg s^(–1), in the rest frame); their cooling rates are very low, qualitatively consistent with expectations for tidal disruption flares. The properties of the TDE candidates are examined using (1) SDSS imaging to compare them to other flares observed in the search, (2) UV emission measured by GALEX, and (3) spectra of the hosts and of one of the flares. Our pipeline excludes optically identifiable AGN hosts, and our variability monitoring over nine 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 these 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 an SN or an AGN-flare explanation, although an entirely new process specific to the inner few hundred parsecs of non-active galaxies cannot be excluded. Based on our observed rate, we infer that hundreds or thousands of TDEs will be present in current and next-generation optical synoptic surveys. Using the approach outlined here, a TDE candidate sample with O(1) purity can be selected using geometric resolution and host and flare color alone, demonstrating that a campaign to create a large sample of TDEs, with immediate and detailed multi-wavelength follow-up, is feasible. A by-product of this work is quantification of the power spectrum of extreme flares in AGNs

Extreme Variability in a Broad Absorption Line Quasar

CRTS J084133.15+200525.8 is an optically bright quasar at z=2.345 that has shown extreme spectral variability over the past decade. Photometrically, the source had a visual magnitude of V~17.3 between 2002 and 2008. Then, over the following five years, the source slowly brightened by approximately one magnitude, to V~16.2. Only ~1 in 10,000 quasars show such extreme variability, as quantified by the extreme parameters derived for this quasar assuming a damped random walk model. A combination of archival and newly acquired spectra reveal the source to be an iron low-ionization broad absorption line (FeLoBAL) quasar with extreme changes in its absorption spectrum. Some absorption features completely disappear over the 9 years of optical spectra, while other features remain essentially unchanged. We report the first definitive redshift for this source, based on the detection of broad H-alpha in a Keck/MOSFIRE spectrum. Absorption systems separated by several 1000 km/s in velocity show coordinated weakening in the depths of their troughs as the continuum flux increases. We interpret the broad absorption line variability to be due to changes in photoionization, rather than due to motion of material along our line of sight. This source highlights one sort of rare transition object that astronomy will now be finding through dedicated time-domain surveys.Comment: 6 pages, 4 figures; accepted for publication in Ap

Machine-assisted discovery of relationships in astronomy

High-volume feature-rich data sets are becoming the bread-and-butter of 21st century astronomy but present significant challenges to scientific discovery. In particular, identifying scientifically significant relationships between sets of parameters is non-trivial. Similar problems in biological and geosciences have led to the development of systems which can explore large parameter spaces and identify potentially interesting sets of associations. In this paper, we describe the application of automated discovery systems of relationships to astronomical data sets, focusing on an evolutionary programming technique and an information-theory technique. We demonstrate their use with classical astronomical relationships – the Hertzsprung–Russell diagram and the Fundamental Plane of elliptical galaxies. We also show how they work with the issue of binary classification which is relevant to the next generation of large synoptic sky surveys, such as the Large Synoptic Survey Telescope (LSST). We find that comparable results to more familiar techniques, such as decision trees, are achievable. Finally, we consider the reality of the relationships discovered and how this can be used for feature selection and extraction

Connecting the time domain community with the Virtual Astronomical Observatory

The time domain has been identified as one of the most important areas of astronomical research for the next decade. The Virtual Observatory is in the vanguard with dedicated tools and services that enable and facilitate the discovery, dissemination and analysis of time domain data. These range in scope from rapid notifications of time-critical astronomical transients to annotating long-term variables with the latest modeling results. In this paper, we will review the prior art in these areas and focus on the capabilities that the VAO is bringing to bear in support of time domain science. In particular, we will focus on the issues involved with the heterogeneous collections of (ancillary) data associated with astronomical transients, and the time series characterization and classification tools required by the next generation of sky surveys, such as LSST and SKA.Comment: Submitted to Proceedings of SPIE Observatory Operations: Strategies, Processes and Systems IV, Amsterdam, 2012 July 2-

Eclipsing binary stars in the era of massive surveys: First results and future prospects

Our thinking about eclipsing binary stars has undergone a tremendous change in the last decade. Eclipsing binary stars are one of nature’s best laboratories for determining the fundamental physical properties of stars and thus for testing the predictions of theoretical models. Some of the largest ongoing variable star surveys include the Catalina Real-time Transient Survey (CRTS) and the VISTA Variables in the Vía Láctea survey (VVV). They both contain a large amount of photometric data and plenty of information about eclipsing binaries that wait to be extracted and exploited. Here we briefly describe our efforts in this direction

A systematic search for close supermassive black hole binaries in the Catalina Real-Time Transient Survey

Hierarchical assembly models predict a population of supermassive black hole (SMBH) binaries. These are not resolvable by direct imaging but may be detectable via periodic variability (or nanohertz frequency gravitational waves). Following our detection of a 5.2 year periodic signal in the quasar PG 1302-102 (Graham et al. 2015), we present a novel analysis of the optical variability of 243,500 known spectroscopically confirmed quasars using data from the Catalina Real-time Transient Survey (CRTS) to look for close (< 0.1 pc) SMBH systems. Looking for a strong Keplerian periodic signal with at least 1.5 cycles over a baseline of nine years, we find a sample of 111 candidate objects. This is in conservative agreement with theoretical predictions from models of binary SMBH populations. Simulated data sets, assuming stochastic variability, also produce no equivalent candidates implying a low likelihood of spurious detections. The periodicity seen is likely attributable to either jet precession, warped accretion disks or periodic accretion associated with a close SMBH binary system. We also consider how other SMBH binary candidates in the literature appear in CRTS data and show that none of these are equivalent to the identified objects. Finally, the distribution of objects found is consistent with that expected from a gravitational wave-driven population. This implies that circumbinary gas is present at small orbital radii and is being perturbed by the black holes. None of the sources is expected to merge within at least the next century. This study opens a new unique window to study a population of close SMBH binaries that must exist according to our current understanding of galaxy and SMBH evolution.Comment: 29 pages, 10 figures, accepted for publication in MNRAS - this version contains extended table and figur

Understanding extreme quasar optical variability with CRTS: I. Major AGN flares

There is a large degree of variety in the optical variability of quasars and it is unclear whether this is all attributable to a single (set of) physical mechanism(s). We present the results of a systematic search for major flares in AGN in the Catalina Real-time Transient Survey as part of a broader study into extreme quasar variability. Such flares are defined in a quantitative manner as being atop of the normal, stochastic variability of quasars. We have identified 51 events from over 900,000 known quasars and high probability quasar candidates, typically lasting 900 days and with a median peak amplitude of $\Delta m = 1.25$ mag. Characterizing the flare profile with a Weibull distribution, we find that nine of the sources are well described by a single-point single-lens model. This supports the proposal by Lawrence et al. (2016) that microlensing is a plausible physical mechanism for extreme variability. However, we attribute the majority of our events to explosive stellar-related activity in the accretion disk: superluminous supernovae, tidal disruption events, and mergers of stellar mass black holes.Comment: 25 pages, 18 figures, accepted for publication by MNRA

Pan-chromatic observations of the remarkable nova LMC 2012

We present the results of an intensive multiwavelength campaign on nova LMC 2012. This nova evolved very rapidly in all observed wavelengths. The time to fall two magnitudes in the V band was only 2 days. In X-rays the super soft phase began 13$\pm$5 days after discovery and ended around day 50 after discovery. During the super soft phase, the \Swift/XRT and \Chandra\ spectra were consistent with the underlying white dwarf being very hot, $\sim$ 1 MK, and luminous, $\sim$ 10$^{38}$ erg s$^{-1}$. The UV, optical, and near-IR photometry showed a periodic variation after the initial and rapid fading had ended. Timing analysis revealed a consistent 19.24$\pm$0.03 hr period in all UV, optical, and near-IR bands with amplitudes of $\sim$ 0.3 magnitudes which we associate with the orbital period of the central binary. No periods were detected in the corresponding X-ray data sets. A moderately high inclination system, $i$ = 60$\pm$10^{\arcdeg}, was inferred from the early optical emission lines. The {\it HST}/STIS UV spectra were highly unusual with only the \ion{N}{5} (1240\AA) line present and superposed on a blue continuum. The lack of emission lines and the observed UV and optical continua from four epochs can be fit with a low mass ejection event, $\sim$ 10$^{-6}$ M$_{\odot}$, from a hot and massive white dwarf near the Chandrasekhar limit. The white dwarf, in turn, significantly illuminated its subgiant companion which provided the bulk of the observed UV/optical continuum emission at the later dates. The inferred extreme white dwarf characteristics and low mass ejection event favor nova LMC 2012 being a recurrent nova of the U Sco subclass.Comment: 18 figures, 6 tables (one online only containing all the photometry