We present multi-wavelength observations of the tidal disruption event (TDE) iPTF15af, discovered by the intermediate Palomar Transient Factory (iPTF) survey at redshift $z=0.07897$. The optical and ultraviolet (UV) light curves of the transient show a slow decay over five months, in agreement with previous optically discovered TDEs. It also has a comparable black-body peak luminosity of $L_{\rm{peak}} \approx 1.5 \times 10^{44}$ erg/s. The inferred temperature from the optical and UV data shows a value of (3$-$5) $\times 10^4$ K. The transient is not detected in X-rays up to $L_X < 3 \times 10^{42}$erg/s within the first five months after discovery. The optical spectra exhibit two distinct broad emission lines in the He II region, and at later times also H$\alpha$ emission. Additionally, emission from [N III] and [O III] is detected, likely produced by the Bowen fluorescence effect. UV spectra reveal broad emission and absorption lines associated with high-ionization states of N V, C IV, Si IV, and possibly P V. These features, analogous to those of broad absorption line quasars (BAL QSOs), require an absorber with column densities $N_{\rm{H}} > 10^{23}$ cm$^{-2}$. This optically thick gas would also explain the non-detection in soft X-rays. The profile of the absorption lines with the highest column density material at the largest velocity is opposite that of BAL QSOs. We suggest that radiation pressure generated by the TDE flare at early times could have provided the initial acceleration mechanism for this gas. Spectral UV line monitoring of future TDEs could test this proposal.

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)

Despite the utility of neural networks (NNs) for astronomical time-series classification, the proliferation of learning architectures applied to diverse datasets has thus far hampered a direct intercomparison of different approaches. Here we perform the first comprehensive study of variants of NN-based learning and inference for astronomical time-series, aiming to provide the community with an overview on relative performance and, hopefully, a set of best-in-class choices for practical implementations. In both supervised and self-supervised contexts, we study the effects of different time-series-compatible layer choices, namely the dilated temporal convolutional neural network (dTCNs), Long-Short Term Memory (LSTM) NNs, Gated Recurrent Units (GRUs) and temporal convolutional NNs (tCNNs). We also study the efficacy and performance of encoder-decoder (i.e., autoencoder) networks compared to direct classification networks, different pathways to include auxiliary (non-time-series) metadata, and different approaches to incorporate multi-passband data (i.e., multiple time-series per source). Performance---applied to a sample of 17,604 variable stars from the MACHO survey across 10 imbalanced classes---is measured in training convergence time, classification accuracy, reconstruction error, and generated latent variables. We find that networks with Recurrent NN (RNNs) generally outperform dTCNs and, in many scenarios, yield to similar accuracy as tCNNs. In learning time and memory requirements, convolution-based layers are more performant. We conclude by discussing the advantages and limitations of deep architectures for variable star classification, with a particular eye towards next-generation surveys such as LSST, WFIRST and ZTF2.

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 $\Delta 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--$\Delta m_{15}$ space, Model~II predicts phases with a root-mean-square error (RMSE) of 1.00\,d and Model~III predicts $\Delta m_{15}$ values with an RMSE of 0.068\,mag.

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.

Synoptic searches for the optical counterpart to a binary neutron star (BNS) or neutron star-black hole (NSBH) merger can pose significant challenges towards the discovery of kilonovae and performing multi-messenger science. In this work, we describe the advantage of a global multi-telescope network towards this end, with a particular focus on the key and complementary role the Dark Energy Camera (DECam) plays in multi-facility follow-up. We describe the Global Relay of Observatories Watching Transients Happen (GROWTH) Target-of-Opportunity (ToO) Marshal, a common web application we built to ingest events, plan observations, search for transient candidates, and retrieve performance summary statistics for all of the telescopes in our network. Our infrastructure enabled us to conduct observations of two events during O3a, S190426c and S190510g. Furthermore, our analysis of deep DECam observations of S190814bv conducted by the DESGW team, and access to a variety of global follow-up facilities allowed us to place meaningful constraints on the parameters of the kilonova and the merging binary. We emphasize the importance of a global telescope network in conjunction with a power telescope like DECam in performing searches for the counterparts to gravitational-wave sources.

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 $\kappa < 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 $ \chi < 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.

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 $\alpha_\mathrm{O,decay}=1.32\pm 0.01$, and a $\Delta 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.

The LIGO/Virgo gravitational wave events S190828j and S190828l were detected only 21 minutes apart, from nearby regions of sky, and with the same source classifications (binary black hole mergers). It is therefore natural to speculate that the two signals are actually strongly lensed images of the same merger. However, an estimate of the separation of the (unknown) positions of the two events requires them to be >10 deg apart, much wider than the arcsecond-scale separations that usually arise in extragalactic lensing. The large separation is much more consistent with two independent, unrelated events that occurred close in time by chance. We quantify the overlap between simulated pairs of lensed events, and use frequentist hypothesis testing to reject S190828j/l as a lensed pair at 99.8% confidence.

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.

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.

The first two months of the third Advanced LIGO and Virgo observing run (2019 April-May) showed that distant gravitational wave (GW) events can now be readily detected. Three candidate mergers containing neutron stars (NS) were reported in a span of 15 days, all likely located more than 100 Mpc away. However, distant events such as the three new NS mergers are likely to be coarsely localized, which highlights the importance of facilities and scheduling systems that enable deep observations over hundreds to thousands of square degrees to detect the electromagnetic counterparts. On 2019-05-10 02:59:39.292 UT the GW candidate S190510g was discovered and initially classified as a BNS merger with 98% probability. The GW event was localized within an area of 3462 deg2, later refined to 1166 deg2 (90%) at a distance of 227 +- 92 Mpc. We triggered Target of Opportunity observations with the Dark Energy Camera (DECam), a wide-field optical imager mounted at the prime focus of the 4m Blanco Telescope at CTIO in Chile. This Letter describes our DECam observations and our real-time analysis results, focusing in particular on the design and implementation of the observing strategy. Within 24 hours of the merger time, we observed 65% of the total enclosed probability of the final skymap with an observing efficiency of 94%. We identified and publicly announced 13 candidate counterparts. S190510g was re-classified 1.7 days after the merger, after our observations were completed, with a "binary neutron star merger" probability reduced from 98% to 42% in favor of a "terrestrial" classification.

We present observations of the possible short GRB 180418A in $\gamma$-rays, X-rays, and in the optical. Early optical photometry with the TAROT and RATIR instruments show a bright peak ($\approx$ 14.2 AB mag) between $T+28$ and $T+90$ seconds that we interpret as the signature of a reversal shock. Later observations can be modeled by a standard forward shock model and show no evidence of jet break, allowing us to constrain the jet collimation to $\theta_j> 7^\circ$. Using deep late-time optical observations we place an upper limit of $r>24$ AB mag on any underlying host galaxy. The detection of the afterglow in the \textit{Swift} UV filters constrains the GRB redshift to $z<1.3$ and places an upper bound on the $\gamma$-ray isotropic equivalent energy $E_{\rm{\gamma,iso}} < 3 \times 10^{51}$ erg. The properties of this GRB (e.g. duration, hardness ratio, energetic, and environment) lie at the intersection between short and long bursts, and we can not conclusively identify its type. We estimate that the probability that it is drawn from the population of short GRBs is 10\%-30\%.

The number of long gamma-ray bursts (GRBs) known to have occurred in the distant Universe (z greater than 5) is small (approx 15), however these events provide a powerful way of probing star formation at the onset of galaxy evolution. In this paper, we present the case for GRB100205A being a largely overlooked high-redshift event. While initially noted as a high-z candidate, this event and its host galaxy have not been explored in detail. By combining optical and near-infrared Gemini afterglow imaging (at t less than 1.3 days since burst) with deep late-time limits on host emission from the Hubble Space Telescope, we show that the most likely scenario is that GRB100205A arose in the redshift range 4-8. GRB100205A is an example of a burst whose afterglow, even at 1 hour post-burst, could only be identified by 8m class IR observations, and suggests that such observations of all optically dark bursts may be necessary to significantly enhance the number of high-redshift GRBs known.

The discovery of a transient kilonova following the gravitational-wave event GW170817 highlighted the critical need for coordinated rapid and wide-field observations, inference, and follow-up across the electromagnetic spectrum. In the Southern hemisphere, the Dark Energy Camera (DECam) on the Blanco 4-m telescope is well-suited to this task, as it is able to cover wide-fields quickly while still achieving the depths required to find kilonovae like the one accompanying GW170817 to $\sim$500 Mpc, the binary neutron star horizon distance for current generation of LIGO/Virgo collaboration (LVC) interferometers. Here, as part of the multi-facility followup by the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration, we describe the observations and automated data movement, data reduction, candidate discovery, and vetting pipeline of our target-of-opportunity DECam observations of S190426c, the first possible neutron star--black hole merger detected via gravitational waves. Starting 7.5hr after S190426c, over 11.28\,hr of observations, we imaged an area of 525\,deg$^2$ ($r$-band) and 437\,deg$^2$ ($z$-band); this was 16.3\% of the total original localization probability and nearly all of the probability density visible from the Southern hemisphere. The machine-learning based pipeline was optimized for fast turnaround, delivering transient candidates for human vetting within 17 minutes, on average, of shutter closure. We reported nine promising counterpart candidates 2.5 hours before the end of our observations. Our observations yielded no detection of a bona fide counterpart to $m_z = 22.5$ and $m_r = 22.9$ at the 5$\sigma$ level of significance, consistent with the refined LVC positioning. We view these observations and rapid inferencing as an important real-world test for this novel end-to-end wide-field pipeline.

Whether supernovae are a significant source of dust has been a long-standing debate. The large quantities of dust observed in high-redshift galaxies raise a fundamental question as to the origin of dust in the Universe since stars cannot have evolved to the AGB dust-producing phase in high-redshift galaxies. In contrast, supernovae occur within several millions of years after the onset of star formation. This white paper focuses on dust formation in supernova ejecta with US-Extremely Large Telescope (ELT) perspective during the era of JWST and LSST.

Supernovae are among the most powerful and influential explosions in the universe. They are also ideal multi-messenger laboratories to study extreme astrophysics. However, many fundamental properties of supernovae related to their diverse progenitor systems and explosion mechanisms remain poorly constrained. Here we outline how late-time spectroscopic observations obtained during the nebular phase (several months to years after explosion), made possible with the next generation of Extremely Large Telescopes, will facilitate transformational science opportunities and rapidly accelerate the community towards our goal of achieving a complete understanding of supernova explosions. We highlight specific examples of how complementary GMT and TMT instrumentation will enable high fidelity spectroscopy from which the line profiles and luminosities of elements tracing mass loss and ejecta can be used to extract kinematic and chemical information with unprecedented detail, for hundreds of objects. This will provide uniquely powerful constraints on the evolutionary phases stars may experience approaching a supernova explosion; the subsequent explosion dynamics; their nucleosynthesis yields; and the formation of compact objects that may act as central engines.

We present the volumetric rate of normal type Ia supernovae (SNe Ia) discovered by the Palomar Transient Factory (PTF). Using strict data-quality cuts, and considering only periods when the PTF maintained a regular cadence, PTF discovered 90 SNe Ia at $z\le0.09$ in a well-controlled sample over three years of operation (2010-2012). We use this to calculate the volumetric rate of SN Ia events by comparing this sample to simulations of hundreds of millions of SN Ia light curves produced in statistically representative realisations of the PTF survey. This quantifies the recovery efficiency of each PTF SN Ia event, and thus the relative weighting of each event. From this, the volumetric SN Ia rate was found to be $r_v=2.43\pm0.29\,\text{(stat)}_{-0.19}^{+0.33}\text{(sys)}\times10^{-5}\,\text{SNe yr}^{-1}\,\text{Mpc}^{-3}\, h_{70}^{3}$. This represents the most precise local measurement of the SN Ia rate. We fit a simple SN Ia delay-time distribution model, $\propto\mathrm{t}^{-\beta}$, to our PTF rate measurement combined with a literature sample of rate measurements from surveys at higher-redshifts. We find $\beta{\sim}1$, consistent with a progenitor channel governed by the gravitational in-spiral of binary white dwarfs.

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.

The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up, and analysis across the different messengers will also explosively grow. The cyberinfrastructure requirements to enhance multi-messenger astrophysics will both be a major challenge and opportunity for astronomers, physicists, computer scientists and cyberinfrastructure specialists. Here we outline some of these requirements and argue for a distributed cyberinfrastructure institute for multi-messenger astrophysics to meet these challenges.

Astrophysical transients have been observed for millennia and have shaped our most basic assumptions about the Universe. In the last century, systematic searches have grown from detecting handfuls of transients per year to over 7000 in 2018 alone. As these searches have matured, we have discovered both large samples of "normal" classes and new, rare classes. Recently, a transient was the first object observed in both gravitational waves and light. Ground-based observatories, including LSST, will discover thousands of transients in the optical, but these facilities will not provide the high-fidelity near-infrared (NIR) photometry and high-resolution imaging of a space-based observatory. WFIRST can fill this gap. With its survey designed to measure the expansion history of the Universe with Type Ia supernovae, WFIRST will also discover and monitor thousands of other transients in the NIR, revealing the physics for these high-energy events. Small-scale GO programs, either as a supplement to the planned survey or as specific target-of-opportunity observations, would significantly expand the scope of transient science that can be studied with WFIRST.

We investigate the prompt emission and the afterglow properties of short duration gamma-ray burst (sGRB) 130603B and another eight sGRB events during 2012-2015, observed by several multi-wavelength facilities including the GTC 10.4m telescope. Prompt emission high energy data of the events were obtained by INTEGRAL/SPI/ACS, Swift/BAT and Fermi/GBM satellites. The prompt emission data by INTEGRAL in the energy range of 0.1-10 MeV for sGRB 130603B, sGRB 140606A, sGRB 140930B, sGRB 141212A and sGRB 151228A do not show any signature of the extended emission or precursor activity and their spectral and temporal properties are similar to those seen in case of other short bursts. For sGRB130603B, our new afterglow photometric data constraints the pre jet-break temporal decay due to denser temporal coverage. For sGRB 130603B, the afterglow light curve, containing both our new as well as previously published photometric data is broadly consistent with the ISM afterglow model. Modeling of the host galaxies of sGRB 130603B and sGRB 141212A using the LePHARE software supports a scenario in which the environment of the burst is undergoing moderate star formation activity. From the inclusion of our late-time data for 8 other sGRBs we are able to; place tight constraints on the non-detection of the afterglow, host galaxy or any underlying kilonova emission. Our late-time afterglow observations of the sGRB 170817A/GW170817 are also discussed and compared with the sub-set of sGRBs.

We have obtained repeated images of 6 fields towards the Galactic bulge in 5 passbands (u, g, r, i, z) with the DECam imager on the Blanco 4m telescope at CTIO. From over 1.6 billion individual photometric measurements in the field centered on Baade's window, we have detected 4877 putative variable stars. 474 of these have been confirmed as fundamental mode RR Lyrae stars, whose colors at minimum light yield line-of-sight reddening determinations as well as a reddening law towards the Galactic Bulge which differs significantly from the standard R_V = 3.1 formulation. Assuming that the stellar mix is invariant over the 3 square-degree field, we are able to derive a line-of-sight reddening map with sub-arcminute resolution, enabling us to obtain de-reddened and extinction corrected color-magnitude diagrams (CMD's) of this bulge field using up to 2.5 million well-measured stars. The corrected CMD's show unprecedented detail and expose sparsely populated sequences: e.g., delineation of the very wide red giant branch, structure within the red giant clump, the full extent of the horizontal branch, and a surprising bright feature which is likely due to stars with ages younger than 1 Gyr. We use the RR Lyrae stars to trace the spatial structure of the ancient stars, and find an exponential decline in density with Galactocentric distance. We discuss ways in which our data products can be used to explore the age and metallicity properties of the bulge, and how our larger list of all variables is useful for learning to interpret future LSST alerts.

