Source author record

Joshua S. Bloom

Joshua S. Bloom 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

99works
14topics
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

99 published item(s)

preprint2023arXiv

Science with the Ultraviolet Explorer (UVEX)

UVEX is a proposed medium class Explorer mission designed to provide crucial missing capabilities that will address objectives central to a broad range of modern astrophysics. The UVEX design has two co-aligned wide-field imagers operating in the FUV and NUV and a powerful broadband medium resolution spectrometer. In its two-year baseline mission, UVEX will perform a multi-cadence synoptic all-sky survey 50/100 times deeper than GALEX in the NUV/FUV, cadenced surveys of the Large and Small Magellanic Clouds, rapid target of opportunity followup, as well as spectroscopic followup of samples of stars and galaxies. The science program is built around three pillars. First, UVEX will explore the low-mass, low-metallicity galaxy frontier through imaging and spectroscopic surveys that will probe key aspects of the evolution of galaxies by understanding how star formation and stellar evolution at low metallicities affect the growth and evolution of low-metallicity, low-mass galaxies in the local universe. Such galaxies contain half the mass in the local universe, and are analogs for the first galaxies, but observed at distances that make them accessible to detailed study. Second, UVEX will explore the dynamic universe through time-domain surveys and prompt spectroscopic followup capability will probe the environments, energetics, and emission processes in the early aftermaths of gravitational wave-discovered compact object mergers, discover hot, fast UV transients, and diagnose the early stages of stellar explosions. Finally, UVEX will become a key community resource by leaving a large all-sky legacy data set, enabling a wide range of scientific studies and filling a gap in the new generation of wide-field, sensitive optical and infrared surveys provided by the Rubin, Euclid, and Roman observatories. This paper discusses the scientific potential of UVEX, and the broad scientific program.

preprint2022arXiv

A Ubiquitous Unifying Degeneracy in Two-Body Microlensing Systems

While gravitational microlensing by planetary systems provides unique vistas on the properties of exoplanets, observations of a given 2-body microlensing event can often be interpreted with multiple distinct physical configurations. Such ambiguities are typically attributed to the close-wide and inner-outer types of degeneracies that arise from transformation invariances and symmetries of microlensing caustics. However, there remain unexplained inconsistencies between aforementioned theories and observations. Here, leveraging a fast machine learning inference framework, we present the discovery of the offset degeneracy, which concerns a magnification-matching behaviour on the lens-axis and is formulated independent of caustics. This offset degeneracy unifies the close-wide and inner-outer degeneracies, generalises to resonant topologies, and upon reanalysis, not only appears ubiquitous in previously published planetary events with 2-fold degenerate solutions, but also resolves prior inconsistencies. Our analysis demonstrates that degenerate caustics do not strictly result in degenerate magnifications and that the commonly invoked close-wide degeneracy essentially never arises in actual events. Moreover, it is shown that parameters in offset degenerate configurations are related by a simple expression. This suggests the existence of a deeper symmetry in the equations governing 2-body lenses than previously recognised.

preprint2022arXiv

HEALPix Alchemy: Fast All-Sky Geometry and Image Arithmetic in a Relational Database for Multimessenger Astronomy Brokers

Efficient searches for electromagnetic counterparts to gravitational wave, high-energy neutrino, and gamma-ray burst events demand rapid processing of image arithmetic and geometry set operations in a database to cross-match galaxy catalogs, observation footprints, and all-sky images. Here we introduce HEALPix Alchemy, an open-source, pure Python implementation of a set of methods that enables rapid all-sky geometry calculations. HEALPix Alchemy is built upon HEALPix, a spatial indexing strategy that is widely used in astronomical databases as well as the native format of LIGO-Virgo-KAGRA gravitational-wave sky localization maps. Our approach leverages new multirange types built into the PostgreSQL 14 database engine. This enables fast all-sky queries against probabilistic multimessenger event localizations and telescope survey footprints. Questions such as "What are the galaxies contained within the 90% credible region of an event?" and "What is the rank-ordered list of the fields within an observing footprint with the highest probability of containing the event?" can be performed in less than a few seconds on commodity hardware using off-the-shelf cloud-managed database implementations without server-side database extensions. Common queries scale roughly linearly with the number of telescope pointings. As the number of fields grows into the hundreds or thousands, HEALPix Alchemy is orders of magnitude faster than other implementations. HEALPix Alchemy is now used as the spatial geometry engine within SkyPortal, which forms the basis of the Zwicky Transient Facility transient marshal, called Fritz.

preprint2022arXiv

In search of short gamma-ray burst optical counterpart with the Zwicky Transient Facility

The Fermi Gamma-ray Burst Monitor (GBM) triggers on-board in response to $\sim$ 40 short gamma-ray bursts (SGRBs) per year; however, their large localization regions have made the search for optical counterparts a challenging endeavour. We have developed and executed an extensive program with the wide field of view of the Zwicky Transient Facility (ZTF) camera, mounted on the Palomar 48 inch Oschin telescope (P48), to perform target-of-opportunity (ToO) observations on 10 Fermi-GBM SGRBs during 2018 and 2020-2021. Bridging the large sky areas with small field of view optical telescopes in order to track the evolution of potential candidates, we look for the elusive SGRB afterglows and kilonovae (KNe) associated with these high-energy events. No counterpart has yet been found, even though more than 10 ground based telescopes, part of the Global Relay of Observatories Watching Transients Happen (GROWTH) network, have taken part in these efforts. The candidate selection procedure and the follow-up strategy have shown that ZTF is an efficient instrument for searching for poorly localized SGRBs, retrieving a reasonable number of candidates to follow-up and showing promising capabilities as the community approaches the multi-messenger era. Based on the median limiting magnitude of ZTF, our searches would have been able to retrieve a GW170817-like event up to $\sim$ 200 Mpc and SGRB afterglows to z = 0.16 or 0.4, depending on the assumed underlying energy model. Future ToOs will expand the horizon to z = 0.2 and 0.7 respectively.

preprint2022arXiv

Target of Opportunity Observations of Gravitational Wave Events with Vera C. Rubin Observatory

The discovery of the electromagnetic counterpart to the binary neutron star merger GW170817 has opened the era of gravitational-wave multi-messenger astronomy. Rapid identification of the optical/infrared kilonova enabled a precise localization of the source, which paved the way to deep multi-wavelength follow-up and its myriad of related science results. Fully exploiting this new territory of exploration requires the acquisition of electromagnetic data from samples of neutron star mergers and other gravitational wave sources. After GW170817, the frontier is now to map the diversity of kilonova properties and provide more stringent constraints on the Hubble constant, and enable new tests of fundamental physics. The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) can play a key role in this field in the 2020s, when an improved network of gravitational-wave detectors is expected to reach a sensitivity that will enable the discovery of a high rate of merger events involving neutron stars (about tens per year) out to distances of several hundred Mpc. We design comprehensive target-of-opportunity observing strategies for follow-up of gravitational-wave triggers that will make the Rubin Observatory the premier instrument for discovery and early characterization of neutron star and other compact object mergers, and yet unknown classes of gravitational wave events.

preprint2021arXiv

Automating Inference of Binary Microlensing Events with Neural Density Estimation

Automated inference of binary microlensing events with traditional sampling-based algorithms such as MCMC has been hampered by the slowness of the physical forward model and the pathological likelihood surface. Current analysis of such events requires both expert knowledge and large-scale grid searches to locate the approximate solution as a prerequisite to MCMC posterior sampling. As the next generation, space-based microlensing survey with the Roman Space Observatory is expected to yield thousands of binary microlensing events, a new scalable and automated approach is desired. Here, we present an automated inference method based on neural density estimation (NDE). We show that the NDE trained on simulated Roman data not only produces fast, accurate, and precise posteriors but also captures expected posterior degeneracies. A hybrid NDE-MCMC framework can further be applied to produce the exact posterior.

preprint2021arXiv

The Gravity Collective: A Search for the Electromagnetic Counterpart to the Neutron Star-Black Hole Merger GW190814

We present optical follow-up imaging obtained with the Katzman Automatic Imaging Telescope, Las Cumbres Observatory Global Telescope Network, Nickel Telescope, Swope Telescope, and Thacher Telescope of the LIGO/Virgo gravitational wave (GW) signal from the neutron star-black hole (NSBH) merger GW190814. We searched the GW190814 localization region (19 deg$^{2}$ for the 90th percentile best localization), covering a total of 51 deg$^{2}$ and 94.6% of the two-dimensional localization region. Analyzing the properties of 189 transients that we consider as candidate counterparts to the NSBH merger, including their localizations, discovery times from merger, optical spectra, likely host-galaxy redshifts, and photometric evolution, we conclude that none of these objects are likely to be associated with GW190814. Based on this finding, we consider the likely optical properties of an electromagnetic counterpart to GW190814, including possible kilonovae and short gamma-ray burst afterglows. Using the joint limits from our follow-up imaging, we conclude that a counterpart with an $r$-band decline rate of 0.68 mag day$^{-1}$, similar to the kilonova AT 2017gfo, could peak at an absolute magnitude of at most $-17.8$ mag (50% confidence). Our data are not constraining for ''red'' kilonovae and rule out ''blue'' kilonovae with $M>0.5 M_{\odot}$ (30% confidence). We strongly rule out all known types of short gamma-ray burst afterglows with viewing angles $<$17$^{\circ}$ assuming an initial jet opening angle of $\sim$$5.2^{\circ}$ and explosion energies and circumburst densities similar to afterglows explored in the literature. Finally, we explore the possibility that GW190814 merged in the disk of an active galactic nucleus, of which we find four in the localization region, but we do not find any candidate counterparts among these sources.

preprint2020arXiv

Deep Generative Modeling of Periodic Variable Stars Using Physical Parameters

The ability to generate physically plausible ensembles of variable sources is critical to the optimization of time-domain survey cadences and the training of classification models on datasets with few to no labels. Traditional data augmentation techniques expand training sets by reenvisioning observed exemplars, seeking to simulate observations of specific training sources under different (exogenous) conditions. Unlike fully theory-driven models, these approaches do not typically allow principled interpolation nor extrapolation. Moreover, the principal drawback of theory-driven models lies in the prohibitive computational cost of simulating source observables from {\it ab initio} parameters. In this work, we propose a computationally tractable machine learning approach to generate realistic light curves of periodic variables capable of integrating physical parameters and variability classes as inputs. Our deep generative model, inspired by the Transparent Latent Space Generative Adversarial Networks (TL-GANs), uses a Variational Autoencoder (VAE) architecture with Temporal Convolutional Network (TCN) layers, trained using the \hbox{OGLE-III} optical light curves and physical characteristics (e.g., effective temperature and absolute magnitude) from Gaia DR2. A test using the temperature-shape relationship of RR\,Lyrae demonstrates the efficacy of our generative "Physics-Enhanced Latent Space VAE" (PELS-VAE) model. Such deep generative models, serving as non-linear non-parametric emulators, present a novel tool for astronomers to create synthetic time series over arbitrary cadences.

preprint2020arXiv

deepSIP: Linking Type Ia Supernova Spectra to Photometric Quantities with Deep Learning

We present {\tt deepSIP} (deep learning of Supernova Ia Parameters), a software package for measuring the phase and -- for the first time using deep learning -- the light-curve shape of a Type Ia supernova (SN~Ia) from an optical spectrum. At its core, {\tt deepSIP} consists of three convolutional neural networks trained on a substantial fraction of all publicly-available low-redshift SN~Ia optical spectra, onto which we have carefully coupled photometrically-derived quantities. We describe the accumulation of our spectroscopic and photometric datasets, the cuts taken to ensure quality, and our standardised technique for fitting light curves. These considerations yield a compilation of 2754 spectra with photometrically characterised phases and light-curve shapes. Though such a sample is significant in the SN community, it is small by deep-learning standards where networks routinely have millions or even billions of free parameters. We therefore introduce a data-augmentation strategy that meaningfully increases the size of the subset we allocate for training while prioritising model robustness and telescope agnosticism. We demonstrate the effectiveness of our models by deploying them on a sample unseen during training and hyperparameter selection, finding that Model~I identifies spectra that have a phase between $-10$ and 18\,d and light-curve shape, parameterised by $Δm_{15}$, between 0.85 and 1.55\,mag with an accuracy of 94.6\%. For those spectra that do fall within the aforementioned region in phase--$Δm_{15}$ space, Model~II predicts phases with a root-mean-square error (RMSE) of 1.00\,d and Model~III predicts $Δm_{15}$ values with an RMSE of 0.068\,mag.

preprint2020arXiv

Kilonova Luminosity Function Constraints based on Zwicky Transient Facility Searches for 13 Neutron Star Mergers

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)

preprint2020arXiv

Modelling the prompt optical emission of GRB 180325A: the evolution of a spike from the optical to gamma-rays

The transition from prompt to the afterglow emission is one of the most exciting and least understood phases in gamma-ray bursts (GRBs). Correlations among optical, X-ray and gamma-ray emission in GRBs have been explored, to attempt to answer whether the earliest optical emission comes from internal and/or external shocks. We present optical photometric observations of GRB 180325A collected with the TAROT and RATIR ground-based telescopes. These observations show two strong optical flashes with separate peaks at $\sim50\;$s and $\sim120\;$s, followed by a temporally extended optical emission. We also present X-rays and gamma-ray observations of GRB 180325A, detected by the Burst Alert Telescope (BAT) and X-ray Telescope (XRT), on the Neil Gehrels Swift observatory, which both observed a narrow flash at $\sim80\;$s. We show that the prompt gamma- and X-ray early emission shares similar temporal and spectral features consistent with internal dissipation within the relativistic outflow (e.g. by internal shocks or magnetic reconnection), while the early optical flashes are likely generated by the reverse shock that decelerates the ejecta as it sweeps up the external medium.

preprint2019arXiv

deepCR: Cosmic Ray Rejection with Deep Learning

Cosmic ray (CR) identification and replacement are critical components of imaging and spectroscopic reduction pipelines involving solid-state detectors. We present deepCR, a deep learning based framework for CR identification and subsequent image inpainting based on the predicted CR mask. To demonstrate the effectiveness of this framework, we train and evaluate models on Hubble Space Telescope ACS/WFC images of sparse extragalactic fields, globular clusters, and resolved galaxies. We demonstrate that at a false positive rate of 0.5%, deepCR achieves close to 100% detection rates in both extragalactic and globular cluster fields, and 91% in resolved galaxy fields, which is a significant improvement over the current state-of-the-art method LACosmic. Compared to a multicore CPU implementation of LACosmic, deepCR CR mask predictions run up to 6.5 times faster on CPU and 90 times faster on a single GPU. For image inpainting, the mean squared errors of deepCR predictions are 20 times lower in globular cluster fields, 5 times lower in resolved galaxy fields, and 2.5 times lower in extragalactic fields, compared to the best performing non-neural technique tested. We present our framework and the trained models as an open-source Python project, with a simple-to-use API. To facilitate reproducibility of the results we also provide a benchmarking codebase.

preprint2019arXiv

GRB 180620A: Evidence for late-time energy injection

