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.

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

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.

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.

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.

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.

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

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

In this work we present the first results of our imaging campaign at Keck Observatory to identify the host galaxies of "dark" gamma-ray bursts (GRBs), events with no detected optical afterglow or with detected optical flux significantly fainter than expected from the observed X-ray afterglow. We find that out of a uniform sample of 29 Swift bursts rapidly observed by the Palomar 60-inch telescope through March 2008 (14 of which we classify as dark), all events have either a detected optical afterglow, a probable optical host-galaxy detection, or both. Our results constrain the fraction of Swift GRBs coming from very high redshift (z > 7), such as the recent GRB 090423, to between 0.2-7 percent at 80% confidence. In contrast, a significant fraction of the sample requires large extinction columns (host-frame A_V > 1 mag, with several events showing A_V > 2-6 mag), identifying dust extinction as the dominant cause of the dark GRB phenomenon. We infer that a significant fraction of GRBs (and, by association, of high-mass star formation) occurs in highly obscured regions. However, the host galaxies of dark GRBs seem to have normal optical colors, suggesting that the source of obscuring dust is local to the vicinity of the GRB progenitor or highly unevenly distributed within the host galaxy.

Swift captured for the first time a smoothly rising X-ray re-brightening of clear non-flaring origin after the steep decay in a long gamma-ray burst (GRB): GRB081028. A rising phase is likely present in all GRBs but is usually hidden by the prompt tail emission and constitutes the first manifestation of what is later to give rise to the shallow decay phase. Contemporaneous optical observations reveal a rapid evolution of the injection frequency of a fast cooling synchrotron spectrum through the optical band, which disfavours the afterglow onset (start of the forward shock emission along our line of sight when the outflow is decelerated) as the origin of the observed re-brightening. We investigate alternative scenarios and find that the observations are consistent with the predictions for a narrow jet viewed off-axis. The high on-axis energy budget implied by this interpretation suggests different physical origins of the prompt and (late) afterglow emission. Strong spectral softening takes place from the prompt to the steep decay phase: we track the evolution of the spectral peak energy from the gamma-rays to the X-rays and highlight the problems of the high latitude and adiabatic cooling interpretations. Notably, a softening of both the high and low spectral slopes with time is also observed. We discuss the low on-axis radiative efficiency of GRB081028 comparing its properties against a sample of Swift long GRBs with secure E_gamma,iso measurements.

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

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

We present extensive early photometric (ultraviolet through near-infrared) and spectroscopic (optical and near-infrared) data on supernova (SN) 2008D as well as X-ray data analysis on the associated Swift/X-ray transient (XRT) 080109. Our data span a time range of 5 hours before the detection of the X-ray transient to 150 days after its detection, and detailed analysis allowed us to derive constraints on the nature of the SN and its progenitor; throughout we draw comparisons with results presented in the literature and find several key aspects that differ. We show that the X-ray spectrum of XRT 080109 can be fit equally well by an absorbed power law or a superposition of about equal parts of both power law and blackbody. Our data first established that SN 2008D is a spectroscopically normal SN Ib (i.e., showing conspicuous He lines), and show that SN 2008D had a relatively long rise time of 18 days and a modest optical peak luminosity. The early-time light curves of the SN are dominated by a cooling stellar envelope (for \Delta t~0.1- 4 day, most pronounced in the blue bands) followed by 56^Ni decay. We construct a reliable measurement of the bolometric output for this stripped-envelope SN, and, combined with estimates of E_K and M_ej from the literature, estimate the stellar radius R_star of its probable Wolf-Rayet progenitor. According to the model of Waxman et al. and of Chevalier & Fransson, we derive R_star^{W07}= 1.2+/-0.7 R_sun and R_star^{CF08}= 12+/-7 R_sun, respectively; the latter being more in line with typical WN stars. Spectra obtained at 3 and 4 months after maximum light show double-peaked oxygen lines that we associate with departures from spherical symmetry, as has been suggested for the inner ejecta of a number of SN Ib cores.

