Source author record

S. G. Djorgovski

S. G. Djorgovski appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

103works
17topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

103 published item(s)

preprint2022arXiv

The WISE-2MASS Survey: Red Quasars Into the Radio Quiet Regime

We present a highly complete sample of broad-line (Type 1) QSOs out to z ~ 3 selected by their mid-infrared colors, a method that is minimally affected by dust reddening. We remove host galaxy emission from the spectra and fit for excess reddening in the residual QSOs, resulting in a Gaussian distribution of colors for unreddened (blue) QSOs, with a tail extending toward heavily reddened (red) QSOs, defined as having E(B - V) > 0.25. This radio-independent selection method enables us to compare red and blue QSO radio properties in both the FIRST (1.4 GHz) and VLASS (2 - 4 GHz) surveys. Consistent with recent results from optically-selected QSOs from SDSS, we find that red QSOs have a significantly higher detection fraction and a higher fraction of compact radio morphologies at both frequencies. We employ radio stacking to investigate the median radio properties of the QSOs including those that are undetected in FIRST and VLASS, finding that red QSOs have significantly brighter radio emission and steeper radio spectral slopes compared with blue QSOs. Finally, we find that the incidence of red QSOs is strongly luminosity dependent, where red QSOs make up > 40% of all QSOs at the highest luminosities. Overall, red QSOs comprise ~ 40% of higher luminosity QSOs, dropping to only a few percent at lower luminosities. Furthermore, red QSOs make up a larger percentage of the radio-detected QSO population. We argue that dusty AGN-driven winds are responsible for both the obscuration as well as excess radio emission seen in red QSOs.

preprint2021arXiv

Gaia GraL: Gaia DR2 Gravitational Lens Systems. VII. XMM-Newton Observations of Lensed Quasars

We present XMM-Newton X-ray observations of nine confirmed lensed quasars at $1 \lesssim z \lesssim 3$ identified by the Gaia Gravitational Lens program. Eight systems are strongly detected, with 0.3--8.0 keV fluxes $F_{0.3-8.0} \gtrsim 5 \times 10^{-14}\ {\rm erg}\ {\rm cm}^{-2}\ {\rm s}^{-1}$. Modeling the X-ray spectra with an absorbed power law, we derive power law photon indices and 2--10 keV luminosities for the eight detected quasars. In addition to presenting sample properties for larger quasar population studies and for use in planning for future caustic crossing events, we also identify three quasars of interest: a quasar that shows evidence of flux variability from previous ROSAT observations, the most closely-separated individual lensed sources resolved by XMM-Newton, and one of the X-ray brightest quasars known at $z>3$. These sources represent the tip of discovery that will be enabled by SRG/eROSITA.

preprint2020arXiv

Candidate Electromagnetic Counterpart to the Binary Black Hole Merger Gravitational Wave Event S190521g

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

preprint2020arXiv

Chandra Observations of Candidate Sub-Parsec Binary Supermassive Black Holes

We present analysis of Chandra X-ray observations of seven quasars that were identified as candidate sub-parsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey (CRTS) based on apparent periodicity in their optical light curves. Simulations predict close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs, including harder or softer X-ray spectra, ripple-like profiles in the Fe K-$α$ line, and distinct peaks in the spectrum due to the separation of the accretion disk into a circumbinary disk and mini-disks around each SMBH. We obtained Chandra observations to test these models and assess whether these quasars could contain binary SMBHs. We instead find that the quasar spectra are all well fit by simple absorbed power law models, with the rest frame 2-10 keV photon indices, $Γ$, and the X-ray-to-optical power slopes, $α_{\rm OX}$, indistinguishable from the larger quasar population. This may indicate that these seven quasars are not truly sub-parsec binary SMBH systems, or it may simply reflect that our sample size was too small to robustly detect any differences. Alternatively, the X-ray spectral changes might only be evident at higher energies than probed by Chandra. Given the available models and current data, no firm conclusions are drawn. These observations will help motivate and direct further work on theoretical models of binary SMBH systems, such as modeling systems with thinner accretion disks and larger binary separations.

preprint2020arXiv

Multiple Chemodynamic Stellar Populations of the Ursa Minor Dwarf Spheroidal Galaxy

We present a Bayesian method to identify multiple (chemodynamic) stellar populations in dwarf spheroidal galaxies (dSphs) using velocity, metallicity, and positional stellar data without the assumption of spherical symmetry. We apply this method to a new Keck/DEIMOS spectroscopic survey of the Ursa Minor (UMi) dSph. We identify 892 likely members, making this the largest UMi sample with line-of-sight velocity and metallicity measurements. Our Bayesian method detects two distinct chemodynamic populations with high significance ($\ln{B}\sim33$). The metal-rich ($[{\rm Fe/H}]=-2.05\pm0.03$) population is kinematically colder (radial velocity dispersion of $σ_v=4.9\pm0.8 \, {\rm km \, s^{-1}}$) and more centrally concentrated than the metal-poor ($[{\rm Fe/H}]=-2.29\pm0.05$) and kinematically hotter population ($σ_v =11.5\pm0.9\, {\rm km \, s^{-1}}$). Furthermore, we apply the same analysis to an independent MMT/Hectochelle data set and confirm the existence of two chemodynamic populations in UMi. In both data sets, the metal-rich population is significantly flattened ($ε=0.75\pm0.03$) and the metal-poor population is closer to spherical ($ε=0.33_{-0.09}^{+0.12}$). Despite the presence of two populations, we are unable to robustly estimate the slope of the dynamical mass profile. We found hints for prolate rotation of order $\sim 2 \, {\rm km \, s^{-1}}$ in the MMT data set, but further observations are required to verify this. The flattened metal-rich population invalidates assumptions built into simple dynamical mass estimators, so we computed new astrophysical dark matter annihilation (J) and decay profiles based on the rounder, hotter metal-poor population and inferred $\log_{10}{(J(0.5^{\circ})/{\rm GeV^{2} \, cm^{-5}})}\approx19.1$ for the Keck data set. Our results paint a more complex picture of the evolution of Ursa Minor than previously discussed.

preprint2019arXiv

Understanding extreme quasar optical variability with CRTS: II. Changing-state quasars

We present the results of a systematic search for quasars in the Catalina Real-time Transient Survey exhibiting both strong photometric and spectroscopic variability over a decadal baseline. We identify 73 sources with specific patterns of optical and mid-IR photometric behavior and a defined spectroscopic change. These "Changing-State" quasars (CSQs) form a higher luminosity sample to complement existing sets of "Changing-Look" AGN and quasars in the literature. The CSQs (by selection) exhibit larger photometric variability than the CLQs. The spectroscopic variability is marginally stronger in the CSQs than CLQs as defined by the change in H$β$/[OIII] ratio. We find 36 sources with declining H$β$ flux, 37 sources with increasing H$β$ flux and discover seven sources with $z > 0.8$, further extending the redshift arm. Our CSQ sample compares to the literature CLQ objects in similar distributions of H$β$ flux ratios and differential Eddington ratios between high (bright) and low (dim) states. Taken as a whole, we find that this population of extreme varying quasars is associated with changes in the Eddington ratio and the timescales imply cooling/heating fronts propagating through the disk.

preprint2016arXiv

Interacting supernovae and supernova impostors. LSQ13zm: an outburst heralds the death of a massive star