The early optical emission of gamma-ray bursts gives an opportunity to understand the central engine and first stages of these events. About 30\% of GRBs present flares whose origin is still a subject of discussion. We present optical photometry of GRB 180620A with the COATLI telescope and RATIR instrument. COATLI started to observe from the end of prompt emission at $T+39.3$~s and RATIR from $T+121.4$~s. We supplement the optical data with the X-ray light curve from \emph{Swift}/XRT. %The optical and X-ray light curves show very unusual behavior with features clearly beyond the standard fireball model. We observe an optical flare from $T+110$ to $T+550$~s, with a temporal index decay $α_\mathrm{O,decay}=1.32\pm 0.01$, and a $Δt/t=1.63$, which we interpret as the signature of a reverse shock component. After the initial normal decay the light curves show a long plateau from $T+500$ to $T+7800$~s both in X-rays and the optical before decaying again after an achromatic jet break at $T+7800$~s. Fluctuations are seen during the plateau phase in the optical. Adding to the complexity of GRB afterglows, the plateau phase (typically associated with the coasting phase of the jet) is seen in this object after the ``normal'' decay phase (emitted during the deceleration phase of the jet) and the jet break phase occurs directly after the plateau. We suggest that this sequence of events can be explained by a rapid deceleration of the jet with $t_d\lesssim 40$ s due to the high density of the environment ($\approx 100$ cm$^{-3}$) followed by reactivation of the central engine which causes the flare and powers the plateau phase.

preprint2019arXiv

GROWTH on S190425z: Searching thousands of square degrees to identify an optical or infrared counterpart to a binary neutron star merger with the Zwicky Transient Facility and Palomar Gattini IR

The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.

preprint2019arXiv

GROWTH on S190814bv: Deep Synoptic Limits on the Optical/Near-Infrared Counterpart to a Neutron Star-Black Hole Merger

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $κ< 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ χ< 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv.

preprint2016arXiv

cesium: Open-Source Platform for Time-Series Inference

Inference on time series data is a common requirement in many scientific disciplines and internet of things (IoT) applications, yet there are few resources available to domain scientists to easily, robustly, and repeatably build such complex inference workflows: traditional statistical models of time series are often too rigid to explain complex time domain behavior, while popular machine learning packages require already-featurized dataset inputs. Moreover, the software engineering tasks required to instantiate the computational platform are daunting. cesium is an end-to-end time series analysis framework, consisting of a Python library as well as a web front-end interface, that allows researchers to featurize raw data and apply modern machine learning techniques in a simple, reproducible, and extensible way. Users can apply out-of-the-box feature engineering workflows as well as save and replay their own analyses. Any steps taken in the front end can also be exported to a Jupyter notebook, so users can iterate between possible models within the front end and then fine-tune their analysis using the additional capabilities of the back-end library. The open-source packages make us of many use modern Python toolkits, including xarray, dask, Celery, Flask, and scikit-learn.

preprint2015arXiv

CfAIR2: Near Infrared Light Curves of 94 Type Ia Supernovae

CfAIR2 is a large homogeneously reduced set of near-infrared (NIR) light curves for Type Ia supernovae (SN Ia) obtained with the 1.3m Peters Automated InfraRed Imaging TELescope (PAIRITEL). This data set includes 4607 measurements of 94 SN Ia and 4 additional SN Iax observed from 2005-2011 at the Fred Lawrence Whipple Observatory on Mount Hopkins, Arizona. CfAIR2 includes JHKs photometric measurements for 88 normal and 6 spectroscopically peculiar SN Ia in the nearby universe, with a median redshift of z~0.021 for the normal SN Ia. CfAIR2 data span the range from -13 days to +127 days from B-band maximum. More than half of the light curves begin before the time of maximum and the coverage typically contains ~13-18 epochs of observation, depending on the filter. We present extensive tests that verify the fidelity of the CfAIR2 data pipeline, including comparison to the excellent data of the Carnegie Supernova Project. CfAIR2 contributes to a firm local anchor for supernova cosmology studies in the NIR. Because SN Ia are more nearly standard candles in the NIR and are less vulnerable to the vexing problems of extinction by dust, CfAIR2 will help the supernova cosmology community develop more precise and accurate extragalactic distance probes to improve our knowledge of cosmological parameters, including dark energy and its potential time variation.

preprint2015arXiv

iPTF14yb: The First Discovery of a GRB Afterglow Independent of a High-Energy Trigger

We report here the discovery by the Intermediate Palomar Transient Factory (iPTF) of iPTF14yb, a luminous ($M_{r}\approx-27.8$ mag), cosmological (redshift 1.9733), rapidly fading optical transient. We demonstrate, based on probabilistic arguments and a comparison with the broader population, that iPTF14yb is the optical afterglow of the long-duration gamma-ray burst GRB 140226A. This marks the first unambiguous discovery of a GRB afterglow prior to (and thus entirely independent of) an associated high-energy trigger. We estimate the rate of iPTF14yb-like sources (i.e., cosmologically distant relativistic explosions) based on iPTF observations, inferring an all-sky value of $\Re_{\mathrm{rel}}=610$ yr$^{-1}$ (68% confidence interval of 110-2000 yr$^{-1}$). Our derived rate is consistent (within the large uncertainty) with the all-sky rate of on-axis GRBs derived by the Swift satellite. Finally, we briefly discuss the implications of the nondetection to date of bona fide "orphan" afterglows (i.e., those lacking detectable high-energy emission) on GRB beaming and the degree of baryon loading in these relativistic jets.

preprint2015arXiv

The detection rate of early UV emission from supernovae: A dedicated GALEX/PTF survey and calibrated theoretical estimates

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

preprint2014arXiv

A Continuum of H- to He-Rich Tidal Disruption Candidates With a Preference for E+A Galaxies

We present the results of a Palomar Transient Factory (PTF) archival search for blue transients which lie in the magnitude range between "normal" core-collapse and superluminous supernovae (i.e. with $-21\,{\leq}M_{R\,(peak)}\,{\leq}-19$). Of the six events found after excluding all interacting Type~IIn and Ia-CSM supernovae, three (PTF09ge, 09axc and 09djl) are coincident with the centers of their hosts, one (10iam) is offset from the center, and for two (10nuj and 11glr) a precise offset can not be determined. All the central events have similar rise times to the He-rich tidal disruption candidate PS1-10jh, and the event with the best-sampled light curve also has similar colors and power-law decay. Spectroscopically, PTF09ge is He-rich, while PTF09axc and 09djl display broad hydrogen features around peak magnitude. All three central events are in low star-formation hosts, two of which are E+A galaxies. Our spectrum of the host of PS1-10jh displays similar properties. PTF10iam, the one offset event, is different photometrically and spectroscopically from the central events and its host displays a higher star formation rate. Finding no obvious evidence for ongoing galactic nuclei activity or recent star formation, we conclude that the three central transients likely arise from the tidal disruption of a star by a super-massive black hole. We compare the spectra of these events to tidal disruption candidates from the literature and find that all of these objects can be unified on a continuous scale of spectral properties. The accumulated evidence of this expanded sample strongly supports a tidal disruption origin for this class of nuclear transients.

preprint2014arXiv

Mid-Infrared Period--Luminosity Relations of RR Lyrae Stars Derived from the AllWISE Data Release

We use photometry from the recent AllWISE Data Release of the Wide-field In-frared Survey Explorer (WISE) of 129 calibration stars, combined with prior distances obtained from the established $M_V-$[Fe/H] relation and Hubble Space Telescope trigonometric parallax, to derive mid-infrared period--luminosity relations for RR Lyrae pulsating variable stars. We derive relations in the W1, W2, and W3 wavebands (3.4, 4.6, and 12 μm, respectively), and for each of the two main RR Lyrae sub-types (RRab and RRc). We report an error on the period--luminosity relation slope for RRab stars of 0.2. We also fit posterior distances for the calibration catalog and find a median fractional distance error of 0.8 per cent.

preprint2014arXiv

Probing the distance and morphology of the Large Magellanic Cloud with RR Lyrae stars

We present a Bayesian analysis of the distances to 15,040 Large Magellanic Cloud (LMC) RR Lyrae stars using $V$- and $I$-band light curves from the Optical Gravitational Lensing Experiment, in combination with new $z$-band observations from the Dark Energy Camera. Our median individual RR Lyrae distance statistical error is 1.89 kpc (fractional distance error of 3.76 per cent). We present three-dimensional contour plots of the number density of LMC RR Lyrae stars and measure a distance to the core LMC RR Lyrae centre of ${50.2482\pm0.0546 {\rm(statistical)} \pm0.4628 {\rm(systematic)} {\rm kpc}}$, equivalently ${μ_{\rm LMC}=18.5056\pm0.0024 {\rm(statistical)} \pm0.02 {\rm(systematic)}}$. This finding is statistically consistent with and four times more precise than the canonical value determined by a recent meta-analysis of 233 separate LMC distance determinations. We also measure a maximum tilt angle of $11.84^{\circ}\pm0.80^{\circ}$ at a position angle of $62^\circ$, and report highly precise constraints on the $V$, $I$, and $z$ RR Lyrae period--magnitude relations. The full dataset of observed mean-flux magnitudes, derived colour excess ${E(V-I)}$ values, and fitted distances for the 15,040 RR Lyrae stars produced through this work is made available through the publication's associated online data.

preprint2014arXiv

Slow-Speed Supernovae from the Palomar Transient Factory: Two Channels

Since the discovery of the unusual prototype SN 2002cx, the eponymous class of low-velocity, hydrogen-poor supernovae has grown to include at most another two dozen members identified from several heterogeneous surveys, in some cases ambiguously. Here we present the results of a systematic study of 1077 hydrogen-poor supernovae discovered by the Palomar Transient Factory, leading to nine new members of this peculiar class. Moreover we find there are two distinct subclasses based on their spectroscopic, photometric, and host galaxy properties: The "SN 2002cx-like" supernovae tend to be in later-type or more irregular hosts, have more varied and generally dimmer luminosities, have longer rise times, and lack a Ti II trough when compared to the "SN 2002es-like" supernovae. None of our objects show helium, and we counter a previous claim of two such events. We also find that these transients comprise 5.6+17-3.7% (90% confidence) of all SNe Ia, lower compared to earlier estimates. Combining our objects with the literature sample, we propose that these subclasses have two distinct physical origins.

preprint2014arXiv

The Highly-Eccentric Detached Eclipsing Binaries in ACVS and MACC

Next-generation synoptic photometric surveys will yield unprecedented (for the astronomical community) volumes of data and the processes of discovery and rare-object identification are, by necessity, becoming more autonomous. Such autonomous searches can be used to find objects of interest applicable to a wide range of outstanding problems in astronomy, and in this paper we present the methods and results of a largely autonomous search for highly eccentric detached eclipsing binary systems in the Machine-learned ASAS Classification Catalog. 106 detached eclipsing binaries with eccentricities greater than 0.1 are presented, most of which are identified here for the first time. We also present new radial-velocity curves and absolute parameters for 6 of those systems with the long-term goal of increasing the number of highly eccentric systems with orbital solutions, thereby facilitating further studies of the tidal circularization process in binary stars.

preprint2014arXiv

Towards precision distances and 3D dust maps using broadband Period--Magnitude relations of RR Lyrae stars

We determine the period-magnitude relations of RR Lyrae stars in 13 photometric bandpasses from 0.4 to 12 μm using timeseries observations of 134 stars. The Bayesian formalism, extended from our previous work to include the effects of line-of-sight dust extinction, allows for the simultaneous inference of the posterior distribution of the mean absolute magnitude, slope of the period-magnitude power-law, and intrinsic scatter about a perfect power-law for each bandpass. In addition, the distance modulus and line-of-sight dust extinction to each RR Lyrae star in the calibration sample is determined, yielding a sample median fractional distance error of 0.66%. The intrinsic scatter in all bands appears to be larger than the photometric errors, except in WISE W1 (3.4 μm) and W2 (4.6 μm) where the photometric error ($σ\approx 0.05$ mag) is to be comparable or larger than the intrinsic scatter. Additional observations at these wavelengths could improve the inferred distances to these sources further. As an application of the methodology, we infer the distance to the RRc-type star RZCep at low Galactic latitude ($b = 5.5^\circ$) to be $μ=8.0397\pm0.0123$ mag ($405.4\pm2.3$ pc) with colour excess $E(B-V)=0.2461\pm0.0089$ mag. This distance, equivalent to a parallax of $2467\pm14$ microarcsec, is consistent with the published HST parallax measurement but with an uncertainty that is 13 times smaller than the HST measurement. If our measurements (and methodology) hold up to scrutiny, the distances to these stars have been determined to an accuracy comparable to those expected with Gaia. As RR Lyrae are one of the primary components of the cosmic distance ladder, the achievement of sub-1% distance errors within a formalism that accounts for dust extinction may be considered a strong buttressing of the path to eventual 1% uncertainties in Hubble's constant.

preprint2013arXiv

A Mid-infrared Study of RR Lyrae Stars with the WISE All-Sky Data Release

We present a group of 3740 previously identified RR Lyrae variables well-observed with the Wide-field Infrared Survey Explorer (WISE). We explore how the shape of the generic RR Lyrae mid-infrared light curve evolves in period-space, comparing light curves in mid-infrared and optical bands. We find that optical light curves exhibit high amplitudes and a large spectrum of light curve shapes, while mid-infrared light curves have low amplitudes and uniform light curve shapes. From the period-space analysis, we hope to improve the classification methods of RR Lyrae variables and enable reliable discovery of these pulsators in the WISE catalog and future mid-infrared surveys such as the James Webb Space Telescope (JWST). We provide mid-infrared templates for typical RR Lyrae stars and demonstrate how these templates can be applied to improve estimates of mid-infrared RR Lyrae mean magnitude, which is used for distance measurement. This method of template fitting is particularly beneficial for improving observational efficiency. For example, using light curves with observational noise of 0.05 mag, we obtain the same level of accuracy in mean magnitude estimates for light curves randomly sampled at 12 data points with template fitting as with light curves randomly sampled at 20 data points with harmonic modelling.

preprint2013arXiv

A Study of the Diverse T Dwarf Population Revealed by WISE

We report the discovery of 87 new T dwarfs uncovered with the Wide-field Infrared Survey Explorer (WISE) and three brown dwarfs with extremely red near-infrared colors that exhibit characteristics of both L and T dwarfs. Two of the new T dwarfs are likely binaries with L7+/-1 primaries and mid-type T secondaries. In addition, our follow-up program has confirmed 10 previously identified T dwarfs and four photometrically-selected L and T dwarf candidates in the literature. This sample, along with the previous WISE discoveries, triples the number of known brown dwarfs with spectral types later than T5. Using the WISE All-Sky Source Catalog we present updated color-color and color-type diagrams for all the WISE-discovered T and Y dwarfs. Near-infrared spectra of the new discoveries are presented, along with spectral classifications. To accommodate later T dwarfs we have modified the integrated flux method of determining spectral indices to instead use the median flux. Furthermore, a newly defined J-narrow index differentiates the early-type Y dwarfs from late-type T dwarfs based on the J-band continuum slope. The K/J indices for this expanded sample show that 32% of late-type T dwarfs have suppressed K-band flux and are blue relative to the spectral standards, while only 11% are redder than the standards. Comparison of the Y/J and K/J index to models suggests diverse atmospheric conditions and supports the possible re-emergence of clouds after the L/T transition. We also discuss peculiar brown dwarfs and candidates that were found not to be substellar, including two Young Stellar Objects and two Active Galactic Nuclei. The coolest WISE-discovered brown dwarfs are the closest of their type and will remain the only sample of their kind for many years to come.

preprint2013arXiv

Direct evidence for a supernova interacting with a large amount of hydrogen-free circumstellar material