We explore the radial distribution of star formation in galaxies in the SAMI Galaxy Survey as a function of their local group environment. Using a sample of galaxies in groups (with halo masses less than $ \simeq 10^{14} \, \mathrm{M_{\odot}}$) from the Galaxy And Mass Assembly Survey, we find signatures of environmental quenching in high-mass groups ($M_{G} > 10^{12.5} \, \mathrm{M_{\odot}}$). The mean integrated specific star formation rate of star-forming galaxies in high-mass groups is lower than for galaxies in low-mass groups or that are ungrouped, with $\Delta \log(sSFR/\mathrm{yr^{-1}}) = 0.45 \pm 0.07$. This difference is seen at all galaxy stellar masses. In high-mass groups, star-forming galaxies more massive than $M_{*} \sim 10^{10} \, \mathrm{M_{\odot}}$ have centrally-concentrated star formation. These galaxies also lie below the star-formation main sequence, suggesting they may be undergoing outside-in quenching. Lower mass galaxies in high-mass groups do not show evidence of concentrated star formation. In groups less massive than $M_{G} = 10^{12.5} \, \mathrm{M_{\odot}}$ we do not observe these trends. In this regime we find a modest correlation between centrally-concentrated star formation and an enhancement in total star formation rate, consistent with triggered star formation in these galaxies.

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

Misalignment of gas and stellar rotation in galaxies can give clues to the origin and processing of accreted gas. Integral field spectroscopic observations of 1213 galaxies from the SAMI Galaxy Survey show that 11% of galaxies with fitted gas and stellar rotation are misaligned by more than 30 degrees in both field/group and cluster environments. Using SAMI morphological classifications and Sersic indices, the misalignment fraction is 45+/-6% in early-type galaxies, but only 5+/-1% in late-type galaxies. The distribution of position angle offsets is used to test the physical drivers of this difference. Slower dynamical settling time of the gas in elliptical stellar mass distributions accounts for a small increase in misalignment in early-type galaxies. However, gravitational dynamical settling time is insufficient to fully explain the observed differences between early- and late-type galaxies in the distributions of the gas/stellar position angle offsets. LTGs have primarily accreted gas close to aligned rather than settled from misaligned based on analysis of the skewed distribution of PA offsets compared to a dynamical settling model. Local environment density is less important in setting the misalignment fractions than morphology, suggesting that mergers are not the main source of accreted gas in these disks. Cluster environments are found to have gas misalignment driven primarily by cluster processes not by gas accretion.

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

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

We present year-long, near-infrared Hubble Space Telescope WFC3 observations of Mira variables in the water megamaser host galaxy NGC 4258. Miras are AGB variables that can be divided into oxygen- (O-) and carbon- (C-) rich subclasses. Oxygen-rich Miras follow a tight (scatter $\sim 0.14$ mag) Period-Luminosity Relation (PLR) in the near-infrared and can be used to measure extragalactic distances. The water megamaser in NGC 4258 gives a geometric distance to the galaxy accurate to 2.6% that can serve to calibrate the Mira PLR. We develop criteria for detecting and classifying O-rich Miras with optical and NIR data as well as NIR data alone. In total, we discover 438 Mira candidates that we classify with high confidence as O-rich. Our most stringent criteria produce a sample of 139 Mira candidates that we use to measure a PLR. We use the OGLE-III sample of O-rich Miras in the LMC to obtain a relative distance modulus, $\mu_{4258} - \mu_{LMC} = 10.95 \pm 0.01 $ (statistical) $\pm 0.06 $ (systematic) mag in good agreement with the relative distance determined using Cepheids. These results demonstrate the feasibility of discovering and characterizing Miras using the near-infrared with the Hubble Space Telescope and the upcoming James Webb Space Telescope and using them to measure extragalactic distances and determine the Hubble constant.

Most Type I superluminous supernovae (SLSNe-I) reported to date have been identified by their high peak luminosities and spectra lacking obvious signs of hydrogen. We demonstrate that these events can be distinguished from normal-luminosity SNe (including Type Ic events) solely from their spectra over a wide range of light-curve phases. We use this distinction to select 19 SLSNe-I and 4 possible SLSNe-I from the Palomar Transient Factory archive (including 7 previously published objects). We present 127 new spectra of these objects and combine these with 39 previously published spectra, and we use these to discuss the average spectral properties of SLSNe-I at different spectral phases. We find that Mn II most probably contributes to the ultraviolet spectral features after maximum light, and we give a detailed study of the O II features that often characterize the early-time optical spectra of SLSNe-I. We discuss the velocity distribution of O II, finding that for some SLSNe-I this can be confined to a narrow range compared to relatively large systematic velocity shifts. Mg II and Fe II favor higher velocities than O II and C II, and we briefly discuss how this may constrain power-source models. We tentatively group objects by how well they match either SN 2011ke or PTF12dam and discuss the possibility that physically distinct events may have been previously grouped together under the SLSN-I label.

In order to determine the causes of kinematic asymmetry in the H$\alpha$ gas in the SAMI Galaxy Survey sample, we investigate the comparative influences of environment and intrinsic properties of galaxies on perturbation. We use spatially resolved H$\alpha$ velocity fields from the SAMI Galaxy Survey to quantify kinematic asymmetry ($\overline{v_{asym}}$) in nearby galaxies and environmental and stellar mass data from the GAMA survey. {We find that local environment, measured as distance to nearest neighbour, is inversely correlated with kinematic asymmetry for galaxies with $\mathrm{\log(M_*/M_\odot)}>10.0$, but there is no significant correlation for galaxies with $\mathrm{\log(M_*/M_\odot)}<10.0$. Moreover, low mass galaxies ($\mathrm{\log(M_*/M_\odot)}<9.0$) have greater kinematic asymmetry at all separations, suggesting a different physical source of asymmetry is important in low mass galaxies.} We propose that secular effects derived from gas fraction and gas mass may be the primary causes of asymmetry in low mass galaxies. High gas fraction is linked to high $\frac{\sigma_{m}}{V}$ (where $\sigma_m$ is H$\alpha$ velocity dispersion and $V$ the rotation velocity), which is strongly correlated with $\overline{v_{asym}}$, and galaxies with $\log(M_*/M_\odot)<9.0$ have offset $\overline{\frac{\sigma_{m}}{V}}$ from the rest of the sample. Further, asymmetry as a fraction of dispersion decreases for galaxies with $\log(M_*/M_\odot)<9.0$. Gas mass and asymmetry are also inversely correlated in our sample. We propose that low gas masses in dwarf galaxies may lead to asymmetric distribution of gas clouds, leading to increased relative turbulence.

We present the ~800 star formation rate maps for the SAMI Galaxy Survey based on H{\alpha} emission maps, corrected for dust attenuation via the Balmer decrement, that are included in the SAMI Public Data Release 1. We mask out spaxels contaminated by non-stellar emission using the [O III]/H{\beta}, [N II]/H{\alpha}, [S II]/H{\alpha}, and [O I]/H{\alpha} line ratios. Using these maps, we examine the global and resolved star-forming main sequences of SAMI galaxies as a function of morphology, environmental density, and stellar mass. Galaxies further below the star-forming main sequence are more likely to have flatter star formation profiles. Early-type galaxies split into two populations with similar stellar masses and central stellar mass surface densities. The main sequence population has centrally-concentrated star formation similar to late-type galaxies, while galaxies >3{\sigma} below the main sequence show significantly reduced star formation most strikingly in the nuclear regions. The split populations support a two-step quenching mechanism, wherein halo mass first cuts off the gas supply and remaining gas continues to form stars until the local stellar mass surface density can stabilize the reduced remaining fuel against further star formation. Across all morphologies, galaxies in denser environments show a decreased specific star formation rate from the outside in, supporting an environmental cause for quenching, such as ram-pressure stripping or galaxy interactions.

We investigate the Tully-Fisher Relation (TFR) for a morphologically and kine- matically diverse sample of galaxies from the SAMI Galaxy Survey using 2 dimensional spatially resolved Halpha velocity maps and find a well defined relation across the stellar mass range of 8.0 < log(M*) < 11.5. We use an adaptation of kinemetry to parametrise the kinematic Halpha asymmetry of all galaxies in the sample, and find a correlation between scatter (i.e. residuals off the TFR) and asymmetry. This effect is pronounced at low stellar mass, corresponding to the inverse relationship between stellar mass and kinematic asymmetry found in previous work. For galaxies with log(M*) < 9.5, 25 +/- 3% are scattered below the root mean square (RMS) of the TFR, whereas for galaxies with log(M*) > 9.5 the fraction is 10 +/- 1% We use 'simulated slits' to directly compare our results with those from long slit spectroscopy and find that aligning slits with the photometric, rather than the kinematic, position angle, increases global scatter below the TFR. Further, kinematic asymmetry is correlated with misalignment between the photometric and kinematic position angles. This work demonstrates the value of 2D spatially resolved kinematics for accurate TFR studies; integral field spectroscopy reduces the underestimation of rotation velocity that can occur from slit positioning off the kinematic axis.

Astronomical surveys of celestial sources produce streams of noisy time series measuring flux versus time ("light curves"). Unlike in many other physical domains, however, large (and source-specific) temporal gaps in data arise naturally due to intranight cadence choices as well as diurnal and seasonal constraints. With nightly observations of millions of variable stars and transients from upcoming surveys, efficient and accurate discovery and classification techniques on noisy, irregularly sampled data must be employed with minimal human-in-the-loop involvement. Machine learning for inference tasks on such data traditionally requires the laborious hand-coding of domain-specific numerical summaries of raw data ("features"). Here we present a novel unsupervised autoencoding recurrent neural network (RNN) that makes explicit use of sampling times and known heteroskedastic noise properties. When trained on optical variable star catalogs, this network produces supervised classification models that rival other best-in-class approaches. We find that autoencoded features learned on one time-domain survey perform nearly as well when applied to another survey. These networks can continue to learn from new unlabeled observations and may be used in other unsupervised tasks such as forecasting and anomaly detection.

Merging neutron stars offer an exquisite laboratory for simultaneously studying strong-field gravity and matter in extreme environments. We establish the physical association of an electromagnetic counterpart EM170817 to gravitational waves (GW170817) detected from merging neutron stars. By synthesizing a panchromatic dataset, we demonstrate that merging neutron stars are a long-sought production site forging heavy elements by r-process nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to classical short gamma-ray bursts with ultra-relativistic jets. Instead, we suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the jet elegantly explains the low-luminosity gamma-rays, the high-luminosity ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit that all merging neutron stars may lead to a wide-angle cocoon breakout; sometimes accompanied by a successful jet and sometimes a choked jet.

We observed the afterglow of GRB 130427A with the RATIR instrument on the 1.5-m Harold L. Johnson telescope of the Observatorio Astron\'omico Nacional on Sierra San Pedro M\'artir. Our homogenous $griZYJH$ photometry extends from the night of burst to three years later. We fit a model for the afterglow. There is a significant positive residual which matches the behavior of SN 1998bw in the $griZ$ filters; we suggest that this is a photometric signature of the supernova SN 2013cq associated with the GRB. The peak absolute magnitude of the supernova is $M_r=-18.43\pm0.11$.

We present the results of a Hubble Space Telescope WFC3/F160W SNAPSHOT sur- vey of the host galaxies of 39 long-duration gamma-ray bursts (LGRBs) at z < 3. We have non-detections of hosts at the locations of 4 bursts. Sufficient accuracy to as- trometrically align optical afterglow images and determine the location of the LGRB within its host was possible for 31/35 detected hosts. In agreement with other work, we find the luminosity distribution of LGRB hosts is significantly fainter than that of a star formation rate-weighted field galaxy sample over the same redshift range, indicating LGRBs are not unbiasedly tracing the star formation rate. Morphologi- cally, the sample of LGRB hosts are dominated by spiral-like or irregular galaxies. We find evidence for evolution of the population of LGRB hosts towards lower-luminosity, higher concentrated hosts at lower redshifts. Their half-light radii are consistent with other LGRB host samples where measurements were made on rest-frame UV obser- vations. In agreement with recent work, we find their 80 per cent enclosed flux radii distribution to be more extended than previously thought, making them intermedi- ate between core-collapse supernova (CCSN) and super-luminous supernova (SLSN) hosts. The galactocentric projected-offset distribution confirms LGRBs as centrally concentrated, much more so than CCSNe and similar to SLSNe. LGRBs are strongly biased towards the brighter regions in their host light distributions, regardless of their offset. We find a correlation between the luminosity of the LGRB explosion site and the intrinsic column density, N_H , towards the burst.

We report on the first open-use based Atacama Large Millimeter/submm Array (ALMA) 345-GHz observation for the late afterglow phase of GRB131030A. The ALMA observation constrained a deep limit at 17.1 d for the afterglow and host galaxy. We also identified a faint submillimeter source (ALMAJ2300-0522) near the GRB131030A position. The deep limit at 345 GHz and multifrequency observations obtained using {\it Swift} and RATIR yielded forward shock modeling with a two-dimensional relativistic hydrodynamic jet simulation and described X-ray excess in the afterglow. The excess was inconsistent with the synchrotron self-inverse Compton radiation from the forward shock. The host galaxy of GRB131030A and optical counterpart of ALMAJ2300-0522 were also identified in the SUBARU image. Based on the deep ALMA limit for the host galaxy, the 3-$\sigma$ upper limits of IR luminosity and the star formation rate (SFR) is estimated as $L_{IR}<1.11\times10^{11} L_{\odot}$ and SFR$<18.7$ ($M_{\odot}$~yr$^{-1}$), respectively. Although the separation angle from the burst location (3.5 arcsec) was rather large, ALMAJ2300-0522 may be one component of the GRB131030A host galaxy, according to previous host galaxy cases.

We study the properties of kinematically disturbed galaxies in the SAMI Galaxy Survey using a quantitative criterion, based on kinemetry (Krajnovic et al.). The approach, similar to the application of kinemetry by Shapiro et al. uses ionised gas kinematics, probed by H{\alpha} emission. By this method 23+/-7% of our 360-galaxy sub-sample of the SAMI Galaxy Survey are kinematically asymmetric. Visual classifications agree with our kinemetric results for 90% of asymmetric and 95% of normal galaxies. We find stellar mass and kinematic asymmetry are inversely correlated and that kinematic asymmetry is both more frequent and stronger in low-mass galaxies. This builds on previous studies that found high fractions of kinematic asymmetry in low mass galaxies using a variety of different methods. Concentration of star forma- tion and kinematic disturbance are found to be correlated, confirming results found in previous work. This effect is stronger for high mass galaxies (log(M*) > 10) and indicates that kinematic disturbance is linked to centrally concentrated star formation. Comparison of the inner (within 0.5Re) and outer H{\alpha} equivalent widths of asymmetric and normal galaxies shows a small but significant increase in inner equivalent width for asymmetric galaxies.