Long-duration gamma-ray bursts (GRBs) are widely believed to be highly-collimated explosions (opening angle theta ~ 1-10 deg). As a result of this beaming factor, the true energy release from a GRB is usually several orders of magnitude smaller than the observed isotropic value. Measuring this opening angle, typically inferred from an achromatic steepening in the afterglow light curve (a "jet" break), has proven exceedingly difficult in the Swift era. Here we undertake a study of five of the brightest (in terms of the isotropic prompt gamma-ray energy release, E(gamma, iso)) GRBs in the Swift era to search for jet breaks and hence constrain the collimation-corrected energy release. We present multi-wavelength (radio through X-ray) observations of GRBs 050820A, 060418, and 080319B, and construct afterglow models to extract the opening angle and beaming-corrected energy release for all three events. Together with results from previous analyses of GRBs 050904 and 070125, we find evidence for an achromatic jet break in all five events, strongly supporting the canonical picture of GRBs as collimated explosions. The most natural explanation for the lack of observed jet breaks from most Swift GRBs is therefore selection effects. However, the opening angles for the events in our sample are larger than would be expected if all GRBs had a canonical energy release of ~ 10e51 erg. The total energy release we measure for those "hyper-energetic" (E(total) >~ 10e52 erg) events in our sample is large enough to start challenging models with a magnetar as the compact central remnant.

We have obtained 1087 NIR (JHKs) measurements of 21 SNe Ia using PAIRITEL, nearly doubling the number of well-sampled NIR SN Ia light curves. These data strengthen the evidence that SNe Ia are excellent standard candles in the NIR, even without correction for optical light-curve shape. We construct fiducial NIR templates for normal SNe Ia from our sample, excluding only the three known peculiar SNe Ia: SN 2005bl, SN 2005hk, and SN 2005ke. The H-band absolute magnitudes in this sample of 18 SNe Ia have an intrinsic rms of only 0.15 mag with no correction for light-curve shape. We found a relationship between the H-band extinction and optical color excess of AH=0.2E(B-V). This variation is as small as the scatter in distance modulus measurements currently used for cosmology based on optical light curves after corrections for light-curve shape. Combining the homogeneous PAIRITEL measurements with 23 SNe Ia from the literature, these 41 SNe Ia have standard H-band magnitudes with an rms scatter of 0.16 mag. The good match of our sample with the literature sample suggests there are few systematic problems with the photometry. We present a nearby NIR Hubble diagram that shows no correlation of the residuals from the Hubble line with light-curve properties. Future samples that account for optical and NIR light-curve shapes, absorption, spectroscopic variation, or host-galaxy properties may reveal effective ways to improve the use of SNe Ia as distance indicators. Since systematic errors due to dust absorption in optical bands remain the leading difficulty in the cosmological use of supernovae, the good behavior of SN Ia NIR light curves and their relative insensitivity to reddening make these objects attractive candidates for future cosmological work.

We report on strong H2 and CO absorption from gas within the host galaxy of gamma-ray burst (GRB) 080607. Analysis of our Keck/LRIS afterglow spectrum reveals a very large HI column density (NHI = 10^22.70 cm^-2) and strong metal-line absorption at z_GRB = 3.0363 with a roughly solar metallicity. We detect a series of A-X bandheads from CO and estimate N(CO) = 10^16.5 cm^-2 and T_ex^CO > 100K. We argue that the high excitation temperature results from UV pumping of the CO gas by the GRB afterglow. Similarly, we observe H2 absorption via the Lyman-Werner bands and estimate N(H2) = 10^21.2 cm^-2 with T_ex^H2 = 10--300K. The afterglow photometry suggests an extinction law with R_V=4 and A_V=3.2 mag and requires the presence of a modest 2175A bump. Additionally, modeling of the Swift/XRT X-ray spectrum confirms a large column density with N(H) = 10^22.58 cm^-2. Remarkably, this molecular gas has extinction properties, metallicity, and a CO/H2 ratio comparable to those of translucent molecular clouds of the Milky Way, suggesting that star formation at high z proceeds in similar environments as today. However, the integrated dust-to-metals ratio is sub-Galactic, suggesting the dust is primarily associated with the molecular phase while the atomic gas has a much lower dust-to-gas ratio. Sightlines like GRB 080607 serve as powerful probes of nucleosynthesis and star-forming regions in the young universe and contribute to the population of "dark" GRB afterglows.