We present our observations of SN 2010mb, a Type Ic SN lacking spectroscopic signatures of H and He. SN 2010mb has a slowly-declining light curve ($\sim600\,$days) that cannot be powered by $^{56}$Ni/$^{56}$Co radioactivity, the common energy source for Type Ic SNe. We detect signatures of interaction with hydrogen-free CSM including a blue quasi-continuum and, uniquely, narrow oxygen emission lines that require high densities ($\sim10^9$cm$^{-3}$). From the observed spectra and light curve we estimate that the amount of material involved in the interaction was $\sim3$M$_{\odot}$. Our observations are in agreement with models of pulsational pair-instability SNe described in the literature.

preprint2013arXiv

Discovery and Early Multi-Wavelength Measurements of the Energetic Type Ic Supernova PTF12gzk: A Massive-Star Explosion in a Dwarf Host Galaxy

We present the discovery and extensive early-time observations of the Type Ic supernova (SN) PTF12gzk. Our finely sampled light curves show a rise of 0.8mag within 2.5hr. Power-law fits [f(t)\sim(t-t_0)^n] to these data constrain the explosion date to within one day. We cannot rule out the expected quadratic fireball model, but higher values of n are possible as well for larger areas in the fit parameter space. Our bolometric light curve and a dense spectral sequence are used to estimate the physical parameters of the exploding star and of the explosion. We show that the photometric evolution of PTF12gzk is slower than that of most SNe Ic, and its high ejecta velocities (~30,000km/s four days after explosion) are closer to the observed velocities of broad-lined SNe Ic associated with gamma-ray bursts (GRBs) than to the observed velocities in normal Type Ic SNe. The high velocities are sustained through the SN early evolution, and are similar to those of GRB-SNe when the SN reach peak magnitude. By comparison with the spectroscopically similar SN 2004aw, we suggest that the observed properties of PTF12gzk indicate an initial progenitor mass of 25-35 solar mass and a large (5-10E51 erg) kinetic energy, close to the regime of GRB-SN properties. The host-galaxy characteristics are consistent with GRB-SN hosts, and not with normal SN Ic hosts as well, yet this SN does not show the broad lines over extended periods of time that are typical of broad-line Type Ic SNe.

preprint2013arXiv

Discovery and redshift of an optical afterglow in 71 square degrees: iPTF13bxl and GRB 130702A

We report the discovery of the optical afterglow of the gamma-ray burst (GRB) 130702A, identified upon searching 71 square degrees surrounding the Fermi Gamma-ray Burst Monitor (GBM) localization. Discovered and characterized by the intermediate Palomar Transient Factory (iPTF), iPTF13bxl is the first afterglow discovered solely based on a GBM localization. Real-time image subtraction, machine learning, human vetting, and rapid response multi-wavelength follow-up enabled us to quickly narrow a list of 27,004 optical transient candidates to a single afterglow-like source. Detection of a new, fading X-ray source by Swift and a radio counterpart by CARMA and the VLA confirmed the association between iPTF13bxl and GRB 130702A. Spectroscopy with the Magellan and Palomar 200-inch telescopes showed the afterglow to be at a redshift of z=0.145, placing GRB 130702A among the lowest redshift GRBs detected to date. The prompt gamma-ray energy release and afterglow luminosity are intermediate between typical cosmological GRBs and nearby sub-luminous events such as GRB 980425 and GRB 060218. The bright afterglow and emerging supernova offer an opportunity for extensive panchromatic follow-up. Our discovery of iPTF13bxl demonstrates the first observational proof-of-principle for ~10 Fermi-iPTF localizations annually. Furthermore, it represents an important step towards overcoming the challenges inherent in uncovering faint optical counterparts to comparably localized gravitational wave events in the Advanced LIGO and Virgo era.

preprint2013arXiv

Discovery of a Cosmological, Relativistic Outburst via its Rapidly Fading Optical Emission

We report the discovery by the Palomar Transient Factory (PTF) of the transient source PTF11agg, which is distinguished by three primary characteristics: (1) bright, rapidly fading optical transient emission; (2) a faint, blue quiescent optical counterpart; and (3) an associated year-long, scintillating radio transient. We argue that these observed properties are inconsistent with any known class of Galactic transients, and instead suggest a cosmological origin. The detection of incoherent radio emission at such distances implies a large emitting region, from which we infer the presence of relativistic ejecta. The observed properties are all consistent with the population of long-duration gamma-ray bursts (GRBs), marking the first time such an outburst has been discovered in the distant universe independent of a high-energy trigger. We searched for possible high-energy counterparts to PTF11agg, but found no evidence for associated prompt emission. We therefore consider three possible scenarios to account for a GRB-like afterglow without a high-energy counterpart: an "untriggered" GRB (lack of satellite coverage), an "orphan" afterglow (viewing-angle effects), and a "dirty fireball" (suppressed high-energy emission). The observed optical and radio light curves appear inconsistent with even the most basic predictions for off-axis afterglow models. The simplest explanation, then, is that PTF11agg is a normal, on-axis long-duration GRB for which the associated high-energy emission was simply missed. However, we have calculated the likelihood of such a serendipitous discovery by PTF and find that it is quite small (~ 2.6%). While not definitive, we nonetheless speculate that PTF11agg may represent a new, more common (> 4 times the on-axis GRB rate at 90% confidence) class of relativistic outbursts lacking associated high-energy emission.

preprint2013arXiv

Genetics of single-cell protein abundance variation in large yeast populations

Many DNA sequence variants influence phenotypes by altering gene expression. Our understanding of these variants is limited by sample sizes of current studies and by measurements of mRNA rather than protein abundance. We developed a powerful method for identifying genetic loci that influence protein expression in very large populations of the yeast Saccharomyes cerevisiae. The method measures single-cell protein abundance through the use of green-fluorescent-protein tags. We applied this method to 160 genes and detected many more loci per gene than previous studies. We also observed closer correspondence between loci that influence protein abundance and loci that influence mRNA abundance of a given gene. Most loci cluster at hotspot locations that influence multiple proteins - in some cases, more than half of those examined. The variants that underlie these hotspots have profound effects on the gene regulatory network and provide insights into genetic variation in cell physiology between yeast strains.

preprint2013arXiv

Keck Observations of 160 Gamma-Ray Burst Host Galaxies

We present a preliminary data release from our multi-year campaign at Keck Observatory to study the host galaxies of a large sample of Swift-era gamma-ray bursts via multi-color ground-based optical imaging and spectroscopy. With over 160 targets observed to date (and almost 100 host detections, most of which have not previously been reported in the literature) our effort represents the broadest GRB host survey to date. While targeting was heterogeneous, our observations span the known diversity of GRBs including short bursts, long bursts, spectrally soft GRBs (XRFs), ultra-energetic GRBs, X-ray faint GRBs, dark GRBs, SN-GRBs, and other sub-classes. We also present a preview of our database (currently available online via a convenient web interface) including a catalog of multi-color photometry, redshifts and line ID's. Final photometry and reduced imaging and spectra will be available in the near future.

preprint2013arXiv

Nearby M, L, and T Dwarfs Discovered by the Wide-field Infrared Survey Explorer (WISE)

In our effort to complete the census of low-mass stars and brown dwarfs in the immediate Solar Neighborhood, we present spectra, photometry, proper motions, and distance estimates for forty-two low-mass star and brown dwarf candidates discovered by the Wide-field Infrared Survey Explorer (WISE). We also present additional follow-up information on twelve candidates selected using WISE data but previously published elsewhere. The new discoveries include fifteen M dwarfs, seventeen L dwarfs, five T dwarfs, and five objects of other type. Among these discoveries is a newly identified "unusually red L dwarf" (WISE J223527.07+451140.9), four peculiar L dwarfs whose spectra are most readily explained as unresolved L+T binary systems, and a T9 dwarf (WISE J124309.61+844547.8). We also show that the recently discovered red L dwarf WISEP J004701.06+680352.1 (Gizis et al. 2012) may be a low-gravity object and hence young and potentially low mass (< 25 MJup).

preprint2013arXiv

SN 2000cx and SN 2013bh: Extremely Rare, Nearly Twin Type Ia Supernovae

The Type Ia supernova (SN Ia) SN 2000cx was one of the most peculiar transients ever discovered, with a rise to maximum brightness typical of a SN Ia, but a slower decline and a higher photospheric temperature. Thirteen years later SN 2013bh (aka iPTF13abc), a near identical twin, was discovered and we obtained optical and near-IR photometry and low-resolution optical spectroscopy from discovery until about 1 month past r-band maximum brightness. The spectra of both objects show iron-group elements (Co II, Ni II, Fe II, Fe III, and high-velocity features [HVFs] of Ti II), intermediate-mass elements (Si II, Si III, and S II), and separate normal velocity features (~12000 km/s) and HVFs (~24000 km/s) of Ca II. Persistent absorption from Fe III and Si III, along with the colour evolution, imply high blackbody temperatures for SNe 2013bh and 2000cx (~12000 K). Both objects lack narrow Na I D absorption and exploded in the outskirts of their hosts, indicating that the SN environments were relatively free of interstellar or circumstellar material and may imply that the progenitors came from a relatively old and low-metallicity stellar population. Models of SN 2000cx, seemingly applicable to SN 2013bh, imply the production of up to ~1 M_Sun of Ni-56 and (4.3-5.5)e-3 M_Sun of fast-moving Ca ejecta.

preprint2013arXiv

The PTF Orion Project: a Possible Planet Transiting a T-Tauri Star

We report observations of a possible young transiting planet orbiting a previously known weak-lined T-Tauri star in the 7-10 Myr old Orion-OB1a/25-Ori region. The candidate was found as part of the Palomar Transient Factory (PTF) Orion project. It has a photometric transit period of 0.448413 +- 0.000040 days, and appears in both 2009 and 2010 PTF data. Follow-up low-precision radial velocity (RV) observations and adaptive optics imaging suggest that the star is not an eclipsing binary, and that it is unlikely that a background source is blended with the target and mimicking the observed transit. RV observations with the Hobby-Eberly and Keck telescopes yield an RV that has the same period as the photometric event, but is offset in phase from the transit center by approximately -0.22 periods. The amplitude (half range) of the RV variations is 2.4 km/s and is comparable with the expected RV amplitude that stellar spots could induce. The RV curve is likely dominated by stellar spot modulation and provides an upper limit to the projected companion mass of M_p sin i_orb < 4.8 +- 1.2 M_Jup; when combined with the orbital inclination, i orb, of the candidate planet from modeling of the transit light curve, we find an upper limit on the mass of the planetary candidate of M_p < 5.5 +- 1.4 M_Jup. This limit implies that the planet is orbiting close to, if not inside, its Roche limiting orbital radius, so that it may be undergoing active mass loss and evaporation.

preprint2013arXiv

The Swift/BAT Hard X-ray Transient Monitor

The Swift/Burst Alert Telescope (BAT) hard X-ray transient monitor provides near real-time coverage of the X-ray sky in the energy range 15-50 keV. The BAT observes 88% of the sky each day with a detection sensitivity of 5.3 mCrab for a full-day observation and a time resolution as fine as 64 seconds. The three main purposes of the monitor are (1) the discovery of new transient X-ray sources, (2) the detection of outbursts or other changes in the flux of known X-ray sources, and (3) the generation of light curves of more than 900 sources spanning over eight years. The primary interface for the BAT transient monitor is a public web page. Between 2005 February 12 and 2013 April 30, 245 sources have been detected in the monitor, 146 of them persistent and 99 detected only in outburst. Among these sources, 17 were previously unknown and were discovered in the transient monitor. In this paper, we discuss the methodology and the data processing and filtering for the BAT transient monitor and review its sensitivity and exposure. We provide a summary of the source detections and classify them according to the variability of their light curves. Finally, we review all new BAT monitor discoveries; for the new sources that are previously unpublished, we present basic data analysis and interpretations.

preprint2013arXiv

Type Ia Supernovae Strongly Interacting with Their Circumstellar Medium

Owing to their utility for measurements of cosmic acceleration, Type Ia supernovae (SNe) are perhaps the best-studied class of SNe, yet the progenitor systems of these explosions largely remain a mystery. A rare subclass of SNe Ia show evidence of strong interaction with their circumstellar medium (CSM), and in particular, a hydrogen-rich CSM; we refer to them as SNe Ia-CSM. In the first systematic search for such systems, we have identified 16 SNe Ia-CSM, and here we present new spectra of 13 of them. Six SNe Ia-CSM have been well-studied previously, three were previously known but are analyzed in-depth for the first time here, and seven are new discoveries from the Palomar Transient Factory. The spectra of all SNe Ia-CSM are dominated by Hα emission (with widths of ~2000 km/s) and exhibit large Hα/Hβ intensity ratios (perhaps due to collisional excitation of hydrogen via the SN ejecta overtaking slower-moving CSM shells); moreover, they have an almost complete lack of He I emission. They also show possible evidence of dust formation through a decrease in the red wing of Hα 75-100 d past maximum brightness, and nearly all SNe Ia-CSM exhibit strong Na I D absorption from the host galaxy. The absolute magnitudes (uncorrected for host-galaxy extinction) of SNe Ia-CSM are found to be -21.3 <= M_R <= -19 mag, and they also seem to show ultraviolet emission at early times and strong infrared emission at late times (but no detected radio or X-ray emission). Finally, the host galaxies of SNe Ia-CSM are all late-type spirals similar to the Milky Way, or dwarf irregulars like the Large Magellanic Cloud, which implies that these objects come from a relatively young stellar population. This work represents the most detailed analysis of the SN Ia-CSM class to date.

preprint2013arXiv

Type IIb Supernova SN 2011dh: Spectra and Photometry from the Ultraviolet to the Near-Infrared

We report spectroscopic and photometric observations of the Type IIb SN 2011dh obtained between 4 and 34 days after the estimated date of explosion (May 31.5 UT). The data cover a wide wavelength range from 2,000 Angstroms in the UV to 2.4 microns in the NIR. Optical spectra provide line profiles and velocity measurements of HI, HeI, CaII and FeII that trace the composition and kinematics of the SN. NIR spectra show that helium is present in the atmosphere as early as 11 days after the explosion. A UV spectrum obtained with the STIS reveals that the UV flux for SN 2011dh is low compared to other SN IIb. The HI and HeI velocities in SN 2011dh are separated by about 4,000 km/s at all phases. We estimate that the H-shell of SN 2011dh is about 8 times less massive than the shell of SN 1993J and about 3 times more massive than the shell of SN 2008ax. Light curves (LC) for twelve passbands are presented. The maximum bolometric luminosity of $1.8 \pm 0.2 \times 10^{42}$ erg s$^{-1}$ occurred about 22 days after the explosion. NIR emission provides more than 30% of the total bolometric flux at the beginning of our observations and increases to nearly 50% of the total by day 34. The UV produces 16% of the total flux on day 4, 5% on day 9 and 1% on day 34. We compare the bolometric light curves of SN 2011dh, SN 2008ax and SN 1993J. The LC are very different for the first twelve days after the explosions but all three SN IIb display similar peak luminosities, times of peak, decline rates and colors after maximum. This suggests that the progenitors of these SN IIb may have had similar compositions and masses but they exploded inside hydrogen shells that that have a wide range of masses. The detailed observations presented here will help evaluate theoretical models for this supernova and lead to a better understanding of SN IIb.

preprint2012arXiv

A Bayesian Approach to Calibrating Period-Luminosity Relations of RR Lyrae Stars in the Mid-Infrared