Galaxy mergers are important events that can determine the fate of a galaxy by changing its morphology, star-formation activity and mass growth. Merger systems have commonly been identified from their disturbed morphologies, and we now can employ Integral Field Spectroscopy to detect and analyze the impact of mergers on stellar kinematics as well. We visually classified galaxy morphology using deep images ($\mu_{\rm r} = 28\,\rm mag\,\, arcsec^{-2}$) taken by the Blanco 4-m telescope at the Cerro Tololo Inter-American Observatory. In this paper we investigate 63 bright ($M_{\rm r}<-19.3$) spectroscopically-selected galaxies in Abell 119; of which 53 are early type and 20 galaxies show a disturbed morphology by visual inspection. A misalignment between the major axes in the photometric image and the kinematic map is conspicuous in morphologically-disturbed galaxies. Our sample is dominated by early-type galaxies, yet it shows a surprisingly tight Tully-Fisher relation except for the morphologically-disturbed galaxies which show large deviations. Three out of the eight slow rotators in our sample are morphology disturbed. The visually-selected morphologically-disturbed galaxies are generally more asymmetric, visually as well as kinematically. Our findings suggest that galaxy interactions, including mergers and perhaps fly-bys, play an important role in determining the orientation and magnitude of galaxy's angular momentum.

We use data from the Sydney-AAO Multi-Object Integral Field Spectrograph (SAMI) Galaxy Survey and the Galaxy And Mass Assembly (GAMA) survey to investigate the spatially-resolved signatures of the environmental quenching of star formation in galaxies. Using dust-corrected measurements of the distribution of H$\alpha$ emission we measure the radial profiles of star formation in a sample of 201 star-forming galaxies covering three orders of magnitude in stellar mass (M$_{*}$; $10^{8.1}$-$10^{10.95}\, $M$_{\odot}$) and in $5^{th}$ nearest neighbour local environment density ($\Sigma_{5}$; $10^{-1.3}$-$10^{2.1}\,$Mpc$^{-2}$). We show that star formation rate gradients in galaxies are steeper in dense ($\log_{10}(\Sigma_{5}/$Mpc$^{2})>0.5$) environments by $0.58\pm 0.29\, dex\, $r$_{e}^{-1}$ in galaxies with stellar masses in the range $10^{10}<$M$_{*}/$M$_{\odot}<10^{11}$ and that this steepening is accompanied by a reduction in the integrated star formation rate. However, for any given stellar mass or environment density the star-formation morphology of galaxies shows large scatter. We also measure the degree to which the star formation is centrally concentrated using the unitless scale-radius ratio ($r_{50,H\alpha}/r_{50,cont}$), which compares the extent of ongoing star formation to previous star formation. With this metric we find that the fraction of galaxies with centrally concentrated star formation increases with environment density, from $\sim 5\pm 4\%$ in low-density environments ($\log_{10}(\Sigma_{5}/$Mpc$^{2})<0.0$) to $30\pm 15\%$ in the highest density environments ($\log_{10}(\Sigma_{5}/$Mpc$^{2})>1.0$). These lines of evidence strongly suggest that with increasing local environment density the star formation in galaxies is suppressed, and that this starts in their outskirts such that quenching occurs in an outside-in fashion in dense environments and is not instantaneous.

