Source author record

Nathaniel R. Butler

Nathaniel R. Butler appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

28 published item(s)

preprint2020arXiv

Limits on the Electromagnetic Counterpart to S190814bv

We derive limits on any electromagnetic counterpart to the compact binary merger S190814bv, whose parameters are consistent with the merger of a black hole and a neutron star. We present observations with the new wide-field optical imager DDOTI and also consider Swift/BAT observations reported by Palmer et al. (2019). We show that Swift/BAT would have detected a counterpart with similar properties to a typical on-axis short GRB at the 98 per cent confidence level, whereas our DDOTI observations only rule out such a counterpart at the 27 per cent confidence level. Neither have sufficient sensitivity to rule out an off-axis counterpart like GW 170817. We compare the efficiency of Swift/BAT and DDOTI for future observations, and show that DDOTI is likely to be about twice as efficient as Swift/BAT for off-axis events up to about 100 Mpc.

preprint2020arXiv

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

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

preprint2020arXiv

The Importance Of Star Formation Intensity In LYα Escape From Green Pea Galaxies And Lyman Break Galaxy Analogs

We have studied ultraviolet images of 40 Green Pea galaxies and 15 local Lyman Break Galaxy Analogs to understand the relation between Ly$α$ photon escape and central UV photometric properties. We measured star formation intensity (SFI, star formation rate per unit area) from the central 250 pc region ($S_{\rm 250pc}$) using COS/NUV images from the \textit{Hubble Space Telescope}. The measured $S_{\rm 250pc}$ of our sample Green Peas ranges from 2.3--46 $M_{\odot} \ \rm{year}^{-1} \ \rm{kpc^{-2}}$, with a geometric mean of $15 M_{\odot} \ \rm{year}^{-1} \ \rm{kpc^{-2}}$ and a standard deviation of 0.266 dex, forming a relatively narrow distribution. The Lyman Break Galaxy Analogs show a similarly narrow distribution of $S_{\rm 250pc}$ (0.271 dex), though with a larger mean of 28 $M_{\odot} \ \rm{year}^{-1} \ \rm{kpc^{-2}}$. We show that while the Ly$α$ equivalent width (EW(Ly$α$)) and the Ly$α$ escape fraction ($f^{Lyα}_{esc}$) are not significantly correlated with the central SFI ($S_{\rm 250pc}$), both are positively correlated with the ratio of surface brightness to galaxy stellar mass ($S_{\rm 250pc}/M_{\rm star}$), with correlation coefficients ($p$-values) of 0.702 ($1\times 10^{-8}$) and 0.529 ($5\times 10^{-4}$) with EW(Ly$α$) and $f^{Lyα}_{esc}$, respectively. These correlations suggest a scenario where intense central star formation can drive a galactic wind in galaxies with relatively shallow gravitational potential wells, thus clearing channels for the escape of Ly$α$ photons.

preprint2019arXiv

GRB 180620A: Evidence for late-time energy injection

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

preprint2016arXiv

An Ultraviolet Spectrum of the Tidal Disruption Flare ASASSN-14li

We present a Hubble Space Telescope STIS spectrum of ASASSN-14li, the first rest-frame UV spectrum of a tidal disruption flare (TDF). The underlying continuum is well fit by a blackbody with $T_{\mathrm{UV}} = 3.5 \times 10^{4}$ K, an order of magnitude smaller than the temperature inferred from X-ray spectra (and significantly more precise than previous efforts based on optical and near-UV photometry). Super-imposed on this blue continuum, we detect three classes of features: narrow absorption from the Milky Way (probably a high-velocity cloud), and narrow absorption and broad (FWHM $\approx 2000$-8000 km s$^{-1}$) emission lines at/near the systemic host velocity. The absorption lines are blueshifted with respect to the emission lines by $Δv = -(250$-400) km s$^{-1}$. Due both to this velocity offset and the lack of common low-ionization features (Mg II, Fe II), we argue these arise from the same absorbing material responsible for the low-velocity outflow discovered at X-ray wavelengths. The broad nuclear emission lines display a remarkable abundance pattern: N III], N IV], He II are quite prominent, while the common quasar emission lines of C III] and Mg II are weak or entirely absent. Detailed modeling of this spectrum will help elucidate fundamental questions regarding the nature of the emission process(es) at work in TDFs, while future UV spectroscopy of ASASSN-14li would help to confirm (or refute) the previously proposed connection between TDFs and "N-rich" quasars.