The first gamma-ray burst (GRB) confirmed to be bright enough to be seen with the naked eye, GRB 080319B at redshift z = 0.937, allowed for exquisite follow-up observations across the electromagnetic spectrum. We present our detailed optical and infrared observations of the afterglow, consisting of over 5000 images starting 51 s after the GRB trigger, in concert with our own analysis of the Swift data. The event is extreme not only in observed properties but intrinsically: it was the most luminous event ever recorded at optical and infrared wavelengths and had an exceedingly high isotropic-equivalent energy release in gamma-rays. At early times, the afterglow evolution is broadly consistent with being reverse-shock dominated, but then is subsumed by a forward shock at around 1000 s. The overall spectral energy distribution, spanning from ultraviolet through near-infrared wavelengths, shows no evidence for a significant amount of dust extinction in the host frame. The afterglow evolution, however, is highly chromatic: starting at about 1000 s the index shifts blueward before shifting back to the red at late times. In our deepest late-time observations, we find tentative evidence for an optical jet break and a luminous supernova. Finally, we examine the detectability of such events with current and future facilities and find that such an event could be detected in gamma-rays by BAT out to z = 10.7 (8 sigma), while the nominal EXIST sensitivity would allow detection to z ~ 32. At K band, this source would have been easily detected with meter-class telescopes to z ~ 17.

Stars of late-M and L spectral types, collectively known as Ultracool Dwarfs (UCDs), may be excellent targets for searches for extrasolar planets. Owing to their small radii, the signal from an Earth-size planet transiting a UCD is, in principle, readily detectable. We present results from a study designed to evaluate the feasibility of using precise near infrared (NIR) photometry to detect terrestrial extrasolar planets orbiting UCDs. We used the Peters Automated InfRared Imaging TELescope (PAIRITEL) to observe a sample of 13 UCDs over a period of 10 months. We consider several important systematic effects in NIR differential photometry and develop techniques for generating photometry with a precision of 0.01 mag and long-term stability. We simulate the planet detection efficiency of an extended campaign to monitor a large sample of UCDs with PAIRITEL. We find that both a targeted campaign with a single telescope lasting several years and a campaign making use of a network of telescopes distributed in longitude could provide significant sensitivity to terrestrial planets orbiting UCDs, potentially in the habitable zone.

We report on multi-wavelength observations, ranging from the X-ray to radio wave bands, of the IPN-localized gamma-ray burst GRB 070125. Spectroscopic observations reveal the presence of absorption lines due to O I, Si II, and C IV, implying a likely redshift of z = 1.547. The well-sampled light curves, in particular from 0.5 to 4 days after the burst, suggest a jet break at 3.7 days, corresponding to a jet opening angle of ~7.0 degrees, and implying an intrinsic GRB energy in the 1 - 10,000 keV band of around E = (6.3 - 6.9)x 10^(51) erg (based on the fluences measured by the gamma-ray detectors of the IPN network). GRB 070125 is among the brightest afterglows observed to date. The spectral energy distribution implies a host extinction of Av < 0.9 mag. Two rebrightening episodes are observed, one with excellent time coverage, showing an increase in flux of 56% in ~8000 seconds. The evolution of the afterglow light curve is achromatic at all times. Late-time observations of the afterglow do not show evidence for emission from an underlying host galaxy or supernova. Any host galaxy would be subluminous, consistent with current GRB host-galaxy samples. Evidence for strong Mg II absorption features is not found, which is perhaps surprising in view of the relatively high redshift of this burst and the high likelihood for such features along GRB-selected lines of sight.

Temporal sampling does more than add another axis to the vector of observables. Instead, under the recognition that how objects change (and move) in time speaks directly to the physics underlying astronomical phenomena, next-generation wide-field synoptic surveys are poised to revolutionize our understanding of just about anything that goes bump in the night (which is just about everything at some level). Still, even the most ambitious surveys will require targeted spectroscopic follow-up to fill in the physical details of newly discovered transients. We are now building a new system intended to ingest and classify transient phenomena in near real-time from high-throughput imaging data streams. Described herein, the Transient Classification Project at Berkeley will be making use of classification techniques operating on ``features'' extracted from time series and contextual (static) information. We also highlight the need for a community adoption of a standard representation of astronomical time series data (i.e., ``VOTimeseries'').