We report photometric and spectroscopic observations of the optical transient LSQ13zm. Historical data reveal the presence of an eruptive episode (that we label as `2013a') followed by a much brighter outburst (`2013b') three weeks later, that we argue to be the genuine supernova explosion. This sequence of events closely resemble those observed for SN2010mc and (in 2012) SN2009ip. The absolute magnitude reached by LSQ13zm during 2013a ($M_R=-14.87\pm0.25\,\rm{mag}$) is comparable with those of supernova impostors, while that of the 2013b event ($M_R=-18.46\pm0.21\,\rm{mag}$) is consistent with those of interacting supernovae. Our spectra reveal the presence of a dense and structured circumstellar medium, probably produced through numerous pre-supernova mass-loss events. In addition, we find evidence for high-velocity ejecta, with a fraction of gas expelled at more than 20000\kms. The spectra of LSQ13zm show remarkable similarity with those of well-studied core-collapse supernovae. From the analysis of the available photometric and spectroscopic data, we conclude that we first observed the last event of an eruptive sequence from a massive star, likely a Luminous Blue Variable, which a short time later exploded as a core-collapse supernova. The detailed analysis of archival images suggest that the host galaxy is a star-forming Blue Dwarf Compact Galaxy.

preprint2016arXiv

Real-Time Data Mining of Massive Data Streams from Synoptic Sky Surveys

The nature of scientific and technological data collection is evolving rapidly: data volumes and rates grow exponentially, with increasing complexity and information content, and there has been a transition from static data sets to data streams that must be analyzed in real time. Interesting or anomalous phenomena must be quickly characterized and followed up with additional measurements via optimal deployment of limited assets. Modern astronomy presents a variety of such phenomena in the form of transient events in digital synoptic sky surveys, including cosmic explosions (supernovae, gamma ray bursts), relativistic phenomena (black hole formation, jets), potentially hazardous asteroids, etc. We have been developing a set of machine learning tools to detect, classify and plan a response to transient events for astronomy applications, using the Catalina Real-time Transient Survey (CRTS) as a scientific and methodological testbed. The ability to respond rapidly to the potentially most interesting events is a key bottleneck that limits the scientific returns from the current and anticipated synoptic sky surveys. Similar challenge arise in other contexts, from environmental monitoring using sensor networks to autonomous spacecraft systems. Given the exponential growth of data rates, and the time-critical response, we need a fully automated and robust approach. We describe the results obtained to date, and the possible future developments.

preprint2016arXiv

The velocity dispersion and mass function of the outer halo globular cluster Palomar 4

We obtained precise line-of-sight radial velocities of 23 member stars of the remote halo globular cluster Palomar 4 (Pal 4) using the High Resolution Echelle Spectrograph (HIRES) at the Keck I telescope. We also measured the mass function of the cluster down to a limiting magnitude of V~28 mag using archival HST/WFPC2 imaging. We derived the cluster's surface brightness profile based on the WFPC2 data and on broad-band imaging with the Low-Resolution Imaging Spectrometer (LRIS) at the Keck II telescope. We find a mean cluster velocity of 72.55+/-0.22 km/s and a velocity dispersion of 0.87+/-0.18 km/s. The global mass function of the cluster, in the mass range 0.55<=M<=0.85 M_solar, is shallower than a Kroupa mass function and the cluster is significantly depleted in low-mass stars in its center compared to its outskirts. Since the relaxation time of Pal 4 is of the order of a Hubble time, this points to primordial mass segregation in this cluster. Extrapolating the measured mass function towards lower-mass stars and including the contribution of compact remnants, we derive a total cluster mass of 29800 M_solar. For this mass, the measured velocity dispersion is consistent with the expectations of Newtonian dynamics and below the prediction of Modified Newtonian Dynamics (MOND). Pal 4 adds to the growing body of evidence that the dynamics of star clusters in the outer Galactic halo can hardly be explained by MOND.

preprint2015arXiv

A Megacam Survey of Outer Halo Satellites. IV. Two foreground populations possibly associated with the Monoceros substructure in the direction of NGC2419 and Koposov2

The origin of the Galactic halo stellar structure known as the Monoceros ring is still under debate. In this work, we study that halo substructure using deep CFHT wide-field photometry obtained for the globular clusters NGC2419 and Koposov2, where the presence of Monoceros becomes significant because of their coincident projected position. Using Sloan Digital Sky Survey photometry and spectroscopy in the area surrounding these globulars and beyond, where the same Monoceros population is detected, we conclude that a second feature, not likely to be associated with Milky Way disk stars along the line-of-sight, is present as foreground population. Our analysis suggests that the Monoceros ring might be composed of an old stellar population of age t ~ 9Gyr and a new component ~ 4Gyr younger at the same heliocentric distance. Alternatively, this detection might be associated with a second wrap of Monoceros in that direction of the sky and also indicate a metallicity spread in the ring. The detection of such a low-density feature in other sections of this halo substructure will shed light on its nature.

preprint2015arXiv

A serendipitous all sky survey for bright objects in the outer solar system

We use seven year's worth of observations from the Catalina Sky Survey and the Siding Spring Survey covering most of the northern and southern hemisphere at galactic latitudes higher than 20 degrees to search for serendipitously imaged moving objects in the outer solar system. These slowly moving objects would appear as stationary transients in these fast cadence asteroids surveys, so we develop methods to discover objects in the outer solar system using individual observations spaced by months, rather than spaced by hours, as is typically done. While we independently discover 8 known bright objects in the outer solar system, the faintest having $V=19.8\pm0.1$, no new objects are discovered. We find that the survey is nearly 100% efficient at detecting objects beyond 25 AU for $V\lesssim 19.1$ ($V\lesssim18.6$ in the southern hemisphere) and that the probability that there is one or more remaining outer solar system object of this brightness left to be discovered in the unsurveyed regions of the galactic plane is approximately 32%.

preprint2015arXiv

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.

preprint2015arXiv

Binary Active Galactic Nuclei in Stripe 82: Constraints on Synchronized Black Hole Accretion in Major Mergers

Representing simultaneous black hole accretion during a merger, binary active galactic nuclei (AGNs) could provide valuable observational constraints to models of galaxy mergers and AGN triggering. High-resolution radio interferometer imaging offers a promising method to identify a large and uniform sample of binary AGNs, because it probes a generic feature of nuclear activity and is free from dust obscuration. Our previous search yielded 52 strong candidates of kpc-scale binaries over the 92 deg$^2$ of the Sloan Digital Sky Survey (SDSS) Stripe 82 area with 2"-resolution Very Large Array (VLA) images. Here we present 0.3"-resolution VLA 6 GHz observations for six candidates that have complete optical spectroscopy. The new data confirm the binary nature of four candidates and identify the other two as line-of-sight projections of radio structures from single AGNs. The four binary AGNs at $z \sim 0.1$ reside in major mergers with projected separations of 4.2-12 kpc. Optical spectral modeling shows that their hosts have stellar masses between $10.3 < log(M_\star/M_\odot) < 11.5$ and velocity dispersions between $120 < σ_\star < 320$ km/s. The radio emission is compact ($<$0.4") and show steep spectrum ($-1.8 < α< -0.5$) at 6 GHz. The host galaxy properties and the Eddington-scaled accretion rates broadly correlate with the excitation state, similar to the general radio-AGN population at low redshifts. Our estimated binary AGN fraction indicates that simultaneous accretion occurs $>23^{+15}_{-8}$% of the time when a kpc-scale galaxy pair is detectable as a radio-AGN. The high duty cycle of the binary phase strongly suggests that major mergers can trigger and synchronize black hole accretion.

preprint2015arXiv

From Stars to Patients: Lessons from Space Science and Astrophysics for Health Care Informatics

Big Data are revolutionizing nearly every aspect of the modern society. One area where this can have a profound positive societal impact is the field of Health Care Informatics (HCI), which faces many challenges. The key idea behind this study is: can we use some of the experience and technical and methodological solutions from the fields that have successfully adapted to the Big Data era, namely astronomy and space science, to help accelerate the progress of HCI? We illustrate this with examples from the Virtual Observatory framework, and the NCI EDRN project. An effective sharing and reuse of tools, methods, and experiences from different fields can save a lot of effort, time, and expense. HCI can thus benefit from the proven solutions to big data challenges from other domains.

preprint2015arXiv

Infrared Time Lags for the Periodic Quasar PG 1302-102

The optical light curve of the quasar PG 1302-102 at $z = 0.278$ shows a strong, smooth 5.2 yr periodic signal, detectable over a period of $\sim 20$ yr. Although 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. At this close separation, the nuclear black holes in PG 1302-102 will likely merge within $\sim 10^{5}$ yr due to gravitational wave emission alone. Here we report the rest-frame near-infrared time lags for PG 1302-102. Compiling data from {\it WISE} and {\it Akari}, we confirm that the periodic behavior reported in the optical light curve from Graham et al. (2015) is reproduced at infrared wavelengths, with best-fit observed-frame 3.4 and $4.6 μ$m time lags of $(2219 \pm 153, 2408 \pm 148)$ days for a near face-on orientation of the torus, or $(4103\pm 153, 4292 \pm 148)$ days for an inclined system with relativistic Doppler boosting in effect. The periodicity in the infrared light curves and the light-travel time of the accretion disk photons to reach the dust glowing regions support that a source within the accretion disk is responsible for the optical variability of PG 1302-102, echoed at the further out dusty regions. The implied distance of this dusty, assumed toroidal region is $\sim$ 1.5 pc for a near face-on geometry, or $\sim$1.1 pc for the relativistic Doppler boosted case.

preprint2015arXiv

Properties and Evolution of the Redback Millisecond Pulsar Binary PSR J2129-0429

PSR J2129-0429 is a "redback" eclipsing millisecond pulsar binary with an unusually long 15.2 hour orbit. It was discovered by the Green Bank Telescope in a targeted search of unidentified Fermi gamma-ray sources. The pulsar companion is optically bright (mean $m_R = 16.6$ mag), allowing us to construct the longest baseline photometric dataset available for such a system. We present ten years of archival and new photometry of the companion from LINEAR, CRTS, PTF, the Palomar 60-inch, and LCOGT. Radial velocity spectroscopy using the Double-Beam Spectrograph on the Palomar 200-inch indicates that the pulsar is massive: $1.74\pm0.18 M_\odot$. The G-type pulsar companion has mass $0.44\pm0.04 M_\odot$, one of the heaviest known redback companions. It is currently 95\% Roche-lobe filling and only mildly irradiated by the pulsar. We identify a clear 13.1 mmag yr$^{-1}$ secular decline in the mean magnitude of the companion as well as smaller-scale variations in the optical lightcurve shape. This behavior may indicate that the companion is cooling. Binary evolution calculations indicate that PSR J2129-0429 has an orbital period almost exactly at the bifurcation period between systems that converge into tighter orbits as black widows and redbacks and those that diverge into wider pulsar--white dwarf binaries. Its eventual fate may depend on whether it undergoes future episodes of mass transfer and increased irradiation.

preprint2014arXiv

Automated Real-Time Classification and Decision Making in Massive Data Streams from Synoptic Sky Surveys

The nature of scientific and technological data collection is evolving rapidly: data volumes and rates grow exponentially, with increasing complexity and information content, and there has been a transition from static data sets to data streams that must be analyzed in real time. Interesting or anomalous phenomena must be quickly characterized and followed up with additional measurements via optimal deployment of limited assets. Modern astronomy presents a variety of such phenomena in the form of transient events in digital synoptic sky surveys, including cosmic explosions (supernovae, gamma ray bursts), relativistic phenomena (black hole formation, jets), potentially hazardous asteroids, etc. We have been developing a set of machine learning tools to detect, classify and plan a response to transient events for astronomy applications, using the Catalina Real-time Transient Survey (CRTS) as a scientific and methodological testbed. The ability to respond rapidly to the potentially most interesting events is a key bottleneck that limits the scientific returns from the current and anticipated synoptic sky surveys. Similar challenge arise in other contexts, from environmental monitoring using sensor networks to autonomous spacecraft systems. Given the exponential growth of data rates, and the time-critical response, we need a fully automated and robust approach. We describe the results obtained to date, and the possible future developments.

preprint2014arXiv

Cataclysmic Variables from the Catalina Real-time Transient Survey

We present 855 cataclysmic variable candidates detected by the Catalina Real-time Transient Survey (CRTS) of which at least 137 have been spectroscopically confirmed and 705 are new discoveries. The sources were identified from the analysis of five years of data, and come from an area covering three quarters of the sky. We study the amplitude distribution of the dwarf novae CVs discovered by CRTS during outburst, and find that in quiescence they are typically two magnitudes fainter compared to the spectroscopic CV sample identified by SDSS. However, almost all CRTS CVs in the SDSS footprint have ugriz photometry. We analyse the spatial distribution of the CVs and find evidence that many of the systems lie at scale heights beyond those expected for a Galactic thin disc population. We compare the outburst rates of newly discovered CRTS CVs with the previously known CV population, and find no evidence for a difference between them. However, we find that significant evidence for a systematic difference in orbital period distribution. We discuss the CVs found below the orbital period minimum and argue that many more are yet to be identified among the full CRTS CV sample. We cross-match the CVs with archival X-ray catalogs and find that most of the systems are dwarf novae rather than magnetic CVs.

preprint2014arXiv

Connection between optical and gamma-ray variability in blazars

We use optical data from the Palomar Transient Factory (PTF) and the Catalina Real-Time Transient Survey (CRTS) to study the variability of gamma-ray detected and non-detected objects in a large population of active galactic nuclei (AGN) selected from the Candidate Gamma-Ray Blazar Survey and Fermi Gamma-Ray Space Telescope catalogs. Our samples include 714 sources with PTF data and 1244 sources with CRTS data. We calculate the intrinsic modulation index to quantify the optical variability amplitude in these samples. We find the gamma-ray detected objects to be more variable than the non-detected ones. The flat spectrum radio quasars (FSRQs) are more variable than the BL Lac objects in our sample, but the significance of the difference depends on the sample used. When dividing the objects based on their synchrotron peak frequency, we find the low synchrotron peaked (LSP) objects to be significantly more variable than the high synchrotron peaked (HSP) ones, explaining the difference between the FSRQs and BL Lacs. This could be due to the LSPs being observed near their electron energy peak, while in the HSPs the emission is caused by lower energy electrons, which cool more slowly. We also find a significant correlation between the optical and gamma-ray fluxes that is stronger in the HSP BL Lacs than in the FSRQs. The FSRQs in our sample are also more Compton dominated than the HSP BL Lacs. These findings are consistent with models where the gamma-ray emission of HSP objects is produced by the synchrotron self-Compton mechanism, while the LSP objects need an additional external Compton component that increases the scatter in the flux-flux correlation.

preprint2014arXiv

Data Driven Discovery in Astrophysics

We review some aspects of the current state of data-intensive astronomy, its methods, and some outstanding data analysis challenges. Astronomy is at the forefront of "big data" science, with exponentially growing data volumes and data rates, and an ever-increasing complexity, now entering the Petascale regime. Telescopes and observatories from both ground and space, covering a full range of wavelengths, feed the data via processing pipelines into dedicated archives, where they can be accessed for scientific analysis. Most of the large archives are connected through the Virtual Observatory framework, that provides interoperability standards and services, and effectively constitutes a global data grid of astronomy. Making discoveries in this overabundance of data requires applications of novel, machine learning tools. We describe some of the recent examples of such applications.

preprint2014arXiv

Discovery of $\sim$ 9,000 new RR Lyrae in the Southern Catalina Surveys

We present the results of a deep, wide-area variability survey in the Southern hemisphere, the first of its kind. As part of the Catalina Sky Surveys, the Siding Spring Survey (SSS) has covered $14,800$ square degrees in the declination range of $-75^{\circ}\leqδ\leq-15^{\circ}$. To mine the enormous SSS dataset efficiently we have developed two algorithms: Automatic Period Selection (APS) and Automatic Fourier Decomposition (AFD), which aim to sharpen the period estimation and produce robust lightcurve models. Armed with the APS and AFD outputs we classify $10,540$ ab-type RR Lyrae (RRab) stars ($\sim$90% of which are new) across the Southern sky. As well as the positional information we supply photometric metallicities, and unreddened distances. For the RRab stars in the halo, a study of the photometric metallicity distribution reveals a nearly Gaussian shape with a mean metallicity of ${\rm [Fe/H]}=-1.4$ dex and a dispersion of $0.3$ dex. A spatial study of the RRab metallicities shows no significant radial gradient in the first $\sim7$ kpc from the Galaxy center. However, further out, a small negative gradient is clearly present. This is complemented by a very obvious correlation of the mean RR Lyrae metallicity with distance above the Galactic plane, $z$. We have also carried out an initial substructure search using the discovered RRab, and present the properties of the candidates with significance greater than $2 σ$. Most prominent among these is a southern extension of the Sagittarius dwarf galaxy's stream system, reaching down to declinations $\sim -40°$.

preprint2014arXiv

Immersive and Collaborative Data Visualization Using Virtual Reality Platforms

Effective data visualization is a key part of the discovery process in the era of big data. It is the bridge between the quantitative content of the data and human intuition, and thus an essential component of the scientific path from data into knowledge and understanding. Visualization is also essential in the data mining process, directing the choice of the applicable algorithms, and in helping to identify and remove bad data from the analysis. However, a high complexity or a high dimensionality of modern data sets represents a critical obstacle. How do we visualize interesting structures and patterns that may exist in hyper-dimensional data spaces? A better understanding of how we can perceive and interact with multi dimensional information poses some deep questions in the field of cognition technology and human computer interaction. To this effect, we are exploring the use of immersive virtual reality platforms for scientific data visualization, both as software and inexpensive commodity hardware. These potentially powerful and innovative tools for multi dimensional data visualization can also provide an easy and natural path to a collaborative data visualization and exploration, where scientists can interact with their data and their colleagues in the same visual space. Immersion provides benefits beyond the traditional desktop visualization tools: it leads to a demonstrably better perception of a datascape geometry, more intuitive data understanding, and a better retention of the perceived relationships in the data.

preprint2014arXiv

Multifrequency Studies of the Peculiar Quasar 4C +21.35 During the 2010 Flaring Activity

The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) gamma-ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE; E > 100 MeV) gamma-rays, poses intriguing questions on the location of the gamma-ray emitting region in this flat spectrum radio quasar. We present multifrequency data of 4C +21.35 collected from centimeter to VHE during 2010 to investigate the properties of this source and discuss a possible emission model. The first hint of detection at VHE was observed by MAGIC on 2010 May 3, soon after a gamma-ray flare detected by Fermi-LAT that peaked on April 29. The same emission mechanism may therefore be responsible for both the HE and VHE emission during the 2010 flaring episodes. Two optical peaks were detected on 2010 April 20 and June 30, close in time but not simultaneous with the two gamma-ray peaks, while no clear connection was observed between the X-ray an gamma-ray emission. An increasing flux density was observed in radio and mm bands from the beginning of 2009, in accordance with the increasing gamma-ray activity observed by Fermi-LAT, and peaking on 2011 January 27 in the mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C +21.35 for the two periods of the VHE detection and a quiescent state, using a one-zone model with the emission coming from a very compact region outside the broad line region. The three SEDs can be fit with a combination of synchrotron self-Compton and external Compton emission of seed photons from a dust torus, changing only the electron distribution parameters between the epochs. The fit of the optical/UV part of the spectrum for 2010 April 29 seems to favor an inner disk radius of <6 gravitational radii, as one would expect from a prograde-rotating Kerr black hole.

preprint2014arXiv

One thousand cataclysmic variables from the Catalina Real-time Transient Survey

Over six years of operation, the Catalina Real-time Transient Survey (CRTS) has identified 1043 cataclysmic variable (CV) candidates --- the largest sample of CVs from a single survey to date. Here we provide spectroscopic identification of 85 systems fainter than g<19, including three AMCVn binaries, one helium-enriched CV, one polar and one new eclipsing CV. We analyse the outburst properties of the full sample and show that it contains a large fraction of low accretion rate CVs with long outburst recurrence times. We argue that most of the high accretion rate dwarf novae in the survey footprint have already been found and that future CRTS discoveries will be mostly low accretion rate systems. We find that CVs with white dwarf dominated spectra have significantly fewer outbursts in their CRTS light curves compared to disc-dominated CVs, reflecting the difference in their accretion rates. Comparing the CRTS sample to other samples of CVs, we estimate the overall external completeness to be 23.6 per cent, but show that as much as 56 per cent of CVs have variability amplitudes that are too small to be selected using the transient selection criteria employed by current ground-based surveys. The full table of CRTS CVs, including their outburst and spectroscopic properties examined in this paper, is provided in the online materials.

preprint2014arXiv

Radio-Selected Binary Active Galactic Nuclei from the Very Large Array Stripe 82 Survey