preprint2016arXiv

SDSS J0159+0105: A Radio-Quiet Quasar with a Centi-Parsec Supermassive Black Hole Binary Candidate

We report a candidate centi-parsec supermassive black hole binary (SMBHB) in the radio-quiet quasar SDSS J0159+0105 at z=0.217. With a modified lomb-scargle code GLSdeDRW and the auto-correlation analysis ACF, we detect two significant (at P>99%) periodic signals at ~741 day and ~1500 day from the 8.1-year Catalina V-band light curve of this quasar. The period ratio, which is close to 1:2, is typical of a black-hole binary system with a mass ratio of 0.05<q<0.8 according to recent numerical simulations. SDSS J0159+0105 has two SDSS spectroscopic observations separated by ~10 years. There is a significant change in the broad H-beta profile between the two epochs, which can be explained by a single broad-line region (BLR) around the binary system illuminated by the aforementioned mini-disks, or a stream of gas flowing from the circumbinary disk to one of the SMBHs. From the single BLR assumption and the orbital period t_orb ~1500 day, we estimate the total virial masses of M_SMBHB ~ 1.3x10^8 M_sun, the average distances of BLR of ~0.04pc (~50 light-day, with +/-0.3 dex uncertainty), and a SMBHB separation of d= (0.01pc)M_{8,tot}^{1/3} (T_rest/3.3yr)^{2/3} ~ 0.013 pc (15 light-day). Based on analytical work, the postulated circumbinary disk has an inner radius of 2d = 0.026 pc (30 light-day). SDSS J0159+0105 also displays unusual spectral energy distribution. The unique properties of SDSS J0159+0105 are consistent with it being a centi-parsec SMBHB. The GLSdeDRW code link: http://butler.lab.asu.edu/qso_period/

preprint2015arXiv

The Energy-Dependence of GRB Minimum Variability Timescales

We constrain the minimum variability timescales for 938 GRBs observed by the Fermi/GBM instrument prior to July 11, 2012. The tightest constraints on progenitor radii derived from these timescales are obtained from light curves in the hardest energy channel. In the softer bands -- or from measurements of the same GRBs in the hard X-rays from Swift -- we show that variability timescales tend to be a factor 2--3 longer. Applying a survival analysis to account for detections and upper limits, we find median minimum timescale in the rest frame for long-duration and short-duration GRBs of 45 ms and 10 ms, respectively. Fewer than 10% of GRBs show evidence for variability on timescales below 2 ms. These shortest timescales require Lorentz factors $\gtrsim 400$ and imply typical emission radii $R \approx 1 {\times} 10^{14}$ cm for long-duration GRBs and $R \approx 3 {\times} 10^{13}$ cm for short-duration GRBs. We discuss implications for the GRB fireball model and investigate whether GRB minimum timescales evolve with cosmic time.

preprint2014arXiv

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

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

preprint2014arXiv

The Transiting Exoplanet Survey Satellite

The Transiting Exoplanet Survey Satellite (TESS) will search for planets transiting bright and nearby stars. TESS has been selected by NASA for launch in 2017 as an Astrophysics Explorer mission. The spacecraft will be placed into a highly elliptical 13.7-day orbit around the Earth. During its two-year mission, TESS will employ four wide-field optical CCD cameras to monitor at least 200,000 main-sequence dwarf stars with I = 4-13 for temporary drops in brightness caused by planetary transits. Each star will be observed for an interval ranging from one month to one year, depending mainly on the star's ecliptic latitude. The longest observing intervals will be for stars near the ecliptic poles, which are the optimal locations for follow-up observations with the James Webb Space Telescope. Brightness measurements of preselected target stars will be recorded every 2 min, and full frame images will be recorded every 30 min. TESS stars will be 10-100 times brighter than those surveyed by the pioneering Kepler mission. This will make TESS planets easier to characterize with follow-up observations. TESS is expected to find more than a thousand planets smaller than Neptune, including dozens that are comparable in size to the Earth. Public data releases will occur every four months, inviting immediate community-wide efforts to study the new planets. The TESS legacy will be a catalog of the nearest and brightest stars hosting transiting planets, which will endure as highly favorable targets for detailed investigations.