A Bayesian approach to calibrating period-luminosity (PL) relations has substantial benefits over generic least-squares fits. In particular, the Bayesian approach takes into account the full prior distribution of the model parameters, such as the a priori distances, and refits these parameters as part of the process of settling on the most highly-constrained final fit. Additionally, the Bayesian approach can naturally ingest data from multiple wavebands and simultaneously fit the parameters of PL relations for each waveband in a procedure that constrains the parameter posterior distributions so as to minimize the scatter of the final fits appropriately in all wavebands. Here we describe the generalized approach to Bayesian model fitting and then specialize to a detailed description of applying Bayesian linear model fitting to the mid-infrared PL relations of RR Lyrae variable stars. For this example application we quantify the improvement afforded by using a Bayesian model fit. We also compare distances previously predicted in our example application to recently published parallax distances measured with the Hubble Space Telescope and find their agreement to be a vindication of our methodology. Our intent with this article is to spread awareness of the benefits and applicability of this Bayesian approach and encourage future PL relation investigations to consider employing this powerful analysis method.

preprint2012arXiv

An Early & Comprehensive Millimeter and Centimeter Wave and X-ray Study of Supernova 2011dh: A Non-Equipartition Blastwave Expanding into A Massive Stellar Wind

Only a handful of supernovae (SNe) have been studied in multi-wavelength from radio to X-rays, starting a few days after explosion. The early detection and classification of the nearby type IIb SN2011dh/PTF11eon in M51 provides a unique opportunity to conduct such observations. We present detailed data obtained at the youngest phase ever of a core-collapse supernova (days 3 to 12 after explosion) in the radio, millimeter and X-rays; when combined with optical data, this allows us to explore the early evolution of the SN blast wave and its surroundings. Our analysis shows that the expanding supernova shockwave does not exhibit equipartition (e_e/e_B ~ 1000), and is expanding into circumstellar material that is consistent with a density profile falling like R^-2. Within modeling uncertainties we find an average velocity of the fast parts of the ejecta of 15,000 +/- 1800 km/s, contrary to previous analysis. This velocity places SN 2011dh in an intermediate blast-wave regime between the previously defined compact and extended SN IIb subtypes. Our results highlight the importance of early (~ 1 day) high-frequency observations of future events. Moreover, we show the importance of combined radio/X-ray observations for determining the microphysics ratio e_e/e_B.

preprint2012arXiv

An Empirical Relation between Sodium Absorption and Dust Extinction

Dust extinction and reddening are ubiquitous in astronomical observations and are often a major source of systematic uncertainty. We present here a study of the correlation between extinction in the Milky Way and the equivalent width of the NaI D absorption doublet. Our sample includes more than 100 high resolution spectra from the KECK telescopes and nearly a million low resolution spectra from the Sloan Digital Sky Survey (SDSS). We measure the correlation to unprecedented precision, constrain its shape, and derive an empirical relation between these quantities with a dispersion of order 0.15 magnitude in E(B-V). From the shape of the curve of growth we further show that a typical sight line through the Galaxy, as seen within the SDSS footprint, crosses about three dust clouds. We provide a brief guide on how to best estimate extinction to extragalactic sources such as supernovae, using the NaI D absorption feature, under a variety of circumstances.

preprint2012arXiv

Berkeley Supernova Ia Program I: Observations, Data Reduction, and Spectroscopic Sample of 582 Low-Redshift Type Ia Supernovae

In this first paper in a series we present 1298 low-redshift (z\leq0.2) optical spectra of 582 Type Ia supernovae (SNe Ia) observed from 1989 through 2008 as part of the Berkeley SN Ia Program (BSNIP). 584 spectra of 199 SNe Ia have well-calibrated light curves with measured distance moduli, and many of the spectra have been corrected for host-galaxy contamination. Most of the data were obtained using the Kast double spectrograph mounted on the Shane 3 m telescope at Lick Observatory and have a typical wavelength range of 3300-10,400 Ang., roughly twice as wide as spectra from most previously published datasets. We present our observing and reduction procedures, and we describe the resulting SN Database (SNDB), which will be an online, public, searchable database containing all of our fully reduced spectra and companion photometry. In addition, we discuss our spectral classification scheme (using the SuperNova IDentification code, SNID; Blondin & Tonry 2007), utilising our newly constructed set of SNID spectral templates. These templates allow us to accurately classify our entire dataset, and by doing so we are able to reclassify a handful of objects as bona fide SNe Ia and a few other objects as members of some of the peculiar SN Ia subtypes. In fact, our dataset includes spectra of nearly 90 spectroscopically peculiar SNe Ia. We also present spectroscopic host-galaxy redshifts of some SNe Ia where these values were previously unknown. [Abridged]

preprint2012arXiv

Classical Novae in Andromeda: Light Curves from the Palomar Transient Factory and GALEX

We present optical light curves of twenty-nine novae in M31 during the 2009 and 2010 observing seasons of the Palomar Transient Factory (PTF). The dynamic and rapid cadences in PTF monitoring of M31, from one day to even ten minutes, provide excellent temporal coverage of nova light curves, enabling us to record the photometric evolution of M31 novae in unprecedented detail. We also detect eight of these novae in the near ultraviolet (UV) band with the Galaxy Evolution Explorer (GALEX) satellite. Novae M31N2009-10b and 2010-11a show prominent UV emission peaking a few days prior to their optical maxima, possibly implying aspherical outbursts. Additionally, our blue-shifted spectrum of the recent outburst of PT And (M31N2010-12a) indicates that it is a recurrent nova in M31 and not a dwarf nova in the Milky Way as was previously assumed. Finally, we systematically searched for novae in all confirmed globular clusters of M31 and found only M31N 2010-10f associated with Bol 126. The specific nova rate in the M31 globular cluster system is thus about one per year which is not enhanced relative to the rate outside the globular cluster system.

preprint2012arXiv

Construction of a Calibrated Probabilistic Classification Catalog: Application to 50k Variable Sources in the All-Sky Automated Survey

With growing data volumes from synoptic surveys, astronomers must become more abstracted from the discovery and introspection processes. Given the scarcity of follow-up resources, there is a particularly sharp onus on the frameworks that replace these human roles to provide accurate and well-calibrated probabilistic classification catalogs. Such catalogs inform the subsequent follow-up, allowing consumers to optimize the selection of specific sources for further study and permitting rigorous treatment of purities and efficiencies for population studies. Here, we describe a process to produce a probabilistic classification catalog of variability with machine learning from a multi-epoch photometric survey. In addition to producing accurate classifications, we show how to estimate calibrated class probabilities, and motivate the importance of probability calibration. We also introduce a methodology for feature-based anomaly detection, which allows discovery of objects in the survey that do not fit within the predefined class taxonomy. Finally, we apply these methods to sources observed by the All Sky Automated Survey (ASAS), and unveil the Machine-learned ASAS Classification Catalog (MACC), which is a 28-class probabilistic classification catalog of 50,124 ASAS sources. We estimate that MACC achieves a sub-20% classification error rate, and demonstrate that the class posterior probabilities are reasonably calibrated. MACC classifications compare favorably to the classifications of several previous domain-specific ASAS papers and to the ASAS Catalog of Variable Stars, which had classified only 24% of those sources into one of 12 science classes. The MACC is publicly available at http://www.bigmacc.info.

preprint2012arXiv

Finding the sources of missing heritability in a yeast cross

For many traits, including susceptibility to common diseases in humans, causal loci uncovered by genetic mapping studies explain only a minority of the heritable contribution to trait variation. Multiple explanations for this "missing heritability" have been proposed. Here we use a large cross between two yeast strains to accurately estimate different sources of heritable variation for 46 quantitative traits and to detect underlying loci with high statistical power. We find that the detected loci explain nearly the entire additive contribution to heritable variation for the traits studied. We also show that the contribution to heritability of gene-gene interactions varies among traits, from near zero to 50%. Detected two-locus interactions explain only a minority of this contribution. These results substantially advance our understanding of the missing heritability problem and have important implications for future studies of complex and quantitative traits.

preprint2012arXiv

Five New Outbursting AM CVn Systems Discovered by the Palomar Transient Factory

We present five new outbursting AM CVn systems and one candidate discovered as part of an ongoing search for such systems using the Palomar Transient Factory (PTF). This is the first large-area, systematic search for AM CVn systems using only large-amplitude photometric variability to select candidates. Three of the confirmed systems and the candidate system were discovered as part of the PTF transient search. Two systems were found as part of a search for outbursts through the PTF photometric database. We discuss the observed characteristics of each of these systems, including the orbital periods of two systems. We also consider the position of these systems, selected in a colour-independent survey, in colour-colour space and compare to systems selected solely by their colours. We find that the colours of our newly discovered systems do not differ significantly from those of previously known systems, but significant errors preclude a definitive answer.

preprint2012arXiv

Highly Variable Extinction and Accretion in the Jet-driving Class I Type Young Star PTF 10nvg (V2492 Cyg, IRAS 20496+4354)

We report extensive new photometry and spectroscopy of the highly variable young stellar object PTF 10nvg including optical and near-infrared time series data as well as mid-infrared and millimeter data. Following the previously reported 2010 rise, during 2011 and 2012 the source underwent additional episodes of brightening and dimming events including prolonged faint states. The observed high-amplitude variations are largely consistent with extinction changes having a 220 day quasi-periodic signal. Spectral evolution includes not only changes in the spectral slope but correlated variation in the prominence of TiO/VO/CO bands and atomic line emission, as well as anticorrelated variation in forbidden line emission which, along with H_2, dominates optical and infrared spectra at faint epochs. Neutral and singly-ionized atomic species are likely formed in an accretion flow and/or impact while the origin of zero-velocity atomic LiI 6707 in emission is unknown. Forbidden lines, including several rare species, exhibit blueshifted emission profiles and likely arise from an outflow/jet. Several of these lines are also seen spatially offset from the continuum source position, presumably in a shocked region of an extended jet. CARMA maps resolve on larger scales a spatially extended outflow in mm-wavelength CO. We attribute the observed photometric and spectroscopic behavior in terms of occultation of the central star as well as the bright inner disk and the accretion/outflow zones that renders shocked gas in the inner part of the jet amenable to observation at the faint epochs. We discuss PTF 10nvg as a source exhibiting both accretion-driven (perhaps analogous to V1647 Ori) and extinction-driven (perhaps analogous to UX Ori or GM Cep) high-amplitude variability phenomena.

preprint2012arXiv

Late Time Radio Emission from X-ray Selected Tidal Disruption Events

We present new observations with the Karl G. Jansky Very Large Array of seven X-ray-selected tidal disruption events (TDEs). The radio observations were carried out between 9 and 22 years after the initial X-ray discovery, and, thus, probe the late-time formation of relativistic jets and jet interactions with the interstellar medium in these systems. We detect a compact radio source in the nucleus of the galaxy IC 3599 and a compact radio source that is a possible counterpart to RX J1420.4+5334. We find no radio counterparts for five other sources with flux density upper limits between 51 and 200 microJy (3 sigma). If the detections truly represent late radio emission associated with a TDE, then our results suggest that a fraction >~ 10% of X-ray-detected TDEs are accompanied by relativistic jets. We explore several models for producing late radio emission, including interaction of the jet with gas in the circumnuclear environment (blast wave model), and emission from the core of the jet itself. Upper limits on the radio flux density from archival observations suggest that the jet formation may have been delayed for years after the TDE, possibly triggered by the accretion rate dropping below a critical threshold of ~ 10^{-2} -- 10^{-3} Mdot_Edd. The non-detections are also consistent with this scenario; deeper radio observations can determine whether relativistic jets are present in these systems. The emission from RX J1420.4+5334 is also consistent with the predictions of the blast wave model, however the radio emission from IC 3599 is substantially underluminous, and its spectral slope is too flat, relative to the blast wave model expectations. Future radio monitoring of IC 3599 and RX J1420.4+5334 will help to better constrain the nature of the jets in these systems.

preprint2012arXiv

Millions of Multiples: Detecting and Characterizing Close-Separation Binary Systems in Synoptic Sky Surveys

The direct detection of binary systems in wide-field surveys is limited by the size of the stars' point-spread-functions (PSFs). A search for elongated objects can find closer companions, but is limited by the precision to which the PSF shape can be calibrated for individual stars. We have developed the BinaryFinder algorithm to search for close binaries by using precision measurements of PSF ellipticity across wide-field survey images. We show that the algorithm is capable of reliably detecting binary systems down to approximately 1/5 of the seeing limit, and can directly measure the systems' position angles, separations and contrast ratios. To verify the algorithm's performance we evaluated 100,000 objects in Palomar Transient Factory (PTF) wide-field-survey data for signs of binarity, and then used the Robo-AO robotic laser adaptive optics system to verify the parameters of 44 high-confidence targets. We show that BinaryFinder correctly predicts the presence of close companions with a <5% false-positive rate, measures the detected binaries' position angles within 2 degrees and separations within 25%, and weakly constrains their contrast ratios. When applied to the full PTF dataset, we estimate that BinaryFinder will discover and characterize ~450,000 physically-associated binary systems with separations <2 arcseconds and magnitudes brighter than R=18. New wide-field synoptic surveys with high sensitivity and sub-arcsecond angular resolution, such as LSST, will allow BinaryFinder to reliably detect millions of very faint binary systems with separations as small as 0.1 arcseconds.

preprint2012arXiv

Optimizing Automated Classification of Periodic Variable Stars in New Synoptic Surveys

Efficient and automated classification of periodic variable stars is becoming increasingly important as the scale of astronomical surveys grows. Several recent papers have used methods from machine learning and statistics to construct classifiers on databases of labeled, multi--epoch sources with the intention of using these classifiers to automatically infer the classes of unlabeled sources from new surveys. However, the same source observed with two different synoptic surveys will generally yield different derived metrics (features) from the light curve. Since such features are used in classifiers, this survey-dependent mismatch in feature space will typically lead to degraded classifier performance. In this paper we show how and why feature distributions change using OGLE and \textit{Hipparcos} light curves. To overcome survey systematics, we apply a method, \textit{noisification}, which attempts to empirically match distributions of features between the labeled sources used to construct the classifier and the unlabeled sources we wish to classify. Results from simulated and real--world light curves show that noisification can significantly improve classifier performance. In a three--class problem using light curves from \textit{Hipparcos} and OGLE, noisification reduces the classifier error rate from 27.0% to 7.0%. We recommend that noisification be used for upcoming surveys such as Gaia and LSST and describe some of the promises and challenges of applying noisification to these surveys.

preprint2012arXiv

Swift J2058.4+0516: Discovery of a Possible Second Relativistic Tidal Disruption Flare?

We report the discovery by the Swift hard X-ray monitor of the transient source Swift J2058.4+0516 (Sw J2058+05). Our multi-wavelength follow-up campaign uncovered a long-lived (duration >~ months), luminous X-ray (L_X,iso ~ 3 x 10^47 erg s^-1) and radio (nu L_nu,iso ~ 10^42 erg s^-1) counterpart. The associated optical emission, however, from which we measure a redshift of 1.1853, is relatively faint, and this is not due to a large amount of dust extinction in the host galaxy. Based on numerous similarities with the recently discovered GRB 110328A / Swift J164449.3+573451 (Sw J1644+57), we suggest that Sw J2058+05 may be the second member of a new class of relativistic outbursts resulting from the tidal disruption of a star by a supermassive black hole. If so, the relative rarity of these sources (compared with the expected rate of tidal disruptions) implies that either these outflows are extremely narrowly collimated (theta < 1 degree), or only a small fraction of tidal disruptions generate relativistic ejecta. Analogous to the case of long-duration gamma-ray bursts and core-collapse supernovae, we speculate that rapid spin of the black hole may be a necessary condition to generate the relativistic component. Alternatively, if powered by gas accretion (i.e., an active galactic nucleus [AGN]), Sw J2058+05 would seem to represent a new mode of variability in these sources, as the observed properties appear largely inconsistent with known classes of AGNs capable of generating relativistic jets (blazars, narrow-line Seyfert 1 galaxies).