We investigate the relationship between stellar and gas specific angular momentum $j$, stellar mass $M_{*}$ and optical morphology for a sample of 488 galaxies extracted from the SAMI Galaxy Survey. We find that $j$, measured within one effective radius, monotonically increases with $M_{*}$ and that, for $M_{*}>$10$^{9.5}$ M$_{\odot}$, the scatter in this relation strongly correlates with optical morphology (i.e., visual classification and S\'ersic index). These findings confirm that massive galaxies of all types lie on a plane relating mass, angular momentum and stellar light distribution, and suggest that the large-scale morphology of a galaxy is regulated by its mass and dynamical state. We show that the significant scatter in the $M_{*}-j$ relation is accounted for by the fact that, at fixed stellar mass, the contribution of ordered motions to the dynamical support of galaxies varies by at least a factor of three. Indeed, the stellar spin parameter (quantified via $\lambda_R$) correlates strongly with S\'ersic and concentration indices. This correlation is particularly strong once slow-rotators are removed from the sample, showing that late-type galaxies and early-type fast rotators form a continuous class of objects in terms of their kinematic properties.

We present optical and near-infrared light curves and optical spectra of SN 2013dx, associated with the nearby (redshift 0.145) gamma-ray burst GRB 130702A. The prompt isotropic gamma-ray energy released from GRB 130702A is measured to be $E_{\gamma,iso}=6.4_{-1.0}^{+1.3}\times10^{50}$erg (1keV-10MeV in the rest frame), placing it intermediate between low-luminosity GRBs like GRB 980425/SN 1998bw and the broader cosmological population. We compare the observed $g'r'i'z'$ light curves of SN 2013dx to a SN 1998bw template, finding that SN 2013dx evolves ~20% faster (steeper rise time), with a comparable peak luminosity. Spectroscopically, SN 2013dx resembles other broad-lined Type Ic supernovae, both associated with (SN 2006aj and SN 1998bw) and lacking (SN 1997ef, SN 2007I, and SN 2010ah) gamma-ray emission, with photospheric velocities around peak of ~21,000km s$^{-1}$. We construct a quasi-bolometric ($g'r'i'z'yJ$) light curve for SN 2013dx, only the fifth GRB-associated SN with extensive NIR coverage and the third with a bolometric light curve extending beyond $\Delta t>40$d. Together with the measured photospheric velocity, we derive basic explosion parameters using simple analytic models. We infer a $^{56}$Ni mass of $M_{\mathrm{Ni}}=0.37\pm0.01$M$_{\odot}$, an ejecta mass of $M_{\mathrm{ej}}=3.1\pm0.1$M$_{\odot}$, and a kinetic energy of $E_{\mathrm{K}}=(8.2\pm0.43)\times10^{51}$ erg (statistical uncertainties only), consistent with previous GRB-associated SNe. When considering the ensemble population of GRB-associated SNe, we find no correlation between the mass of synthesized $^{56}$Ni and high-energy properties, despite clear predictions from numerical simulations that $M_{\mathrm{Ni}}$ should correlate with the degree of asymmetry. On the other hand, $M_{\mathrm{Ni}}$ clearly correlates with the kinetic energy of the supernova ejecta across a wide range of core-collapse events.

We present our extensive observational campaign on the Swift-discovered GRB141121A, al- most ten years after its launch. Our observations covers radio through X-rays, and extends for more than 30 days after discovery. The prompt phase of GRB 141121A lasted 1410 s and, at the derived redshift of z = 1.469, the isotropic energy is E{\gamma},iso = 8.0x10^52 erg. Due to the long prompt duration, GRB141121A falls into the recently discovered class of UL-GRBs. Peculiar features of this burst are a flat early-time optical light curve and a radio-to-X-ray rebrightening around 3 days after the burst. The latter is followed by a steep optical-to-X-ray decay and a much shallower radio fading. We analyze GRB 141121A in the context of the standard forward-reverse shock (FS,RS) scenario and we disentangle the FS and RS contributions. Finally, we comment on the puzzling early-time (t ~3 d) behavior of GRB 141121A, and suggest that its interpretation may require a two-component jet model. Overall, our analysis confirms that the class of UL-GRBs represents our best opportunity to firmly establish the prominent emission mechanisms in action during powerful GRB explosions, and future missions (like SVOM, XTiDE, or ISS-Lobster) will provide many more of such objects.

Gamma-ray bursts (GRBs) are the most luminous explosions in the universe, yet the nature and physical properties of their energy sources are far from understood. Very important clues, however, can be inferred by studying the afterglows of these events. We present optical and X-ray observations of GRB 130831A obtained by Swift, Chandra, Skynet, RATIR, Maidanak, ISON, NOT, LT and GTC. This burst shows a steep drop in the X-ray light-curve at $\simeq 10^5$ s after the trigger, with a power-law decay index of $\alpha \sim 6$. Such a rare behaviour cannot be explained by the standard forward shock (FS) model and indicates that the emission, up to the fast decay at $10^5$ s, must be of "internal origin", produced by a dissipation process within an ultrarelativistic outflow. We propose that the source of such an outflow, which must produce the X-ray flux for $\simeq 1$ day in the cosmological rest frame, is a newly born magnetar or black hole. After the drop, the faint X-ray afterglow continues with a much shallower decay. The optical emission, on the other hand, shows no break across the X-ray steep decrease, and the late-time decays of both the X-ray and optical are consistent. Using both the X-ray and optical data, we show that the emission after $\simeq 10^5$ s can be explained well by the FS model. We model our data to derive the kinetic energy of the ejecta and thus measure the efficiency of the central engine of a GRB with emission of internal origin visible for a long time. Furthermore, we break down the energy budget of this GRB into the prompt emission, the late internal dissipation, the kinetic energy of the relativistic ejecta, and compare it with the energy of the associated supernova, SN 2013fu.

We present late-time multi-wavelength observations of Swift J1644+57, suggested to be a relativistic tidal disruption flare (TDF). Our observations extend to >4 years from discovery, and show that 1.4 years after outburst the relativistic jet switched-off on a timescale less than tens of days, corresponding to a power-law decay faster than $t^{-70}$. Beyond this point weak X-rays continue to be detected at an approximately constant luminosity of $L_X \sim 5 \times 10^{42}$ erg s$^{-1}$, and are marginally inconsistent with a continuing decay of $t^{-5/3}$, similar to that seen prior to the switch-off. Host photometry enables us to infer a black hole mass of $M_{BH}=3 \times 10^6$ M$_{\odot}$, consistent with the late time X-ray luminosity arising from sub-Eddington accretion onto the black hole in the form of either an unusually optically faint AGN or a slowly varying phase of the transient. Optical/IR observations show a clear bump in the light curve at timescales of 30-50 days, with a peak magnitude (corrected for host galaxy extinction) of $M_R \sim -22-23$. The luminosity of the bump is significantly higher than seen in other, non-relativistic TDFs and does not match any re-brightening seen at X-ray or radio wavelengths. Its luminosity, light curve shape and spectrum are broadly similar to those seen in superluminous SNe, although subject to large uncertainties in the correction of the significant host extinction. We discuss these observations in the context of both TDF and massive star origins for Swift J1644+5734 and other candidate relativistic tidal flares.

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.

The IceCube neutrino observatory pursues a follow-up program selecting interesting neutrino events in real-time and issuing alerts for electromagnetic follow-up observations. In March 2012, the most significant neutrino alert during the first three years of operation was issued by IceCube. In the follow-up observations performed by the Palomar Transient Factory (PTF), a Type IIn supernova (SN) PTF12csy was found $0.2^\circ$ away from the neutrino alert direction, with an error radius of $0.54^\circ$. It has a redshift of $z=0.0684$, corresponding to a luminosity distance of about $300 \, \mathrm{Mpc}$ and the Pan-STARRS1 survey shows that its explosion time was at least 158 days (in host galaxy rest frame) before the neutrino alert, so that a causal connection is unlikely. The a posteriori significance of the chance detection of both the neutrinos and the SN at any epoch is $2.2 \, \sigma$ within IceCube's 2011/12 data acquisition season. Also, a complementary neutrino analysis reveals no long-term signal over the course of one year. Therefore, we consider the SN detection coincidental and the neutrinos uncorrelated to the SN. However, the SN is unusual and interesting by itself: It is luminous and energetic, bearing strong resemblance to the SN IIn 2010jl, and shows signs of interaction of the SN ejecta with a dense circumstellar medium. High-energy neutrino emission is expected in models of diffusive shock acceleration, but at a low, non-detectable level for this specific SN. In this paper, we describe the SN PTF12csy and present both the neutrino and electromagnetic data, as well as their analysis.

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.

We present near-infrared (NIR) time-series spectroscopy, as well as complementary ultraviolet (UV), optical, and NIR data, of the Type Ia supernova (SN Ia) iPTF13ebh, which was discovered within two days from the estimated time of explosion. The first NIR spectrum was taken merely 2.3 days after explosion and may be the earliest NIR spectrum yet obtained of a SN Ia. The most striking features in the spectrum are several NIR C I lines, and the C I {\lambda}1.0693 {\mu}m line is the strongest ever observed in a SN Ia. Interestingly, no strong optical C II counterparts were found, even though the optical spectroscopic time series began early and is densely-cadenced. Except at the very early epochs, within a few days from the time of explosion, we show that the strong NIR C I compared to the weaker optical C II appears to be general in SNe Ia. iPTF13ebh is a fast decliner with {\Delta}m15(B) = 1.79 $\pm$ 0.01, and its absolute magnitude obeys the linear part of the width-luminosity relation. It is therefore categorized as a "transitional" event, on the fast-declining end of normal SNe Ia as opposed to subluminous/91bg-like objects. iPTF13ebh shows NIR spectroscopic properties that are distinct from both the normal and subluminous/91bg-like classes, bridging the observed characteristics of the two classes. These NIR observations suggest composition and density of the inner core similar to that of 91bg-like events, and a deep reaching carbon burning layer not observed in slower declining SNe Ia. There is also a substantial difference between the explosion times inferred from the early-time light curve and the velocity evolution of the Si II {\lambda}0.6355 {\mu}m line, implying a long dark phase of ~ 4 days.

We present optical and near-infrared (NIR) photometry of 28 gamma-ray bursts (GRBs) detected by the \textit{Swift} satellite and rapidly observed by the Reionization and Transients Infrared/Optical (RATIR) camera. We compare the optical flux at fiducial times of 5.5 and 11 hours after the high-energy trigger to that in the X-ray regime to quantify optical darkness. 46$\pm$9 per cent (13/28) of all bursts in our sample and 55$\pm$10 per cent (13/26) of long GRBs are optically dark, which is statistically consistently with previous studies. Fitting RATIR optical and NIR spectral energy distributions (SEDs) of 19 GRBs, most (6/7) optically dark GRBs either occur at high-redshift ($z>4.5$) or have a high dust content in their host galaxies ($A_{\rm V} > 0.3$). Performing K-S tests, we compare the RATIR sample to those previously presented in the literature, finding our distributions of redshift, optical darkness, host dust extinction and X-ray derived column density to be consistent. The one reported discrepancy is with host galaxy dust content in the BAT6 sample, which appears inconsistent with our sample and other previous literature. Comparing X-ray derived host galaxy hydrogen column densities to host galaxy dust extinction, we find that GRBs tend to occur in host galaxies with a higher metal-to-dust ratio than our own Galaxy, more akin to the Large and Small Magellanic Clouds. Finally, to mitigate time evolution of optical darkness, we measure $\beta_{\rm OX,rest}$ at a fixed rest frame time, $t_{\rm rest}=1.5$ hours and fixed rest frame energies in the X-ray and optical regimes. Choosing to evaluate optical flux at $\lambda_{\rm rest}=0.25~\mu$m, we remove high-redshift as a source of optical darkness, demonstrating that optical darkness must result from either high-redshift, dust content in the host galaxy along the GRB sight line, or a combination of the two.

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.

The SAMI Galaxy Survey will observe 3400 galaxies with the Sydney-AAO Multi-object Integral-field spectrograph (SAMI) on the Anglo-Australian Telescope (AAT) in a 3-year survey which began in 2013. We present the throughput of the SAMI system, the science basis and specifications for the target selection, the survey observation plan and the combined properties of the selected galaxies. The survey includes four volume limited galaxy samples based on cuts in a proxy for stellar mass, along with low-stellar mass dwarf galaxies all selected from the Galaxy And Mass Assembly (GAMA) survey. The GAMA regions were selected because of the vast array of ancillary data available, including ultraviolet through to radio bands. These fields are on the celestial equator at 9, 12, and 14.5 hours, and cover a total of 144 square degrees (in GAMA-I). Higher density environments are also included with the addition of eight clusters. The clusters have spectroscopy from 2dFGRS and SDSS and photometry in regions covered by the Sloan Digital Sky Survey (SDSS) and/or VLT Survey Telescope/ATLAS. The aim is to cover a broad range in stellar mass and environment, and therefore the primary survey targets cover redshifts 0.004 < z < 0.095, magnitudes r$_{pet}$ < 19.4, stellar masses $10^{7} - 10^{12}$ M$_{sol}$, and environments from isolated field galaxies through groups to clusters of $10^{15}$ M$_{sol}$.

We demonstrate the feasibility and potential of using large integral field spectroscopic surveys to investigate the prevalence of galactic-scale outflows in the local Universe. Using integral field data from SAMI and the Wide Field Spectrograph, we study the nature of an isolated disk galaxy, SDSS J090005.05+000446.7 (z = 0.05386). In the integral field datasets, the galaxy presents skewed line profiles changing with position in the galaxy. The skewed line profiles are caused by different kinematic components overlapping in the line-of-sight direction. We perform spectral decomposition to separate the line profiles in each spatial pixel as combinations of (1) a narrow kinematic component consistent with HII regions, (2) a broad kinematic component consistent with shock excitation, and (3) an intermediate component consistent with shock excitation and photoionisation mixing. The three kinematic components have distinctly different velocity fields, velocity dispersions, line ratios, and electron densities. We model the line ratios, velocity dispersions, and electron densities with our MAPPINGS IV shock and photoionisation models, and we reach remarkable agreement between the data and the models. The models demonstrate that the different emission line properties are caused by major galactic outflows that introduce shock excitation in addition to photoionisation by star-forming activities. Interstellar shocks embedded in the outflows shock-excite and compress the gas, causing the elevated line ratios, velocity dispersions, and electron densities observed in the broad kinematic component. We argue from energy considerations that, with the lack of a powerful active galactic nucleus, the outflows are likely to be driven by starburst activities. Our results set a benchmark of the type of analysis that can be achieved by the SAMI Galaxy Survey on large numbers of galaxies.

A fundamental challenge for wide-field imaging surveys is obtaining follow-up spectroscopic observations: there are > $10^9$ photometrically cataloged sources, yet modern spectroscopic surveys are limited to ~few x $10^6$ targets. As we approach the Large Synoptic Survey Telescope (LSST) era, new algorithmic solutions are required to cope with the data deluge. Here we report the development of a machine-learning framework capable of inferring fundamental stellar parameters (Teff, log g, and [Fe/H]) using photometric-brightness variations and color alone. A training set is constructed from a systematic spectroscopic survey of variables with Hectospec/MMT. In sum, the training set includes ~9000 spectra, for which stellar parameters are measured using the SEGUE Stellar Parameters Pipeline (SSPP). We employed the random forest algorithm to perform a non-parametric regression that predicts Teff, log g, and [Fe/H] from photometric time-domain observations. Our final, optimized model produces a cross-validated root-mean-square error (RMSE) of 165 K, 0.39 dex, and 0.33 dex for Teff, log g, and [Fe/H], respectively. Examining the subset of sources for which the SSPP measurements are most reliable, the RMSE reduces to 125 K, 0.37 dex, and 0.27 dex, respectively, comparable to what is achievable via low-resolution spectroscopy. For variable stars this represents a ~12-20% improvement in RMSE relative to models trained with single-epoch photometric colors. As an application of our method, we estimate stellar parameters for ~54,000 known variables. We argue that this method may convert photometric time-domain surveys into pseudo-spectrographic engines, enabling the construction of extremely detailed maps of the Milky Way, its structure, and history.

We present a methodology for the regularisation and combination of sparse sampled and irregularly gridded observations from fibre-optic multi-object integral-field spectroscopy. The approach minimises interpolation and retains image resolution on combining sub-pixel dithered data. We discuss the methodology in the context of the Sydney-AAO Multi-object Integral-field spectrograph (SAMI) Galaxy Survey underway at the Anglo-Australian Telescope. The SAMI instrument uses 13 fibre bundles to perform high-multiplex integral-field spectroscopy across a one degree diameter field of view. The SAMI Galaxy Survey is targeting 3000 galaxies drawn from the full range of galaxy environments. We demonstrate the subcritical sampling of the seeing and incomplete fill factor for the integral-field bundles results in only a 10% degradation in the final image resolution recovered. We also implement a new methodology for tracking covariance between elements of the resulting datacubes which retains 90% of the covariance information while incurring only a modest increase in the survey data volume.

We take advantage of the first data from the Sydney-AAO Multi-object Integral field (SAMI) Galaxy Survey to investigate the relation between the kinematics of gas and stars, and stellar mass in a comprehensive sample of nearby galaxies. We find that all 235 objects in our sample, regardless of their morphology, lie on a tight relation linking stellar mass ($M_{*}$) to internal velocity quantified by the $S_{0.5}$ parameter, which combines the contribution of both dispersion ($\sigma$) and rotational velocity ($V_{rot}$) to the dynamical support of a galaxy ($S_{0.5}=\sqrt{0.5V_{rot}^{2}+\sigma^{2}}$). Our results are independent of the baryonic component from which $\sigma$ and $V_{rot}$ are estimated, as the $S_{0.5}$ of stars and gas agree remarkably well. This represents a significant improvement compared to the canonical $M_{*}$ vs. $V_{rot}$ and $M_{*}$ vs. $\sigma$ relations. Not only is no sample pruning necessary, but also stellar and gas kinematics can be used simultaneously, as the effect of asymmetric drift is taken into account once $V_{rot}$ and $\sigma$ are combined. Our findings illustrate how the combination of dispersion and rotational velocities for both gas and stars can provide us with a single dynamical scaling relation valid for galaxies of all morphologies across at least the stellar mass range 8.5$<log(M_{*}/M_{\odot})<$11. Such relation appears to be more general and at least as tight as any other dynamical scaling relation, representing a unique tool for investigating the link between galaxy kinematics and baryonic content, and a less biased comparison with theoretical models.

We present the Early Data Release of the Sydney-AAO Multi-object Integral field spectrograph (SAMI) Galaxy Survey. The SAMI Galaxy Survey is an ongoing integral field spectroscopic survey of ~3400 low-redshift (z<0.12) galaxies, covering galaxies in the field and in groups within the Galaxy And Mass Assembly (GAMA) survey regions, and a sample of galaxies in clusters. In the Early Data Release, we publicly release the fully calibrated datacubes for a representative selection of 107 galaxies drawn from the GAMA regions, along with information about these galaxies from the GAMA catalogues. All datacubes for the Early Data Release galaxies can be downloaded individually or as a set from the SAMI Galaxy Survey website. In this paper we also assess the quality of the pipeline used to reduce the SAMI data, giving metrics that quantify its performance at all stages in processing the raw data into calibrated datacubes. The pipeline gives excellent results throughout, with typical sky subtraction residuals in the continuum of 0.9-1.2 per cent, a relative flux calibration uncertainty of 4.1 per cent (systematic) plus 4.3 per cent (statistical), and atmospheric dispersion removed with an accuracy of 0."09, less than a fifth of a spaxel.

We present the discovery of a luminous unresolved H II complex on the edge of dwarf galaxy GAMA J141103.98-003242.3 using data from the Sydney-AAO Multi-object Integral field spectrograph (SAMI) Galaxy Survey. This dwarf galaxy is situated at a distance of ~100 Mpc and contains an unresolved region of H II emission that contributes ~70 per cent of the galaxy's H_alpha luminosity, located at the top end of established H II region luminosity functions. For the H II complex, we measure a star-formation rate of 0.147\pm0.041 M_solar yr^-1 and a metallicity of 12+log(O/H) = 8.01\pm0.05 that is lower than the rest of the galaxy by ~0.2 dex. Data from the H I Parkes All-Sky Survey (HIPASS) indicate the likely presence of neutral hydrogen in the galaxy to potentially fuel ongoing and future star-forming events. We discuss various triggering mechanisms for the intense star-formation activity of this H II complex, where the kinematics of the ionised gas are well described by a rotating disc and do not show any features indicative of interactions. We show that SAMI is an ideal instrument to identify similar systems to GAMA J141103.98-003242.3, and the SAMI Galaxy Survey is likely to find many more of these systems to aid in the understanding of their formation and evolution.

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.

We determine the period-magnitude relations of RR Lyrae stars in 13 photometric bandpasses from 0.4 to 12 {\mu}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 {\mu}m) and W2 (4.6 {\mu}m) where the photometric error ($\sigma \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 $\mu=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.

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.

We present optical photometry and spectroscopy of the broad-lined Type Ic supernova (SNIc-BL) PTF10qts, which was discovered as part of the Palomar Transient Factory. The supernova was located in a dwarf galaxy of magnitude $r=21.1$ at a redshift $z=0.0907$. We find that the $R$-band light curve is a poor proxy for bolometric data and use photometric and spectroscopic data to construct and constrain the bolometric light curve. The derived bolometric magnitude at maximum light is $M_{\rm bol} = -18.51\pm0.2$ mag, comparable to that of SN 1998bw ($M_{\rm bol} = -18.7$ mag) which was associated with a gamma-ray burst (GRB). PTF10qts is one of the most luminous SNIc-BL observed without an accompanying GRB. We estimate the physical parameters of the explosion using data from our programme of follow-up observations, finding that it produced a larger mass of radioactive nickel compared to other SNeIc-BL with similar inferred ejecta masses and kinetic energies. The progenitor of the event was likely a $\sim20$M$_{\odot}$ star.

We present optical and near-infrared (NIR) photometry for three gamma-ray burst supernovae (GRB-SNe): GRB 120729A, GRB 130215A / SN 2013ez and GRB 130831A / SN 2013fu. In the case of GRB 130215A / SN 2013ez, we also present optical spectroscopy at t-t0=16.1 d, which covers rest-frame 3000-6250 Angstroms. Based on Fe II (5169) and Si (II) (6355), our spectrum indicates an unusually low expansion velocity of 4000-6350 km/s, the lowest ever measured for a GRB-SN. Additionally, we determined the brightness and shape of each accompanying SN relative to a template supernova (SN 1998bw), which were used to estimate the amount of nickel produced via nucleosynthesis during each explosion. We find that our derived nickel masses are typical of other GRB-SNe, and greater than those of SNe Ibc that are not associated with GRBs. For GRB 130831A / SN 2013fu, we use our well-sampled R-band light curve (LC) to estimate the amount of ejecta mass and the kinetic energy of the SN, finding that these too are similar to other GRB-SNe. For GRB 130215A, we take advantage of contemporaneous optical/NIR observations to construct an optical/NIR bolometric LC of the afterglow. We fit the bolometric LC with the millisecond magnetar model of Zhang & Meszaros (2001), which considers dipole radiation as a source of energy injection to the forward shock powering the optical/NIR afterglow. Using this model we derive an initial spin period of P=12 ms and a magnetic field of B=1.1 x 10^15 G, which are commensurate with those found for proposed magnetar central engines of other long-duration GRBs.

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 ${\mu_{\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.

We present a densely-sampled, homogeneous set of light curves of 64 low redshift (z < 0.05) stripped-envelope supernovae (SN of type IIb, Ib, Ic and Ic-bl). These data were obtained between 2001 and 2009 at the Fred L. Whipple Observatory (FLWO) on Mt. Hopkins in Arizona, with the optical FLWO 1.2-m and the near-infrared PAIRITEL 1.3-m telescopes. Our dataset consists of 4543 optical photometric measurements on 61 SN, including a combination of UBVRI, UBVr'i', and u'BVr'i', and 2142 JHKs near-infrared measurements on 25 SN. This sample constitutes the most extensive multi-color data set of stripped-envelope SN to date. Our photometry is based on template-subtracted images to eliminate any potential host galaxy light contamination. This work presents these photometric data, compares them with data in the literature, and estimates basic statistical quantities: date of maximum, color, and photometric properties. We identify promising color trends that may permit the identification of stripped-envelope SN subtypes from their photometry alone. Many of these SN were observed spectroscopically by the CfA SN group, and the spectra are presented in a companion paper (Modjaz et al. 2014). A thorough exploration that combines the CfA photometry and spectroscopy of stripped-envelope core-collapse SN will be presented in a follow-up paper.

We present a template fitting algorithm for determining photometric redshifts, $z_{\rm phot}$, of candidate high-redshift gamma-ray bursts (GRBs). Using afterglow photometry, obtained by the Reionization And Transients InfraRed (RATIR) camera, this algorithm accounts for the intrinsic GRB afterglow spectral energy distribution (SED), host dust extinction and the effect of neutral hydrogen (local and cosmological) along the line of sight. We present the results obtained by this algorithm and RATIR photometry of GRB 130606A, finding a range of best fit solutions $5.6 < z_{\rm phot} < 6.0$ for models of several host dust extinction laws (none, MW, LMC and SMC), consistent with spectroscopic measurements of the redshift of this GRB. Using simulated RATIR photometry, we find our algorithm provides precise measures of $z_{\rm phot}$ in the ranges $4 < z_{\rm phot} \lesssim 8$ and $9 < z_{\rm phot} < 10$ and can robustly determine when $z_{\rm phot}>4$. Further testing highlights the required caution in cases of highly dust extincted host galaxies. These tests also show that our algorithm does not erroneously find $z_{\rm phot} < 4$ when $z_{\rm sim}>4$, thereby minimizing false negatives and allowing us to rapidly identify all potential high-redshift events.

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.

The spiral host galaxy of GRB 060505 at z=0.089 was the site of a puzzling long duration burst without an accompanying supernova. Studies of the burst environment by Th\"one et al. (2008) suggested that this GRB came from the collapse of a massive star and that the GRB site was a region with properties different from the rest of the galaxy. We reobserved the galaxy in high spatial resolution using the VIMOS integral-field unit (IFU) at the VLT with a spaxel size of 0.67 arcsec. Furthermore, we use long slit high resolution data from HIRES/Keck at two different slit positions covering the GRB site, the center of the galaxy and an HII region next to the GRB region. We compare the properties of different HII regions in the galaxy with the GRB site and study the global and local kinematic properties of this galaxy. The resolved data show that the GRB site has the lowest metallicity in the galaxy with around 1/3 Z_solar, but its specific SFR (SSFR) of 7.4 M_solar/yr/L/L* and age (determined by the Halpha EW) are similar to other HII regions in the host. The galaxy shows a gradient in metallicity and SSFR from the bulge to the outskirts as it is common for spiral galaxies. This gives further support to the theory that GRBs prefer regions of higher star-formation and lower metallicity, which, in S-type galaxies, are more easily found in the spiral arms than in the centre. Kinematic measurements of the galaxy do not show evidence for large perturbations but a minor merger in the past cannot be excluded. This study confirms the collapsar origin of GRB060505 but reveals that the properties of the HII region surrounding the GRB were not unique to that galaxy. Spatially resolved observations are key to know the implications and interpretations of unresolved GRB hosts observations at higher redshifts.

We report on the discovery of SN2014J in the nearby galaxy M82. Given its proximity, it offers the best opportunity to date to study a thermonuclear supernova over a wide range of the electromagnetic spectrum. The first set of optical, near-IR and mid-IR observations of SN2014J, orchestrated by the intermediate Palomar Transient Factory (iPTF), show that SN2014J is a spectroscopically normal Type Ia supernova, albeit exhibiting high-velocity features in its spectrum and heavily reddened by dust in the host galaxy. Our earliest detections start just hours after the fitted time of explosion. We use high-resolution optical spectroscopy to analyze the dense intervening material and do not detect any evolution in the resolved absorption features during the lightcurve rise. Similarly to other highly reddened Type Ia supernovae, a low value of total-to-selective extinction, Rv < 2, provides the best match to our observations. We also study pre-explosion optical and near-IR images from HST with special emphasis on the sources nearest to the SN location.

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 {\mu}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.

A significant proportion ($\sim30\%$) of the short-duration gamma-ray bursts (SGRBs) localised by Swift have no detected host galaxy coincident with the burst location to deep limits, and also no high-likelihood association with proximate galaxies on the sky. These SGRBs may represent a population at moderately high redshifts ($z\gtrsim1$), for which the hosts are faint, or a population where the progenitor has been kicked far from its host or is sited in an outlying globular cluster. We consider the afterglow and host observations of three 'hostless' bursts (GRBs 090305A, 091109B and 111020A), coupled with a new observational diagnostic to aid the association of SGRBs with putative host galaxies to investigate this issue. Considering the well localised SGRB sample, 7/25 SGRBs can be classified as 'hostless' by our diagnostic. Statistically, however, the proximity of these seven SGRBs to nearby galaxies is higher than is seen for random positions on the sky. This suggests that the majority of 'hostless' SGRBs have likely been kicked from proximate galaxies at moderate redshift. Though this result still suggests only a small proportion of SGRBs will be within the AdLIGO horizon for NS-NS or NS-BH inspiral detection ($z\sim0.1$), in the particular case of GRB 111020A a plausible host candidate is at $z=0.02$.

We present the discovery, classification, and extensive panchromatic (from radio to X-ray) follow-up observations of PTF11qcj, a supernova discovered by the Palomar Transient Factory. PTF11qcj is located at a distance of dL ~ 124 Mpc. Our observations with the Karl G. Jansky Very Large Array show that this event is radio-loud: PTF11qcj reached a radio peak luminosity comparable to that of the famous gamma-ray-burst-associated supernova 1998bw (L_{5GHz} ~ 10^{29} erg/s/Hz). PTF11qcj is also detected in X-rays with the Chandra observatory, and in the infrared band with Spitzer. Our multi-wavelength analysis probes the supernova interaction with circumstellar material. The radio observations suggest a progenitor mass-loss rate of ~10^{-4} Msun/yr x (v_w/1000 km/s), and a velocity of ~(0.3-0.5)c for the fastest moving ejecta (at ~10d after explosion). However, these estimates are derived assuming the simplest model of supernova ejecta interacting with a smooth circumstellar material characterized by radial power-law density profile, and do not account for possible inhomogeneities in the medium and asphericity of the explosion. The radio light curve shows deviations from such a simple model, as well as a re-brightening at late times. The X-ray flux from PTF11qcj is compatible with the high-frequency extrapolation of the radio synchrotron emission (within the large uncertainties). An IR light echo from pre-existing dust is in agreement with our infrared data. Our analysis of pre-explosion data from the Palomar Transient Factory suggests that a precursor eruption of absolute magnitude M_r ~ -13 mag may have occurred ~ 2.5 yr prior to the supernova explosion. Based on our panchromatic follow-up campaign, we conclude that PTF11qcj fits the expectations from the explosion of a Wolf-Rayet star. Precursor eruptions may be a feature characterizing the final pre-explosion evolution of such stars.

There is a growing number of supernovae (SNe), mainly of Type IIn, which present an outburst prior to their presumably final explosion. These precursors may affect the SN display, and are likely related to some poorly charted phenomena in the final stages of stellar evolution. Here we present a sample of 16 SNe IIn for which we have Palomar Transient Factory (PTF) observations obtained prior to the SN explosion. By coadding these images taken prior to the explosion in time bins, we search for precursor events. We find five Type IIn SNe that likely have at least one possible precursor event, three of which are reported here for the first time. For each SN we calculate the control time. Based on this analysis we find that precursor events among SNe IIn are common: at the one-sided 99% confidence level, more than 50% of SNe IIn have at least one pre-explosion outburst that is brighter than absolute magnitude -14, taking place up to 1/3 yr prior to the SN explosion. The average rate of such precursor events during the year prior to the SN explosion is likely larger than one per year, and fainter precursors are possibly even more common. We also find possible correlations between the integrated luminosity of the precursor, and the SN total radiated energy, peak luminosity, and rise time. These correlations are expected if the precursors are mass-ejection events, and the early-time light curve of these SNe is powered by interaction of the SN shock and ejecta with optically thick circumstellar material.

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.

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.

We present observations and analysis of the host galaxies of 23 heavily dust-obscured gamma-ray bursts (GRBs) observed by the Swift satellite during the years 2005-2009, representing all GRBs with an unambiguous host-frame extinction of A_V>1 mag from this period. Deep observations with Keck, Gemini, VLT, HST, and Spitzer successfully detect the host galaxies and establish spectroscopic or photometric redshifts for all 23 events, enabling us to provide measurements of the intrinsic host star-formation rates, stellar masses, and mean extinctions. Compared to the hosts of unobscured GRBs at similar redshifts, we find that the hosts of dust-obscured GRBs are (on average) more massive by about an order of magnitude and also more rapidly star-forming and dust-obscured. While this demonstrates that GRBs populate all types of star-forming galaxies including the most massive, luminous systems at z~2, at redshifts below 1.5 the overall GRB population continues to show a highly significant aversion away from massive galaxies and a preference for low-mass systems relative to what would be expected given a purely SFR-selected galaxy sample. This supports the notion that the GRB rate is strongly dependent on metallicity, and may suggest that the most massive galaxies in the Universe underwent a transition in their chemical properties ~9 Gyr ago. We also conclude that, based on the absence of unobscured GRBs in massive galaxies and the absence of obscured GRBs in low-mass galaxies, the dust distributions of the lowest-mass and the highest-mass galaxies are relatively homogeneous, while intermediate-mass galaxies (~10^9 M_sun) have diverse internal properties.

During the LIGO and Virgo joint science runs in 2009-2010, gravitational wave (GW) data from three interferometer detectors were analyzed within minutes to select GW candidate events and infer their apparent sky positions. Target coordinates were transmitted to several telescopes for follow-up observations aimed at the detection of an associated optical transient. Images were obtained for eight such GW candidates. We present the methods used to analyze the image data as well as the transient search results. No optical transient was identified with a convincing association with any of these candidates, and none of the GW triggers showed strong evidence for being astrophysical in nature. We compare the sensitivities of these observations to several model light curves from possible sources of interest, and discuss prospects for future joint GW-optical observations of this type.

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.

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.

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.

The intermediate Palomar Transient Factory reports our discovery of a young supernova, iPTF13bvn, in the nearby galaxy, NGC5806 (22.5Mpc). Our spectral sequence in the optical and infrared suggests a likely Type Ib classification. We identify a single, blue progenitor candidate in deep pre-explosion imaging within a 2{\sigma} error circle of 80 mas (8.7 pc). The candidate has a MB luminosity of -5.2 +/- 0.4 mag and a B-I color of 0.1+/-0.3 mag. If confirmed by future observations, this would be the first direct detection for a progenitor of a Type Ib. Fitting a power law to the early light curve, we find an extrapolated explosion date around 1.1 days before our first detection. We see no evidence of shock cooling. The pre-explosion detection limits constrain the radius of the progenitor to be smaller than a few solar radii. iPTF13bvn is also detected in cm and mm-wavelengths. Fitting a synchrotron self-absorption model to our radio data, we find a mass loading parameter of 1.3*10^12 g/cm. Assuming a wind velocity of 10^3km/s, we derive a progenitor mass loss rate of 3*10^-5Msun/yr. Our observations, taken as a whole, are consistent with a Wolf Rayet progenitor of the supernova iPTF13bvn.

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.

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.

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.

We present early optical photometry and spectroscopy of the afterglow and host galaxy of the bright short-duration gamma-ray burst GRB 130603B discovered by the Swift satellite. Using the Gemini South telescope, our prompt optical spectra reveal a strong trace from the afterglow superimposed on continuum and emission lines from the $z = 0.3568 \pm 0.0005$ host galaxy. The combination of a relatively bright optical afterglow (r' = 21.52 at $\Delta_t $= 8.4hr), together with an observed offset of 0\farcs9 from the host nucleus (4.8kpc projected distance at z=0.3568), allow us to extract a relatively clean spectrum dominated by afterglow light . The spatially resolved spectrum allows us to constrain the properties of the explosion site directly, and compare these with the host galaxy nucleus, as well as other short-duration GRB host galaxies. We find that while the host is a relatively luminous ($L \approx 0.8 L^{*}_{B}$), star-forming galaxy with solar metallicity, the spectrum of the afterglow exhibits weak CaII absorption features but negligible emission features. The explosion site therefore lacks evidence of recent star formation, consistent with the relatively long delay time distribution expected in a compact binary merger scenario. The star formation rate (both in an absolute sense and normalized to the luminosity) and metallicity of the host are both consistent with the known sample of short-duration GRB hosts and with recent results which suggest GRB130603B emission to be the product of the decay of radioactive species produced during the merging process of a NS-NS binary ("kilonova"). Ultimately, the discovery of more events similar to GRB130603B and their rapid follow-up from 8-m class telescopes will open new opportunities for our understanding of the final stages of compact-objects binary systems.

(Abridged) Some supernovae (SNe) may be powered by the interaction of the SN ejecta with a large amount of circumstellar matter (CSM). Here we outline a method to measure the mass of the optically thick CSM around such SNe. We present observations of SN2010jl, including the first detection of a SN using NuSTAR. The total radiated luminosity of SN2010jl is extreme, at least 9e50 erg. By modeling the visible-light data, we robustly show that the mass of the circumstellar material within ~1e16 cm of the progenitor was in excess of 10 solar masses, likely ejected tens of years prior to the SN explosion. Our modeling suggests that the shock velocity during shock breakout was ~6000 km/s, decelerating to ~2600 km/s about two years after maximum light. Our late-time NuSTAR+XMM spectra of the SN presumably provide the first direct measurement of SN shock velocity two years after the SN maximum light -- measured to be in the range of 2000 to 4500 km/s if the ions and electrons are in equilibrium, and >~2000 km/s if they are not in equilibrium. This measurement is in agreement with the shock velocity predicted by our modeling of the optical data. We also show that the mean radial density distribution of the CSM roughly follows an r^-2 law. A possible explanation for the massive CSM with a wind-like profile is that they are the result of multiple pulsational pair instability events prior to the SN explosion, separated from each other by years.

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.

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

We present broadband observations and analysis of Swift gamma-ray burst (GRB) 120119A. Our early-time afterglow detections began under 15 s after the burst in the host frame (redshift z = 1.73), and they yield constraints on the burst energetics and local environment. Late-time afterglow observations of the burst show evidence for a moderate column of dust (A_V ~ 1.1 mag) similar to, but statistically distinct from, dust seen along Small Magellanic Cloud sightlines. Deep late-time observations reveal a dusty, rapidly star-forming host galaxy. Most notably, our early-time observations exhibit a significant red-to-blue colour change in the first ~200 s after the trigger at levels heretofore unseen in GRB afterglows. This colour change, which is coincident with the final phases of the prompt emission, is a hallmark prediction of the photodestruction of dust in GRB afterglows. We test whether dust-destruction signatures are significantly distinct from other sources of colour change, namely a change in the intrinsic spectral index {\beta}. We find that a time-varying power-law spectrum alone cannot adequately describe the observed colour change, and allowing for dust destruction (via a time-varying A_V) significantly improves the fit. While not definitively ruling out other possibilities, this event provides the best support yet for the direct detection of dust destruction in the local environment of a GRB.

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.

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{\alpha} emission (with widths of ~2000 km/s) and exhibit large H{\alpha}/H{\beta} 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{\alpha} 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.

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.

Various lines of evidence suggest that very massive stars experience extreme mass-loss episodes shortly before they explode as a supernova. Interestingly, several models predict such pre-explosion outbursts. Establishing a causal connection between these mass-loss episodes and the final supernova explosion will provide a novel way to study pre-supernova massive-star evolution. Here we report on observations of a remarkable mass-loss event detected 40 days prior to the explosion of the Type IIn supernova SN 2010mc (PTF 10tel). Our photometric and spectroscopic data suggest that this event is a result of an energetic outburst, radiating at least 6x10^47 erg of energy, and releasing about 0.01 Solar mass at typical velocities of 2000 km/s. We show that the temporal proximity of the mass-loss outburst and the supernova explosion implies a causal connection between them. Moreover, we find that the outburst luminosity and velocity are consistent with the predictions of the wave-driven pulsation model and disfavor alternative suggestions.

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.

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.

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.

In 2006, Prochter et al. reported a statistically significant enhancement of very strong Mg II absorption systems intervening the sightlines to gamma-ray bursts (GRBs) relative to the in- cidence of such absorption along quasar sightlines. This counterintuitive result, has inspired a diverse set of astrophysical explanations (e.g. dust, gravitational lensing) but none of these has obviously resolved the puzzle. Using the largest set of GRB afterglow spectra available, we reexamine the purported enhancement. In an independent sample of GRB spectra with a survey path 3 times larger than Prochter et al., we measure the incidence per unit redshift of $\geq 1$\AA rest-frame equivalent width Mg II absorbers at $z \approx 1$ to be l(z)= 0.18 $\pm$ 0.06. This is fully consistent with current estimates for the incidence of such absorbers along quasar sightlines. Therefore, we do not confirm the original enhancement and suggest those results suffered from a statistical fluke. Signatures of the original result do remain in our full sample (l(z) shows an $\approx 1.5$ enhancement over l(z)QSO), but the statistical significance now lies at $\approx 90%$ c.l. Restricting our analysis to the subset of high-resolution spectra of GRB afterglows (which overlaps substantially with Prochter et al.), we still reproduce a statistically significant enhancement of Mg II absorption. The reason for this excess, if real, is still unclear since there is no connection between the rapid afterglow follow-up process with echelle (or echellette) spectrographs and the detectability of strong Mg II doublets. Only a larger sample of such high-resolution data will shed some light on this matter.

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.

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.

We present near- and mid-infrared photometry and spectroscopy from PAIRITEL, IRTF, and Spitzer of a metallicity-unbiased sample of 117 cool, hydrogen-atmosphere white dwarfs from the Palomar-Green survey and find five with excess radiation in the infrared, translating to a 4.3+2.7-1.2% frequency of debris disks. This is slightly higher than, but consistent with the results of previous surveys. Using an initial-final mass relation, we apply this result to the progenitor stars of our sample and conclude that 1-7Msol stars have at least a 4.3% chance of hosting planets; an indirect probe of the intermediate-mass regime eluding conventional exoplanetary detection methods. Alternatively, we interpret this result as a limit on accretion timescales as a fraction of white dwarf cooling ages; white dwarfs accrete debris from several generations of disks for ~10Myr. The average total mass accreted by these stars ranges from that of 200km asteroids to Ceres-sized objects, indicating that white dwarfs accrete moons and dwarf planets as well as Solar System asteroid analogues.

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.

We present the results of a machine-learning (ML) based search for new R Coronae Borealis (RCB) stars and DY Persei-like stars (DYPers) in the Galaxy using cataloged light curves from the All-Sky Automated Survey (ASAS) Catalog of Variable Stars (ACVS). RCB stars - a rare class of hydrogen-deficient carbon-rich supergiants - are of great interest owing to the insights they can provide on the late stages of stellar evolution. DYPers are possibly the low-temperature, low-luminosity analogs to the RCB phenomenon, though additional examples are needed to fully establish this connection. While RCB stars and DYPers are traditionally identified by epochs of extreme dimming that occur without regularity, the ML search framework more fully captures the richness and diversity of their photometric behavior. We demonstrate that our ML method can use newly discovered RCB stars to identify additional candidates within the same data set. Our search yields 15 candidates that we consider likely RCB stars/DYPers: new spectroscopic observations confirm that four of these candidates are RCB stars and four are DYPers. Our discovery of four new DYPers increases the number of known Galactic DYPers from two to six; noteworthy is that one of the new DYPers has a measured parallax and is m ~ 7 mag, making it the brightest known DYPer to date. Future observations of these new DYPers should prove instrumental in establishing the RCB connection. We consider these results, derived from a machine-learned probabilistic classification catalog, as an important proof-of-concept for the efficient discovery of rare sources with time-domain surveys.

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.

There is a consensus that Type-Ia supernovae (SNe Ia) arise from the thermonuclear explosion of white dwarf stars that accrete matter from a binary companion. However, direct observation of SN Ia progenitors is lacking, and the precise nature of the binary companion remains uncertain. A temporal series of high-resolution optical spectra of the SN Ia PTF 11kx reveals a complex circumstellar environment that provides an unprecedentedly detailed view of the progenitor system. Multiple shells of circumsteller are detected and the SN ejecta are seen to interact with circumstellar material (CSM) starting 59 days after the explosion. These features are best described by a symbiotic nova progenitor, similar to RS Ophiuchi.

We construct a photometrically calibrated catalog of non-variable sources from the Palomar Transient Factory (PTF) observations. The first version of this catalog presented here, the PTF photometric catalog 1.0, contains calibrated R_PTF-filter magnitudes for about 21 million sources brighter than magnitude 19, over an area of about 11233 deg^2. The magnitudes are provided in the PTF photometric system, and the color of a source is required in order to convert these magnitudes into other magnitude systems. We estimate that the magnitudes in this catalog have typical accuracy of about 0.02 mag with respect to magnitudes from the Sloan Digital Sky Survey. The median repeatability of our catalog's magnitudes for stars between 15 and 16 mag, is about 0.01 mag, and it is better than 0.03 mag for 95% of the sources in this magnitude range. The main goal of this catalog is to provide reference magnitudes for photometric calibration of visible light observations. Subsequent versions of this catalog, which will be published incrementally online, will be extended to a larger sky area and will also include g_PTF-filter magnitudes, as well as variability and proper motion information.

(Abridged). The optical light curve of some SNe may be powered by the outward diffusion of the energy deposited by the explosion shock in optically thick circumstellar matter (CSM). Recently, it was shown that the radiation-mediated and -dominated shock in an optically thick wind must transform into a collisionless shock and can produce hard X-rays. The X-rays are expected to peak at late times, relative to maximum visible light. Here we report on a search, using Swift and Chandra, for X-ray emission from 28 SNe that belong to classes whose progenitors are suspected to be embedded in dense CSM (IIn/Ibn/SLSN-I). Two SNe in our sample have X-ray properties that are roughly consistent with the expectation for X-rays from a collisionless shock in optically thick CSM. Therefore, we suggest that their optical light curves are powered by shock breakout in CSM. We show that two other events were too X-ray bright during the SN maximum optical light to be explained by the shock breakout model. We conclude that the light curves of some, but not all, type-IIn/Ibn SNe are powered by shock breakout in CSM. For the rest of the SNe in our sample, including all the SLSN-I events, our X-ray limits are not deep enough and were typically obtained at too early times to conclude about their nature. We argue that the optical light curves of SNe, for which the X-ray emission peaks at late times, are likely powered by the diffusion of shock energy from a dense CSM. We comment about the possibility to detect some of these events in radio.

The nearby Type Ia supernova SN 2011fe in M101 (cz=241 km s^-1) provides a unique opportunity to study the early evolution of a "normal" Type Ia supernova, its compositional structure, and its elusive progenitor system. We present 18 high signal-to-noise spectra of SN 2011fe during its first month beginning 1.2 days post-explosion and with an average cadence of 1.8 days. This gives a clear picture of how various line-forming species are distributed within the outer layers of the ejecta, including that of unburned material (C+O). We follow the evolution of C II absorption features until they diminish near maximum light, showing overlapping regions of burned and unburned material between ejection velocities of 10,000 and 16,000 km s^-1. This supports the notion that incomplete burning, in addition to progenitor scenarios, is a relevant source of spectroscopic diversity among SNe Ia. The observed evolution of the highly Doppler-shifted O I 7774 absorption features detected within five days post-explosion indicate the presence of O I with expansion velocities from 11,500 to 21,000 km s^-1. The fact that some O I is present above C II suggests that SN 2011fe may have had an appreciable amount of unburned oxygen within the outer layers of the ejecta.

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

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]

The unequivocal, spectroscopic detection of the 2175 bump in extinction curves outside the Local Group is rare. To date, the properties of the bump have been examined in only two GRB afterglows (GRB 070802 and GRB 080607). In this work we analyse in detail the detections of the 2175 extinction bump in the optical spectra of the two further GRB afterglows: GRB 080605 and 080805. We gather all available optical/NIR photometric, spectroscopic and X-ray data to construct multi-epoch SEDs for both GRB afterglows. We fit the SEDs with the Fitzpatrick & Massa (1990) model with a single or broken PL. We also fit a sample of 38 GRB afterglows, known to prefer a SMC-type extinction curve, with the same model. We find that the SEDs of GRB 080605 and GRB 080805 at two epochs are fit well with a single PL with a derived extinction of A_V = 0.52(+0.13 -0.16) and 0.50 (+0.13 -0.10), and 2.1(+0.7-0.6) and 1.5+/-0.2 respectively. While the slope of the extinction curve of GRB 080805 is not well-constrained, the extinction curve of GRB 080605 has an unusual very steep far-UV rise together with the 2175 bump. Such an extinction curve has previously been found in only a small handful of sightlines in the MW. One possible explanation of such an extinction curve may be dust arising from two different regions with two separate grain populations, however we cannot distinguish the origin of the curve. We finally compare the four 2175 bump sightlines to the larger GRB afterglow sample and to Local Group sightlines. We find that while the width and central positions of the bumps are consistent with what is observed in the Local Group, the relative strength of the detected bump (A_bump) for GRB afterglows is weaker for a given A_V than for almost any Local Group sightline. Such dilution of the bump strength may offer tentative support to a dual dust-population scenario.

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.

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.

We present the discovery of PTF 10vgv, a Type Ic supernova detected by the Palomar Transient Factory, using the Palomar 48-inch telescope (P48). R-band observations of the PTF 10vgv field with P48 probe the supernova emission from its very early phases (about two weeks before R-band maximum), and set limits on its flux in the week prior to the discovery. Our sensitive upper limits and early detections constrain the post-shock-breakout luminosity of this event. Via comparison to numerical (analytical) models, we derive an upper-limit of R \lesssim 4.5 Rsun (R \lesssim 1 Rsun) on the radius of the progenitor star, a direct indication in favor of a compact Wolf-Rayet star. Applying a similar analysis to the historical observations of SN 1994I, yields R \lesssim 1/4 Rsun for the progenitor radius of this supernova.

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.

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.

The Palomar Transient Factory (PTF) is a synoptic survey designed to explore the transient and variable sky in a wide variety of cadences. We use PTF observations of fields that were observed multiple times (>=10) per night, for several nights, to find asteroids, construct their lightcurves and measure their rotation periods. Here we describe the pipeline we use to achieve these goals and present the results from the first four (overlapping) PTF fields analyzed as part of this program. These fields, which cover an area of 21 deg^2, were observed on four nights with a cadence of ~20 min. Our pipeline was able to detect 624 asteroids, of which 145 (~20%) were previously unknown. We present high quality rotation periods for 88 main-belt asteroids and possible period or lower limit on the period for an additional 85 asteroids. For the remaining 451 asteroids, we present lower limits on their photometric amplitudes. Three of the asteroids have lightcurves that are characteristic of binary asteroids. We estimate that implementing our search for all existing high-cadence PTF data will provide rotation periods for about 10,000 asteroids mainly in the magnitude range ~14 to ~20.

The Palomar Transient Factory (PTF) provides multiple epoch imaging for a large fraction of the sky. Here we describe the photometric calibration of the PTF data products that relates the PTF magnitudes to other mag systems. The calibration process utilizes SDSS r~16 mag point source objects as photometric standards. During photometric conditions, this allows us to solve for the extinction coefficients and color terms, and to estimate the camera illumination correction. This also enables the calibration of fields that are outside the SDSS footprint. We test the precision and repeatability of the PTF photometric calibration. Given that PTF is observing in a single filter each night, we define a PTF calibrated magnitude system for R-band and g-band. We show that, in this system, 59% (47%) of the photometrically calibrated PTF R-band (g-band) data achieve a photometric precision of 0.02-0.04 mag. Given the objects' color, the PTF magnitude system can be converted to other systems. A night-by-night comparison of the calibrated magnitudes of individual stars observed on multiple nights shows that they are consistent to a level of ~0.02 mag. Most of the data that were taken under non-photometric conditions can be calibrated relative to other epochs of the same sky footprint obtained during photometric conditions. We provide a guide describing the use of the PTF photometric calibration data products, as well as the transformations between the PTF magnitude system and the SDSS and Johnson-Cousins systems. (abridged).

The known host galaxies of short-hard gamma-ray bursts (GRBs) to date are characterized by low to moderate star-formation rates and a broad range of stellar masses. In this paper, we positionally associate the recent unambiguously short-hard Swift GRB 100206A with a disk galaxy at redshift z=0.4068 that is rapidly forming stars at a rate of ~30 M_sun/yr, almost an order of magnitude higher than any previously identified short GRB host. Using photometry from Gemini, Keck, PAIRITEL, and WISE, we show that the galaxy is very red (g-K = 4.3 AB mag), heavily obscured (A_V ~ 2 mag), and has the highest metallicity of any GRB host to date (12 + log[O/H]_KD02 = 9.2): it is a classical luminous infrared galaxy (LIRG), with L_IR ~ 4 x 10^11 L_sun. While these properties could be interpreted to support an association of this GRB with very recent star formation, modeling of the broadband spectral energy distribution also indicates that a substantial stellar mass of mostly older stars is present. The current specific star-formation rate is modest (specific SFR ~ 0.5 Gyr^-1), the current star-formation rate is not substantially elevated above its long-term average, and the host morphology shows no sign of recent merger activity. Our observations are therefore equally consistent with an older progenitor, similar to what is inferred for other short-hard GRBs. Given the precedent established by previous short GRB hosts and the significant fraction of the Universe's stellar mass in LIRG-like systems at z >~0.3, an older progenitor represents the most likely origin of this event.

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.

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.

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.

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.

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.

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.

VOEvent defines the content and meaning of a standard information packet for representing, transmitting, publishing and archiving information about a transient celestial event, with the implication that timely follow-up is of interest. The objective is to motivate the observation of targets-of-opportunity, to drive robotic telescopes, to trigger archive searches, and to alert the community. VOEvent is focused on the reporting of photon events, but events mediated by disparate phenomena such as neutrinos, gravitational waves, and solar or atmospheric particle bursts may also be reported. Structured data is used, rather than natural language, so that automated systems can effectively interpret VOEvent packets. Each packet may contain zero or more of the "who, what, where, when & how" of a detected event, but in addition, may contain a hypothesis (a "why") regarding the nature of the underlying physical cause of the event. Citations to previous VOEvents may be used to place each event in its correct context. Proper curation is encouraged throughout each event's life cycle from discovery through successive follow-ups. VOEvent packets gain persistent identifiers and are typically stored in databases reached via registries. VOEvent packets may therefore reference other packets in various ways. Packets are encouraged to be small and to be processed quickly. This standard does not define a transport layer or the design of clients, repositories, publishers or brokers; it does not cover policy issues such as who can publish, who can build a registry of events, who can subscribe to a particular registry, nor the intellectual property issues.

[Abridged] We present comprehensive X-ray, optical, near- and mid-infrared, and sub-mm observations of GRB 080207 and its host galaxy. The afterglow was undetected in the optical and near-IR, implying an optical to X-ray index <0.3, identifying GRB 080207 as a dark burst. Swift X-ray observations show extreme absorption in the host, which is confirmed by the unusually large optical extinction found by modelling the X-ray to nIR afterglow spectral energy distribution. Our Chandra observations obtained 8 days post-burst allow us to place the afterglow on the sky to sub-arcsec accuracy, enabling us to pinpoint an extremely red galaxy (ERO). Follow-up host observations with HST, Spitzer, Gemini, Keck and the James Clerk Maxwell Telescope (JCMT) provide a photometric redshift solution of z ~1.74 (+0.05,-0.06) (1 sigma), 1.56 < z < 2.08 at 2 sigma) for the ERO host, and suggest that it is a massive and morphologically disturbed ultra-luminous infrared galaxy (ULIRG) system, with L_FIR ~ 2.4 x 10^12 L_solar. These results add to the growing evidence that GRBs originating in very red hosts always show some evidence of dust extinction in their afterglows (though the converse is not true -- some extinguished afterglows are found in blue hosts). This indicates that a poorly constrained fraction of GRBs occur in very dusty environments. By comparing the inferred stellar masses, and estimates of the gas phase metallicity in both GRB hosts and sub-mm galaxies we suggest that many GRB hosts, even at z>2 are at lower metallicity than the sub-mm galaxy population, offering a likely explanation for the dearth of sub-mm detected GRB hosts. However, we also show that the dark GRB hosts are systematically more massive than those hosting optically bright events, perhaps implying that previous host samples are severely biased by the exclusion of dark events.

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

Understanding the origin and diversity of emission processes responsible for Gamma-ray Bursts (GRBs) remains a pressing challenge. While prompt and contemporaneous panchromatic observations have the potential to test predictions of the internal-external shock model, extensive multiband imaging has been conducted for only a few GRBs. We present rich, early-time, multiband datasets for two \swift\ events, GRB 110205A and GRB 110213A. The former shows optical emission since the early stages of the prompt phase, followed by the steep rising in flux up to ~1000s after the burst ($t^{-\alpha}$ with $\alpha=-6.13 \pm 0.75$). We discuss this feature in the context of the reverse-shock scenario and interpret the following single power-law decay as being forward-shock dominated. Polarization measurements, obtained with the RINGO2 instrument mounted on the Liverpool Telescope, also provide hints on the nature of the emitting ejecta. The latter event, instead, displays a very peculiar optical to near-infrared lightcurve, with two achromatic peaks. In this case, while the first peak is probably due to the onset of the afterglow, we interpret the second peak to be produced by newly injected material, signifying a late-time activity of the central engine.

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.

We present the discovery and follow-up observations of a broad-line type-Ic supernova (SN), PTF 10bzf (SN 2010ah), detected by the Palomar Transient Factory (PTF) on 2010 February 23. The SN distance is \cong 218 Mpc, greater than GRB 980425 / SN 1998bw and GRB 060218 / SN 2006aj, but smaller than the other SNe firmly associated with gamma-ray bursts (GRBs). We conducted a multi-wavelength follow-up campaign with Palomar-48 inch, Palomar 60-inch, Gemini-N, Keck, Wise, Swift, the Allen Telescope Array, CARMA, WSRT, and EVLA. Here we compare the properties of PTF 10bzf with those of SN 1998bw and other broad-line SNe. The optical luminosity and spectral properties of PTF 10bzf suggest that this SN is intermediate, in kinetic energy and amount of 56Ni, between non GRB-associated SNe like 2002ap or 1997ef, and GRB-associated SNe like 1998bw. No X-ray or radio counterpart to PTF 10bzf was detected. X-ray upper-limits allow us to exclude the presence of an underlying X-ray afterglow as luminous as that of other SN-associated GRBs like GRB 030329 or GRB 031203. Early-time radio upper-limits do not show evidence for mildly-relativistic ejecta. Late-time radio upper-limits rule out the presence of an underlying off-axis GRB, with energy and wind density similar to the SN-associated GRB 030329 and GRB 031203. Finally, by performing a search for a GRB in the time window and at the position of PTF 10bzf, we find that no GRB in the IPN catalog could be associated with this SN.

Type Ia supernovae are key tools for measuring distances on a cosmic scale. They are generally thought to be the thermonuclear explosion of an accreting white dwarf in a close binary system. The nature of the mass donor is still uncertain. In the single-degenerate model it is a main-sequence star or an evolved star, whereas in the double-degenerate model it is another white dwarf. We show that the velocity structure of absorbing material along the line of sight to 35 type Ia supernovae tends to be blueshifted. These structures are likely signatures of gas outflows from the supernova progenitor systems. Thus many type Ia supernovae in nearby spiral galaxies may originate in single-degenerate systems.

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$\alpha$ 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.

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.

Stellar rotation periods measured from single-age populations are critical for investigating how stellar angular momentum content evolves over time, how that evolution depends on mass, and how rotation influences the stellar dynamo and the magnetically heated chromosphere and corona. We report rotation periods for 40 late-K to mid-M stars members of the nearby, rich, intermediate-age (~600 Myr) open cluster Praesepe. These rotation periods were derived from ~200 observations taken by the Palomar Transient Factory of four cluster fields from 2010 February to May. Our measurements indicate that Praesepe's mass-period relation transitions from a well-defined singular relation to a more scattered distribution of both fast and slow rotators at ~0.6 Msun. The location of this transition is broadly consistent with expectations based on observations of younger clusters and the assumption that stellar-spin down is the dominant mechanism influencing angular momentum evolution at 600 Myr. However, a comparison to data recently published for the Hyades, assumed to be coeval to Praesepe, indicates that the divergence from a singular mass-period relation occurs at different characteristic masses, strengthening the finding that Praesepe is the younger of the two clusters. We also use previously published relations describing the evolution of rotation periods as a function of color and mass to evolve the sample of Praesepe periods in time. Comparing the resulting predictions to periods measured in M35 and NGC 2516 (~150 Myr) and for kinematically selected young and old field star populations suggests that stellar spin-down may progress more slowly than described by these relations.

We report the detection of a faint optical flash by the 2-m Faulkes Telescope North simultaneously with the second of two prompt gamma-ray pulses in INTEGRAL gamma-ray burst (GRB) 080603A, beginning at t_rest = 37 s after the onset of the GRB. This optical flash appears to be distinct from the subsequent emerging afterglow emission, for which we present comprehensive broadband radio to X-ray light curves to 13 days post-burst and rigorously test the standard fireball model. The intrinsic extinction toward GRB 080603A is high (A_V,z = 0.8 mag), and the well-sampled X-ray-to-near-infrared spectral energy distribution is interesting in requiring an LMC2 extinction profile, in contrast to the majority of GRBs. Comparison of the gamma-ray and extinction-corrected optical flux densities of the flash rules out an inverse-Compton origin for the prompt gamma-rays; instead, we suggest that the optical flash could originate from the inhomogeneity of the relativistic flow. In this scenario, a large velocity irregularity in the flow produces the prompt gamma-rays, followed by a milder internal shock at a larger radius that would cause the optical flash. Flat gamma-ray spectra, roughly F propto nu^-0.1, are observed in many GRBs. If the flat spectrum extends down to the optical band in GRB 080603A, the optical flare could be explained as the low-energy tail of the gamma-ray emission. If this is indeed the case, it provides an important clue to understanding the nature of the emission process in the prompt phase of GRBs and highlights the importance of deep (R> 20 mag), rapid follow-up observations capable of detecting faint, prompt optical emission.

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.

(Abridged) Active galactic nuclei show a wealth of interesting physical processes, some of which are poorly understood. We want to address a number of open questions, including the location and physics of the outflow from AGN, the nature of the continuum emission, the geometry and physical state of the X-ray broad emission line region, the Fe-K line complex, the metal abundances of the nucleus and finally the interstellar medium of our own Galaxy. We study one of the best targets for these aims, the Seyfert 1 galaxy Mrk 509 with a multiwavelength campaign using five satellites (XMM-Newton, INTEGRAL, Chandra, HST and Swift) and two ground-based facilities (WHT and PAIRITEL). Our observations cover more than five decades in frequency, from 2 um to 200 keV. The combination of high-resolution spectroscopy and time variability allows us to disentangle and study the different components. Our campaign covers 100 days from September to December 2009, and is centred on a simultaneous set of deep XMM-Newton and INTEGRAL observations with regular time intervals, spanning seven weeks. We obtain a continuous light curve in the X-ray and UV band, showing a strong, up to 60% flux increase in the soft X-ray band during the three weeks in the middle of our deepest monitoring campaign, and which is correlated with an enhancement of the UV flux. This allows us to study the time evolution of the continuum and the outflow. By stacking the observations, we have also obtained one of the best X-ray and UV spectra of a Seyfert galaxy ever obtained. In this paper we also study the effects of the spectral energy distribution (SED) that we obtained on the photo-ionisation equilibrium. Thanks to our broad-band coverage, uncertainties on the SED do not strongly affect the determination of this equilibrium.

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.

The rate of image acquisition in modern synoptic imaging surveys has already begun to outpace the feasibility of keeping astronomers in the real-time discovery and classification loop. Here we present the inner workings of a framework, based on machine-learning algorithms, that captures expert training and ground-truth knowledge about the variable and transient sky to automate 1) the process of discovery on image differences and, 2) the generation of preliminary science-type classifications of discovered sources. Since follow-up resources for extracting novel science from fast-changing transients are precious, self-calibrating classification probabilities must be couched in terms of efficiencies for discovery and purity of the samples generated. We estimate the purity and efficiency in identifying real sources with a two-epoch image-difference discovery algorithm for the Palomar Transient Factory (PTF) survey. Once given a source discovery, using machine-learned classification trained on PTF data, we distinguish between transients and variable stars with a 3.8% overall error rate (with 1.7% errors for imaging within the Sloan Digital Sky Survey footprint). At >96% classification efficiency, the samples achieve 90% purity. Initial classifications are shown to rely primarily on context-based features, determined from the data itself and external archival databases. In the ~one year since autonomous operations, this discovery and classification framework has led to several significant science results, from outbursting young stars to subluminous Type IIP supernovae to candidate tidal disruption events. We discuss future directions of this approach, including the possible roles of crowdsourcing and the scalability of machine learning to future surveys such a the Large Synoptical Survey Telescope (LSST).

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.

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.

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.

