Cosmic Ray Rejection by Linear Filtering of Single Images

We present a convolution-based algorithm for finding cosmic rays in single well-sampled astronomical images. The spatial filter used is the point spread function (approximated by a Gaussian) minus a scaled delta function, and cosmic rays are identified by thresholding the filtered image. This filter searches for features with significant power at spatial frequencies too high for legitimate objects. Noise properties of the filtered image are readily calculated, which allows us to compute the probability of rejecting a pixel not contaminated by a cosmic ray (the false alarm probability). We demonstrate that the false alarm probability for a pixel containing object flux will never exceed the corresponding probability for a blank sky pixel, provided we choose the convolution kernel appropriately. This allows confident rejection of cosmic rays superposed on real objects. Identification of multiple-pixel cosmic ray hits can be enhanced by running the algorithm iteratively, replacing flagged pixels with the background level at each iteration.Comment: Accepted for publication in PASP (May 2000 issue). An iraf script implementing the algorithm is available from the author, or from http://sol.stsci.edu/~rhoads/ . 16 pages including 3 figures. Uses AASTeX aaspp4 styl

The Dynamics and Light Curves of Beamed Gamma Ray Burst Afterglows

The energy requirements of gamma ray bursts have in past been poorly constrained because of three major uncertainties: The distances to bursts, the degree of burst beaming, and the efficiency of gamma ray production. The first of these has been resolved, with both indirect evidence (the distribution of bursts in flux and position) and direct evidence (redshifted absorption features in the afterglow spectrum of GRB 970508) pointing to cosmological distances. We now wish to address the second uncertainty. Afterglows allow a statistical test of beaming, described in an earlier paper. In this paper, we modify a standard fireball afterglow model to explore the effects of beaming on burst remnant dynamics and afterglow emission. If the burst ejecta are beamed into angle zeta, the burst remnant's evolution changes qualitatively once its bulk Lorentz factor Gamma < 1/zeta: Before this, Gamma declines as a power law of radius, while afterwards, it declines exponentially. This change results in a broken power law light curve whose late-time decay is faster than expected for a purely spherical geometry. These predictions disagree with afterglow observations of GRB 970508. We explored several variations on our model, but none seems able to change this result. We therefore suggest that this burst is unlikely to have been highly beamed, and that its energy requirements were near those of isotropic models. More recent afterglows may offer the first practical applications for our beamed models.Comment: 18 pages, uses emulateapj.sty, four embedded postscript figures. Submitted to The Astrophysical Journal, 199

The Near Infrared and Multiwavelength Afterglow of GRB 000301c