preprint2012arXiv

Using Machine Learning for Discovery in Synoptic Survey Imaging

Modern time-domain surveys continuously monitor large swaths of the sky to look for astronomical variability. Astrophysical discovery in such data sets is complicated by the fact that detections of real transient and variable sources are highly outnumbered by bogus detections caused by imperfect subtractions, atmospheric effects and detector artefacts. In this work we present a machine learning (ML) framework for discovery of variability in time-domain imaging surveys. Our ML methods provide probabilistic statements, in near real time, about the degree to which each newly observed source is astrophysically relevant source of variable brightness. We provide details about each of the analysis steps involved, including compilation of the training and testing sets, construction of descriptive image-based and contextual features, and optimization of the feature subset and model tuning parameters. Using a validation set of nearly 30,000 objects from the Palomar Transient Factory, we demonstrate a missed detection rate of at most 7.7% at our chosen false-positive rate of 1% for an optimized ML classifier of 23 features, selected to avoid feature correlation and over-fitting from an initial library of 42 attributes. Importantly, we show that our classification methodology is insensitive to mis-labelled training data up to a contamination of nearly 10%, making it easier to compile sufficient training sets for accurate performance in future surveys. This ML framework, if so adopted, should enable the maximization of scientific gain from future synoptic survey and enable fast follow-up decisions on the vast amounts of streaming data produced by such experiments.

preprint2011arXiv

A Compact Degenerate Primary-Star Progenitor of SN 2011fe

While a white dwarf is, from a theoretical perspective, the most plausible primary star in Type Ia supernova (SN Ia), many other candidates have not been formally ruled out. Shock energy deposited in the envelope of any exploding primary contributes to the early SN brightness and, since this radiation energy is degraded by expansion after the explosion, the diffusive luminosity depends on the initial primary radius. We present a new non-detection limit of the nearby SN Ia 2011fe, obtained what appears to be just 4 hours after explosion, allowing us to directly constrain the initial primary radius, R_p. Coupled with the non-detection of a quiescent X-ray counterpart and the inferred synthesized Ni mass, we show that R_p <~ 0.02 R_sun (a factor of 5 smaller than previously inferred), that the average density of the primary must be rho_p > 10,000 gm cm^{-3}, and that the effective temperature must be less than a few x 10^5 K. This rules out hydrogen burning main sequence stars and giants. Constructing the helium-burning main sequence and carbon-burning main sequence, we find such objects are also excluded. By process of elimination, we find that only degeneracy-supported compact objects---WDs and neutron stars---are viable as the primary star of SN 2011fe. With few caveats, we also restrict the companion (secondary) star radius to R_ c <~ 0.1 R_sun, excluding Roche-Lobe overflowing red giant and main-sequence companions to high significance.

preprint2011arXiv

A relativistic jetted outburst from a massive black hole fed by a tidally disrupted star

While gas accretion onto some massive black holes (MBHs) at the centers of galaxies actively powers luminous emission, the vast majority of MBHs are considered dormant. Occasionally, a star passing too near a MBH is torn apart by gravitational forces, leading to a bright panchromatic tidal disruption flare (TDF). While the high-energy transient Swift J164449.3+573451 ("Sw 1644+57") initially displayed none of the theoretically anticipated (nor previously observed) TDF characteristics, we show that the observations (Levan et al. 2011) suggest a sudden accretion event onto a central MBH of mass ~10^6-10^7 solar masses. We find evidence for a mildly relativistic outflow, jet collimation, and a spectrum characterized by synchrotron and inverse Compton processes; this leads to a natural analogy of Sw 1644+57 with a smaller-scale blazar. The phenomenologically novel Sw 1644+57 thus connects the study of TDFs and active galaxies, opening a new vista on disk-jet interactions in BHs and magnetic field generation and transport in accretion systems.

preprint2011arXiv

Active Learning to Overcome Sample Selection Bias: Application to Photometric Variable Star Classification

Despite the great promise of machine-learning algorithms to classify and predict astrophysical parameters for the vast numbers of astrophysical sources and transients observed in large-scale surveys, the peculiarities of the training data often manifest as strongly biased predictions on the data of interest. Typically, training sets are derived from historical surveys of brighter, more nearby objects than those from more extensive, deeper surveys (testing data). This sample selection bias can cause catastrophic errors in predictions on the testing data because a) standard assumptions for machine-learned model selection procedures break down and b) dense regions of testing space might be completely devoid of training data. We explore possible remedies to sample selection bias, including importance weighting (IW), co-training (CT), and active learning (AL). We argue that AL---where the data whose inclusion in the training set would most improve predictions on the testing set are queried for manual follow-up---is an effective approach and is appropriate for many astronomical applications. For a variable star classification problem on a well-studied set of stars from Hipparcos and OGLE, AL is the optimal method in terms of error rate on the testing data, beating the off-the-shelf classifier by 3.4% and the other proposed methods by at least 3.0%. To aid with manual labeling of variable stars, we developed a web interface which allows for easy light curve visualization and querying of external databases. Finally, we apply active learning to classify variable stars in the ASAS survey, finding dramatic improvement in our agreement with the ACVS catalog, from 65.5% to 79.5%, and a significant increase in the classifier's average confidence for the testing set, from 14.6% to 42.9%, after a few AL iterations.

preprint2011arXiv

Calcium-rich gap transients in the remote outskirts of galaxies

From the first two seasons of the Palomar Transient Factory, we identify three peculiar transients (PTF09dav, PTF10iuv, PTF11bij) with five distinguishing characteristics: peak luminosity in the gap between novae and supernovae (M_R = 15.5 to -16.5), rapid photometric evolution (rise-time ~12--15 days), large photospheric velocities (~6000 to 11000 km/s), early spectroscopic evolution into nebular phase (~1 to 3 months) and peculiar nebular spectra dominated by Calcium. We also culled the extensive decade-long Lick Observatory Supernova Search database and identified an additional member of this group, SN 2007ke. Our choice of photometric and spectroscopic properties was motivated by SN 2005E (Perets et al. 2010). To our surprise, as in the case of SN 2005E, all four members of this group are also clearly offset from the bulk of their host galaxy. Given the well-sampled early and late-time light curves, we derive ejecta masses in the range of 0.4--0.7 Msun. Spectroscopically, we find that there may be a diversity in the photospheric phase, but the commonality is in the unusual nebular spectra. Our extensive follow-up observations rule out standard thermonuclear and standard core-collapse explosions for this class of "Calcium-rich gap" transients. If the progenitor is a white dwarf, we are likely seeing a detonation of the white dwarf core and perhaps, even shockfront interaction with a previously ejected nova shell. In the less likely scenario of a massive star progenitor, a very non-standard channel specific to a low-metallicity environment needs to be invoked (e.g., ejecta fallback leading to black hole formation). Detection (or lack thereof) of a faint underlying host (dwarf galaxy, cluster) will provide a crucial and decisive diagnostic to choose between these alternatives.

preprint2011arXiv

Constraints on the Progenitor System of the Type Ia Supernova SN 2011fe/PTF11kly

Type Ia supernovae (SNe) serve as a fundamental pillar of modern cosmology, owing to their large luminosity and a well-defined relationship between light-curve shape and peak brightness. The precision distance measurements enabled by SNe Ia first revealed the accelerating expansion of the universe, now widely believed (though hardly understood) to require the presence of a mysterious "dark" energy. General consensus holds that Type Ia SNe result from thermonuclear explosions of a white dwarf (WD) in a binary system; however, little is known of the precise nature of the companion star and the physical properties of the progenitor system. Here we make use of extensive historical imaging obtained at the location of SN 2011fe/PTF11kly, the closest SN Ia discovered in the digital imaging era, to constrain the visible-light luminosity of the progenitor to be 10-100 times fainter than previous limits on other SN Ia progenitors. This directly rules out luminous red giants and the vast majority of helium stars as the mass-donating companion to the exploding white dwarf. Any evolved red companion must have been born with mass less than 3.5 times the mass of the Sun. These observations favour a scenario where the exploding WD of SN 2011fe/PTF11kly, accreted matter either from another WD, or by Roche-lobe overflow from a subgiant or main-sequence companion star.

preprint2011arXiv

Data Mining and Machine-Learning in Time-Domain Discovery & Classification

The changing heavens have played a central role in the scientific effort of astronomers for centuries. Galileo's synoptic observations of the moons of Jupiter and the phases of Venus starting in 1610, provided strong refutation of Ptolemaic cosmology. In more modern times, the discovery of a relationship between period and luminosity in some pulsational variable stars led to the inference of the size of the Milky Way, the distance scale to the nearest galaxies, and the expansion of the Universe. Distant explosions of supernovae were used to uncover the existence of dark energy and provide a precise numerical account of dark matter. Indeed, time-domain observations of transient events and variable stars, as a technique, influences a broad diversity of pursuits in the entire astronomy endeavor. While, at a fundamental level, the nature of the scientific pursuit remains unchanged, the advent of astronomy as a data-driven discipline presents fundamental challenges to the way in which the scientific process must now be conducted. Digital images (and data cubes) are not only getting larger, there are more of them. On logistical grounds, this taxes storage and transport systems. But it also implies that the intimate connection that astronomers have always enjoyed with their data---from collection to processing to analysis to inference---necessarily must evolve. The pathway to scientific inference is now influenced (if not driven by) modern automation processes, computing, data-mining and machine learning. The emerging reliance on computation and machine learning is a general one, but the time-domain aspect of the data and the objects of interest presents some unique challenges, which we describe and explore in this chapter.

preprint2011arXiv

Evidence for an FU Orionis-like Outburst from a Classical T Tauri Star

We present pre- and post-outburst observations of the new FU Orionis-like young stellar object PTF 10qpf (also known as LkHa 188-G4 and HBC 722). Prior to this outburst, LkHa 188-G4 was classified as a classical T Tauri star on the basis of its optical emission-line spectrum superposed on a K8-type photosphere, and its photometric variability. The mid-infrared spectral index of LkHa 188-G4 indicates a Class II-type object. LkHa 188-G4 exhibited a steady rise by ~1 mag over ~11 months starting in Aug. 2009, before a subsequent more abrupt rise of > 3 mag on a time scale of ~2 months. Observations taken during the eruption exhibit the defining characteristics of FU Orionis variables: (i) an increase in brightness by > 4 mag, (ii) a bright optical/near-infrared reflection nebula appeared, (iii) optical spectra are consistent with a G supergiant and dominated by absorption lines, the only exception being Halpha which is characterized by a P Cygni profile, (iv) near-infrared spectra resemble those of late K--M giants/supergiants with enhanced absorption seen in the molecular bands of CO and H_2O, and (v) outflow signatures in H and He are seen in the form of blueshifted absorption profiles. LkHa 188-G4 is the first member of the FU Orionis-like class with a well-sampled optical to mid-infrared spectral energy distribution in the pre-outburst phase. The association of the PTF 10qpf outburst with the previously identified classical T Tauri star LkHa 188-G4 (HBC 722) provides strong evidence that FU Orionis-like eruptions represent periods of enhanced disk accretion and outflow, likely triggered by instabilities in the disk. The early identification of PTF 10qpf as an FU Orionis-like variable will enable detailed photometric and spectroscopic observations during its post-outburst evolution for comparison with other known outbursting objects.

preprint2011arXiv

Mid-infrared Period-Luminosity Relations of RR Lyrae Stars Derived from the WISE Preliminary Data Release

Interstellar dust presents a significant challenge to extending parallax-determined distances of optically observed pulsational variables to larger volumes. Distance ladder work at mid-infrared wavebands, where dust effects are negligible and metallicity correlations are minimized, have been largely focused on few-epoch Cepheid studies. Here we present the first determination of mid-infrared period-luminosity (PL) relations of RR Lyrae stars from phase-resolved imaging using the preliminary data release of the Wide-Field Infrared Survey Explorer (WISE). We present a novel statistical framework to predict posterior distances of 76 well-observed RR Lyrae that uses the optically constructed prior distance moduli while simultaneously imposing a power-law PL relation to WISE-determined mean magnitudes. We find that the absolute magnitude in the bluest WISE filter is M_W1 = (-0.421+-0.014) - (1.681+-0.147)*log(P/0.50118 day), with no evidence for a correlation with metallicity. Combining the results from the three bluest WISE filters, we find that a typical star in our sample has a distance measurement uncertainty of 0.97% (statistical) plus 1.17% (systematic). We do not fundamentalize the periods of RRc stars to improve their fit to the relations. Taking the Hipparcos-derived mean V-band magnitudes, we use the distance posteriors to determine a new optical metallicity-luminosity relation which we present in Section 5. The results of this analysis will soon be tested by HST parallax measurements and, eventually, with the Gaia astrometric mission.

preprint2011arXiv

On Machine-Learned Classification of Variable Stars with Sparse and Noisy Time-Series Data

With the coming data deluge from synoptic surveys, there is a growing need for frameworks that can quickly and automatically produce calibrated classification probabilities for newly-observed variables based on a small number of time-series measurements. In this paper, we introduce a methodology for variable-star classification, drawing from modern machine-learning techniques. We describe how to homogenize the information gleaned from light curves by selection and computation of real-numbered metrics ("feature"), detail methods to robustly estimate periodic light-curve features, introduce tree-ensemble methods for accurate variable star classification, and show how to rigorously evaluate the classification results using cross validation. On a 25-class data set of 1542 well-studied variable stars, we achieve a 22.8% overall classification error using the random forest classifier; this represents a 24% improvement over the best previous classifier on these data. This methodology is effective for identifying samples of specific science classes: for pulsational variables used in Milky Way tomography we obtain a discovery efficiency of 98.2% and for eclipsing systems we find an efficiency of 99.1%, both at 95% purity. We show that the random forest (RF) classifier is superior to other machine-learned methods in terms of accuracy, speed, and relative immunity to features with no useful class information; the RF classifier can also be used to estimate the importance of each feature in classification. Additionally, we present the first astronomical use of hierarchical classification methods to incorporate a known class taxonomy in the classifier, which further reduces the catastrophic error rate to 7.8%. Excluding low-amplitude sources, our overall error rate improves to 14%, with a catastrophic error rate of 3.5%.

preprint2011arXiv

Optimal Time-Series Selection of Quasars

We present a novel method for the optimal selection of quasars using time-series observations in a single photometric bandpass. Utilizing the damped random walk model of Kelly et al. (2009), we parameterize the ensemble quasar structure function in Sloan Stripe 82 as a function of observed brightness. The ensemble model fit can then be evaluated rigorously for and calibrated with individual light curves with no parameter fitting. This yields a classification in two statistics --- one describing the fit confidence and one describing the probability of a false alarm --- which can be tuned, a priori, to achieve high quasar detection fractions (99% completeness with default cuts), given an acceptable rate of false alarms. We establish the typical rate of false alarms due to known variable stars as <3% (high purity). Applying the classification, we increase the sample of potential quasars relative to those known in Stripe 82 by as much as 29%, and by nearly a factor of two in the redshift range 2.5<z<3, where selection by color is extremeley inefficient. This represents 1875 new quasars in a 290 deg^2 field. The observed rates of both quasars and stars agree well with the model predictions, with >99% of quasars exhibiting the expected variability profile. We discus the utility of the method at high-redshift and in the regime of noisy and sparse data. Our time series selection complements well independent selection based on quasar colors and has strong potential for identifying high redshift quasars for BAO and other cosmology studies in the LSST era.

preprint2011arXiv

PTF1 J071912.13+485834.0: An outbursting AM CVn system discovered by a synoptic survey