Galaxy mergers play an important role in the growth of galaxies and their supermassive black holes. Simulations suggest that tidal interactions could enhance black hole accretion, which can be tested by the fraction of binary active galactic nuclei (AGNs) among galaxy mergers. But determining the fraction requires a statistical sample of binaries. We have identified kpc-scale binary AGNs directly from high-resolution radio imaging. Inside the 92 square deg covered by the high-resolution Very Large Array survey of the Sloan Digital Sky Survey (SDSS) Stripe 82 field, we identified 22 grade A and 30 grade B candidates of binary radio AGNs with angular separations less than 5" (10 kpc at z = 0.1). Eight of the candidates have optical spectra for both components from the SDSS spectroscopic surveys and our Keck program. Two grade B candidates are projected pairs, but the remaining six candidates are all compelling cases of binary AGNs based on either emission line ratios or the excess in radio power compared to the H-alpha-traced star formation rate. Only two of the six binaries were previously discovered by an optical spectroscopic search. Based on these results, we estimate that ~60% of our binary candidates would be confirmed once we obtain complete spectroscopic information. We conclude that wide-area high-resolution radio surveys offer an efficient method to identify large samples of binary AGNs. These radio-selected binary AGNs complement binaries identified at other wavelengths and are useful for understanding the triggering mechanisms of black hole accretion.

preprint2014arXiv

The Catalina Surveys Periodic Variable Star Catalog

We present ~47,000 periodic variables found during the analysis of 5.4 million variable star candidates within a 20,000 square degree region covered by the Catalina Surveys Data Release-1 (CSDR1). Combining these variables with type-ab RR Lyrae from our previous work, we produce an on-line catalog containing periods, amplitudes, and classifications for ~61,000 periodic variables. By cross-matching these variables with those from prior surveys, we find that > 90% of the ~8,000 known periodic variables in the survey region are recovered. For these sources we find excellent agreement between our catalog and prior values of luminosity, period and amplitude, as well as classification. We investigate the rate of confusion between objects classified as contact binaries and type-c RR Lyrae (RRc's) based on periods, colours, amplitudes, metalicities, radial velocities and surface gravities. We find that no more than few percent of these variables in these classes are misidentified. By deriving distances for this clean sample of ~5,500 RRc's, we trace the path of the Sagittarius tidal streams within the Galactic halo. Selecting 146 outer-halo RRc's with SDSS radial velocities, we confirm the presence of a coherent halo structure that is inconsistent with current N-body simulations of the Sagittarius tidal stream. We also find numerous long-period variables that are very likely associated within the Sagittarius tidal streams system. Based on the examination of 31,000 contact binary light curves we find evidence for two subgroups exhibiting irregular lightcurves. One subgroup presents significant variations in mean brightness that are likely due to chromospheric activity. The other subgroup shows stable modulations over more than a thousand days and thereby provides evidence that the O'Connell effect is not due to stellar spots.

preprint2014arXiv

The supernova CSS121015:004244+132827: a clue for understanding super-luminous supernovae

We present optical photometry and spectra of the super luminous type II/IIn supernova CSS121015:004244+132827 (z=0.2868) spanning epochs from -30 days (rest frame) to more than 200 days after maximum. CSS121015 is one of the more luminous supernova ever found and one of the best observed. The photometric evolution is characterized by a relatively fast rise to maximum (~40 days in the SN rest frame), and by a linear post-maximum decline. The light curve shows no sign of a break to an exponential tail. A broad Halpha is first detected at ~ +40 days (rest-frame). Narrow, barely-resolved Balmer and [O III] 5007 A lines, with decreasing strength, are visible along the entire spectral evolution. The spectra are very similar to other super luminous supernovae (SLSNe) with hydrogen in their spectrum, and also to SN 2005gj, sometimes considered a type Ia interacting with H-rich CSM. The spectra are also similar to a subsample of H-deficient SLSNe. We propose that the properties of CSS121015 are consistent with the interaction of the ejecta with a massive, extended, opaque shell, lost by the progenitor decades before the final explosion, although a magnetar powered model cannot be excluded. Based on the similarity of CSS121015 with other SLSNe (with and without H), we suggest that the shocked-shell scenario should be seriously considered as a plausible model for both types of SLSN.

preprint2014arXiv

Ultra-short Period Binaries from the Catalina Surveys

We investigate the properties of 367 ultra-short period binary candidates selected from 31,000 sources recently identified from Catalina Surveys data. Based on light curve morphology, along with WISE, SDSS and GALEX multi-colour photometry, we identify two distinct groups of binaries with periods below the 0.22 day contact binary minimum. In contrast to most recent work, we spectroscopically confirm the existence of M-dwarf+M-dwarf contact binary systems. By measuring the radial velocity variations for five of the shortest-period systems, we find examples of rare cool-white dwarf+M-dwarf binaries. Only a few such systems are currently known. Unlike warmer white dwarf systems, their UV flux and their optical colours and spectra are dominated by the M-dwarf companion. We contrast our discoveries with previous photometrically-selected ultra-short period contact binary candidates, and highlight the ongoing need for confirmation using spectra and associated radial velocity measurements. Overall, our analysis increases the number of ultra-short period contact binary candidates by more than an order of magnitude.

preprint2013arXiv

A comparison of period finding algorithms

This paper presents a comparison of popular period finding algorithms applied to the light curves of variable stars from the Catalina Real-time Transient Survey (CRTS), MACHO and ASAS data sets. We analyze the accuracy of the methods against magnitude, sampling rates, quoted period, quality measures (signal-to-noise and number of observations), variability, and object classes. We find that measure of dispersion-based techniques - analysis-of-variance with harmonics and conditional entropy - consistently give the best results but there are clear dependencies on object class and light curve quality. Period aliasing and identifying a period harmonic also remain significant issues. We consider the performance of the algorithms and show that a new conditional entropy-based algorithm is the most optimal in terms of completeness and speed. We also consider a simple ensemble approach and find that it performs no better than individual algorithms.

preprint2013arXiv

A Megacam Survey of Outer Halo Satellites. II. Blue Stragglers in the Lowest Stellar Density Systems

We present a homogeneous study of blue straggler stars across ten outer halo globular clusters, three classical dwarf spheroidal and nine ultra-faint galaxies based on deep and wide-field photometric data taken with MegaCam on the Canada-France-Hawaii Telescope. We find blue straggler stars to be ubiquitous among these Milky Way satellites. Based on these data, we can test the importance of primordial binaries or multiple systems on blue straggler star formation in low density environments. For the outer halo globular clusters we find an anti-correlation between the specific frequency of blue straggler and absolute magnitude, similar to that previously observed for inner halo clusters. When plotted against density and encounter rate, the frequency of blue stragglers are well fitted by single trends with smooth transitions between dwarf galaxies and globular clusters, which points to a common origin for their blue stragglers. The fraction of blue stragglers stays constant and high in the low encounter rate regime spanned by our dwarf galaxies, and decreases with density and encounter rate in the range spanned by our globular clusters. We find that young stars can mimic blue stragglers in dwarf galaxies only if their ages are 2.5+/-0.5 Gyr and they represent ~1-7% of the total number of stars, which we deem highly unlikely. These results point to mass-transfer or mergers of primordial binaries or multiple systems as the dominant blue straggler formation mechanism in low density systems.

preprint2013arXiv

A novel variability-based method for quasar selection: evidence for a rest frame ~54 day characteristic timescale

We compare quasar selection techniques based on their optical variability using data from the Catalina Real-time Transient Survey (CRTS). We introduce a new technique based on Slepian wavelet variance (SWV) that shows comparable or better performance to structure functions and damped random walk models but with fewer assumptions. Combining these methods with WISE mid-IR colors produces a highly efficient quasar selection technique which we have validated spectroscopically. The SWV technique also identifies characteristic timescales in a time series and we find a characteristic rest frame timescale of ~54 days, confirmed in the light curves of ~18000 quasars from CRTS, SDSS and MACHO data, and anticorrelated with absolute magnitude. This indicates a transition between a damped random walk and $P(f) \propto f^{-1/3}$ behaviours and is the first strong indication that a damped random walk model may be too simplistic to describe optical quasar variability.

preprint2013arXiv

Dust Reddened Quasars in FIRST and UKIDSS: Beyond the Tip of the Iceberg

We present the results of a pilot survey to find dust-reddened quasars by matching the FIRST radio catalog to the UKIDSS near-infrared survey, and using optical data from SDSS to select objects with very red colors. The deep K-band limit provided by UKIDSS allows for finding more heavily-reddened quasars at higher redshifts as compared with previous work using FIRST and 2MASS. We selected 87 candidates with K<=17.0 from the UKIDSS Large Area Survey (LAS) First Data Release (DR1) which covers 190 deg2. These candidates reach up to ~1.5 magnitudes below the 2MASS limit and obey the color criteria developed to identify dust-reddened quasars. We have obtained 61 spectroscopic observations in the optical and/or near-infrared as well as classifications in the literature and have identified 14 reddened quasars with E(B-V)>0.1, including three at z>2. We study the infrared properties of the sample using photometry from the WISE Observatory and find that infrared colors improve the efficiency of red quasar selection, removing many contaminants in an infrared-to-optical color-selected sample alone. The highest-redshift quasars (z > 2) are only moderately reddened, with E(B-V) ~ 0.2-0.3. We find that the surface density of red quasars rises sharply with faintness, comprising up to 17% of blue quasars at the same apparent K-band flux limit. We estimate that to reach more heavily reddened quasars (i.e., E(B-V) > 0.5) at z>2 and a depth of K=17 we would need to survey at least ~2.5 times more area.

preprint2013arXiv

Evidence for a Milky Way Tidal Stream Reaching Beyond 100 kpc

We present the analysis of 1,207 RR Lyrae found in photometry taken by the Catalina Survey's Mount Lemmon telescope. By combining accurate distances for these stars with measurements for ~14,000 type-AB RR Lyrae from the Catalina Schmid telescope, we reveal an extended association that reaches Galactocentric distances beyond 100 kpc and overlaps the Sagittarius streams system. This result confirms earlier evidence for the existence of an outer halo tidal stream resulting from a disrupted stellar system. By comparing the RR Lyrae source density with that expected based on halo models, we find the detection has ~8 sigma significance. We investigate the distances, radial velocities, metallicities, and period-amplitude distribution of the RR Lyrae. We find that both radial velocities and distances are inconsistent with current models of the Sagittarius stream. We also find tentative evidence for a division in source metallicities for the most distant sources. Following prior analyses, we compare the locations and distances of the RR Lyrae with photometrically selected candidate horizontal branch stars and find supporting evidence that this structure spans at least 60 deg of the sky. We investigate the prospects of an association between the stream and unusual globular cluster NGC 2419.

preprint2013arXiv

Feature Selection Strategies for Classifying High Dimensional Astronomical Data Sets

The amount of collected data in many scientific fields is increasing, all of them requiring a common task: extract knowledge from massive, multi parametric data sets, as rapidly and efficiently possible. This is especially true in astronomy where synoptic sky surveys are enabling new research frontiers in the time domain astronomy and posing several new object classification challenges in multi dimensional spaces; given the high number of parameters available for each object, feature selection is quickly becoming a crucial task in analyzing astronomical data sets. Using data sets extracted from the ongoing Catalina Real-Time Transient Surveys (CRTS) and the Kepler Mission we illustrate a variety of feature selection strategies used to identify the subsets that give the most information and the results achieved applying these techniques to three major astronomical problems.

preprint2013arXiv

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

preprint2013arXiv

Quasars Probing Quasars VI. Excess HI Absorption Within One Proper Mpc of z~2 Quasars

With close pairs of quasars at different redshifts, a background quasar sightline can be used to study a foreground quasar's environment in absorption. We use a sample of 650 projected quasar pairs to study the HI Lya absorption transverse to luminous, z~2 quasars at proper separations of 30kpc < R < 1Mpc. In contrast to measurements along the line-of-sight, regions transverse to quasars exhibit enhanced HI Lya absorption and a larger variance than the ambient intergalactic medium, with increasing absorption and variance toward smaller scales. Analysis of composite spectra reveals excess absorption characterized by a Lya equivalent width profile W = 2.3A (R/100kpc)^-0.46. We also observe a high (~60%) covering factor of strong, optically thick HI absorbers (HI column log NHI > 17.3) at separations R<200kpc, which decreases to ~20% at R~1Mpc, but still represents a significant excess over the cosmic average. This excess of optically thick absorption can be described by a quasar-absorber cross-correlation function xi_QA(r) = (r/r_0)^gamma with a large correlation length r_0 = 12.5+2.7-1.4 Mpc/h (comoving) and gamma = 1.68+0.14-0.30. The HI absorption measured around quasars exceeds that of any previously studied population, consistent with quasars being hosted by massive dark matter halos Mhalo~10^12.5 Msun at z~2.5. The environments of these massive halos are highly biased towards producing optically thick gas, and may even dominate the cosmic abundance of Lyman limit systems and hence the intergalactic opacity to ionizing photons at z~2.5. The anisotropic absorption around quasars implies the transverse direction is much less likely to be illuminated by ionizing radiation than the line-of-sight, which we interpret in terms of the same obscuration effects frequently invoked in unified models of active galactic nuclei.

preprint2013arXiv

Scaling the smoothness of the IGM

We use, for the first time, the Detrend Fluctuation Analysis (DFA) to study the correlation properties of the transmitted flux fluctuations, in the Lyman-$α$ (Ly$α$) Forest along the lines of sight (LOS) to QSOs, at different space scales. We consider in our analysis the transmitted flux in the intergalactic medium over the redshift range 2 $\leq$ z $\leq$ 4.5 from a sample of 45 high-quality medium resolution (R $\sim$ 4300) quasar spectra obtained with Echelle Spectrograph and Imager (ESI) mounted on the Keck II 10-m telescope, and from a sample of 19 high-quality high resolution (R $\sim$ 50000) quasar spectra obtained with Ultra-Violet and Visible Echelle Spectrograh (UVES) mounted on the ESO KUEYEN 8.2 m telescope. The result of the DFA method applied to both datasets, shows that there exists a difference in the correlation properties between the short and long-range regimes: the slopes of the transmitted flux fluctuation function are different on small and large scales. The scaling exponents, $α_{1}$ = 1.635$\pm$ 0.115 and $α_{2}$ = 0.758$\pm$ 0.085 for the ESI/Keck sample and $α_{1}$ = 1.763 $\pm 0.128$ and $α_{2}$ = 0.798 $\pm 0.084$ for the UVES/VLT sample for the short and long range regime respectively. The transition between the two regims is observed at about $\sim 1.4 h^{-1} Mpc$ (comoving). The fact that $α_{1}$ is always larger than $α_{2}$ for each spectrum supports the common view that the Universe is smoother on large scales than on small scales. The non detection of considerable variations in the scaling exponents from LOS to LOS confirms that anisotropies cannot be ubiquitous, at least on these scales.

preprint2013arXiv

Structure and Dynamics of the Globular Cluster Palomar 13

We present Keck/DEIMOS spectroscopy and CFHT/MegaCam photometry for the Milky Way globular cluster Palomar 13. We triple the number of spectroscopically confirmed members, including many repeat velocity measurements. Palomar 13 is the only known globular cluster with possible evidence for dark matter, based on a Keck/HIRES 21 star velocity dispersion of sigma=2.2+/-0.4 km/s. We reproduce this measurement, but demonstrate that it is inflated by unresolved binary stars. For our sample of 61 stars, the velocity dispersion is sigma=0.7(+0.6/-0.5) km/s. Combining our DEIMOS data with literature values, our final velocity dispersion is sigma=0.4(+0.4/-0.3) km/s. We determine a spectroscopic metallicity of [Fe/H]=-1.6+/-0.1 dex, placing a 1-sigma upper limit of sigma_[Fe/H]~0.2 dex on any internal metallicity spread. We determine Palomar 13's total luminosity to be M_V=-2.8+/-0.4, making it among the least luminous known globular clusters. The photometric isophotes are regular out to the half-light radius and mildly irregular outside this radius. The outer surface brightness profile slope is shallower than typical globular clusters (eta=-2.8+/-0.3). Thus at large radius, tidal debris is likely affecting the appearance of Palomar 13. Combining our luminosity with the intrinsic velocity dispersion, we find a dynamical mass of of M_1/2=1.3(+2.7/-1.3)x10^3 M_sun and a mass-to-light ratio of M/L_V=2.4(+5.0/-2.4) M_sun/L_sun. Within our measurement errors, the mass-to-light ratio agrees with the theoretical predictions for a single stellar population. We conclude that, while there is some evidence for tidal stripping at large radius, the dynamical mass of Palomar 13 is consistent with its stellar mass and neither significant dark matter, nor extreme tidal heating, is required to explain the cluster dynamics.

preprint2013arXiv

The MICA Experiment: Astrophysics in Virtual Worlds

We describe the work of the Meta-Institute for Computational Astrophysics (MICA), the first professional scientific organization based in virtual worlds. MICA was an experiment in the use of this technology for science and scholarship, lasting from the early 2008 to June 2012, mainly using the Second Life and OpenSimulator as platforms. We describe its goals and activities, and our future plans. We conducted scientific collaboration meetings, professional seminars, a workshop, classroom instruction, public lectures, informal discussions and gatherings, and experiments in immersive, interactive visualization of high-dimensional scientific data. Perhaps the most successful of these was our program of popular science lectures, illustrating yet again the great potential of immersive VR as an educational and outreach platform. While the members of our research groups and some collaborators found the use of immersive VR as a professional telepresence tool to be very effective, we did not convince a broader astrophysics community to adopt it at this time, despite some efforts; we discuss some possible reasons for this non-uptake. On the whole, we conclude that immersive VR has a great potential as a scientific and educational platform, as the technology matures and becomes more broadly available and accepted.

preprint2013arXiv

Using conditional entropy to identify periodicity

This paper presents a new period finding method based on conditional entropy that is both efficient and accurate. We demonstrate its applicability on simulated and real data. We find that it has comparable performance to other information-based techniques with simulated data but is superior with real data, both for finding periods and just identifying periodic behaviour. In particular, it is robust against common aliasing issues found with other period-finding algorithms.

preprint2012arXiv

A Comprehensive View of a Strongly Lensed Planck-Associated Submillimeter Galaxy

We present high-resolution maps of stars, dust, and molecular gas in a strongly lensed submillimeter galaxy (SMG) at z = 3.259. HATLAS12--00 is selected from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) as a strong lens candidate mainly based on its unusually high 500um flux density (~300 mJy). It is the only high-redshift Planck detection in the 130 deg^2 H-ATLAS Phase 1 area. Keck Adaptive Optics images reveal a quadruply imaged galaxy in the K-band while the Submillimeter Array and the Extended Very Large Array show doubly imaged 880um and CO(1-0) sources, indicating differentiated distributions of the various components in the galaxy. In the source plane, the stars reside in three major kpc-scale clumps extended over ~1.6 kpc, the dust in a compact (~1 kpc) region ~3 kpc north of the stars, and the cold molecular gas in an extended (~7 kpc) disk ~5 kpc northeast of the stars. The emission from the stars, dust, and gas are magnified by ~17, 8, and 7 times, respectively, by four lensing galaxies at z ~ 1. Intrinsically, the lensed galaxy is a warm (T_dust ~ 40-65 K), hyper-luminous (L_IR ~ 1.7e13 Lsun; SFR ~ 2000 Msun/yr), gas-rich (M_gas/M_baryon ~ 70%), young (M_stellar/SFR ~ 20 Myr), and short-lived (M_gas/SFR ~ 40 Myr) starburst, without a significant active galactic nucleus. With physical properties similar to unlensed z > 2 SMGs, HATLAS12--00 offers a detailed view of a typical SMG through a powerful cosmic microscope.