We present near-infrared observations of the counterpart of GRB 000301c. The K' filter (2.1 micron) light curve shows a well-sampled break in the decay slope at t=3.5 days post-burst. The early time slope is very shallow (~ -0.1), while the late time slope is steep (-2.2). Comparison with the optical (R band) light curve shows marginally significant differences, especially in the early time decay slope (which is steeper in the optical) and the break time (which occurs later in the optical). This is contrary to the general expectation that light curve breaks should either be achromatic (e.g., for breaks due to collimation effects) or should occur later at longer wavelengths (for most other breaks). The observed color variations might be intrinsic to the afterglow, or might indicate systematic errors of > 0.08 magnitude in all fluxes. Even if the break is achromatic, we argue that its sharpness poses difficulties for explanations that depend on collimated ejecta. The R light curve shows further signs of fairly rapid variability (a bump, steep drop, and plateau) that are not apparent in the K' light curve. In addition, by combining the IR-optical-UV data with millimeter and radio fluxes, we are able to constrain the locations of the self-absorption break and cooling break and to infer the location of the spectral peak at t=3 days: f_nu = 3.4 mJy at nu=1e12 Hz. Using the multiwavelength spectral energy distribution, we are able to constrain the blast wave energy, which was E > 3e53 erg if the explosion was isotropic. This implies a maximum gamma ray production efficiency of ~ 0.15 for GRB 000301C.Comment: Accepted to The Astrophysical Journal. 24 pages, 4 figures, 3 tables; uses AASTeX 5 macros. This version includes a new figure (R-K' color vs. time), a better sampled R band light curve, and more extensive discussion of the optical data and error analysi

GRB Energetics and the GRB Hubble Diagram: Promises and Limitations

We present a complete sample of 29 GRBs for which it has been possible to determine temporal breaks (or limits) from their afterglow light curves. We interpret these breaks within the framework of the uniform conical jet model, incorporating realistic estimates of the ambient density and propagating error estimates on the measured quantities. In agreement with our previous analysis of a smaller sample, the derived jet opening angles of those 16 bursts with redshifts result in a narrow clustering of geometrically-corrected gamma-ray energies about E_gamma = 1.33e51 erg; the burst-to-burst variance about this value is a factor of 2.2. Despite this rather small scatter, we demonstrate in a series of GRB Hubble diagrams, that the current sample cannot place meaningful constraints upon the fundamental parameters of the Universe. Indeed for GRBs to ever be useful in cosmographic measurements we argue the necessity of two directions. First, GRB Hubble diagrams should be based upon fundamental physical quantities such as energy, rather than empirically-derived and physically ill-understood distance indicators. Second, a more homogeneous set should be constructed by culling sub-classes from the larger sample. These sub-classes, though now first recognizable by deviant energies, ultimately must be identifiable by properties other than those directly related to energy. We identify a new sub-class of GRBs (``f-GRBs'') which appear both underluminous by factors of at least 10 and exhibit a rapid fading at early times. About 10-20% of observed long-duration bursts appear to be f-GRBs.Comment: Accepted to the Astrophysical Journal (20 May 2003). 19 pages, 3 Postscript figure

On the energy of gamma-ray bursts

We show that gamma-ray burst (GRB) afterglow observations strongly suggest, within the fireball model framework, that radiating electrons are shock accelerated to a power-law energy distribution, with universal index p \approx 2.2, and that the fraction of shock energy carried by electrons, \xi_e, is universal and close to equipartition, \xi_e ~ 1/3. For universal p and \xi_e, a single measurement of the X-ray afterglow flux on the time scale of a day provides a robust estimate of the fireball energy per unit solid angle, \epsilon, averaged over a conical section of the fireball of opening angle \theta ~ 0.1. Applying our analysis to BeppoSAX afterglow data we find that: (i) Fireball energies are in the range of 4\pi\epsilon=10^{51.5} to 10^{53.5} erg; (ii) The ratio of observed $\gamma$-ray to total fireball energy per unit solid angle, \epsilon_\gamma / \epsilon, is of order unity, satisfying abs[log10(\epsilon_\gamma/\epsilon)]<0.5; (iii) If fireballs are jet like, their opening angle should satisfy \theta>=0.1. Our results imply that if typical opening angles are \theta ~ 0.1, a value consistent with our analysis, the total energy associated with a GRB event is in the range of 10^{50} erg to 10^{51.5} erg.Comment: 16 pages; Submitted to Ap

First Results from the Large Area Lyman Alpha Survey

We report on a new survey for z=4.5 Lyman alpha sources, the Large Area Lyman Alpha (LALA) survey. Our survey achieves an unprecedented combination of volume and sensitivity by using narrow-band filters on the new 8192x8192 pixel CCD Mosaic Camera at the 4 meter Mayall telescope of Kitt Peak National Observatory. Well-detected sources with flux and equivalent width matching known high redshift Lyman alpha galaxies (i.e., observed equivalent width above 80 Angstroms and line+continuum flux between 2.6e-17 and 5.2e-17 erg/cm^2/sec in an 80 Angstrom filter) have an observed surface density corresponding to 11000 +- 700 per square degree per unit redshift at z=4.5. Spatial variation in this surface density is apparent on comparison between counts in 6561 and 6730 Angstrom filters. Early spectroscopic followup results from the Keck telescope included three sources meeting our criteria for good Lyman alpha candidates. Of these, one is confirmed as a z=4.52 source, while another remains consistent with either z=4.55 or z=0.81. We infer that 30 to 50% of our good candidates are bona fide Lyman alpha emitters, implying a net density of about 4000 Lyman alpha galaxies per square degree per unit redshift.Comment: 10 pages, 2 figures (3 .ps files), uses AASTeX 4. Submitted to The Astrophysical Journal Letter

Stars in the Hubble Ultra Deep Field

We identified 46 unresolved source candidates in the Hubble Ultra Deep Field, down to i775 = 29.5. Unresolved objects were identified using a parameter S, which measures the deviation from the curve-of-growth of a point source. Extensive testing of this parameter was carried out, including the effects of decreasing signal-to-noise and of the apparent motions of stars, which demonstrated that stars brighter than i775 = 27.0 could be robustly identified. Low resolution grism spectra of the 28 objects brighter than i775 = 27.0 identify 18 M and later stellar type dwarfs, 2 candidate L-dwarfs, 2 QSOs, and 4 white dwarfs. Using the observed population of dwarfs with spectral type M4 or later, we derive a Galactic disk scale height of 400 \pm 100 pc for M and L stars. The local white dwarf density is computed to be as high as (1.1 \pm 0.3) x10^(-2) stars/pc^3. Based on observations taken 73 days apart, we determined that no object in the field has a proper motion larger than 0.027"/year (3 sigma detection limit). No high velocity white dwarfs were identified in the HUDF, and all four candidates appear more likely to be part of the Galactic thick disk. The lack of detected halo white dwarfs implies that, if the dark matter halo is 12 Gyr old, white dwarfs account for less than 10% of the dark matter halo mass.Comment: 35 pages, 11 figures, accepted by Ap

How to Tell a Jet from a Balloon: A Proposed Test for Beaming in Gamma Ray Bursts

If gamma ray bursts are highly collimated, the energy requirements of each event may be reduced by several (~ 4-6) orders of magnitude, and the event rate increased correspondingly. Extreme conditions in gamma ray bursters lead to highly relativistic motions (bulk Lorentz factors Gamma > 100). This results in strong forward beaming of the emitted radiation in the observer's rest frame. Thus, all information on gamma ray bursts comes from those ejecta emitted in a narrow cone (opening angle 1/Gamma) pointing towards the observer. We are at present ignorant of whether there are ejecta outside that cone or not. The recent detection of longer wavelength transients following gamma ray bursts allows an empirical test of whether gamma ray bursts are collimated jets or spherical fireballs. The bulk Lorentz factor of the burst ejecta will decrease with time after the event, as the ejecta sweep up the surrounding medium. Thus, radiation from the ejecta is beamed into an ever increasing solid angle as the burst remnant evolves. It follows that if gamma ray bursts are highly collimated, many more optical and radio transients should be observed without associated gamma rays than with them. Published supernova searches may contain enough data to test the most extreme models of gamma ray beaming. We close with a brief discussion of other possible consequences of beaming, including its effect on the evolution of burst remnants.Comment: Original replaced with accepted refereed manuscript. 11 pages, uses AASTeX 4.0 LaTeX macros. To be published in The Astrophysical Journal Letters, vol. 487, p. L1 (20 September 1997

Microlensing of Globular Clusters as a Probe of Galactic Structure

The spatial distribution of compact dark matter in our Galaxy can be determined in a few years of monitoring Galactic globular clusters for microlensing. Globular clusters are the only dense fields of stars distributed throughout the three-dimensional halo and hence are uniquely suited to probe its structure. The microlensing optical depths towards different clusters have varying contributions from the thin disk, thick disk, bulge, and halo of the Galaxy. Although measuring individual optical depths to all the clusters is a daunting task, we show that interesting Galactic structure information can be extracted with as few as $40$--$120$ events in total for the entire globular cluster system (observable with 2--5 years of monitoring). The globular cluster microlensing is particularly sensitive to the core radius of the halo mass distribution and to the scale length, surface mass density, and radial scale height variations of the thin disk.Comment: 14 pages, 1 figure. Submitted to ApJ Letters. Uses aastex macro

Variable polarization in the optical afterglow of GRB 021004

We present polarimetric observations of the afterglow of gamma-ray burst (GRB) 021004, obtained with the Nordic Optical Telescope (NOT) and the Very Large Telescope (VLT) between 8 and 17 hours after the burst. Comparison among the observations shows a 45 degree change in the position angle from 9 hours after the burst to 16 hours after the burst, and comparison with published data from later epochs even shows a 90 degree change between 9 and 89 hours after the burst. The degree of linear polarization shows a marginal change, but is also consistent with being constant in time. In the context of currently available models for changes in the polarization of GRBs, a homogeneous jet with an early break time of t_b ~ 1 day provides a good explanation of our data. The break time is a factor 2 to 6 earlier than has been found from the analysis of the optical light curve. The change in the position angle of the polarization rules out a structured jet model for the GRB.Comment: 5 pages, 2 figures. Published in A&A letter