We present extensive photometric and spectroscopic observations of PTF1 J071912.13+485834.0, an outbursting AM CVn system discovered by the Palomar Transient Factory (PTF). AM CVn systems are stellar binaries with some of the smallest separations known and orbital periods ranging from 5 to 65 minutes. They are believed to be composed of a white dwarf accretor and a (semi)-degenerate He-rich donor and are considered to be the helium equivalents of Cataclysmic Variables. We have spectroscopically and photometrically identified an orbital period of 26.77 \pm 0.02 minutes for PTF1 J071912.13+485834.0 and found a super-outburst recurrence time of greater than 65 days along with the presence of "normal" outbursts - rarely seen in AM CVn systems but well known in super-outbursting Cataclysmic Variables. We present a long-term light curve over two super-cycles as well as high cadence photometry of both outburst and quiescent stages, both of which show clear variability. We also compare both the outburst and quiescent spectra of PTF1 J071912.13+485834.0 to other known AM CVn systems, and use the quiescent phase-resolved spectroscopy to determine the origin of the photometric variability. Finally, we draw parallels between the different subclasses of SU UMa-type Cataclysmic Variables and outbursting AM CVn systems. We conclude by predicting that the Palomar Transient Factory may more than double the number of outbursting AM CVn systems known, which would greatly increase our understanding of AM CVn systems.

preprint2011arXiv

PTF10fqs: A Luminous Red Nova in the Spiral Galaxy Messier 99

The Palomar Transient Factory (PTF) is systematically charting the optical transient and variable sky. A primary science driver of PTF is building a complete inventory of transients in the local Universe (distance less than 200 Mpc). Here, we report the discovery of PTF10fqs, a transient in the luminosity "gap" between novae and supernovae. Located on a spiral arm of Messier 99, PTF 10fqs has a peak luminosity of Mr = -12.3, red color (g-r = 1.0) and is slowly evolving (decayed by 1 mag in 68 days). It has a spectrum dominated by intermediate-width H (930 km/s) and narrow calcium emission lines. The explosion signature (the light curve and spectra) is overall similar to thatof M85OT2006-1, SN2008S, and NGC300OT. The origin of these events is shrouded in mystery and controversy (and in some cases, in dust). PTF10fqs shows some evidence of a broad feature (around 8600A) that may suggest very large velocities (10,000 km/s) in this explosion. Ongoing surveys can be expected to find a few such events per year. Sensitive spectroscopy, infrared monitoring and statistics (e.g. disk versus bulge) will eventually make it possible for astronomers to unravel the nature of these mysterious explosions.

preprint2011arXiv

PTF10iya: A short-lived, luminous flare from the nuclear region of a star-forming galaxy

We present the discovery and characterisation of PTF10iya, a short-lived (dt ~ 10 d, with an optical decay rate of ~ 0.3 mag per d), luminous (M_g ~ -21 mag) transient source found by the Palomar Transient Factory. The ultraviolet/optical spectral energy distribution is reasonably well fit by a blackbody with T ~ 1-2 x 10^4 K and peak bolometric luminosity L_BB ~ 1-5 x 10^44 erg per s (depending on the details of the extinction correction). A comparable amount of energy is radiated in the X-ray band that appears to result from a distinct physical process. The location of PTF10iya is consistent with the nucleus of a star-forming galaxy (z = 0.22405 +/- 0.00006) to within 350 mas (99.7 per cent confidence radius), or a projected distance of less than 1.2 kpc. At first glance, these properties appear reminiscent of the characteristic "big blue bump" seen in the near-ultraviolet spectra of many active galactic nuclei (AGNs). However, emission-line diagnostics of the host galaxy, along with a historical light curve extending back to 2007, show no evidence for AGN-like activity. We therefore consider whether the tidal disruption of a star by an otherwise quiescent supermassive black hole may account for our observations. Though with limited temporal information, PTF10iya appears broadly consistent with the predictions for the early "super-Eddington" phase of a solar-type star disrupted by a ~ 10^7 M_sun black hole. Regardless of the precise physical origin of the accreting material, the large luminosity and short duration suggest that otherwise quiescent galaxies can transition extremely rapidly to radiate near the Eddington limit; many such outbursts may have been missed by previous surveys lacking sufficient cadence.

preprint2011arXiv

PTF10ops - a subluminous, normal-width lightcurve Type Ia supernova in the middle of nowhere

PTF10ops is a Type Ia supernova (SN Ia), whose lightcurve and spectral properties place it outside the current SN Ia subtype classifications. Its spectra display the characteristic lines of subluminous SNe Ia, but it has a normal-width lightcurve with a long rise-time, typical of normal luminosity SNe Ia. The early-time optical spectra of PTF10ops were modelled using a spectral fitting code and found to have all the lines typically seen in subluminous SNe Ia, without the need to invoke more uncommon elements. The host galaxy environment of PTF10ops is also unusual with no galaxy detected at the position of the SN down to an absolute limiting magnitude of r \geq -12.0 mag, but a very massive galaxy is present at a separation of ~148 kpc and at the same redshift as suggested by the SN spectral features. The progenitor of PTF10ops is most likely a very old star, possibly in a low metallicity environment, which affects its explosion mechanism and observational characteristics. PTF10ops does not easily fit into any of the current models of either subluminous or normal SN Ia progenitor channels.

preprint2011arXiv

PTF11eon/SN2011dh: Discovery of a Type IIb Supernova From a Compact Progenitor in the Nearby Galaxy M51

On May 31, 2011 UT a supernova (SN) exploded in the nearby galaxy M51 (the Whirlpool Galaxy). We discovered this event using small telescopes equipped with CCD cameras, as well as by the Palomar Transient Factory (PTF) survey, and rapidly confirmed it to be a Type II supernova. Our early light curve and spectroscopy indicates that PTF11eon resulted from the explosion of a relatively compact progenitor star as evidenced by the rapid shock-breakout cooling seen in the light curve, the relatively low temperature in early-time spectra and the prompt appearance of low-ionization spectral features. The spectra of PTF11eon are dominated by H lines out to day 10 after explosion, but initial signs of He appear to be present. Assuming that He lines continue to develop in the near future, this SN is likely a member of the cIIb (compact IIb; Chevalier and Soderberg 2010) class, with progenitor radius larger than that of SN 2008ax and smaller than the eIIb (extended IIb) SN 1993J progenitor. Our data imply that the object identified in pre-explosion Hubble Space Telescope images at the SN location is possibly a companion to the progenitor or a blended source, and not the progenitor star itself, as its radius (~10^13 cm) would be highly inconsistent with constraints from our post-explosion photometric and spectroscopic data.

preprint2011arXiv

Rapid, Machine-Learned Resource Allocation: Application to High-redshift GRB Follow-up

As the number of observed Gamma-Ray Bursts (GRBs) continues to grow, follow-up resources need to be used more efficiently in order to maximize science output from limited telescope time. As such, it is becoming increasingly important to rapidly identify bursts of interest as soon as possible after the event, before the afterglows fade beyond detectability. Studying the most distant (highest redshift) events, for instance, remains a primary goal for many in the field. Here we present our Random forest Automated Triage Estimator for GRB redshifts (RATE GRB-z) for rapid identification of high-redshift candidates using early-time metrics from the three telescopes onboard Swift. While the basic RATE methodology is generalizable to a number of resource allocation problems, here we demonstrate its utility for telescope-constrained follow-up efforts with the primary goal to identify and study high-z GRBs. For each new GRB, RATE GRB-z provides a recommendation - based on the available telescope time - of whether the event warrants additional follow-up resources. We train RATE GRB-z using a set consisting of 135 Swift bursts with known redshifts, only 18 of which are z > 4. Cross-validated performance metrics on this training data suggest that ~56% of high-z bursts can be captured from following up the top 20% of the ranked candidates, and ~84% of high-z bursts are identified after following up the top ~40% of candidates. We further use the method to rank 200+ Swift bursts with unknown redshifts according to their likelihood of being high-z.

preprint2011arXiv

Real-Time Detection and Rapid Multiwavelength Follow-up Observations of a Highly Subluminous Type II-P Supernova from the Palomar Transient Factory Survey

The Palomar Transient Factory (PTF) is an optical wide-field variability survey carried out using a camera with a 7.8 square degree field of view mounted on the 48-in Oschin Schmidt telescope at Palomar Observatory. One of the key goals of this survey is to conduct high-cadence monitoring of the sky in order to detect optical transient sources shortly after they occur. Here, we describe the real-time capabilities of the PTF and our related rapid multiwavelength follow-up programs, extending from the radio to the gamma-ray bands. We present as a case study observations of the optical transient PTF10vdl (SN 2010id), revealed to be a very young core-collapse (Type II-P) supernova having a remarkably low luminosity. Our results demonstrate that the PTF now provides for optical transients the real-time discovery and rapid-response follow-up capabilities previously reserved only for high-energy transients like gamma-ray bursts.

preprint2011arXiv

SN2010jp (PTF10aaxi): A Jet-Driven Type II Supernova

We present photometry and spectroscopy of the peculiar TypeII supernova (SN) 2010jp, also named PTF10aaxi. The light curve exhibits a linear decline with a relatively low peak absolute magnitude of only -15.9, and a low radioactive decay luminosity at late times that suggests a nickel mass below 0.003 $M_{\odot}$. Spectra of SN2010jp display an unprecedented triple-peaked H$α$ line profile, showing: (1) a narrow (800 km/s) central component that suggests shock interaction with dense CSM; (2) high-velocity blue and red emission features centered at -12600 and +15400 km/s; and (3) broad wings extending from -22000 to +25000 km/s. These features persist during 100 days after explosion. We propose that this line profile indicates a bipolar jet-driven explosion, with the central component produced by normal SN ejecta and CSM interaction at mid latitudes, while the high-velocity bumps and broad line wings arise in a nonrelativistic bipolar jet. Two variations of the jet interpretation seem plausible: (1) A fast jet mixes 56Ni to high velocities in polar zones of the H-rich envelope, or (2) the reverse shock in the jet produces blue and red bumps in Balmer lines when a jet interacts with dense CSM. Jet-driven SNeII are predicted for collapsars resulting from a wide range of initial masses above 25 $M_{\odot}$ at sub-solar metallicity. This seems consistent with the SN host environment, which is either an extremely low-luminosity dwarf galaxy or very remote parts of an interacting pair of star-forming galaxies. It also seems consistent with the low 56Ni mass that may accompany black hole formation. We speculate that the jet survives to produce observable signatures because the star's H envelope was mostly stripped away by previous eruptive mass loss.

preprint2011arXiv

Supernova 2011fe from an Exploding Carbon-Oxygen White Dwarf Star

Type Ia supernovae (SNe Ia) have been used empirically as standardized candles to reveal the accelerating universe even though fundamental details, such as the nature of the progenitor system and how the star explodes, remained a mystery. There is consensus that a white dwarf star explodes after accreting matter in a binary system, but the secondary could be anything from a main sequence star to a red giant, or even another white dwarf. The uncertainty stems from the fact that no recent SN Ia has been discovered close enough to detect the stars before explosion. Here we report early observations of SN 2011fe (PTF11kly) in M101 at a distance of 6.4 Mpc, the closest SN Ia in the past 25 years. We find that the exploding star was likely a carbon-oxygen white dwarf, and from the lack of an early shock we conclude that the companion was most likely a main sequence star. Early spectroscopy shows high-velocity oxygen that varies on a time scale of hours and extensive mixing of newly synthesized intermediate mass elements in the outermost layers of the supernova. A companion paper uses pre-explosion images to rule out luminous red giants and most helium stars as companions.

preprint2011arXiv

The Gamma-Ray Burst - Supernova Connection

A preponderance of evidence links long-duration, soft-spectrum gamma-ray bursts (GRBs) with the death of massive stars. The observations of the GRB-supernova (SN) connection present the most direct evidence of this physical link. We summarize 30 GRB-SN associations and focus on five ironclad cases, highlighting the subsequent insight into the progenitors enabled by detailed observations. We also address the SN association (or lack thereof) with several sub-classes of GRBs, finding that the X-ray Flash (XRF) population is likely associated with massive stellar death whereas short-duration events likely arise from an older population not readily capable of producing a SN concurrent with a GRB. Interestingly, a minority population of seemingly long-duration, soft-spectrum GRBs show no evidence for SN-like activity; this may be a natural consequence of the range of Ni-56 production expected in stellar deaths.

preprint2011arXiv

The Palomar Transient Factory Orion Project: Eclipsing Binaries and Young Stellar Objects

The Palomar Transient Factory (PTF) Orion project is an experiment within the broader PTF survey, a systematic automated exploration of the sky for optical transients. Taking advantage of the wide field of view available using the PTF camera at the Palomar 48" telescope, 40 nights were dedicated in December 2009-January 2010 to perform continuous high-cadence differential photometry on a single field containing the young (7-10Myr) 25 Ori association. The primary motivation for the project is to search for planets around young stars in this region. The unique data set also provides for much ancillary science. In this first paper we describe the survey and data reduction pipeline, and present initial results from an inspection of the most clearly varying stars relating to two of the ancillary science objectives: detection of eclipsing binaries and young stellar objects. We find 82 new eclipsing binary systems, 9 of which we are candidate 25 Ori- or Orion OB1a-association members. Of these, 2 are potential young W UMa type systems. We report on the possible low-mass (M-dwarf primary) eclipsing systems in the sample, which include 6 of the candidate young systems. 45 of the binary systems are close (mainly contact) systems; one shows an orbital period among the shortest known for W UMa binaries, at 0.2156509 \pm 0.0000071d, with flat-bottomed primary eclipses, and a derived distance consistent with membership in the general Orion association. One of the candidate young systems presents an unusual light curve, perhaps representing a semi-detached binary system with an inflated low-mass primary or a star with a warped disk, and may represent an additional young Orion member. Finally, we identify 14 probable new classical T-Tauri stars in our data, along with one previously known (CVSO 35) and one previously reported as a candidate weak-line T-Tauri star (SDSS J052700.12+010136.8).

preprint2011arXiv

Three New Eclipsing White-dwarf - M-dwarf Binaries Discovered in a Search for Transiting Planets Around M-dwarfs

We present three new eclipsing white-dwarf / M-dwarf binary systems discovered during a search for transiting planets around M-dwarfs. Unlike most known eclipsing systems of this type, the optical and infrared emission is dominated by the M-dwarf components, and the systems have optical colors and discovery light curves consistent with being Jupiter-radius transiting planets around early M-dwarfs. We detail the PTF/M-dwarf transiting planet survey, part of the Palomar Transient Factory (PTF). We present a Graphics Processing Unit (GPU)-based box-least-squares search for transits that runs approximately 8X faster than similar algorithms implemented on general purpose systems. For the discovered systems, we decompose low-resolution spectra of the systems into white-dwarf and M-dwarf components, and use radial velocity measurements and cooling models to estimate masses and radii for the white dwarfs. The systems are compact, with periods between 0.35 and 0.45 days and semimajor axes of approximately 2 solar radii (0.01 AU). We use the Robo-AO laser guide star adaptive optics system to tentatively identify one of the objects as a triple system. We also use high-cadence photometry to put an upper limit on the white dwarf radius of 0.025 solar radii (95% confidence) in one of the systems. We estimate that 0.08% (90% confidence) of M-dwarfs are in these short-period, post-common-envelope white-dwarf / M-dwarf binaries where the optical light is dominated by the M-dwarf. Similar eclipsing binary systems can have arbitrarily small eclipse depths in red bands and generate plausible small-planet-transit light curves. As such, these systems are a source of false positives for M-dwarf transiting planet searches. We present several ways to rapidly distinguish these binaries from transiting planet systems.