preprint2012arXiv

Classification by Boosting Differences in Input Vectors: An application to datasets from Astronomy

There are many occasions when one does not have complete information in order to classify objects into different classes, and yet it is important to do the best one can since other decisions depend on that. In astronomy, especially time-domain astronomy, this situation is common when a transient is detected and one wishes to determine what it is in order to decide if one must follow it. We propose to use the Difference Boosting Neural Network (DBNN) which can boost differences between feature vectors of different objects in order to differentiate between them. We apply it to the publicly available data of the Catalina Real-Time Transient Survey (CRTS) and present preliminary results. We also describe another use with a stellar spectral library to identify spectra based on a few features. The technique itself is more general and can be applied to a varied class of problems.

preprint2012arXiv

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.

preprint2012arXiv

Data challenges of time domain astronomy

Astronomy has been at the forefront of the development of the techniques and methodologies of data intensive science for over a decade with large sky surveys and distributed efforts such as the Virtual Observatory. However, it faces a new data deluge with the next generation of synoptic sky surveys which are opening up the time domain for discovery and exploration. This brings both new scientific opportunities and fresh challenges, in terms of data rates from robotic telescopes and exponential complexity in linked data, but also for data mining algorithms used in classification and decision making. In this paper, we describe how an informatics-based approach-part of the so-called "fourth paradigm" of scientific discovery-is emerging to deal with these. We review our experiences with the Palomar-Quest and Catalina Real-Time Transient Sky Surveys; in particular, addressing the issue of the heterogeneity of data associated with transient astronomical events (and other sensor networks) and how to manage and analyze it.

preprint2012arXiv

Flashes in a Star Stream: Automated Classification of Astronomical Transient Events

An automated, rapid classification of transient events detected in the modern synoptic sky surveys is essential for their scientific utility and effective follow-up using scarce resources. This presents some unusual challenges: the data are sparse, heterogeneous and incomplete; evolving in time; and most of the relevant information comes not from the data stream itself, but from a variety of archival data and contextual information (spatial, temporal, and multi-wavelength). We are exploring a variety of novel techniques, mostly Bayesian, to respond to these challenges, using the ongoing CRTS sky survey as a testbed. The current surveys are already overwhelming our ability to effectively follow all of the potentially interesting events, and these challenges will grow by orders of magnitude over the next decade as the more ambitious sky surveys get under way. While we focus on an application in a specific domain (astrophysics), these challenges are more broadly relevant for event or anomaly detection and knowledge discovery in massive data streams.

preprint2012arXiv

Probing the Outer Galactic halo with RR Lyrae from the Catalina Surveys

We present the analysis of 12227 type-ab RR Lyrae found among the 200 million public lightcurves in the Catalina Surveys Data Release 1 (CSDR1). These stars span the largest volume of the Milky Way ever surveyed with RR Lyrae, covering ~20,000 square degrees of the sky (0 < RA < 360, -22 < Dec < 65 deg) to heliocentric distances of up to 60kpc. Each of the RR Lyrae are observed between 60 and 419 times over a six-year period. Using period finding and Fourier fitting techniques we determine periods and apparent magnitudes for each source. We find that the periods at generally accurate to sigma = 0.002% by comparison with 2842 previously known RR Lyrae and 100 RR Lyrae observed in overlapping survey fields. We photometrically calibrate the light curves using 445 Landolt standard stars and show that the resulting magnitudes are accurate to ~0.05 mags using SDSS data for ~1000 blue horizontal branch stars and 7788 of the RR Lyrae. By combining Catalina photometry with SDSS spectroscopy, we analyze the radial velocity and metallicity distributions for > 1500 of the RR Lyrae. Using the accurate distances derived for the RR Lyrae, we show the paths of the Sagittarius tidal streams crossing the sky at heliocentric distances from 20 to 60 kpc. By selecting samples of Galactic halo RR Lyrae, we compare their velocity, metallicity, and distance with predictions from a recent detailed N-body model of the Sagittarius system. We find that there are some significant differences between the distances and structures predicted and our observations.

preprint2012arXiv

Sky Surveys

Sky surveys represent a fundamental data basis for astronomy. We use them to map in a systematic way the universe and its constituents, and to discover new types of objects or phenomena. We review the subject, with an emphasis on the wide-field imaging surveys, placing them in a broader scientific and historical context. Surveys are the largest data generators in astronomy, propelled by the advances in information and computation technology, and have transformed the ways in which astronomy is done. We describe the variety and the general properties of surveys, the ways in which they may be quantified and compared, and offer some figures of merit that can be used to compare their scientific discovery potential. Surveys enable a very wide range of science; that is perhaps their key unifying characteristic. As new domains of the observable parameter space open up thanks to the advances in technology, surveys are often the initial step in their exploration. Science can be done with the survey data alone or a combination of different surveys, or with a targeted follow-up of potentially interesting selected sources. Surveys can be used to generate large, statistical samples of objects that can be studied as populations, or as tracers of larger structures. They can be also used to discover or generate samples of rare or unusual objects, and may lead to discoveries of some previously unknown types. We discuss a general framework of parameter spaces that can be used for an assessment and comparison of different surveys, and the strategies for their scientific exploration. As we move into the Petascale regime, an effective processing and scientific exploitation of such large data sets and data streams poses many challenges, some of which may be addressed in the framework of Virtual Observatory and Astroinformatics, with a broader application of data mining and knowledge discovery technologies.

preprint2012arXiv

VAST: An ASKAP Survey for Variables and Slow Transients

The Australian Square Kilometre Array Pathfinder (ASKAP) will give us an unprecedented opportunity to investigate the transient sky at radio wavelengths. In this paper we present VAST, an ASKAP survey for Variables and Slow Transients. VAST will exploit the wide-field survey capabilities of ASKAP to enable the discovery and investigation of variable and transient phenomena from the local to the cosmological, including flare stars, intermittent pulsars, X-ray binaries, magnetars, extreme scattering events, interstellar scintillation, radio supernovae and orphan afterglows of gamma ray bursts. In addition, it will allow us to probe unexplored regions of parameter space where new classes of transient sources may be detected. In this paper we review the known radio transient and variable populations and the current results from blind radio surveys. We outline a comprehensive program based on a multi-tiered survey strategy to characterise the radio transient sky through detection and monitoring of transient and variable sources on the ASKAP imaging timescales of five seconds and greater. We also present an analysis of the expected source populations that we will be able to detect with VAST.

preprint2011arXiv

A Kiloparsec-Scale Binary Active Galactic Nucleus Confirmed by the Expanded Very Large Array

We report the confirmation of a kpc-scale binary active galactic nucleus (AGN) with high-resolution radio images from the Expanded Very Large Array (EVLA). SDSS J150243.1+111557 is a double-peaked [O III] AGN at z = 0.39 selected from the Sloan Digital Sky Survey. Our previous near-infrared adaptive optics imaging reveals two nuclei separated by 1.4" (7.4 kpc), and our optical integral-field spectroscopy suggests that they are a type-1--type-2 AGN pair. However, these data alone cannot rule out the single AGN scenario where the narrow emission-line region associated with the secondary is photoionized by the broad-line AGN in the primary. Our new EVLA images at 1.4, 5.0, and 8.5 GHz show two steep-spectrum compact radio sources spatially coincident with the optical nuclei. The radio power of the type-2 AGN is an order-of-magnitude in excess of star-forming galaxies with similar extinction-corrected [O II] 3727 luminosities, indicating that the radio emission is powered by accretion. Therefore, SDSS J150243.1+111557 is one of the few confirmed kpc-scale binary AGN systems. Spectral-energy-distribution modeling shows that SDSS J150243.1+111557 is a merger of two ~10^{11} M_sun galaxies. With both black hole masses around 10^8 Msun, the AGNs are accreting at ~10 times below the Eddington limit.

preprint2011arXiv

Discovery of a Multiply-Lensed Submillimeter Galaxy in Early HerMES Herschel/SPIRE Data