Variable X-ray and gamma-ray emission is characteristic of the most extreme physical processes in the Universe, and studying the sources of these energetic photons has been a major driver in astronomy for the past 50 years. Here we present multiwavelength observations of a unique gamma-ray selected transient, discovered by Swift, which was accompanied by bright emission across the electromagnetic spectrum, and whose properties are unlike any previously observed source. We pinpoint the event to the center of a small, star-forming galaxy at redshift z=0.3534. Its high-energy emission has lasted much longer than any gamma-ray burst, while its peak luminosity was about 100 times higher than the brightest active galactic nuclei. The association of the outburst with the center of its host galaxy suggests that this phenomenon has its origin in a new, rare mechanism associated with a massive black hole in the nucleus of a galaxy.

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.

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.

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.

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.

Metallicity is expected to influence not only the lives of massive stars but also the outcome of their deaths as supernovae (SNe) and as gamma-ray bursts (GRBs). However, there are surprisingly few direct measurements of the local metallicities of different flavors of core-collapse SNe. Here we present the largest existing set of host-galaxy spectra with H II region emission lines at the sites of 35 stripped-envelope core-collapse SNe. We derive local oxygen abundances in a robust manner in order to constrain the SN Ib/c progenitor population. We obtain spectra at the SN sites, include SNe from targeted and untargeted surveys, and perform the abundance determinatinos using three different oxygen-abundance calibrations. The sites of SNe Ic (the demise of the most heavily stripped stars having lost both the H and He layers) are systematically more metal rich than those of SNe Ib (arising from stars that retained their He layer) in all calibrations. A Kolmogorov-Smirnov test yields the very low probability of 1% that SN Ib and SN Ic environment abundances, which are different on average by ~0.2 dex (in the Pettini & Pagel scale), are drawn from the same parent population. Broad-lined SNe Ic (without GRBs) occur at metallicities between those of SNe Ib and SNe Ic. Lastly, we find that the host-galaxy central oxygen abundance is not a good indicator of the local SN metallicity; hence, large-scale SN surveys need to obtain local abundance measurements in order to quantify the impact of metallicity on stellar death.