preprint2010arXiv

A Mature Dusty Star-forming Galaxy Hosting GRB080607 at z=3.036

We report the discovery of the host galaxy of dark burst GRB080607 at z_GRB=3.036. GRB080607 is a unique case of a highly extinguished (A_V~3 mag) afterglow that was yet sufficiently bright for high-quality absorption-line spectroscopy. The host galaxy is clearly resolved in deep HST WF3/IR F160W images and well detected in the Spitzer IRAC 3.5 micron and 4.5 micron channels, while displaying little/no fluxes in deep optical images from Keck and Magellan. The extremely red optical-infrared colors are consistent with the large extinction seen in the afterglow light, suggesting that the large amount of dust and gas surface mass density seen along the afterglow sightline is not merely local but likely reflects the global dust content across the entire host galaxy. Adopting the dust properties and metallicity of the host ISM derived from studies of early-time afterglow light and absorption-line spectroscopy, we perform a stellar population synthesis analysis of the observed spectral energy distribution to constrain the intrinsic luminosity and stellar population of this dark burst host. The host galaxy is best described by an exponentially declining star formation rate of e-folding time tau=2 Gyr and an age of ~2 Gyr. We also derive an extinction corrected star formation rate of SFR 125 h^{-2} M_sun/yr and a total stellar mass of M_* ~ 4x10^11 h^{-2} M_sun. Our study provides an example of massive, dusty star-forming galaxies contributing to the GRB host galaxy population, supporting the notion that long-duration GRBs trace the bulk of cosmic star formation.

preprint2010arXiv

Core-Collapse Supernovae from the Palomar Transient Factory: Indications for a Different Population in Dwarf Galaxies

We use the first compilation of 72 core-collapse supernovae (SNe) from the Palomar Transient Factory (PTF) to study their observed subtype distribution in dwarf galaxies compared to giant galaxies. Our sample is the largest single-survey, untargeted, spectroscopically classified, homogeneous collection of core-collapse events ever assembled, spanning a wide host-galaxy luminosity range (down to M_r ~ -14 mag) and including a substantial fraction (>20%) of dwarf (M_r >= -18 mag) hosts. We find more core-collapse SNe in dwarf galaxies than expected and several interesting trends emerge. We use detailed subclassifications of stripped-envelope core-collapse SNe and find that all Type I core-collapse events occurring in dwarf galaxies are either SNe Ib or broad-lined SNe Ic (SNe Ic-BL), while "normal" SNe Ic dominate in giant galaxies. We also see a significant excess of SNe IIb in dwarf hosts. We hypothesize that in lower metallicity hosts, metallicity-driven mass loss is reduced, allowing massive stars that would have appeared as "normal" SNe Ic in metal-rich galaxies to retain some He and H, exploding as Ib/IIb events. At the same time, another mechanism allows some stars to undergo extensive stripping and explode as SNe Ic-BL (and presumably also as long-duration gamma-ray bursts). Our results are still limited by small-number statistics, and our measurements of the observed N(Ib/c)/N(II) number ratio in dwarf and giant hosts (0.25_{-0.15}^{+0.3} and 0.23_{-0.08}^{+0.11}, respectively; 1 sigma uncertainties) are consistent with previous studies and theoretical predictions. As additional PTF data accumulate, more robust statistical analyses will be possible, allowing the evolution of massive stars to be probed via the dwarf-galaxy SN population.

preprint2010arXiv

Discovery of Precursor LBV Outbursts in Two Recent Optical Transients: The Fitfully Variable Missing Links UGC 2773-OT and SN 2009ip

We present progenitor-star detections, light curves, and optical spectra of SN2009ip and the 2009 optical transient in UGC2773 (U2773-OT), which were not genuine SNe. Precursor variability in the decade before outburst indicates that both of the progenitor stars were LBVs. Their pre-outburst light curves resemble the S Doradus phases that preceded giant eruptions of eta Carinae and SN1954J (V12 in NGC2403), with intermediate progenitor luminosities. HST detections a decade before discovery indicate that the SN2009ip and U2773-OT progenitors were supergiants with likely initial masses of 50-80 Msun and $\ga$20 Msun, respectively. Both outbursts had spectra befitting known LBVs, although in different physical states. SN 2009ip exhibited a hot LBV spectrum with characteristic speeds of 550 km/s, plus faster material up to 5000 km/s, resembling the slow Homunculus and fast blast wave of eta Carinae. U2773-OT shows a forest of narrow absorption and emission lines comparable to that of S Dor in its cool state, plus [CaII] emission and an IR excess indicative of dust, similar to SN2008S and N300-OT. [CaII] emission is probably tied to a dusty pre-outburst environment, and not the outburst mechanism. SN2009ip and U2773-OT may provide a critical link between historical LBV eruptions, while U2773-OT may provide a link between LBVs and SN2008S and N300-OT. Future searches will uncover more examples of precursor LBV variability of this kind, providing key clues that may help unravel the instability driving LBVs.

preprint2010arXiv

Hubble Space Telescope Studies of Nearby Type Ia Supernovae: The Mean Maximum Light Ultraviolet Spectrum and its Dispersion

We present the first results of an ongoing campaign using the STIS spectrograph on-board the Hubble Space Telescope (HST) whose primary goal is the study of near ultraviolet (UV) spectra of local Type Ia supernovae (SNe Ia). Using events identified by the Palomar Transient Factory and subsequently verified by ground-based spectroscopy, we demonstrate the ability to locate and classify SNe Ia as early as 16 days prior to maximum light. This enables us to trigger HST in a non-disruptive mode to obtain near UV spectra within a few days of maximum light for comparison with earlier equivalent ground-based spectroscopic campaigns conducted at intermediate redshifts, z ~ 0.5. We analyze the spectra of 12 Type Ia supernovae located in the Hubble flow with 0.01 < z < 0.08. Although a fraction of our eventual sample, these data, together with archival data, already provide a substantial advance over that previously available. Restricting samples to those of similar phase and stretch, the mean UV spectrum agrees reasonably closely with that at intermediate redshift, although some differences are found in the metallic absorption features. A larger sample will determine whether these differences reflect possible sample biases or are a genuine evolutionary effect. Significantly, the wavelength-dependent dispersion, which is larger in the UV, follows similar trends to that observed at intermediate redshift and is driven, in part, by differences in the various metallic features. While the origin of the UV dispersion remains uncertain, our comparison suggests that it may reflect compositional variations amongst our sample rather than being predominantly an evolutionary effect.

preprint2010arXiv

PTF10nvg: An Outbursting Class I Protostar in the Pelican/North American Nebula

During a synoptic survey of the North American Nebula region, the Palomar Transient Factory (PTF) detected an optical outburst (dubbed PTF10nvg) associated with the previously unstudied flat or rising spectrum infrared source IRAS 20496+4354. The PTF R-band light curve reveals that PTF10nvg brightened by more than 5 mag during the current outburst, rising to a peak magnitude of R~13.5 in 2010 Sep. Follow-up observations indicate PTF10nvg has undergone a similar ~5 mag brightening in the K band, and possesses a rich emission-line spectrum, including numerous lines commonly assumed to trace mass accretion and outflows. Many of these lines are blueshifted by ~175 km/s from the North American Nebula's rest velocity, suggesting that PTF10nvg is driving an outflow. Optical spectra of PTF10nvg show several TiO/VO bandheads fully in emission, indicating the presence of an unusual amount of dense (> 10^10 cm^-3), warm (1500-4000 K) circumstellar material. Near-infrared spectra of PTF10nvg appear quite similar to a spectrum of McNeil's Nebula/V1647 Ori, a young star which has undergone several brightenings in recent decades, and 06297+1021W, a Class I protostar with a similarly rich near--infrared emission line spectrum. While further monitoring is required to fully understand this event, we conclude that the brightening of PTF10nvg is indicative of enhanced accretion and outflow in this Class-I-type protostellar object, similar to the behavior of V1647 Ori in 2004-2005.

preprint2010arXiv

Rapidly Decaying Supernova 2010X: A Candidate ".Ia" Explosion

We present the discovery, photometric and spectroscopic follow-up observations of SN 2010X (PTF 10bhp). This supernova decays exponentially with tau_d=5 days, and rivals the current recordholder in speed, SN 2002bj. SN 2010X peaks at M_r=-17mag and has mean velocities of 10,000 km/s. Our light curve modeling suggests a radioactivity powered event and an ejecta mass of 0.16 Msun. If powered by Nickel, we show that the Nickel mass must be very small (0.02 Msun) and that the supernova quickly becomes optically thin to gamma-rays. Our spectral modeling suggests that SN 2010X and SN 2002bj have similar chemical compositions and that one of Aluminum or Helium is present. If Aluminum is present, we speculate that this may be an accretion induced collapse of an O-Ne-Mg white dwarf. If Helium is present, all observables of SN 2010X are consistent with being a thermonuclear Helium shell detonation on a white dwarf, a ".Ia" explosion. With the 1-day dynamic-cadence experiment on the Palomar Transient Factory, we expect to annually discover a few such events.

preprint2010arXiv

The Cosmic Rate, Luminosity Function and Intrinsic Correlations of Long GRBs

We calculate durations and spectral parameters for 207 Swift bursts detected by the BAT instrument from April 2007 to August 2009, including 67 events with measured redshifts. This is the first supplement to our catalog of 425 Swift GRBs (147 with redshifts) starting from GRB041220. This complete and extensive data set, analyzed with a unified methodology, allows us to conduct an accurate census of intrinsic GRB energetics, hardnesses, durations, and redshifts. The GRB world model we derive reproduces well the observables from both Swift and pre-Swift satellites. Comparing to the cosmic star formation rate, we estimate that only about 0.1% of massive stars explode as bright GRBs. There is strong evidence for evolution in the Swift population at intermediate and high-z, and we can rule out (at the 5-sigma level) that this is due to evolution in the luminosity function of GRBs. Instead, the Swift sample suggests a modest propensity for low-metallicity, evidenced by an increase in the rate density with redshift. Treating the multivariate data and selection effects rigorously, we find a real, intrinsic correlation between E_iso and E_pk (and possibly also duration); however, the correlation {\it is not} a narrow log-log relation and its observed appearance is strongly detector-dependent. We also estimate the high-z rate (3--9% of GRBs at z beyond 5) and discuss the extent of a large missing population of low-E_pk XRFs as well as a potentially large missing population of short-duration GRBs that will be probed by EXIST.

preprint2009arXiv

An Unusually Fast-Evolving Supernova

Analyses of supernovae (SNe) have revealed two main types of progenitors: exploding white dwarfs and collapsing massive stars. We present SN2002bj, which stands out as different from any SN reported to date. Its light curve rises and declines very rapidly, yet reaches a peak intrinsic brightness greater than -18 mag. A spectrum obtained 7 days after discovery shows the presence of helium and intermediate-mass elements, yet no clear hydrogen or iron-peak elements. The spectrum only barely resembles that of a Type Ia supernova, with added carbon and helium. Its properties suggest that SN2002bj may be representative of a class of progenitors that previously has been only hypothesized: a helium detonation on a white dwarf, ejecting a small envelope of material. New surveys should find many such objects, despite their scarcity.

preprint2009arXiv

Astro2010 Decadal Survey Whitepaper: Coordinated Science in the Gravitational and Electromagnetic Skies

It is widely expected that the coming decade will witness the first direct detection of gravitational waves (GWs). The ground-based LIGO and Virgo GW observatories are being upgraded to advanced sensitivity, and are expected to observe a significant binary merger rate. The launch of The Laser Interferometer Space Antenna (LISA) would extend the GW window to low frequencies, opening new vistas on dynamical processes involving massive (M >~ 10^5 M_Sun) black holes. GW events are likely to be accompanied by electromagnetic (EM) counterparts and, since information carried electromagnetically is complementary to that carried gravitationally, a great deal can be learned about an event and its environment if it becomes possible to measure both forms of radiation in concert. Measurements of this kind will mark the dawn of trans-spectral astrophysics, bridging two distinct spectral bands of information. The aim of this whitepaper is to articulate future directions in both theory and observation that are likely to impact broad astrophysical inquiries of general interest. What will EM observations reflect on the nature and diversity of GW sources? Can GW sources be exploited as complementary probes of cosmology? What cross-facility coordination will expand the science returns of gravitational and electromagnetic observations?

preprint2009arXiv

Exploring the Optical Transient Sky with the Palomar Transient Factory

The Palomar Transient Factory (PTF) is a wide-field experiment designed to investigate the optical transient and variable sky on time scales from minutes to years. PTF uses the CFH12k mosaic camera, with a field of view of 7.9 deg^2 and a plate scale of 1 asec/pixel, mounted on the the Palomar Observatory 48-inch Samuel Oschin Telescope. The PTF operation strategy is devised to probe the existing gaps in the transient phase space and to search for theoretically predicted, but not yet detected, phenomena, such as fallback supernovae, macronovae, .Ia supernovae and the orphan afterglows of gamma-ray bursts. PTF will also discover many new members of known source classes, from cataclysmic variables in their various avatars to supernovae and active galactic nuclei, and will provide important insights into understanding galactic dynamics (through RR Lyrae stars) and the Solar system (asteroids and near-Earth objects). The lessons that can be learned from PTF will be essential for the preparation of future large synoptic sky surveys like the Large Synoptic Survey Telescope. In this paper we present the scientific motivation for PTF and describe in detail the goals and expectations for this experiment.

preprint2009arXiv

GJ 3236: a new bright, very low-mass eclipsing binary system discovered by the MEarth observatory

We report the detection of eclipses in GJ 3236, a bright (I = 11.6) very low mass binary system with an orbital period of 0.77 days. Analysis of light- and radial velocity curves of the system yielded component masses of 0.38 +/- 0.02 and 0.28 +/- 0.02 Msol. The central values for the stellar radii are larger than the theoretical models predict for these masses, in agreement with the results for existing eclipsing binaries, although the present 5% observational uncertainties limit the significance of the larger radii to approximately 1 sigma. Degeneracies in the light curve models resulting from the unknown configuration of surface spots on the components of GJ 3236 currently dominate the uncertainties in the radii, and could be reduced by obtaining precise, multi-band photometry covering the full orbital period. The system appears to be tidally synchronized and shows signs of high activity levels as expected for such a short orbital period, evidenced by strong Halpha emission lines in the spectra of both components. These observations probe an important region of mass-radius parameter space around the predicted transition to fully-convective stellar interiors, where there are a limited number of precise measurements available in the literature.

preprint2009arXiv

GRB 090426: The Environment of a Rest-Frame 0.35-second Gamma-Ray Burst at Redshift z=2.609

We present the discovery of an absorption-line redshift of z = 2.609 for GRB 090426, establishing the first firm lower limit to a redshift for a gamma-ray burst with an observed duration of <2 s. With a rest-frame burst duration of T_90z = 0.35 s and a detailed examination of the peak energy of the event, we suggest that this is likely (at >90% confidence) a member of the short/hard phenomenological class of GRBs. From analysis of the optical-afterglow spectrum we find that the burst originated along a very low HI column density sightline, with N_HI < 3.2 x 10^19 cm^-2. Our GRB 090426 afterglow spectrum also appears to have weaker low-ionisation absorption (Si II, C II) than ~95% of previous afterglow spectra. Finally, we also report the discovery of a blue, very luminous, star-forming putative host galaxy (~2 L*) at a small angular offset from the location of the optical afterglow. We consider the implications of this unique GRB in the context of burst duration classification and our understanding of GRB progenitor scenarios.