We report the discovery of a bright ($f(250\mum) > 400$ mJy), multiply-lensed submillimeter galaxy \obj\ in {\it Herschel}/SPIRE Science Demonstration Phase data from the HerMES project. Interferometric 880\mum\ Submillimeter Array observations resolve at least four images with a large separation of $\sim 9\arcsec$. A high-resolution adaptive optics $K_p$ image with Keck/NIRC2 clearly shows strong lensing arcs. Follow-up spectroscopy gives a redshift of $z=2.9575$, and the lensing model gives a total magnification of $μ\sim 11 \pm 1$. The large image separation allows us to study the multi-wavelength spectral energy distribution (SED) of the lensed source unobscured by the central lensing mass. The far-IR/millimeter-wave SED is well described by a modified blackbody fit with an unusually warm dust temperature, $88 \pm 3$ K. We derive a lensing-corrected total IR luminosity of $(1.43 \pm 0.09) \times 10^{13}\, \mathrm{L}_{\odot}$, implying a star formation rate of $\sim 2500\, \mathrm{M}_{\odot}\, \mathrm{yr}^{-1}$. However, models primarily developed from brighter galaxies selected at longer wavelengths are a poor fit to the full optical-to-millimeter SED. A number of other strongly lensed systems have already been discovered in early {\it Herschel} data, and many more are expected as additional data are collected.

preprint2011arXiv

Discovery, classification, and scientific exploration of transient events from the Catalina Real-time Transient Survey

Exploration of the time domain - variable and transient objects and phenomena - is rapidly becoming a vibrant research frontier, touching on essentially every field of astronomy and astrophysics, from the Solar system to cosmology. Time domain astronomy is being enabled by the advent of the new generation of synoptic sky surveys that cover large areas on the sky repeatedly, and generating massive data streams. Their scientific exploration poses many challenges, driven mainly by the need for a real-time discovery, classification, and follow-up of the interesting events. Here we describe the Catalina Real-Time Transient Survey (CRTS), that discovers and publishes transient events at optical wavelengths in real time, thus benefiting the entire community. We describe some of the scientific results to date, and then focus on the challenges of the automated classification and prioritization of transient events. CRTS represents a scientific and a technological testbed and precursor for the larger surveys in the future, including the Large Synoptic Survey Telescope (LSST) and the Square Kilometer Array (SKA).

preprint2011arXiv

Exploring the Time Domain With Synoptic Sky Surveys

Synoptic sky surveys are becoming the largest data generators in astronomy, and they are opening a new research frontier, that touches essentially every field of astronomy. Opening of the time domain to a systematic exploration will strengthen our understanding of a number of interesting known phenomena, and may lead to the discoveries of as yet unknown ones. We describe some lessons learned over the past decade, and offer some ideas that may guide strategic considerations in planning and execution of the future synoptic sky surveys.

preprint2011arXiv

Extracting Knowledge From Massive Astronomical Data Sets

The exponential growth of astronomical data collected by both ground based and space borne instruments has fostered the growth of Astroinformatics: a new discipline laying at the intersection between astronomy, applied computer science, and information and computation (ICT) technologies. At the very heart of Astroinformatics is a complex set of methodologies usually called Data Mining (DM) or Knowledge Discovery in Data Bases (KDD). In the astronomical domain, DM/KDD are still in a very early usage stage, even though new methods and tools are being continuously deployed in order to cope with the Massive Data Sets (MDS) that can only grow in the future. In this paper, we briefly outline some general problems encountered when applying DM/KDD methods to astrophysical problems, and describe the DAME (DAta Mining & Exploration) web application. While specifically tailored to work on MDS, DAME can be effectively applied also to smaller data sets. As an illustration, we describe two application of DAME to two different problems: the identification of candidate globular clusters in external galaxies, and the classification of active galactic nuclei (AGN). We believe that tools and services of this nature will become increasingly necessary for the data-intensive astronomy (and indeed all sciences) in the 21st century.

preprint2011arXiv

Mergers in Double-Peaked [O III] Active Galactic Nuclei

As a natural consequence of galaxy mergers, binary active galactic nuclei (AGNs) should be commonplace. Nevertheless, observational confirmations are rare, especially for binaries with separations less than ten kpc. Such a system may show two sets of narrow emission lines in a single spectrum owing to the orbital motion of the binary. We have obtained high-resolution near-infrared images of 50 double-peaked [O III] 5007 AGNs with the Keck II laser guide star adaptive optics system. The Sloan Digital Sky Survey sample is compiled from the literature and consists of 17 type-1 AGNs between 0.18 < z < 0.56 and 33 type-2 AGNs between 0.03 < z < 0.24. The new images reveal eight type-1 and eight type-2 sources that are apparently undergoing mergers. These are strong candidates of kpc-scale binary AGNs, because they show multiple components separated between 0.6 and 12 kpc and often disturbed morphologies. Because most of the type-1s are at higher redshifts than the type-2s, the higher merger fraction of type-1s (47+/-20%) compared to that of type-2s (24+/-10%) can be attributed to the general evolution of galaxy merger fraction with redshift. Furthermore, we show that AGN mergers are outliers of the M_BH-sigma relation because of over-estimated stellar velocity dispersions, illustrating the importance of removing mergers from the samples defining the M_BH-sigma relations. Finally, we find that the emission-line properties are indistinguishable for spatially resolved and unresolved sources, emphasizing that scenarios involving a single AGN can produce the same double-peaked line profiles and they account for at least 70% of the double-peaked [O III] AGNs.

preprint2011arXiv

Modeling of the HERMES J105751.1+573027 submillimeter source lensed by a dark matter dominated foreground group of galaxies

We present the results of a gravitational lensing analysis of the bright $\zs=2.957$ sub-millimeter galaxy (SMG), HERMES J105751.1+573027 found in {\it Herschel}/SPIRE Science Demonstration Phase data from the Herschel Multi-tiered Extragalactic Survey (HerMES) project. The high resolution imaging available in optical and Near-IR channels, along with CO emission obtained with the Plateau de Bure Interferometer, allow us to precisely estimate the intrinsic source extension and hence estimate the total lensing magnification to be $μ=10.9\pm 0.7$. We measure the half-light radius $R_{\rm eff}$ of the source in the rest-frame Near-UV and $V$ bands that characterize the unobscured light coming from stars and find $R_{\rm eff,*}= [2.0 \pm 0.1]$ kpc, in good agreement with recent studies on the Submillimeter Galaxy population. This lens model is also used to estimate the size of the gas distribution ($R_{\rm eff,gas}= [1.1\pm0.5]$) kpc by mapping back in the source plane the CO (J=5-4) transition line emission. The lens modeling yields a relatively large Einstein radius $R_{\rm Ein}= 4\farcs10 \pm 0\farcs02$, corresponding to a deflector velocity dispersion of [$483\pm 16] \,\kms$. This shows that HERMES J105751.1+573027 is lensed by a {\it galaxy group-size} dark matter halo at redshift $\zl\sim 0.6$. The projected dark matter contribution largely dominates the mass budget within the Einstein radius with $f_{\rm dm}(<R_{\rm Ein})\sim 80%$. This fraction reduces to $f_{\rm dm}(<R_{\rm eff,G1}\simeq 4.5\kpc)\sim 47%$ within the effective radius of the main deflecting galaxy of stellar mass $M_{\rm *,G1}=[8.5\pm 1.6] \times 10^{11}\msun$. At this smaller scale the dark matter fraction is consistent with results already found for massive lensing ellipticals at $z\sim0.2$ from the SLACS survey.

preprint2011arXiv

New Cataclysmic Variable 1RXS J073346.0+261933 in Gemini

In course of the search for the optical identifications associated with ROSAT X-ray sources we have found a highly variable object with the very unusual long-term behavior, color indices and high X-ray-to-optical flux ratio. We report the archival photometric light curve from the Catalina Sky Survey, optical spectroscopy from RTT150 and time-resolved photometry from Astrotel-Caucasus telescope. The object appears to be the magnetic cataclysmic variable (polar) with orbital period of P=3.20 hr.

preprint2011arXiv

Real Time Classification of Transient Events in Synoptic Sky Surveys

An automated, rapid classification of transient events detected in the modern synoptic sky surveys is essential for their scientific utility and effective follow-up using scarce resources. This problem will grow by orders of magnitude with the next generation of surveys. We are exploring a variety of novel automated classification techniques, mostly Bayesian, to respond to these challenges, using the ongoing CRTS sky survey as a testbed. We describe briefly some of the methods used.

preprint2011arXiv

Recommendations of the Virtual Astronomical Observatory (VAO) Science Council for the VAO second year activity

The VAO (Virtual Astronomical Observatory) Science Council (VAO-SC) met on July 27-28, 2011 at the Harvard-Smithsonian Center for Astrophysics in Cambridge MA, to review the VAO performance during its first year of operations. In this meeting the VAO demonstrated the new tools for astronomers that are being released in September 2011 and presented plans for the second year of activities, resulting from studies conducted during the first year. This document contains the recommendations of the VAO-SC for the second year of activity of the VAO.

preprint2011arXiv

SN 2008jb: A "Lost" Core-Collapse Supernova in a Star-Forming Dwarf Galaxy at ~10 Mpc

We present the discovery and follow-up observations of SN 2008jb, a core-collapse supernova in the dwarf irregular galaxy ESO 302-14 at 9.6 Mpc. This transient was missed by galaxy-targeted surveys and was only found in archival optical images obtained by CRTS and ASAS. It was detected shortly after explosion and reached a bright optical maximum, Vmax = 13.6 mag (M_Vmax ~ -16.5). The shape of the light curve shows a plateau of 100 days, followed by a drop of 1.4 mag in V-band to a decline with the approximate Co 56 decay slope, consistent with 0.04 Msun of Ni 56 synthesized in the explosion. A spectrum obtained 2 years after explosion shows a broad, boxy Halpha emission line, which is unusual for type IIP supernovae. We detect the supernova in archival Spitzer and WISE images obtained 8-14 months after explosion, which show clear signs of warm dust emission. The dwarf irregular host galaxy has a low gas-phase oxygen abundance, 12 + log(O/H) = 8.2 (~1/5 Solar), similar to those of the SMC and the hosts of long gamma-ray bursts and luminous core-collapse supernovae. We study the host environment using GALEX FUV, R-band, and Halpha images and find that the supernova occurred in a large star-formation complex. The morphology of the Halpha emission appears as a large shell (R = 350 pc) surrounding the FUV and optical emission. We estimate an age of ~9 Myr and a total mass of ~2 x 10^5 Msun for the star-formation complex. These properties are consistent with the expanding Halpha supershells observed in well-studied nearby dwarf galaxies, which are tell-tale signs of feedback from the cumulative effect of massive star winds and supernovae. The age estimated for the star-forming region suggests a relatively high-mass progenitor star with initial mass of ~20 Msun. We discuss the implications of these findings in the study of core-collapse supernova progenitors. (Abridged)

preprint2011arXiv

The Catalina Real-time Transient Survey

The Catalina Real-time Transient Survey (CRTS) currently covers 33,000 deg^2 of the sky in search of transient astrophysical events, with time baselines ranging from 10 minutes to ~7 years. Data provided by the Catalina Sky Survey provides an unequaled baseline against which >4,000 unique optical transient events have been discovered and openly published in real-time. Here we highlight some of the discoveries of CRTS.

preprint2011arXiv

The Catalina Real-Time Transient Survey (CRTS)

Catalina Real-Time Transient Survey (CRTS) is a synoptic sky survey uses data streams from 3 wide-field telescopes in Arizona and Australia, covering the total area of ~30,000 deg2, down to the limiting magnitudes ~ 20 - 21 mag per exposure, with time baselines from 10 min to 6 years (and growing); there are now typically ~ 200 - 300 exposures per pointing, and coadded images reach deeper than 23 mag. The basic goal of CRTS is a systematic exploration and characterization of the faint, variable sky. The survey has detected ~ 3,000 high-amplitude transients to date, including ~ 1,000 supernovae, hundreds of CVs (the majority of them previously uncatalogued), and hundreds of blazars / OVV AGN, highly variable and flare stars, etc. CRTS has a complete open data philosophy: all transients are published immediately electronically, with no proprietary period at all, and all of the data (images, light curves) will be publicly available in the near future, thus benefiting the entire astronomical community. CRTS is a scientific and technological testbed and precursor for the grander synoptic sky surveys to come.

preprint2011arXiv

The DAME/VO-Neural Infrastructure: an Integrated Data Mining System Support for the Science Community

Astronomical data are gathered through a very large number of heterogeneous techniques and stored in very diversified and often incompatible data repositories. Moreover in the e-science environment, it is needed to integrate services across distributed, heterogeneous, dynamic "virtual organizations" formed by different resources within a single enterprise and/or external resource sharing and service provider relationships. The DAME/VONeural project, run jointly by the University Federico II, INAF (National Institute of Astrophysics) Astronomical Observatories of Napoli and the California Institute of Technology, aims at creating a single, sustainable, distributed e-infrastructure for data mining and exploration in massive data sets, to be offered to the astronomical (but not only) community as a web application. The framework makes use of distributed computing environments (e.g. S.Co.P.E.) and matches the international IVOA standards and requirements. The integration process is technically challenging due to the need of achieving a specific quality of service when running on top of different native platforms. In these terms, the result of the DAME/VO-Neural project effort will be a service-oriented architecture, obtained by using appropriate standards and incorporating Grid paradigms and restful Web services frameworks where needed, that will have as main target the integration of interdisciplinary distributed systems within and across organizational domains.

preprint2011arXiv

The Discovery and Nature of Optical Transient CSS100217:102913+404220

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.

preprint2011arXiv

The Faint End of the Quasar Luminosity Function at z~4: Implications for Ionization of the Intergalactic Medium and Cosmic Downsizing

We present an updated determination of the z ~ 4 QSO luminosity function (QLF), improving the quality of the determination of the faint end of the QLF presented in Glikman et al. (2010). We have observed an additional 43 candidates from our survey sample, yielding one additional QSO at z = 4.23 and increasing the completeness of our spectroscopic follow-up to 48% for candidates brighter than R = 24 over our survey area of 3.76 deg2. We study the effect of using K-corrections to compute the rest-frame absolute magnitude at 1450A compared with measuring M1450 directly from the object spectra. We find a luminosity-dependent bias: template-based K-corrections overestimate the luminosity of low-luminosity QSOs, likely due to their reliance on templates derived from higher luminosity QSOs. Combining our sample with bright quasars from the Sloan Digital Sky Survey and using spectrum-based M1450 for all the quasars, we fit a double-power-law to the binned QLF. Our best fit has a bright-end slope, α = 3.3\pm0.2, and faint-end slope, β = 1.6(+0.8/-0.6). Our new data revise the faint-end slope of the QLF down to flatter values similar to those measured at z ~ 3. The break luminosity, though poorly constrained, is at M* = -24.1(+0.7/-1.9), approximately 1 - 1.5 mag fainter than at z ~ 3. This QLF implies that QSOs account for about half the radiation needed to ionize the IGM at these redshifts.