PTF09dav is a peculiar subluminous type Ia supernova (SN) discovered by the Palomar Transient Factory (PTF). Spectroscopically, it appears superficially similar to the class of subluminous SN1991bg-like SNe, but it has several unusual features which make it stand out from this population. Its peak luminosity is fainter than any previously discovered SN1991bg-like SN Ia (M_B -15.5), but without the unusually red optical colors expected if the faint luminosity were due to extinction. The photospheric optical spectra have very unusual strong lines of Sc II and Mg I, with possible Sr II, together with stronger than average Ti II and low velocities of ~6000 km/s. The host galaxy of PTF09dav is ambiguous. The SN lies either on the extreme outskirts (~41kpc) of a spiral galaxy, or in an very faint (M_R>-12.8) dwarf galaxy, unlike other 1991bg-like SNe which are invariably associated with massive, old stellar populations. PTF09dav is also an outlier on the light-curve-width--luminosity and color--luminosity relations derived for other sub-luminous SNe Ia. The inferred 56Ni mass is small (0.019+/-0.003Msun), as is the estimated ejecta mass of 0.36Msun. Taken together, these properties make PTF09dav a remarkable event. We discuss various physical models that could explain PTF09dav. Helium shell detonation or deflagration on the surface of a CO white-dwarf can explain some of the features of PTF09dav, including the presence of Sc and the low photospheric velocities, but the observed Si and Mg are not predicted to be very abundant in these models. We conclude that no single model is currently capable of explaining all of the observed signatures of PTF09dav.

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.