preprint2009arXiv

Improved Standardization of Type II-P Supernovae: Application to an Expanded Sample

In the epoch of precise and accurate cosmology, cross-confirmation using a variety of cosmographic methods is paramount to circumvent systematic uncertainties. Owing to progenitor histories and explosion physics differing from those of Type Ia SNe (SNe Ia), Type II-plateau supernovae (SNe II-P) are unlikely to be affected by evolution in the same way. Based on a new analysis of 17 SNe II-P, and on an improved methodology, we find that SNe II-P are good standardizable candles, almost comparable to SNe Ia. We derive a tight Hubble diagram with a dispersion of 10% in distance, using the simple correlation between luminosity and photospheric velocity introduced by Hamuy & Pinto 2002. We show that the descendent method of Nugent et al. 2006 can be further simplified and that the correction for dust extinction has low statistical impact. We find that our SN sample favors, on average, a very steep dust law with total to selective extinction R_V<2. Such an extinction law has been recently inferred for many SNe Ia. Our results indicate that a distance measurement can be obtained with a single spectrum of a SN II-P during the plateau phase combined with sparse photometric measurements.

preprint2009arXiv

Limits on Radioactive-Powered Emission Associated With a Short-Hard GRB 070724A in a Star-Forming Galaxy

We present results of an extensive observing campaign of the short duration, hard spectrum gamma-ray burst (GRB) 070724A, aimed at detecting the radioactively-powered emission that might follow from a binary merger or collapse involving compact objects. Our multi-band observations span the range in time over which this so-called Li-Paczynski mini-supernova could be active, beginning within 3 hours of the GRB trigger, and represent some of the deepest and most comprehensive searches for such emission. We find no evidence for such activity and place limits on the abundances and the lifetimes of the possible radioactive nuclides that could form in the rapid decompression of nuclear-density matter. Furthermore, our limits are significantly fainter than the peak magnitude of any previously detected broad-lined Type Ic supernova (SN) associated with other GRBs, effectively ruling out a long GRB-like SN for with this event. Given the unambiguous redshift of the host galaxy (z=0.456), GRB 070724A represents one of a small, but growing, number of short-hard GRBs for which firm physical/restframe quantities currently exist. The host of GRB 070724A is a moderately star-forming galaxy with an older stellar population component and a relatively high metallicity of 12+log(O/H)_KD02=9.1. We find no significant evidence for large amounts of extinction along the line of sight that could mask the presence of a SN explosion and estimate a small probability for chance alignment with the putative host. We discuss how our derived constraints fit into the evolving picture of short-hard GRBs, their potential progenitors, and the host environments in which they are thought to be produced.

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.

preprint2008arXiv

A New Low-Mass Eclipsing Binary from SDSS-II

We present observations of a new low-mass double-lined eclipsing binary system discovered using repeat observations of the celestial equator from the Sloan Digital Sky Survey II. Using near-infrared photometry and optical spectroscopy we have measured the properties of this short-period [P=0.407037(14) d] system and its two components. We find the following parameters for the two components: M_1=0.272+/-0.020 M_sun, R_1=0.268+/-0.010 R_sun, M_2=0.240+/-0.022 M_sun, R_2=0.248+/-0.0090 R_sun, T_1=3320+/-130 K, T_2=3300+/-130 K. The masses and radii of the two components of this system agree well with theoretical expectations based on models of low-mass stars, within the admittedly large errors. Future synoptic surveys like Pan-STARRS and LSST will produce a wealth of information about low-mass eclipsing systems and should make it possible, with an increased reliance on follow-up observations, to detect many systems with low-mass and sub-stellar companions. With the large numbers of objects for which these surveys will produce high-quality photometry, we suggest that it becomes possible to identify such systems even with sparse time sampling and a relatively small number of individual observations.

preprint2007arXiv

Keck and ESO-VLT View of the Symmetry of the Ejecta of the XRF/SN 2006aj

Nebular-phase spectra of SN 2006aj, which was discovered in coincidence with X-ray flash 060218, were obtained with Keck in 2006 July and the Very Large Telescope in 2006 September. At the latter epoch spectropolarimetry was also attempted, yielding an upper limit of ~ 2% for the polarization. The spectra show strong emission lines of [OI] and MgI], as expected from a Type Ic supernova, but weak CaII lines. The [FeII] lines that were strong in the spectra of SN 1998bw are much weaker in SN 2006aj, consistent with the lower luminosity of this SN. The outer velocity of the line-emitting ejecta is ~ 8000 km/s in July and ~ 7400 km/s in September, consistent with the relatively low kinetic energy of expansion of SN 2006aj. All emission lines have similar width, and the profiles are symmetric, indicating that no major asymmetries are present in the ejecta at the velocities sampled by the nebular lines (v < 8000 km/s), except perhaps in the innermost part. The spectra were modelled with a non-LTE code. The mass of 56Ni required to power the emission spectrum is ~ 0.20 Msun, in excellent agreement with the results of early light curve modelling. The oxygen mass is ~ 1.5 Msun, again much less than in SN 1998bw but larger by ~ 0.7 Msun than the value derived from the early-time modelling. The total ejected mass is ~ 2 Msun below 8000 km/s. This confirms that SN 2006aj was only slightly more massive and energetic than the prototypical Type Ic SN 1994I, but also indicates the presence of a dense inner core, containing ~ 1 Msun of mostly oxygen and carbon. The presence of such a core is inferred for all broad-lined SNe Ic. This core may have the form of an equatorial oxygen-dominated region, but it is too deep to affect the early light curve and too small to affect the late polarization spectrum.

preprint2007arXiv

Multicolor Infrared Observations of SN 2006aj, the Supernova Associated with XRF 060218 - Paper I

We report simultaneous multicolor near-infrared (NIR) observations of the supernova associated with x-ray Flash 060218 during the first 16 days after the high energy event. We find that the light curve rises and peaks relatively fast compared to other SN Ic, with the characteristic broad NIR peak seen in all three bands. We find that the rise profile before the peak is largely independent of NIR wavelength, each band appearing to transition into a plateau phase around day 10--13. Since the light curve is in the plateau phase when our observations end at day 16, we can only place limits on the peak absolute magnitudes, but we estimate that SN 2006aj is one of the lowest NIR luminosity XRF/GRB associated SNe observed to date. The broad peaks observed in the {\em JHK$_s$} bands point to a large increase in the NIR contribution of the total flux output from days 10--16. This evolution can be seen in the broad color and SED diagrams constructed using {\em UBVRIJHK$_s$} monochromatic flux measurements for the first 16 days of the event. Ultimately, a 10-day rise time would make SN 2006aj an extremely fast rise SN Ic event, faster than SN 1998bw and SN 2003dh, which combined with its underluminous nature, indicates a lower amount of $^{56}$Ni ejected by the progenitor compared to other XRF/GRB-SNe. Furthermore, the lack of significant color change during the rise portion of the burst points to little or no spectral evolution over the first 10 days of activity in the NIR.

preprint2007arXiv

Probing the ISM Near Star Forming Regions with GRB Afterglow Spectroscopy: Gas, Metals, and Dust

We study the chemical abundances of the interstellar medium surrounding high z gamma-ray bursts (GRBs) through analysis of the damped Lya systems (DLAs) identified in afterglow spectra. These GRB-DLAs are characterized by large HI column densities N(HI) and metallicities [M/H] spanning 1/100 to nearly solar, with median [M/H]>-1. The majority of GRB-DLAs have [M/H] values exceeding the cosmic mean metallicity of atomic gas at z>2, i.e. if anything, the GRB-DLAs are biased to larger metallicity. We also observe (i) large [Zn/Fe] values (>+0.6) and sub-solar Ti/Fe ratios which imply substantial differential depletion, (ii) large a/Fe ratios suggesting nucleosynthetic enrichment by massive stars, and (iii) low C^0/C^+ ratios (<10^{-4}). Quantitatively, the observed depletion levels and C^0/C^+ ratios of the gas are not characteristic of cold, dense HI clouds in the Galactic ISM. We argue that the GRB-DLAs represent the ISM near the GRB but not gas directly local to the GRB (e.g. its molecular cloud or circumstellar material). We compare these observations with DLAs intervening background quasars (QSO-DLAs). The GRB-DLAs exhibit larger N(HI) values, higher a/Fe and Zn/Fe ratios, and have higher metallicity than the QSO-DLAs. We argue that the differences primarily result from galactocentric radius-dependent differences in the ISM: GRB-DLAs preferentially probe denser, more depleted, higher metallicity gaslocated in the inner few kpc whereas QSO-DLAs are more likely to intersect the less dense, less enriched, outer regions of the galaxy. Finally, we investigate whether dust obscuration may exclude GRB-DLA sightlines from QSO-DLA samples; we find that the majority of GRB-DLAs would be recovered which implies little observational bias against large N(HI) systems.

preprint2007arXiv

SN 2006gy: Discovery of the most luminous supernova ever recorded, powered by the death of an extremely massive star like Eta Carinae

(abridged) We report our discovery and observations of the peculiar Type IIn supernova SN2006gy in NGC1260, revealing that it reached a peak magnitude of -22, making it the most luminous supernova ever recorded. It is not yet clear what powers the total radiated energy of 1e51 erg, but we argue that any mechanism -- thermal emission, circumstellar interaction, or 56Ni decay -- requires a very massive progenitor star. The circumstellar interaction hypothesis would require truly exceptional conditions around the star probably experienced an LBV eruption like the 19th century eruption of eta Carinae. Alternatively, radioactive decay of 56Ni may be a less objectionable hypothesis. That power source would imply a large Ni mass of 22 Msun, requiring that SN2006gy was a pair-instability supernova where the star's core was obliterated. SN2006gy is the first supernova for which we have good reason to suspect a pair-instability explosion. Based on a number of lines of evidence, we rule out the hypothesis that SN 2006gy was a ``Type IIa'' event. Instead, we propose that the progenitor may have been a very massive evolved object like eta Carinae that, contrary to expectations, failed to completely shed its massive hydrogen envelope before it died. Our interpretation of SN2006gy implies that the most massive stars can explode earlier than expected, during the LBV phase, preventing them from ever becoming Wolf-Rayet stars. SN2006gy also suggests that the most massive stars can create brilliant supernovae instead of dying ignominious deaths through direct collapse to a black hole.

preprint2006arXiv

Temporal Variation in Excited Fe^+ Near a GRB Afterglow

Excited Si^+ and Fe^+ species are routinely observed in the host environment of gamma-ray burst (GRB) afterglows, but are not commonly seen in other extragalactic locations. Their presence signals unusual properties in the gaseous environment of these GRB hosts that arise either as a result of the intense ionizing radiation of the afterglow or through collision excitation in a dense cloud. In particular, the photon pumping scenario has explicit expectations for temporal variation in the strength of the excited lines, owing to the decline in the ionizing flux of the GRB afterglow. We analyze afterglow spectra of GRB 020813 obtained in two epochs of ~16 hours apart, and examine transitions from the first excited state of Fe^+ at J=7/2 in these two sets of data. We report a significant decline by at least a factor of five in the equivalent width of the FeII 2396 transition, the strongest from the J=7/2 state. We perform a Monte-Carlo analysis and determine that this temporal variation is present at more than 3sigma level of significance. This observation represents the first detection in the temporal variation of the excited Fe+ states in the GRB host ISM, a direct influence of the burst itself on its environment. We further estimate that the Fe^+ gas resides in 50-100 pc from the afterglow, based on the afterglow lightcurve and the presence and absence of the excited FeII 2396 in the two-epoch observations.

preprint2005arXiv

Present and Future Prospects for GRB Standard Candles

Following our previous work, we conclude that a GRB standard candle constructed from the Ghirlanda et al. power-law relation between the geometry-corrected energy (E_gamma) and the peak of the rest-frame prompt burst spectrum (E_p) is not yet cosmographically useful, despite holding some potential advantages over SNe Ia. This is due largely to the small sample of \~20 GRBs with the required measured redshifts, jet-breaks, and peak energies, and to the strong sensitivity of the goodness-of-fit of the power-law to input assumptions. The most important such finding concerns the sensitivity to the generally unknown density (and density profile), of the circumburst medium. Although the E_p-E_gamma relation is a highly significant correlation over many cosmologies, until the sample expands to include many low-z events, it will be most sensitive to Omega_M but essentially insensitive to Omega_Lambda and w, with some hope of constraining dw/dt with high-z GRB data alone. The relation clearly represents a significant improvement in the search for an empirical GRB standard candle, but is further hindered by an unknown physical basis for the relation, the lack of a low-z training set to calibrate the relation in a cosmology-independent way, and several major potential systematic uncertainties and selection effects. Until these concerns are addressed, a larger sample is acquired, and attempts are made to marginalize or perform Monte Carlo simulations over the unknown density distribution, we urge caution concerning claims of the utility of GRBs for cosmography and especially the attempts to combine GRBs with SNe Ia.

preprint2004arXiv

The Host Galaxy of GRB 031203: Implications of its low metallicity, low redshift, and starburst nature

We present Keck/NIRSPEC near-IR images and Magellan/IMACS optical spectroscopy of the host galaxy of GRB 031203. The host is an actively star-forming galaxy at z=0.1055 +/- 0.0001. This is the lowest redshift GRB to-date, aside from GRB 980425. From the hydrogen Balmer lines, we infer an extinction of A_V = 3.62 +/- 0.25 or a total reddening E_T(B-V) = 1.17 +/- 0.1 toward the sightline to the nebular regions. After correcting for reddening, we perform an emission-line analysis and derive an ISM temperature of T=13400+/-2000K and electron density of n_e = 300 cm^(-3). These imply a metallicity [O/H]=-0.72+/-0.15 dex and a roughly solar abundance pattern for N, Ne, S, and Ar. Integrating Ha, we infer a dust-corrected star formation rate (SFR) of > 11 Msol/yr. These observations have the following implications: (1) the galaxy has a low K'-band luminosity L ~ L^*/5, typical of GRB host galaxies; (2) the low redshift indicates GRB 031203 had an isotropic-equivalent gamma-ray energy release smaller than all previous confirmed GRB events. The burst discovery raises the likelihood of identifying many additional low z, low flux events with Swift; (3) the large SFR, low metallicity, and the inferred hard radiation field is suggestive of massive star formation, supporting the collapsar model; (4) several lines of evidence argue against the identification of GRB 031203 as an X-ray flash event.

preprint1997arXiv

Using STIS to find Gamma-Ray Burst Redshifts

A recent spectrum of the optical afterglow of GRB 970508 suggests that gamma-ray bursts (GRBs) are cosmological in origin and it is of crucial importance to derive an accurate distance to each burst. If GRBs occur near their host galaxies then Lyman limit absorption should be observable in roughly half the GRB afterglow spectra. Here we outline the methodology to obtain a redshift from the GRB afterglow spectrum using the recently installed Space Telescope Imaging Spectrograph (STIS) instrument onboard the Hubble Space Telescope. A low-resolution spectrum with the Multi-Anode Microchannel Array (MAMA) detector gives complete spectral coverage over the wavelength range 157.0-318.0 nm (Near UV) and 115.0-174.0 nm (Far UV). Assuming a Target of Opportunity observation is conducted soon (< 3 weeks) after a bright burst, a relatively small integration time (about 3 orbits) would be sufficient to detect the Lyman limit over a wide redshift range (0.3 < z < 2.2). Detection (or non-detection) of the Lyman limit, in concert with ground-based observations of nearby galaxies and Mg II and C IV absorption lines, should provide meaningful constraints on the relationship of GRBs to galaxies.