preprint2011arXiv

The Nature of Double-Peaked [O III] Active Galactic Nuclei

Active galactic nuclei (AGNs) with double-peaked [O III] lines are suspected to be sub-kpc or kpc-scale binary AGNs. However, pure gas kinematics can produce the same double-peaked line profile in spatially integrated spectra. Here we combine integral-field spectroscopy and high-resolution imaging of 42 double-peaked [O III] AGNs from Sloan Digital Sky Survey to investigate the constituents of the population. We find two binary AGNs where the line-splitting is driven by the orbital motion of the merging nuclei. Such objects account for only ~2% of the double-peaked AGNs. Almost all (~98%) of the double-peaked AGNs were selected because of gas kinematics; and half of those show spatially resolved narrow-line regions that extend 4-20 kpc from the nuclei. Serendipitously, we find two spectrally unresolved binary AGNs where gas kinematics produced the double-peaked [O III] lines. The relatively frequent serendipitous discoveries indicate that only ~1% of binary AGNs would appear double-peaked in Sloan spectra and 2.2_{-0.8}^{+2.5}% of all Sloan AGNs are binary AGNs. Therefore, the double-peaked sample does not offer much advantage over any other AGN samples in finding binary AGNs. The binary AGN fraction implies an elevated AGN duty cycle (8_{-3}^{+8}%), suggesting galaxy interactions enhance nuclear accretion. We illustrate that integral-field spectroscopy is crucial for identifying binary AGNs: several objects previously classified as "binary AGNs" with long-slit spectra are most likely single AGNs with extended narrow-line regions. The formation of extended narrow-line regions driven by radiation pressure is also discussed.

preprint2011arXiv

The VAO Transient Facility

The time domain community wants robust and reliable tools to enable production of and subscription to community-endorsed event notification packets (VOEvent). The VAO Transient Facility (VTF) is being designed to be the premier brokering service for the community, both collecting and disseminating observations about time-critical astronomical transients but also supporting annotations and the application of intelligent machine-learning to those observations. This distinguishes two types of activity associated with the facility: core infrastructure and user services. In this paper, we will review the prior art in both areas and describe the planned capabilities of the VTF. In particular, we will focus on scalability and quality-of-service issues required by the next generation of sky surveys, such as LSST and SKA.

preprint2011arXiv

Three QSOs acting as strong gravitational lenses

We report the discovery of three new cases of QSOs acting as strong gravitational lenses on background emission line galaxies: SDSS J0827+5224 (zQSO = 0.293, zs = 0.412), SDSS J0919+2720 (zQSO = 0.209, zs = 0.558), SDSS J1005+4016 (zQSO = 0.230, zs = 0.441). The selection was carried out using a sample of 22,298 SDSS spectra displaying at least four emission lines at a redshift beyond that of the foreground QSO. The lensing nature is confirmed from Keck imaging and spectroscopy, as well as from HST/WFC3 imaging in the F475W and F814W filters. Two of the QSOs have face-on spiral host galaxies and the third is a QSO+galaxy pair. The velocity dispersion of the host galaxies, inferred from simple lens modeling, is between σ_v = 210 and 285 km/s, making these host galaxies comparable in mass with the SLACS sample of early-type strong lenses.

preprint2011arXiv

Towards an Automated Classification of Transient Events in Synoptic Sky Surveys

We describe the development of a system for an automated, iterative, real-time classification of transient events discovered in synoptic sky surveys. The system under development incorporates a number of Machine Learning techniques, mostly using Bayesian approaches, due to the sparse nature, heterogeneity, and variable incompleteness of the available data. The classifications are improved iteratively as the new measurements are obtained. One novel feature is the development of an automated follow-up recommendation engine, that suggest those measurements that would be the most advantageous in terms of resolving classification ambiguities and/or characterization of the astrophysically most interesting objects, given a set of available follow-up assets and their cost functions. This illustrates the symbiotic relationship of astronomy and applied computer science through the emerging discipline of AstroInformatics.

preprint2010arXiv

Discovery of eclipsing white dwarf systems in a search for Earth-size companions

Although white dwarfs are believed to be the end point of most stellar evolution, unlike main sequence stars, they have not yet been the subject of dedicated time-domain surveys for exoplanets. We discuss how their size and distinctive colour make them excellent targets for wide-field searches for exoplanets. In particular, we note that planets of Earth-size can give rise to multi-magnitude eclipses of massive white dwarfs. Such a large signal is almost unmistakable and would be detectable even with very low-precision photometry. For objects of smaller size, the high accuracy photometry currently being used to detect Super-Earth and smaller planets transiting Sun-sized stars, is capable of revealing minor planets down to R~100km as they transit white dwarfs. Such observations can be used to test current evidence for asteroid-size objects being the cause for dust rings which have recently been observed for a number of white dwarfs. No other current exoplanet search method is capable of detecting such exo-asteroids. As an initial test of this search strategy, we combine synoptic data from the Catalina Sky Survey with multi-colour photometry and spectra from the Sloan Digital Sky Survey to search ~12,000 white dwarf lightcurves for eclipsing events. We find 20 new eclipsing white dwarf binary systems with low-mass companions. This doubles the number of known eclipsing white dwarfs and is expected to enable the determination of accurate white dwarf radii. Three of the discoveries have radii consistent with substellar systems and show no evidence of flux from the eclipsing object in their SDSS optical spectra, or near-IR data.

preprint2010arXiv

First case of strong gravitational lensing by a QSO : SDSS J0013+1523 at z = 0.120

We present the first case of strong gravitational lensing by a QSO : SDSS J0013+1523, at z = 0.120. The discovery is the result of a systematic search for emission lines redshifted behind QSOs, among 22298 spectra of the SDSS data release 7. Apart from the z = 0.120 spectral features of the foreground QSO, the spectrum of SDSS J0013+1523 also displays the OII and Hbeta emission lines and the OIII doublet, all at the same redshift, z = 0.640. Using sharp Keck adaptive optics K-band images obtained using laser guide stars, we unveil two objects within a radius of 2 arcsec from the QSO. Deep Keck optical spectroscopy clearly confirms one of these objects at z = 0.640 and shows traces of the OIII, emission line of the second object, also at z = 0.640. Lens modeling suggests that they represent two images of the same z = 0.640 emission-line galaxy. Our Keck spectra also allow us to measure the redshift of an intervening galaxy at z = 0.394, located 3.2 arcsec away from the line of sight to the QSO. If the z = 0.120 QSO host galaxy is modeled as a singular isothermal sphere, its mass within the Einstein radius is M_E(r < 1 kpc) = 2.16e10 M_Sun and its velocity dispersion is sigma_SIS = 169 km/s. This is about 1 sigma away from the velocity dispersion estimated from the width of the QSO Hbeta emission line, sigma_*(M_BH) = 124 +/- 47 km/s. Deep optical HST imaging will be necessary to constrain the total radial mass profile of the QSO host galaxy using the detailed shape of the lensed source. This first case of a QSO acting as a strong lens on a more distant object opens new directions in the study of QSO host galaxies.

preprint2010arXiv

First Earth-based Detection of a Superbolide on Jupiter

Cosmic collisions on planets cause detectable optical flashes that range from terrestrial shooting stars to bright fireballs. On June 3, 2010 a bolide in Jupiter's atmosphere was simultaneously observed from the Earth by two amateur astronomers observing Jupiter in red and blue wavelengths. The bolide appeared as a flash of 2 s duration in video recording data of the planet. The analysis of the light curve of the observations results in an estimated energy of the impact of 0.9-4.0x10^{15} J which corresponds to a colliding body of 8-13 m diameter assuming a mean density of 2 g cm^{-3}. Images acquired a few days later by the Hubble Space Telescope and other large ground-based facilities did not show any signature of aerosol debris, temperature or chemical composition anomaly, confirming that the body was small and destroyed in Jupiter's upper atmosphere. Several collisions of this size may happen on Jupiter on a yearly basis. A systematic study of the impact rate and size of these bolides can enable an empirical determination of the flux of meteoroids in Jupiter with implications for the populations of small bodies in the outer Solar System and may allow a better quantification of the threat of impacting bodies to Earth. The serendipitous recording of this optical flash opens a new window in the observation of Jupiter with small telescopes.

preprint2010arXiv

Transverse Sizes of CIV Absorption Systems Measured from Multiple QSO Sightlines

We present tomography of the circum-galactic metal distribution at redshift 1.7 to 4.5 derived from echellete spectroscopy of binary quasars. We find CIV systems at similar redshifts in paired sightlines more often than expected for sightline-independent redshifts. As the separation of the sightlines increases from 36 kpc to 907 kpc, the amplitude of this clustering decreases. At the largest separations, the CIV systems cluster similar to Lyman-break galaxies (Adelberger et al. 2005a). The CIV systems are significantly less correlated than these galaxies, however, at separations less than R_1 ~ 0.42 +/- 0.15 h-1 comoving Mpc. Measured in real space, i.e., transverse to the sightlines, this length scale is significantly smaller than the break scale estimated from the line-of-sight correlation function in redshift space (Scannapieco et al. 2006a). Using a simple model, we interpret the new real-space measurement as an indication of the typical physical size of enriched regions. We adopt this size for enriched regions and fit the redshift-space distortion in the line-of-sight correlation function. The fitted velocity kick is consistent with the peculiar velocity of galaxies as determined by the underlying mass distribution and places an upper limit on the outflow (or inflow) speed of metals. The implied time scale for dispersing metals is larger than the typical stellar ages of Lyman-break galaxies (Shapley et al. 2001), and we argue that enrichment by galaxies at z > 4.3 played a greater role in dispersing metals. To further constrain the growth of enriched regions, we discuss empirical constraints on the evolution of the CIV correlation function with cosmic time. This study demonstrates the potential of tomography for measuring the metal enrichment history of the circum-galactic medium.

preprint2009arXiv

Binary Quasars at High Redshift I: 24 New Quasar Pairs at z ~ 3-4

The clustering of quasars on small scales yields fundamental constraints on models of quasar evolution and the buildup of supermassive black holes. This paper describes the first systematic survey to discover high redshift binary quasars. Using color-selection and photometric redshift techniques, we searched 8142 deg^2 of SDSS imaging data for binary quasar candidates, and confirmed them with follow-up spectroscopy. Our sample of 27 high redshift binaries (24 of them new discoveries) at redshifts 2.9 < z < 4.3 with proper transverse separations 10 kpc < R_{\perp} < 650 kpc increases the number of such objects known by an order of magnitude. Eight members of this sample are very close pairs with R_{\perp} < 100 kpc, and of these close systems four are at z > 3.5. The completeness and efficiency of our well-defined selection algorithm are quantified using simulated photometry and we find that our sample is ~ 50% complete. Our companion paper uses this knowledge to make the first measurement of the small scale clustering (R < 1 Mpc/h comoving) of high-redshift quasars. High redshift binaries constitute exponentially rare coincidences of two extreme (M >~ 10^9 Msun) supermassive black holes. At z ~ 4 there is about one close binary per 10 Gpc^3, thus these could be the highest sigma peaks, the analogs of superclusters, in the early Universe.

preprint2009arXiv

Binary Quasars at High Redshift II: Sub-Mpc Clustering at z ~ 3-4

We present measurements of the small-scale (0.1<~ r <~ 1 Mpc/h) quasar two-point correlation function at z>2.9, for a flux-limited (i<21) sample of 15 binary quasars compiled by Hennawi et al. (2009). The amplitude of the small-scale clustering increases from z ~ 3 to z ~ 4. The small-scale clustering amplitude is comparable to or lower than power-law extrapolations (with slope gamma=2) from the large-scale correlation function of the i<20.2 quasar sample from the Sloan Digital Sky Survey. Using simple prescriptions relating quasars to dark matter halos, we model the observed small-scale clustering with halo occupation models. Reproducing the large-scale clustering amplitude requires that the active fraction of the black holes in the central galaxies of halos is near unity, but the level of small-scale clustering favors an active fraction of black holes in satellite galaxies 0.1 <~ f_s <~ 0.5 at z >~ 3.

preprint2009arXiv

Discovery of the Extremely Energetic Supernova 2008fz

We report on the discovery and initial observations of the energetic type IIn supernova (SN), 2008fz. The optical energy emitted by SN 2008fz (based on the light curve over a 88 day period), is possibly the most ever observed for a supernova (1.4 x 10^51 erg). The event was more luminous than the type IIn SN 2006gy, but exhibited same smooth, slowly evolving light curve. As is characteristic of type IIn SN, the early spectra of 2008fz initially exhibited narrow Balmer lines which were replaced by a broader component at later times. The spectra also show a blue continuum with no signs of Ca or Na absorption, suggesting that there is little extinction due to intragalatic dust in the host or circumstellar material. No host galaxy is identified in prior coadded images reaching R ~ 22. From the supernova's redshift, z=0.133, we place an upper limit on the host of M_R=-17. The presence of the SN within such a faint host follows the majority of recently discovered highly luminous SN. A possible reason for this occurrence is the very high star formation rate occurring in low-mass galaxies in combination with the low metallicity environment, which makes the production of very massive stars possible. We determine the peak absolute magnitude of the event to be M_V = -22.3 from the initial photometry and the redshift distance, placing it among the most luminous supernovae discovered.

preprint2009arXiv

Exploring the Use of Virtual Worlds as a Scientific Research Platform: The Meta-Institute for Computational Astrophysics (MICA)

We describe the Meta-Institute for Computational Astrophysics (MICA), the first professional scientific organization based exclusively in virtual worlds (VWs). The goals of MICA are to explore the utility of the emerging VR and VWs technologies for scientific and scholarly work in general, and to facilitate and accelerate their adoption by the scientific research community. MICA itself is an experiment in academic and scientific practices enabled by the immersive VR technologies. We describe the current and planned activities and research directions of MICA, and offer some thoughts as to what the future developments in this arena may be.