Supernova (SN) 2006gy was a hydrogen-rich core-collapse SN that remains one of the most luminous optical supernovae ever observed. The total energy budget (> 2 x 10^51 erg radiated in the optical alone) poses many challenges for standard SN theory. We present new ground-based near-infrared (NIR) observations of SN 2006gy, as well as a single epoch of Hubble Space Telescope (HST) imaging obtained more than two years after the explosion. Our NIR data taken around peak optical emission show an evolution that is largely consistent with a cooling blackbody, with tentative evidence for a growing NIR excess starting at day ~100. Our late-time Keck adaptive optics (AO) NIR image, taken on day 723, shows little change from previous NIR observations taken around day 400. Furthermore, the optical HST observations show a reduced decline rate after day 400, and the SN is bluer on day 810 than it was at peak. This late-time decline is inconsistent with Co56 decay, and thus is problematic for the various pair-instability SN models used to explain the nature of SN 2006gy. The slow decline of the NIR emission can be explained with a light echo, and we confirm that the late-time NIR excess is the result of a massive (>10 Msun) dusty shell heated by the SN peak luminosity. The late-time optical observations require the existence of a scattered light echo, which may be generated by the same dust that contributes to the NIR echo. Both the NIR and optical echoes originate in the proximity of the progenitor, ~10^18 cm for the NIR echo and <~10-40 pc for the optical echo, which provides further evidence that the progenitor of SN 2006gy was a very massive star.

We present spectroscopic and photometric observations of the Type IIn supernova (SN) 2008iy. SN 2008iy showed an unprecedentedly long rise time of ~400 days, making it the first SN to take significantly longer than 100 days to reach peak optical luminosity. The peak absolute magnitude of SN 2008iy was M_r ~ -19.1 mag, and the total radiated energy over the first ~700 days was ~2 x 10^50 erg. Spectroscopically, SN 2008iy is very similar to the Type IIn SN 1988Z at late times, and, like SN 1988Z, it is a luminous X-ray source (both supernovae had an X-ray luminosity L_ X > 10^41 erg/s). The Halpha emission profile of SN 2008iy shows a narrow P Cygni absorption component, implying a pre-SN wind speed of ~100 km/s. We argue that the luminosity of SN 2008iy is powered via the interaction of the SN ejecta with a dense, clumpy circumstellar medium. The ~400 day rise time can be understood if the number density of clumps increases with distance over a radius ~1.7 x 10^16 cm from the progenitor. This scenario is possible if the progenitor experienced an episodic phase of enhanced mass-loss < 1 century prior to explosion or the progenitor wind speed increased during the decades before core collapse. We favour the former scenario, which is reminiscent of the eruptive mass-loss episodes observed for luminous blue variable (LBV) stars. The progenitor wind speed and increased mass-loss rates serve as further evidence that at least some, and perhaps all, Type IIn supernovae experience LBV-like eruptions shortly before core collapse. We also discuss the host galaxy of SN 2008iy, a subluminous dwarf galaxy, and offer a few reasons why the recent suggestion that unusual, luminous supernovae preferentially occur in dwarf galaxies may be the result of observational biases.

We present observations and interpretation of the Type IIn supernova SN 2008am discovered by the ROTSE Supernova Verification Project (RSVP). SN 2008am peaked at approximately -22.3 mag at a redshift of z=0.2338, giving it a peak luminosity of 3 x 10^{44}erg/s and making it one of the most luminous supernovae ever observed. The total radiated energy is ~ 2 x 10^{51} erg. Photometric observations in the ultraviolet, optical and infrared bands (J,H,Ks) constrain the SED evolution. We obtained six optical spectra of the supernova, five on the early decline from maximum light and a sixth nearly a year later plus a very late-time spectrum (~2 yr) of the host galaxy. The spectra of SN 2008am show strong Balmer-line and He I lambda 5876A emission with intermediate widths (~25A) in the first ~40 days after optical maximum. We examine a variety of models for the line wings and conclude that multiple scattering is most likely, implying that our spectra contain no specific information on the bulk flow velocity. We examine a variety of models for the ROTSE light curve subject to the rise time and the nature of the spectra, including radioactive decay, shocks in optically-thick and optically-thin circumstellar media (CSM) and a magnetar. The most successful model is one for which the CSM is optically-thick and in which diffusion of forward shock-deposited luminosity gives rise to the observed light curve. Diffusion of the shock-deposited energy from the forward shock is found to be important to account for the rising part of the light curve. Although there are differences in detail, SN 2008am appears to be closely related to other super-luminous Type IIn supernovae, SN 2006gy, SN 2006tf and perhaps SN 2008iy, that may represent the deaths of very massive LBV-type progenitors and for which the luminosity is powered by the interaction of the ejecta with a dense circumstellar medium.

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%.

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.

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.

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.

This paper presents the first results from a new citizen science project: Galaxy Zoo Supernovae. This proof of concept project uses members of the public to identify supernova candidates from the latest generation of wide-field imaging transient surveys. We describe the Galaxy Zoo Supernovae operations and scoring model, and demonstrate the effectiveness of this novel method using imaging data and transients from the Palomar Transient Factory (PTF). We examine the results collected over the period April-July 2010, during which nearly 14,000 supernova candidates from PTF were classified by more than 2,500 individuals within a few hours of data collection. We compare the transients selected by the citizen scientists to those identified by experienced PTF scanners, and find the agreement to be remarkable - Galaxy Zoo Supernovae performs comparably to the PTF scanners, and identified as transients 93% of the ~130 spectroscopically confirmed SNe that PTF located during the trial period (with no false positive identifications). Further analysis shows that only a small fraction of the lowest signal-to-noise SN detections (r > 19.5) are given low scores: Galaxy Zoo Supernovae correctly identifies all SNe with > 8{\sigma} detections in the PTF imaging data. The Galaxy Zoo Supernovae project has direct applicability to future transient searches such as the Large Synoptic Survey Telescope, by both rapidly identifying candidate transient events, and via the training and improvement of existing machine classifier algorithms.

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.

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.