We report on observations of a gamma-ray burst (GRB 061126) with an extremely bright (R ~ 12 mag at peak) early-time optical afterglow. The optical afterglow is already fading as a power law 22 seconds after the trigger, with no detectable prompt contribution in our first exposure, which was coincident with a large prompt-emission gamma-ray pulse. The optical--infrared photometric spectral energy distribution is an excellent fit to a power law, but it exhibits a moderate red-to-blue evolution in the spectral index at about 500 s after the burst. This color change is contemporaneous with a switch from a relatively fast decay to slower decay. The rapidly decaying early afterglow is broadly consistent with synchrotron emission from a reverse shock, but a bright forward-shock component predicted by the intermediate- to late-time X-ray observations under the assumptions of standard afterglow models is not observed. Indeed, despite its remarkable early-time brightness, this burst would qualify as a dark burst at later times on the basis of its nearly flat optical-to-X-ray spectral index. Our photometric spectral energy distribution provides no evidence of host-galaxy extinction, requiring either large quantities of grey dust in the host system (at redshift 1.1588 +/- 0.0006, based upon our late-time Keck spectroscopy) or separate physical origins for the X-ray and optical afterglows.

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

The Peters Automated Infrared Imaging Telescope (PAIRITEL) is the first meter-class telescope operating as a fully robotic IR imaging system. Dedicated in October 2004, PAIRITEL began regular observations in mid-December 2004 as part of a 1.5 year commissioning period. The system was designed to respond without human intervention to new gamma-ray burst transients: this milestone was finally reached on November 9, 2005 but the telescope had a number of semi-automated sub-10 minute responses throughout early commissioning. When not operating in Target of Opportunity mode, PAIRITEL performs a number of queue scheduled transient monitoring campaigns. To achieve this level of automation, we have developed communicating tools to connect the various sub-systems: an intelligent queue scheduling database, run-time configurable observation sequence software, a data reduction pipeline, and a master state machine which monitors and controls all functions within and affecting the observatory.

We present lightcurves of the afterglow of GRB050502A, including very early data at t-t_{GRB} < 60s. The lightcurve is composed of unfiltered ROTSE-IIIb optical observations from 44s to 6h post-burst, R-band MDM observations from 1.6 to 8.4h post-burst, and PAIRITEL J H K_s observations from 0.6 to 2.6h post-burst. The optical lightcurve is fit by a broken power law, where t^{alpha} steepens from alpha = -1.13 +- 0.02 to alpha = -1.44 +- 0.02 at \~5700s. This steepening is consistent with the evolution expected for the passage of the cooling frequency nu_c through the optical band. Even in our earliest observation at 44s post-burst, there is no evidence that the optical flux is brighter than a backward extrapolation of the later power law would suggest. The observed decay indices and spectral index are consistent with either an ISM or a Wind fireball model, but slightly favor the ISM interpretation. The expected spectral index in the ISM interpretation is consistent within 1 sigma with the observed spectral index beta = -0.8 +- 0.1; the Wind interpretation would imply a slightly (~2 sigma) shallower spectral index than observed. A small amount of dust extinction at the source redshift could steepen an intrinsic spectrum sufficiently to account for the observed value of beta. In this picture, the early optical decay, with the peak at or below 4.7e14 Hz at 44s, requires very small electron and magnetic energy partitions from the fireball.

The explosion that results in a cosmic gamma-ray burst (GRB) is thought to produce emission from two physical processes -- the activity of the central engine gives rise to the high-energy emission of the burst through internal shocking and the subsequent interaction of the flow with the external environment produces long-wavelength afterglow. While afterglow observations continue to refine our understanding of GRB progenitors and relativistic shocks, gamma-ray observations alone have not yielded a clear picture of the origin of the prompt emission nor details of the central engine. Only one concurrent visible-light transient has been found and was associated with emission from an external shock. Here we report the discovery of infrared (IR) emission contemporaneous with a GRB, beginning 7.2 minutes after the onset of GRB 041219a. Our robotic telescope acquired 21 images during the active phase of the burst, yielding the earliest multi-colour observations of any long-wavelength emission associated with a GRB. Analysis of an initial IR pulse suggests an origin consistent with internal shocks. This opens a new possibility to study the central engine of GRBs with ground-based observations at long wavelengths.