preprint2009arXiv

Highly Variable Objects in the Palomar-QUEST Survey: A Blazar Search using Optical Variability

We identify 3,113 highly variable objects in 7,200 square degrees of the Palomar-QUEST Survey, which each varied by more than 0.4 magnitudes simultaneously in two broadband optical filters on timescales from hours to roughly 3.5 years. The primary goal of the selection is to find blazars by their well-known violent optical variability. Because most known blazars have been found in radio and/or X-ray wavelengths, a sample discovered through optical variability may have very different selection effects, elucidating the range of behavior possible in these systems. A set of blazars selected in this unusual manner will improve our understanding of the physics behind this extremely variable and diverse class of AGN. The object positions, variability statistics, and color information are available using the Palomar-QUEST CasJobs server. The time domain is just beginning to be explored over large sky areas; we do not know exactly what a violently variable sample will hold. About 20% of the sample has been classified in the literature; over 70% of those objects are known or likely AGN. The remainder largely consists of a variety of variable stars, including a number of RR Lyrae and cataclysmic variables.

preprint2009arXiv

LSST Science Book, Version 2.0

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

preprint2009arXiv

The Faint End of the Quasar Luminosity Function at z~4

We have conducted a spectroscopic survey to find faint quasars (-26.0 < M_{1450} < -22.0) at redshifts z=3.8-5.2 in order to measure the faint end of the quasar luminosity function at these early times. Using available optical imaging data from portions of the NOAO Deep Wide-Field Survey and the Deep Lens Survey, we have color-selected quasar candidates in a total area of 3.76 deg^2. Thirty candidates have R <= 23 mags. We conducted spectroscopic followup for 28 of our candidates and found 23 QSOs, 21 of which are reported here for the first time, in the 3.74 < z <5.06 redshift range. We estimate our survey completeness through detailed Monte Carlo simulations and derive the first measurement of the density of quasars in this magnitude and redshift interval. We find that the binned luminosity function is somewhat affected by the K-correction used to compute the rest-frame absolute magnitude at 1450A. Considering only our R <= 23 sample, the best-fit single power-law (Phi \propto L^beta) gives a faint-end slope beta = -1.6+/-0.2. If we consider our larger, but highly incomplete sample going one magnitude fainter, we measure a steeper faint-end slope -2 < beta < -2.5. In all cases, we consistently find faint-end slopes that are steeper than expected based on measurements at z ~ 3. We combine our sample with bright quasars from the Sloan Digital Sky Survey to derive parameters for a double-power-law luminosity function. Our best fit finds a bright-end slope, alpha = -2.4+/-0.2, and faint-end slope, beta = -2.3+/-0.2, without a well-constrained break luminosity. This is effectively a single power-law, with beta = -2.7+/-0.1. We use these results to place limits on the amount of ultraviolet radiation produced by quasars and find that quasars are able to ionize the intergalactic medium at these redshifts.

preprint2009arXiv

The Palomar Transient Factory: System Overview, Performance and First Results

The Palomar Transient Factory (PTF) is a fully-automated, wide-field survey aimed at a systematic exploration of the optical transient sky. The transient survey is performed using a new 8.1 square degree camera installed on the 48-inch Samuel Oschin telescope at Palomar Observatory; colors and light curves for detected transients are obtained with the automated Palomar 60-inch telescope. PTF uses eighty percent of the 1.2-m and fifty percent of the 1.5-m telescope time. With an exposure of 60-s the survey reaches a depth of approximately 21.3 in g' and 20.6 in R (5 sigma, median seeing). Four major experiments are planned for the five-year project: 1) a 5-day cadence supernova search; 2) a rapid transient search with cadences between 90 seconds and 1 day; 3) a search for eclipsing binaries and transiting planets in Orion; and 4) a 3-pi sr deep H-alpha survey. PTF provides automatic, realtime transient classification and follow up, as well as a database including every source detected in each frame. This paper summarizes the PTF project, including several months of on-sky performance tests of the new survey camera, the observing plans and the data reduction strategy. We conclude by detailing the first 51 PTF optical transient detections, found in commissioning data.

preprint2009arXiv

The velocity dispersion and mass-to-light ratio of the remote halo globular cluster NGC 2419

Precise radial velocity measurements from HIRES on the Keck I telescope are presented for 40 stars in the outer halo globular cluster NGC 2419. These data are used to probe the cluster's stellar mass function and search for the presence of dark matter in this cluster. NGC 2419 is one of the best Galactic globular clusters for such a study due to its long relaxation time (T_{r0} ~ 10^{10} yr) and large Galactocentric distance (R_{GC} ~ 90 kpc) -- properties that make significant evolutionary changes in the low-mass end of the cluster mass function unlikely. We find a mean cluster velocity of <v_r>=-20.3 +- 0.7 km/sec and an internal velocity dispersion of σ= 4.14 +- 0.48 km/sec, leading to a total mass of (9.0 +- 2.2) * 10^5 Msun and a global mass-to-light ratio of M/L_V = 2.05 +- 0.50 in solar units. This mass-to-light ratio is in good agreement with what one would expect for a pure stellar system following a standard mass function at the metallicity of NGC 2419. In addition, the mass-to-light ratio does not appear to rise towards the outer parts of the cluster. Our measurements therefore rule out the presence of a dark matter halo with mass larger than ~10^7 Msun inside the central 500 pc, which is lower than what is found for the central dark matter densities of dSph galaxies. We also discuss the relevance of our measurements for alternative gravitational theories such as MOND, and for possible formation scenarios of ultra-compact dwarf galaxies.

preprint2009arXiv

Thirty Meter Telescope Site Testing I: Overview

As part of the conceptual and preliminary design processes of the Thirty Meter Telescope (TMT), the TMT site testing team has spent the last five years measuring the atmospheric properties of five candidate mountains in North and South America with an unprecedented array of instrumentation. The site testing period was preceded by several years of analyses selecting the five candidates, Cerros Tolar, Armazones and Tolonchar in northern Chile; San Pedro Martir in Baja California, Mexico and the 13 North (13N) site on Mauna Kea, Hawaii. Site testing was concluded by the selection of two remaining sites for further consideration, Armazones and Mauna Kea 13N. It showed that all five candidates are excellent sites for an extremely large astronomical observatory and that none of the sites stands out as the obvious and only logical choice based on its combined properties. This is the first article in a series discussing the TMT site testing project.

preprint2008arXiv

First Results from the Catalina Real-time Transient Survey

We report on the results from the first six months of the Catalina Real-time Transient Survey (CRTS). In order to search for optical transients with timescales of minutes to years, the CRTS analyses data from the Catalina Sky Survey which repeatedly covers twenty six thousand of square degrees on the sky. The CRTS provides a public stream of transients that are bright enough to be followed up using small telescopes. Since the beginning of the survey, all CRTS transients have been made available to astronomers around the world in real-time using HTML tables, RSS feeds and VOEvents. As part of our public outreach program the detections are now also available in KML through Google Sky. The initial discoveries include over 350 unique optical transients rising more than two magnitudes from past measurements. Sixty two of these are classified as supernovae, based on light curves, prior deep imaging and spectroscopic data. Seventy seven are due to cataclysmic variables (only 13 previously known), while an additional 100 transients were too infrequently sampled to distinguish between faint CVs and SNe. The remaining optical transients include AGN, Blazars, high proper motions stars, highly variable stars (such as UV Ceti stars) and transients of an unknown nature. Our results suggest that there is a large population of SNe missed by many current supernova surveys because of selection biases. These objects appear to be associated with faint host galaxies. We also discuss the unexpected discovery of white dwarf binary systems through dramatic eclipses.

preprint2008arXiv

The Environments of High Redshift QSOs

We present a sample of $i_{775}$-dropout candidates identified in five Hubble Advanced Camera for Surveys fields centered on Sloan Digital Sky Survey QSOs at redshift $z\sim 6$. Our fields are as deep as the Great Observatory Origins Deep Survey (GOODS) ACS images which are used as a reference field sample. We find them to be overdense in two fields, underdense in two fields, and as dense as the average density of GOODS in one field. The two excess fields show significantly different color distributions from that of GOODS at the 99% confidence level, strengthening the idea that the excess objects are indeed associated with the QSO. The distribution of $i_{775}$-dropout counts in the five fields is broader than that derived from GOODS at the 80% to 96% confidence level, depending on which selection criteria were adopted to identify $i_{775}$-dropouts; its width cannot be explained by cosmic variance alone. Thus, QSOs seem to affect their environments in complex ways. We suggest the picture where the highest redshift QSOs are located in very massive overdensities and are therefore surrounded by an overdensity of lower mass halos. Radiative feedback by the QSO can in some cases prevent halos from becoming galaxies, thereby generating in extreme cases an underdensity of galaxies. The presence of both enhancement and suppression is compatible with the expected differences between lines of sight at the end of reionization as the presence of residual diffuse neutral hydrogen would provide young galaxies with shielding from the radiative effects of the QSO.

preprint2008arXiv

The Origins and the Early Evolution of Quasars and Supermassive Black Holes

The relationship between galaxies and supermassive black holes (SMBH) found in their cores plays a key role in the formation and evolution of both of these major constituents of the universe, as well as the evolution of the intergalactic medium. Neither can be fully understood on their own, and studies of galaxy and SMBH co-formation and co-evolution are now among the central topics of research in cosmology. Yet the very origins, and the early growth phases of the SMBH are still not firmly established. We review our current understanding of the relevant processes and their astrophysical and cosmological context, with an emphasis on the observability of the SMBH growth mechanisms at high redshifts, and their leftover progeny at low redshifts.

preprint2006arXiv

HST and Spitzer Observations of the Host Galaxy of GRB 050904: A Metal-Enriched, Dusty Starburst at z=6.295

We present deep Hubble Space Telescope and Spitzer Space Telescope observations of the host galaxy of GRB 050904 at z=6.295. The host is detected in the H-band and marginally at 3.6 micron. From these detections, and limits in the z'-band and 4.5 micron, we infer an extinction-corrected absolute magnitude, M(UV)=-20.7 mag, or ~L*, a substantial star formation rate of 15 solar masses per year, and a stellar mass of a few 10^9 solar masses. A comparison to the published sample of spectroscopically-confirmed galaxies at z>5.5 reveals that the host of GRB 050904 would evade detection and/or confirmation in any of the current surveys due to the lack of detectable Ly-alpha emission, which is likely the result of dust extinction (A[1200]~1.5 mag). This suggests that not all luminous starburst galaxies at z~6 are currently being accounted for. Most importantly, using the metallicity of Z~0.05 solar inferred from the afterglow absorption spectrum, our observations indicate for the first time that the observed evolution in the mass- and luminosity-metallicity relations from z=0 to z~2 continues on to z>6. The ease of measuring redshifts and metallicities from the afterglow emission suggests that in tandem with the next generation ground- and space-based telescopes, a GRB mission with dedicated near-IR follow-up can provide unique information on the evolution of stars and galaxies through the epoch of re-ionization.

preprint2006arXiv

Some Pattern Recognition Challenges in Data-Intensive Astronomy

We review some of the recent developments and challenges posed by the data analysis in modern digital sky surveys, which are representative of the information-rich astronomy in the context of Virtual Observatory. Illustrative examples include the problems of an automated star-galaxy classification in complex and heterogeneous panoramic imaging data sets, and an automated, iterative, dynamical classification of transient events detected in synoptic sky surveys. These problems offer good opportunities for productive collaborations between astronomers and applied computer scientists and statisticians, and are representative of the kind of challenges now present in all data-intensive fields. We discuss briefly some emergent types of scalable scientific data analysis systems with a broad applicability.

preprint2006arXiv

The ERO Host Galaxy of GRB 020127: Implications for the Metallicity of GRB Progenitors

We present optical and near-IR observations of the host galaxy of GRB 020127, for which we measure R-K=6.2 mag. This is only the second GRB host to date, which is classified as an ERO. The spectral energy distribution is typical of a dusty starburst galaxy, with a redshift, z~1.9, a luminosity, L~5L*, and an inferred stellar mass of ~10^11-10^12 Msun, two orders of magnitude more massive than typical GRB hosts. A comparison to the z~2 mass-metallicity (M-Z) relation indicates that the host metallicity is about 0.5 Zsun. This result shows that at least some GRBs occur in massive, metal-enriched galaxies, and that the proposed low metallicity bias of GRB progenitors is not as severe as previously claimed. Instead we propose that the blue colors and sub-L* luminosities of most GRB hosts reflect their young starburst populations. This, along with the locally increased fraction of starbursts at lower stellar mass, may in fact be the underlying reason for the claimed low metallicity bias of z<0.2 GRB hosts. Consequently, GRBs and their host galaxies may serve as a reliable tracer of cosmic star formation, particularly at z>1 where the M-Z relation is systematically lower, and the star formation rate is dominated by sub-L* galaxies.

preprint2005arXiv

The ACS Virgo Cluster Survey. VII. Resolving the Connection Between Globular Clusters and Ultra-Compact Dwarf Galaxies

We investigate the connection between globular clusters and ultra-compact dwarf galaxies (UCDs) by examining the properties of several compact objects associated with M87, all of which were previously classified as globular clusters. Combining imaging from the Hubble Space Telescope with ground-based Keck spectroscopy, we find two objects to have half-light radii, velocity dispersions and mass-to-light ratios that are consistent with the predictions of population synthesis models for old, metal-rich, luminous globular clusters. Three other objects are much larger, with half-light radii of approximately 20pc, and have V-band mass-to-light ratios in the range 6-9. These objects, which we consider to be UCDs, resemble the nuclei of nucleated dwarf elliptical galaxies in Virgo, having similar mass-to-light ratios, luminosities and colors. These UCDs are found to obey the extrapolated scaling relations of galaxies more closely than those of Galactic globular clusters. There appears to be a transition between the two types of stellar systems at a mass of about two million solar masses. If the UCDs are gravitationally bound, then we suggest that the presence of dark matter is the fundamental property distinguishing globular clusters from UCDs. More than half of the UCD candidates uncovered in the ACS Virgo Cluster Survey are associated with a single galaxy -- M87 -- which suggests that proximity to the Virgo center may be of critical importance for the formation of these objects. These results show that distinguishing bonafide UCDs from bright globular clusters requires a careful analysis of their detailed structural and dynamical properties, particularly their mass-to-light ratios. (ABRIDGED)