We report SMARTS, Gemini and Swift-UVOT observations of the optical transient (OT) associated with gamma-ray burst (GRB) 091127, at redshift 0.49, taken between 0.9 hr and 102 days following the Swift trigger. In our early-time observations, the OT fades in a manner consistent with previously observed GRB afterglows. However, after 9 days post-burst, the OT is observed to brighten for a period of ~2 weeks, after which the source resumes fading. A comparison of this late-time "bump" to SN 1998bw (the broad-lined Type Ic supernova associated with GRB 980425), and several other GRB supernovae (SNe), indicates that the most straightforward explanation is that GRB 091127 was accompanied by a contemporaneous SN (SN 2009nz) that peaked at a magnitude of M_V=-19.0+/-0.2. SN 2009nz is globally similar to other GRB supernovae, but evolves slightly faster than SN 1998bw and reaches a slightly dimmer peak magnitude. We also analyze the early-time UV-optical-IR spectral energy distribution of the afterglow of GRB 091127 and find that there is little to no reddening in the host galaxy along the line-of-slight to this burst.

Type-IIn supernovae (SNe), which are characterized by strong interaction of their ejecta with the surrounding circumstellar matter (CSM), provide a unique opportunity to study the mass-loss history of massive stars shortly before their explosive death. We present the discovery and follow-up observations of a Type IIn SN, PTF 09uj, detected by the Palomar Transient Factory (PTF). Serendipitous observations by GALEX at ultraviolet (UV) wavelengths detected the rise of the SN light curve prior to the PTF discovery. The UV light curve of the SN rose fast, with a time scale of a few days, to a UV absolute AB magnitude of about -19.5. Modeling our observations, we suggest that the fast rise of the UV light curve is due to the breakout of the SN shock through the dense CSM (n~10^10 cm^-3). Furthermore, we find that prior to the explosion the progenitor went through a phase of high mass-loss rate (~0.1 solar mass per year) that lasted for a few years. The decay rate of this SN was fast relative to that of other SNe IIn.

We present early-time optical through infrared photometry of the bright gamma-ray burst GRB 080607, starting only 6 s following the initial trigger in the rest frame. Complemented by our previously published spectroscopy, this high-quality photometric dataset allows us to solve for the extinction properties of the redshift 3.036 sightline, giving perhaps the most detailed information on the ultraviolet continuum absorption properties of any sightline outside our Local Group to date. The extinction properties are not adequately modeled by any ordinary extinction template (including the average Milky Way, Large Magellanic Cloud, and Small Magellanic Cloud curves), partially because the 2175-Angstrom feature (while present) is weaker by about a factor of two than when seen under similar circumstances locally. However, the spectral energy distribution is exquisitely fitted by the more general Fitzpatrick & Massa (1990) parameterization of Local-Group extinction, putting it in the same family as some peculiar Milky Way extinction curves. After correcting for this (considerable, A_V = 3.3 +/- 0.4 mag) extinction, GRB 080607 is revealed to have been among the most optically luminous events ever observed, comparable to the naked-eye burst GRB 080319B. Its early peak time (t_rest < 6 s) indicates a high initial Lorentz factor (Gamma > 600), while the extreme luminosity may be explained in part by a large circumburst density. Only because of its early high luminosity could the afterglow of GRB 080607 be studied in such detail in spite of the large attenuation and great distance, making this burst an excellent prototype for the understanding of other highly obscured extragalactic objects, and of the class of "dark" GRBs in particular.

The Energetic X-ray Imaging Survey Telescope (EXIST) is designed to i) use the birth of stellar mass black holes, as revealed by cosmic Gamma-Ray Bursts (GRBs), as probes of the very first stars and galaxies to exist in the Universe. Both their extreme luminosity (~104 times larger than the most luminous quasars) and their hard X-ray detectability over the full sky with wide-field imaging make them ideal "back-lights" to measure cosmic structure with X-ray, optical and near-IR (nIR) spectra over many sight lines to high redshift. The full-sky imaging detection and rapid followup narrow-field imaging and spectroscopy allow two additional primary science objectives: ii) novel surveys of supermassive black holes (SMBHs) accreting as very luminous but rare quasars, which can trace the birth and growth of the first SMBHs as well as quiescent SMBHs (non-accreting) which reveal their presence by X-ray flares from the tidal disruption of passing field stars; and iii) a multiwavelength Time Domain Astrophysics (TDA) survey to measure the temporal variability and physics of a wide range of objects, from birth to death of stars and from the thermal to non-thermal Universe. These science objectives are achieved with the telescopes and mission as proposed for EXIST described here.

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.

Aims: With this paper we want to investigate the highly variable afterglow light curve and environment of gamma-ray burst (GRB) 060526 at $z=3.221$. Methods: We present one of the largest photometric datasets ever obtained for a GRB afterglow, consisting of multi-color photometric data from the ultraviolet to the near infrared. The data set contains 412 data points in total to which we add additional data from the literature. Furthermore, we present low-resolution high signal-to-noise spectra of the afterglow. The afterglow light curve is modeled with both an analytical model using broken power law fits and with a broad-band numerical model which includes energy injections. The absorption lines detected in the spectra are used to derive column densities using a multi-ion single-component curve-of-growth analysis from which we derive the metallicity of the host of GRB 060526. Results: The temporal behaviour of the afterglow follows a double broken power law with breaks at $t=0.090\pm0.005$ and $t=2.401\pm0.061$ days. It shows deviations from the smooth set of power laws that can be modeled by additional energy injections from the central engine, although some significant microvariability remains. The broadband spectral-energy distribution of the afterglow shows no significant extinction along the line of sight. The metallicity derived from \ion{S}{II} and \ion{Fe}{II} of [S/H] = --0.57 $\pm$0.25 and [Fe/H] = --1.09$\pm$0.24 is relatively high for a galaxy at that redshift but comparable to the metallicity of other GRB hosts at similar redshifts. At the position of the afterglow, no host is detected to F775W(AB) = 28.5 mag with the HST, implying an absolute magnitude of the host M(1500 \AA{})$>$--18.3 mag which is fainter than most long-duration hosts, although the GRB may be associated with a faint galaxy at a distance of 11 kpc.

In August 2008, the accreting milli-second X-ray pulsar (AMXP), IGR J00291+5934, underwent an outburst lasting ~ 100 days, the first since its discovery in 2004. We present data from the double-peaked outburst from Faulkes Telescope North, the INT, the Keck Telescope, PAIRITEL, the Westerbork Synthesis Radio Telescope and the Swift, XMM-Newton and RXTE X-ray missions. We study the outburst's evolution at various wavelengths. We study the light curve morphology, presenting the first radio-X-ray Spectral Energy Distributions (SEDs) for this source and the most detailed UV-IR SEDs for any outbursting AMXP. We show simple models that attempt to identify the emission mechanisms responsible. We analyse short-timescale optical variability, and compare a medium resolution optical spectrum with those from 2004. The outburst morphology is unusual for an AMXP, comprising two peaks, the second containing a 'plateau' of ~ 10 days at maximum brightness within 30 days of the initial activity. This has implications on duty cycles of short-period X-ray transients. The X-ray spectrum can be fitted by a single, hard power-law. We detect optical variability of ~ 0.05 magnitudes, on timescales of minutes, but find no periodic modulation. In the optical, the SEDs contain a blue component, indicative of an irradiated disc, and a transient near-infrared (NIR) excess. This excess is consistent with a simple model of an optically thick synchrotron jet (as seen in other outbursting AMXPs). The optical spectrum shows a double-peaked H alpha profile, a diagnostic of an accretion disc, but we do not clearly see other lines (e.g. He I, II) reported in 2004. Optical/IR observations of AMXPs are excellent for studying the evolution of both the outer accretion disc and the inner jet, and may eventually provide us with tight constraints to model disc-jet coupling in accreting neutron stars.

We present broadband (radio, optical, and X-ray) light curves and spectra of the afterglows of four long-duration gamma-ray bursts (GRBs 090323, 090328, 090902B, and 090926A) detected by the Gamma-Ray Burst Monitor (GBM) and Large Area Telescope (LAT) instruments on the Fermi satellite. With its wide spectral bandpass, extending to GeV energies, Fermi is sensitive to GRBs with very large isotropic energy releases (10e54 erg). Although rare, these events are particularly important for testing GRB central-engine models. When combined with spectroscopic redshifts, our afterglow data for these four events are able to constrain jet collimation angles, the density structure of the circumburst medium, and both the true radiated energy release and the kinetic energy of the outflows. In agreement with our earlier work, we find that the relativistic energy budget of at least one of these events (GRB 090926A) exceeds the canonical value of 10e51 erg by an order of magnitude. Such energies pose a severe challenge for models in which the GRB is powered by a magnetar or neutrino-driven collapsar, but remain compatible with theoretical expectations for magneto-hydrodynamical (MHD) collapsar models. Our jet opening angles (theta) are similar to those found for pre-Fermi GRBs, but the large initial Lorentz factors (Gamma_0) inferred from the detection of GeV photons imply theta Gamma_0 ~ 70-90, values which are above those predicted in MHD models of jet acceleration. Finally, we find that these Fermi-LAT events preferentially occur in a low-density circumburst environment, and we speculate that this might result from the lower mass-loss rates of their lower-metallicity progenitor stars. Future studies of Fermi-LAT afterglows in the radio with the order-of-magnitude improvement in sensitivity offered by the EVLA should definitively establish the relativistic energy budgets of these events.

We present early-time photometric and spectroscopic observations of supernova (SN) 2009kr in NGC 1832. We find that its properties to date support its classification as Type II-linear (SN II-L), a relatively rare subclass of core-collapse supernovae (SNe). We have also identified a candidate for the SN progenitor star through comparison of pre-explosion, archival images taken with WFPC2 onboard the Hubble Space Telescope (HST) with SN images obtained using adaptive optics (AO) plus NIRC2 on the 10-m Keck-II telescope. Although the host galaxy's substantial distance (~26 Mpc) results in large uncertainties in the relative astrometry, we find that if this candidate is indeed the progenitor, it is a highly luminous (M_V = -7.8 mag) yellow supergiant with initial mass ~18-24 M_sun. This would be the first time that a SN II-L progenitor has been directly identified. Its mass may be a bridge between the upper initial mass limit for the more common Type II-plateau SNe (SNe II-P) and the inferred initial mass estimate for one Type II-narrow SN (SN IIn).

The epochs of origin of the first stars and galaxies, and subsequent growth of the first supermassive black holes, are among the most fundamental questions. Observations of the highest redshift Gamma-Ray Bursts (GRBs) will be the most compelling in situ probe of the history of initial star formation and consequent epoch of reionization if their prompt and precise detection can be followed immediately by sensitive near-IR imaging and spectroscopy. Blazars are the persistent analogs of GRBs and for the same reason (beaming) can be observed at highest redshifts where they might best trace the high accretion rate-driven jets and growth of supermassive black holes in galaxies. The proposed EXIST mission can uniquely probe these questions, and many others, given its unparalled combination of sensitivity and spatial-spectral-temporal coverage and resolution. Here we provide a brief summary of the mission design, key science objectives, mission plan and readiness for EXIST, as proposed to Astro2010.

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.

We present observations and analysis of the broadband afterglow of Swift GRB 071025. Using optical and infrared (RIYJHK) photometry, we derive a photometric redshift of 4.4 < z < 5.2; at this redshift our simultaneous multicolour observations begin at ~30 s after the GRB trigger in the host frame and during the initial rising phase of the afterglow. We associate the light curve peak at 580 s in the observer frame with the formation of the forward shock, giving an estimate of the initial Lorentz factor Gamma_0 ~ 200. The red spectral energy distribution (even in regions not affected by the Lyman-alpha break) provides secure evidence of a large dust column. However, the inferred extinction curve shows a prominent flat component between 2000-3000 Angstroms in the rest-frame, inconsistent with any locally observed template but well-fit by models of dust formed by supernovae. Time-dependent fits to the extinction profile reveal no evidence of dust destruction and limit the decrease in the extinction column to Delta A_3000 < 0.54 mag after t = 50 s in the rest frame. Our observations provide evidence of a transition in dust properties at z~5, in agreement with studies of high-z quasars, and suggest that SN-formed dust continues to dominate the opacity of typical galaxies at this redshift.

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.

Stars with initial masses 10 M_{solar} < M_{initial} < 100 M_{solar} fuse progressively heavier elements in their centres, up to inert iron. The core then gravitationally collapses to a neutron star or a black hole, leading to an explosion -- an iron-core-collapse supernova (SN). In contrast, extremely massive stars (M_{initial} > 140 M_{solar}), if such exist, have oxygen cores which exceed M_{core} = 50 M_{solar}. There, high temperatures are reached at relatively low densities. Conversion of energetic, pressure-supporting photons into electron-positron pairs occurs prior to oxygen ignition, and leads to a violent contraction that triggers a catastrophic nuclear explosion. Tremendous energies (>~ 10^{52} erg) are released, completely unbinding the star in a pair-instability SN (PISN), with no compact remnant. Transitional objects with 100 M_{solar} < M_{initial} < 140 M_{solar}, which end up as iron-core-collapse supernovae following violent mass ejections, perhaps due to short instances of the pair instability, may have been identified. However, genuine PISNe, perhaps common in the early Universe, have not been observed to date. Here, we present our discovery of SN 2007bi, a luminous, slowly evolving supernova located within a dwarf galaxy (~1% the size of the Milky Way). We measure the exploding core mass to be likely ~100 M_{solar}, in which case theory unambiguously predicts a PISN outcome. We show that >3 M_{solar} of radioactive 56Ni were synthesized, and that our observations are well fit by PISN models. A PISN explosion in the local Universe indicates that nearby dwarf galaxies probably host extremely massive stars, above the apparent Galactic limit, perhaps resulting from star formation processes similar to those that created the first stars in the Universe.

The quasar SDSS J153636.22+044127.0, exhibiting peculiar broad emission-line profiles with multiple components, was proposed as a candidate sub-parsec binary supermassive black hole system. More recently, imaging revealed two spatially distinct sources, leading some to suggest the system to be a quasar pair separated by ~5 kpc. We present Palomar and Keck optical spectra of this system from which we identify a third velocity component to the emission lines. We argue that the system is more likely an unusual member of the class of active galactic nuclei (AGNs) known as "double-peaked emitters" than a sub-parsec black hole binary or quasar pair. We find no significant velocity evolution of the two main peaks over the course of 0.95 yr, with a 3-sigma upper limit on any secular change of 70 km/s/yr. We also find that the three velocity components of the emission lines are spatially coincident to within 0.015" along the slit, apparently ruling out the double-quasar hypothesis.

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.

We present broadband (gamma-ray, X-ray, near-infrared, optical, and radio) observations of the gamma-ray burst (GRB) 090709A and its afterglow in an effort to ascertain the origin of this high-energy transient. Previous analyses suggested that GRB 090709A exhibited quasi-periodic oscillations with a period of 8.06 s, a trait unknown in long-duration GRBs but typical of flares from soft gamma-ray repeaters. When properly accounting for the underlying shape of the power-density spectrum of GRB 090709A, we find no conclusive (> 3 sigma) evidence for the reported periodicity. In conjunction with the location of the transient (far from the Galactic plane and absent any nearby host galaxy in the local universe) and the evidence for extinction in excess of the Galactic value, we consider a magnetar origin relatively unlikely. A long-duration GRB, however, can account for the majority of the observed properties of this source. GRB 090709A is distinguished from other long-duration GRBs primarily by the large amount of obscuration from its host galaxy (A_K,obs >~ 2 mag).