The internet has brought about great change in the astronomical community, but this interconnectivity is just starting to be exploited for use in instrumentation. Utilizing the internet for communicating between distributed astronomical systems is still in its infancy, but it already shows great potential. Here we present an example of a distributed network of telescopes that performs more efficiently in synchronous operation than as individual instruments. RAPid Telescopes for Optical Response (RAPTOR) is a system of telescopes at LANL that has intelligent intercommunication, combined with wide-field optics, temporal monitoring software, and deep-field follow-up capability all working in closed-loop real-time operation. The Telescope ALert Operations Network (TALON) is a network server that allows intercommunication of alert triggers from external and internal resources and controls the distribution of these to each of the telescopes on the network. TALON is designed to grow, allowing any number of telescopes to be linked together and communicate. Coupled with an intelligent alert client at each telescope, it can analyze and respond to each distributed TALON alert based on the telescopes needs and schedule.

The mining of Virtual Observatories (VOs) is becoming a powerful new method for discovery in astronomy. Here we report on the development of SkyDOT (Sky Database for Objects in the Time domain), a new Virtual Observatory, which is dedicated to the study of sky variability. The site will confederate a number of massive variability surveys and enable exploration of the time domain in astronomy. We discuss the architecture of the database and the functionality of the user interface. An important aspect of SkyDOT is that it is continuously updated in near real time so that users can access new observations in a timely manner. The site will also utilize high level machine learning tools that will allow sophisticated mining of the archive. Another key feature is the real time data stream provided by RAPTOR (RAPid Telescopes for Optical Response), a new sky monitoring experiment under construction at Los Alamos National Laboratory (LANL).

Fast variability of optical objects is an interesting though poorly explored subject in modern astronomy. Real-time data processing and identification of transient celestial events in the images is very important for such study as it allows rapid follow-up with more sensitive instruments. We discuss an approach which we have developed for the RAPTOR project, a pioneering closed-loop system combining real-time transient detection with rapid follow-up. RAPTOR's data processing pipeline is able to identify and localize an optical transient within seconds after the observation. The testing we performed so far have been confirming the effectiveness of our method for the optical transient detection. The software pipeline we have developed for RAPTOR can easily be applied to the data from other experiments.

The Rapid Telescope for Optical Response (RAPTOR) program consists of a network of robotic telescopes dedicated to the search for fast optical transients. The pilot project is composed of three observatories separated by approximately 38 kilometers located near Los Alamos, New Mexico. Each of these observatories is composed of a telescope, mount, enclosure, and weather station, all operating robotically to perform individual or coordinated transient searches. The telescopes employ rapidly slewing mounts capable of slewing a 250 pound load 180 degrees in under 2 seconds with arcsecond precision. Each telescope consists of wide-field cameras for transient detection and a narrow-field camera with greater resolution and sensitivity. The telescopes work together by employing a closed-loop system for transient detection and follow-up. Using the combined data from simultaneous observations, transient alerts are generated and distributed via the Internet. Each RAPTOR telescope also has the capability of rapidly responding to external transient alerts received over the Internet from a variety of ground-based and satellite sources. Each observatory may be controlled directly, remotely, or robotically while providing state-of-health and observational results to the client and the other RAPTOR observatories. We discuss the design and implementation of the spatially distributed RAPTOR system.

The Rapid Telescopes for Optical Response (RAPTOR) experiment is a spatially distributed system of autonomous robotic telescopes that is designed to monitor the sky for optical transients. The core of the system is composed of two telescope arrays, separated by 38 kilometers, that stereoscopically view the same 1500 square-degree field with a wide-field imaging array and a central 4 square-degree field with a more sensitive narrow-field "fovea" imager. Coupled to each telescope array is a real-time data analysis pipeline that is designed to identify interesting transients on timescales of seconds and, when a celestial transient is identified, to command the rapidly slewing robotic mounts to point the narrow-field ``fovea'' imagers at the transient. The two narrow-field telescopes then image the transient with higher spatial resolution and at a faster cadence to gather light curve information. Each "fovea" camera also images the transient through a different filter to provide color information. This stereoscopic monitoring array is supplemented by a rapidly slewing telescope with a low resolution spectrograph for follow-up observations of transients and a sky patrol telescope that nightly monitors about 10,000 square-degrees for variations, with timescales of a day or longer, to a depth about 100 times fainter. In addition to searching for fast transients, we will use the data stream from RAPTOR as a real-time sentinel for recognizing important variations in known sources. Altogether, the RAPTOR project aims to construct a new type of system for discovery in optical astronomy--one that explores the time domain by "mining the sky in real time".