preprint2005arXiv

Virtual Astronomy, Information Technology, and the New Scientific Methodology

All sciences, including astronomy, are now entering the era of information abundance. The exponentially increasing volume and complexity of modern data sets promises to transform the scientific practice, but also poses a number of common technological challenges. The Virtual Observatory concept is the astronomical community's response to these challenges: it aims to harness the progress in information technology in the service of astronomy, and at the same time provide a valuable testbed for information technology and applied computer science. Challenges broadly fall into two categories: data handling (or "data farming"), including issues such as archives, intelligent storage, databases, interoperability, fast networks, etc., and data mining, data understanding, and knowledge discovery, which include issues such as automated clustering and classification, multivariate correlation searches, pattern recognition, visualization in highly hyperdimensional parameter spaces, etc., as well as various applications of machine learning in these contexts. Such techniques are forming a methodological foundation for science with massive and complex data sets in general, and are likely to have a much broather impact on the modern society, commerce, information economy, security, etc. There is a powerful emerging synergy between the computationally enabled science and the science-driven computing, which will drive the progress in science, scholarship, and many other venues in the 21st century.

preprint2004arXiv

Black Holes From the Dark Ages: Exploring the Reionization Era and Early Structure Formation With Quasars and Gamma-Ray Bursts

The cosmic reionization era, which includes formation of the first stars, galaxies, and AGN, is now one of the most active frontiers of cosmological research. We review briefly our current understanding of the early structure formation, and use the ideas about a joint formation of massive black holes (which power the early QSOs) and their host galaxies to employ high-redshift QSOs as probes of the early galaxy formation and primordial large-scale structure. There is a growing evidence for a strong biasing in the formation of the first luminous sources, which would lead to a clumpy reionization. Absorption spectroscopy of QSOs at z > 6 indicates the end of the reionization era at z ~ 6; yet measurements from the WMAP satellite suggest and early reionization at z ~ 10 - 20. The first generation of massive stars, perhaps aided by the early mini-quasars, may have reionized the universe at such high redshifts, but their feedback may have disrupted the subsequent star and galaxy formation, leading to an extended and perhaps multimodal reionization history ending by z ~ 6. Observations of gamma-ray bursts from the death events of these putative Population III stars may provide essential insight into the primordial structure formation, reionization, early chemical enrichment, and formation of seed black holes which may grow to become central engines of luminous quasars.

preprint2003arXiv

The bright optical afterglow of the nearby gamma-ray burst of 29 March 2003

Many past studies of cosmological gamma-ray bursts (GRBs) have been limited because of the large distance to typical GRBs, resulting in faint afterglows. There has long been a recognition that a nearby GRB would shed light on the origin of these mysterious cosmic explosions, as well as the physics of their fireballs. However, GRBs nearer than z=0.2 are extremely rare, with an estimated rate of localisation of one every decade. Here, we report the discovery of bright optical afterglow emission from GRB 030329. Our prompt dissemination and the brilliance of the afterglow resulted in extensive followup (more than 65 telescopes) from radio through X-ray bands, as well as measurement of the redshift, z=0.169. The gamma-ray and afterglow properties of GRB 030329 are similar to those of cosmological GRBs (after accounting for the small distance), making this the nearest known cosmological GRB. Observations have already securely identified the progenitor as a massive star that exploded as a supernova, and we anticipate futher revelations of the GRB phenomenon from studies of this source.

preprint2003arXiv

The Northern Sky Optical Cluster Survey II: An Objective Cluster Catalog for 5800 Square Degrees

We present a new, objectively defined catalog of candidate galaxy clusters based on the galaxy catalogs from the Digitized Second Palomar Observatory Sky Survey (DPOSS). This cluster catalog, derived from the best calibrated plates in the high latitude (|b|>30) Northern Galactic Cap region, covers 5,800 square degrees, and contains 8,155 candidate clusters. A simple adaptive kernel density mapping technique, combined with the SExtractor object detection algorithm, is used to detect galaxy overdensities, which we identify as clusters. Simulations of the background galaxy distribution and clusters of varying richnesses and redshifts allow us to optimize detection parameters, and measure the completeness and contamination rates for our catalog. Cluster richnesses and photometric redshifts are measured, using integrated colors and magnitudes for each cluster. An extensive spectroscopic survey is used to confirm the photometric results. This catalog, with well-characterized sample properties, provides a sound basis for future studies of cluster physics and large scale structure.

preprint2002arXiv

Detection of a supernova signature associated with GRB 011121

Using observations from an extensive monitoring campaign with the Hubble Space Telescope we present the detection of an intermediate-time flux excess which is redder in color relative to the afterglow of GRB 011121, currently distinguished as the gamma-ray burst with the lowest known redshift. The red ``bump,'' which exhibits a spectral roll-over at ~7200 Angstrom, is well described by a redshifted Type Ic supernova that occurred approximately at the same time as the gamma-ray burst event. The inferred luminosity is about half that of the bright supernova 1998bw. These results serve as compelling evidence for a massive star origin of long-duration gamma-ray bursts. Models that posit a supernova explosion weeks to months preceding the gamma-ray burst event are excluded by these observations. Finally, we discuss the relationship between spherical core-collapse supernovae and gamma-ray bursts.

preprint2002arXiv

GRB 010921: Discovery of the First HETE Afterglow

We report the discovery of the optical and radio afterglow of GRB 010921, the first gamma-ray burst afterglow to be found from a localization by the High Energy Transient Explorer (HETE) satellite. We present optical spectroscopy of the host galaxy which we find to be a dusty and apparently normal star-forming galaxy at z = 0.451. The unusually steep optical spectral slope of the afterglow can be explained by heavy extinction, A_V > 0.5 mag, along the line of sight to the GRB. Dust with similar A_V for the the host galaxy as a whole appears to be required by the measurement of a Balmer decrement in the spectrum of the host galaxy. Thanks to the low redshift, continued observations of the afterglow will enable the strongest constraints, to date, on the existence of a possible underlying supernova.

preprint2002arXiv

GRB 010921: Strong Limits on an Underlying Supernova from HST

GRB 010921 was the first HETE-2 GRB to be localized via its afterglow emission. The low-redshift of the host galaxy, z=0.451, prompted us to undertake intensive multi-color observations with the Hubble Space Telescope with the goal of searching for an underlying supernova component. We do not detect any coincident supernova to a limit 1.34 mag fainter than SN 1998bw at 99.7% confidence, making this one of the most sensitive searches for an underlying SN. Analysis of the afterglow data allow us to infer that the GRB was situated behind a net extinction (Milky Way and the host galaxy) of A_V ~ 1.8 mag in the observer frame. Thus, had it not been for such heavy extinction our data would have allowed us to probe for an underlying SN with brightness approaching those of more typical Type Ib/c supernovae.

preprint2002arXiv

The Broadband Afterglow of GRB980329

We present radio observations of the afterglow of the bright gamma-ray burst GRB980329 made between one month and several years after the burst, a re-analysis of previously published submillimeter data, and late-time optical and near-infrared (NIR) observations of the host galaxy. From the absence of a spectral break in the optical/NIR colors of the host galaxy, we exclude the earlier suggestion that GRB980329 lies at a redshift of z >~5. We combine our data with the numerous multi-wavelength observations of the early afterglow, fit a comprehensive afterglow model to the entire broadband dataset, and derive fundamental physical parameters of the blast-wave and its host environment. Models for which the ejecta expand isotropically require both a high circumburst density and extreme radiative losses from the shock. No low density model (n << 10 cm^{-3}) fits the data. A burst with a total energy of ~ 10^{51} erg, with the ejecta narrowly collimated to an opening angle of a few degrees, driven into a surrounding medium with density ~ 20 cm^{-3}, provides a satisfactory fit to the lightcurves over a range of redshifts.

preprint2002arXiv

The Faint Optical Afterglow and Host Galaxy of GRB 020124: Implications for the Nature of Dark Gamma-Ray Bursts

We present ground-based optical observations of GRB 020124 starting 1.6 hours after the burst, as well as subsequent Very Large Array (VLA) and Hubble Space Telescope (HST) observations. The optical afterglow of GRB 020124 is one of the faintest afterglows detected to date, and it exhibits a relatively rapid decay, $F_ν\propto t^{-1.60\pm 0.04}$, followed by further steepening. In addition, a weak radio source was found coincident with the optical afterglow. The HST observations reveal that a positionally coincident host galaxy must be the faintest host to date, R>29.5 mag. The afterglow observations can be explained by several models requiring little or no extinction within the host galaxy, A_V~0-0.9 mag. These observations have significant implications for the statistics of the so-called dark bursts (bursts for which no optical afterglow is detected), which are usually attributed to dust extinction within the host galaxy. The faintness and relatively rapid decay of the afterglow of GRB 020124, combined with the low inferred extinction indicate that some dark bursts are intrinsically dim and not dust obscured. Thus, the diversity in the underlying properties of optical afterglows must be observationally determined before substantive inferences can be drawn from the statistics of dark bursts.

preprint2001arXiv

The GRB Host Galaxies and Redshifts

Observations of GRB host galaxies and their environments in general can provide valuable clues about the nature of progenitors. Bursts are associated with faint, <R> ~ 25 mag, galaxies at cosmological redshifts, <z> ~ 1. The host galaxies span a range of luminosities and morphologies, and appear to be broadly typical for the normal, evolving, actively star-forming galaxy populations at comparable redshifts and magnitudes, but may have somewhat elevated SFR per unit luminosity. There are also spectroscopic hints of massive star formation, from the ratios of [Ne III] and [O II] lines. The observed, unobscured star formation rates are typically a few M_sun/yr, but a considerable fraction of the total star formation in the hosts may be obscured by dust. A census of detected optical afterglows provides a powerful new handle on the obscured fraction of star formation in the universe; the current results suggest that at most a half of the massive star formation was hidden by dust.

preprint1999arXiv

The Northern Sky Optical Cluster Survey I: Detection of Galaxy Clusters in DPOSS

The Northern Sky Optical Cluster Survey is a project to create an objective catalog of galaxy clusters over the entire high-galactic-latitude Northern sky, with well understood selection criteria. We use the object catalogs generated from the Digitized Second Palomar Sky Survey (DPOSS, Djorgovski et al. 1999) as the basis for this survey. We apply a color criterion to select against field galaxies, and use a simple adaptive kernel technique to create galaxy density maps, combined with the bootstrap technique to make significance maps, from which density peaks are selected. This survey attempts to eliminate some of the subjective criteria and assumptions of past surveys, including detection by eye (Abell 1958, ACO 1989) and assumed luminosity functions and cluster profiles (PDCS, Postman et al. 1995). We also utilize more information (especially colors) than the most similar recent survey, the APM (Dalton et al. 1992). This paper presents the details of our cluster detection technique, as well as some initial results for two small areas totaling ~60 square degrees. We find a mean surface density of ~1.5 clusters per square degree, consistent with the detection of richness class 0 and higher clusters to z~0.3. In addition, we demonstrate an effective photometric redshift estimator for our clusters.

preprint1998arXiv

Spectroscopy of the Host Galaxy of the Gamma--Ray Burst 980703

We present spectroscopic observations of the host galaxy of the gamma-ray burst (GRB) 980703. Several emission and absorption features are detected, making the redshift, z = 0.966, completely unambiguous. This is only the third known redshift for a GRB host. The implied isotropic gamma-ray energy release from the burst is in excess of 10^{53} erg, for a reasonable choice of cosmological parameters. The spectroscopic properties of the host galaxy are typical for a star formation powered object. Using the observed value of the Balmer decrement, we derived the extinction in the galaxy's restframe, A_V = 0.3 +- 0.3 mag. Using three different star formation rate indicators, we estimate SFR ~ 10 Msun/yr, or higher, depending on the extinction, with a lower limit of SFR > 7 Msun/yr. This is the highest value of the star formation rate measured for a GRB galaxy so far, and it gives some support to the idea that GRBs are closely related to massive star formation.

preprint1997arXiv

Data-Mining a Large Digital Sky Survey: From the Challenges to the Scientific Results

The analysis and an efficient scientific exploration of the Digital Palomar Observatory Sky Survey (DPOSS) represents a major technical challenge. The input data set consists of 3 Terabytes of pixel information, and contains a few billion sources. We describe some of the specific scientific problems posed by the data, including searches for distant quasars and clusters of galaxies, and the data-mining techniques we are exploring in addressing them. Machine-assisted discovery methods may become essential for the analysis of such multi-Terabyte data sets. New and future approaches involve unsupervised classification and clustering analysis in the Giga-object data space, including various Bayesian techniques. In addition to the searches for known types of objects in this data base, these techniques may also offer the possibility of discovering previously unknown, rare types of astronomical objects.

preprint1996arXiv

Discovery of a Galaxy Responsible for a DLA System at z=3.15 and a Near-Infrared Search for Primeval Galaxies

We report the detection, both in line and continuum emission, of a galaxy responsible for a damped Ly-alpha absorption (DLA) system at z=3.15 projected 2.3 arcsec from QSO 2233+131. The star formation rate implied by both the Ly-alpha and restframe far-UV continuum is approximately 7 M_sun / yr. The galaxy also shows no sign of an active nucleus and has blue colors, both consistent with it being observed early in its star formation history. The Ly-alpha emission is offset approximately 200 km/s from the DLA velocity, which could be explained as the rotational signature of a disk galaxy. We also report on the progress of an ongoing near-infrared narrow-band search with the Keck telescope for line emission from a population of primeval galaxies at high redshift. We have found several promising candidates which are being followed up with optical spectroscopy. Our survey is beginning to place constraints on a possible population of slightly dusty primeval galaxies at high redshift.