preprint2014arXiv

Uncovering the Intrinsic Variability of Gamma-ray Bursts

We develop a robust technique to determine the minimum variability timescale for Gamma-ray Burst light curves, utilizing Haar wavelets. Our approach averages over the data for a given GRB, providing an aggregate measure of signal variation while also retaining sensitivity to narrow pulses within complicated time-series. In contrast to previous studies using wavelets, which simply define the minimum timescale in reference to the measurement noise floor, our approach identifies the signature of temporally-smooth features in the wavelet scaleogram and then additionally identifies a break in the scaleogram on longer timescales as signature of a true, temporally-unsmooth light curve feature or features. We apply our technique to the large sample of Swift GRB Gamma-ray light curves and for the first time -- due to the presence of a large number of GRBs with measured redshift -- determine the distribution of minimum variability timescales in the source frame. We find a median minimum timescale for long-duration GRBs in the source frame of $Δt_{\rm min}=0.5$ s, with the shortest timescale found being on the order of 10 ms. This short timescale suggests a compact central engine ($3 \times 10^3$ km). We discuss further implications for the GRB fireball model and present a tantalizing correlation between minimum timescale and redshift, which may in part be due to cosmological time-dilation.

preprint2013arXiv

A Tidal Disruption Event in a Nearby Galaxy Hosting an Intermediate Mass Black Hole

We report the serendipitous discovery of a bright point source flare in the Abell cluster 1795 with archival EUVE and Chandra observations. Assuming the EUVE emission is associated with the Chandra source, the X-ray 0.5-7 keV flux declined by a factor of ~2300 over a time span of 6 years, following a power-law decay with index ~2.44+-0.40. The Chandra data alone vary by a factor of ~20. The spectrum is well fit by a blackbody with a constant temperature of kT~0.09 keV (~10^6 K). The flare is spatially coincident with the nuclear region of a faint, inactive galaxy with a photometric redshift consistent at the one sigma level with the cluster (z=0.062476). We argue that these properties are indicative of a tidal disruption of a star by a black hole with log(M_BH/M_sun)~5.5+-0.5. If so, such a discovery indicates that tidal disruption flares may be used to probe black holes in the intermediate mass range, which are very difficult to study by other means.

preprint2012arXiv

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

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

preprint2012arXiv

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

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

preprint2011arXiv

A Compact Degenerate Primary-Star Progenitor of SN 2011fe

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

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2011arXiv

Optimal Time-Series Selection of Quasars

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

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2010arXiv

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

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

preprint2010arXiv

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

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

preprint2010arXiv

The long rapid decay phase of the extended emission from the short GRB 080503

GRB080503 was classified as a short GRB with extended emission (Perley et al. 2009). The origin of such extended emission (found in about a quarter of Swift short GRBs) is still unclear and may provide some clues to the identity of the elusive progenitors of short GRBs. The extended emission from GRB 080503 is followed by a rapid decay phase (RDP) that is detected over an unusually large dynamical range (one decade in time and ~3.5 decades in flux). We model the broad envelope of extended emission and the subsequent RDP using a physical model (Genet & Granot 2009), in which the prompt emission (and its tail) is the sum of its individual pulses (and their tails). For GRB 080503, a single pulse fit is found to be unacceptable. The RDP displays very strong spectral evolution and shows some evidence for the presence of two spectral components with different temporal behaviour, likely arising from distinct physical regions. A two pulse fit provides a much better fit to the data. The shallow gamma-ray and steep hard X-ray decays are hard to account for simultaneously, and require the second pulse to deviate from the simplest version of the model we use. Therefore, while high latitude emission is a viable explanation for the RDP in GRB080503, it is quite plausible that another mechanism is at work here. Finally, we note that the properties of the RDP following the extended emission of short GRBs appear to have different properties than that following the prompt emission of long GRBs. However, a larger sample of short GRBs with extended emission is required before any strong conclusion can be drawn.

preprint2009arXiv

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

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

preprint2009arXiv

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

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

preprint2007arXiv

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

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