Identifying the sources of the highest energy cosmic rays requires understanding how they are deflected by the stochastic, spatially intermittent intergalactic magnetic field. Here we report measurements of energetic charged-particle propagation through a laser-produced magnetized plasma with these properties. We characterize the diffusive transport of the particles experimentally. The results show that the transport is diffusive and that, for the regime of interest for the highest-energy cosmic rays, the diffusion coefficient is unaffected by the spatial intermittency of the magnetic field.

Magnetic fields are ubiquitous in the Universe. Extragalactic disks, halos and clusters have consistently been shown, via diffuse radio-synchrotron emission and Faraday rotation measurements, to exhibit magnetic field strengths ranging from a few nG to tens of $\mu$G. The energy density of these fields is typically comparable to the energy density of the fluid motions of the plasma in which they are embedded, making magnetic fields essential players in the dynamics of the luminous matter. The standard theoretical model for the origin of these strong magnetic fields is through the amplification of tiny seed fields via turbulent dynamo to the level consistent with current observations. Here we demonstrate, using laser-produced colliding plasma flows, that turbulence is indeed capable of rapidly amplifying seed fields to near equipartition with the turbulent fluid motions. These results support the notion that turbulent dynamo is a viable mechanism responsible for the observed present-day magnetization of the Universe.

The universe is permeated by magnetic fields, with strengths ranging from a femtogauss in the voids between the filaments of galaxy clusters to several teragauss in black holes and neutron stars. The standard model behind cosmological magnetic fields is the nonlinear amplification of seed fields via turbulent dynamo to the values observed. We have conceived experiments that aim to demonstrate and study the turbulent dynamo mechanism in the laboratory. Here we describe the design of these experiments through simulation campaigns using FLASH, a highly capable radiation magnetohydrodynamics code that we have developed, and large-scale three-dimensional simulations on the Mira supercomputer at Argonne National Laboratory. The simulation results indicate that the experimental platform may be capable of reaching a turbulent plasma state and study dynamo amplification. We validate and compare our numerical results with a small subset of experimental data using synthetic diagnostics.

Long-lived rotational and meridional flows are important ingredients of the solar cycle. Magnetic field images have typically been used to measure these flows on the solar surface by cross-correlating thin longitudinal strips or square patches across sufficiently long time gaps. Here, I use one month of SDO/HMI line-of-sight magnetic field observations, combined with the SWAMIS magnetic feature tracking algorithm to measure the motion of individual features in these magnetograms. By controlling for perturbations due to short-lived flows and due to false motions from feature interactions, I effectively isolate the long-lived flows traced by the magnetic features. This allows me to produce high-fidelity differential rotation measurements with well-characterized variances and covariances of the fit parameters.I find a sidereal rotational profile of $(14.296\pm0.006)+(-1.847\pm0.056)\sin^{2}b+(-2.615\pm0.093)\sin^{4}b$, with units of $\textrm{ deg d}^{-1}$, and a large covariance $\sigma_{BC}^{2}=-4.87\times10^{-3}(\textrm{ deg d}^{-1})^{2}$. I also produce medium-fidelity measurements of the much weaker meridional flow that is broadly consistent with previous results. This measurement shows a peak flow of $16.7\pm0.6\text{ m s}^{-1}$ at latitude $b=45^\circ$ but is insufficiently characterized at higher latitudes to ascertain whether the chosen functional form $2\cos b\sin b$ is appropriate. This work shows that measuring the motions of individual features in photospheric magnetograms can produce high precision results in relatively short time spans, and suggests that high resolution non-longitudinally averaged photospheric velocity residual measurements could be produced to compare with coronal results, and to provide other diagnostics of the solar dynamo.

We propose a new way of generating magnetized supersonic jets using a ring laser to irradiate a flat surface target. Using 2D FLASH code simulations which include the Biermann Battery term, we demonstrate that strong toroidal fields can be generated and sustained downstream in the collimated jet outflow far from the target surface. The field strength can be controlled by varying the ring laser separation, thereby providing a versatile laboratory platform for studying the effects of magnetic field in a variety of astrophysical settings.

The Biermann Battery effect is frequently invoked in cosmic magnetogenesis and studied in High-Energy Density laboratory physics experiments. Generation of magnetic fields by the Biermann effect due to mis-aligned density and temperature gradients in smooth flow <i>behind</i> shocks is well known. We show that a Biermann-effect magnetic field is also generated <i>within</i> shocks. Direct implementation of the Biermann effect in MHD codes does not capture this physical process, and worse, produces unphysical magnetic fields at shocks whose value does not converge with resolution. We show that this convergence breakdown is due to naive discretization, which fails to account for the fact that discretized irrotational vector fields have spurious solenoidal components that grow without bound near a discontinuity. We show that careful consideration of the kinetics of ion viscous shocks leads to a formulation of the Biermann effect that gives rise to a convergent algorithm. We note two novel physical effects: a <i>resistive magnetic precursor</i> in which Biermann-generated field in the shock "leaks" resistively upstream; and a <i>thermal magnetic precursor</i>, in which field is generated by the Biermann effect ahead of the shock front due to gradients created by the shock's electron thermal conduction precursor. Both effects appear to be potentially observable in experiments at laser facilities. We re-examine published studies of magnetogenesis in galaxy cluster formation, and conclude that the simulations in question had inadequate resolution to reliably estimate the field generation rate. Corrected estimates suggest primordial field values in the range $B\sim 10^{-22}$G --- $10^{-19}$G by $z=3$.

We present a set of tools for detecting small-scale solar magnetic cancellations and the disk counterpart of type II spicules (the so-called Rapid Blueshifted Excursions (RBEs)), using line-of-sight photospheric magnetograms and chromospheric spectroscopic observations, respectively. For tracking magnetic cancellation, we improve the Southwest Automatic Magnetic Identification Suite (SWAMIS) so that it is able to detect certain obscure cancellations that can be easily missed. For detecting RBEs, we use a normalized reference profile to reduce false-positive detections caused by the non-uniform background and seeing condition. Similar to the magnetic feature tracking in SWAMIS, we apply a dual-threshold method to enhance the accuracy of RBE detection. These tools are employed to analyze our coordinated observations using the Interferometric BIdimensional Spectrometer at Dunn Solar Telescope (DST) of the National Solar Observatory (NSO) and Hinode. We present the statistical properties of magnetic cancellations and RBEs, and explore their correlation using this data set.

Chromospheric rapid blueshifted excursions (RBEs) are suggested to be the disk counterparts of type II spicules at the limb and believed to contribute to the coronal heating process. Previous identification of RBEs was mainly based on feature detection using Dopplergrams. In this paper, we study RBEs on 2011 October 21 in a very quiet region at the disk center, which were observed with the high-cadence imaging spectroscopy of the Ca II 8542 A line from the Interferometric Bidimensional Spectrometer (IBIS). By using an automatic spectral analysis algorithm, a total of 98 RBEs are identified during a 11 minute period. Most of these RBEs have either a round or elongated shape, with an average area of 1.2 arcsec^2. The detailed temporal evolution of spectra from IBIS makes possible a quantitative determination of the velocity (~16 km/s) and acceleration (~400 m/s^2) of Ca II 8542 RBEs, and reveal an additional deceleration (~-160 m/s^2) phase that usually follows the initial acceleration. In addition, we also investigate the association of RBEs with the concomitant photospheric magnetic field evolution, using coordinated high-resolution and high-sensitivity magnetograms made by Hinode. Clear examples are found where RBEs appear to be associated with the preceding magnetic flux emergence and/or the subsequent flux cancellation. However, a further analysis with the aid of the Southwest Automatic Magnetic Identification Suite does not yield a significant statistical association between these RBEs and magnetic field evolution. We discuss the implications of our results in the context of understanding the driving mechanism of RBEs.

We present a systematic study of the pure deflagration model of Type Ia supernovae using three-dimensional, high-resolution, full-star hydrodynamical simulations, nucleosynthetic yields calculated using Lagrangian tracer particles, and light curves calculated using radiation transport. We evaluate the simulations by comparing their predicted light curves with many observed SNe Ia using the SALT2 data-driven model and find that the simulations may correspond to under-luminous SNe Iax. We explore the effects of the initial conditions on our results by varying the number of randomly selected ignition points from 63 to 3500, and the radius of the centered sphere they are confined in from 128 to 384 km. We find that the rate of nuclear burning depends on the number of ignition points at early times, the density of ignition points at intermediate times, and the radius of the confining sphere at late times. The results depend primarily on the number of ignition points, but we do not expect this to be the case in general. The simulations with few ignition points release more nuclear energy $E_{\mathrm{nuc}}$, have larger kinetic energies $E_{\mathrm{K}}$, and produce more $^{56}$Ni than those with many ignition points, and differ in the distribution of $^{56}$Ni, Si, and C/O in the ejecta. For these reasons, the simulations with few ignition points exhibit higher peak B-band absolute magnitudes $M_\mathrm{B}$ and light curves that rise and decline more quickly; their $M_\mathrm{B}$ and light curves resemble those of under-luminous SNe Iax, while those for simulations with many ignition points are not.

We explore the nature of the small-scale solar dynamo by tracking magnetic features. We investigate two previously-explored categories of the small-scale solar dynamo: shallow and deep. Recent modeling work on the shallow dynamo has produced a number of scenarios for how a strong network concentration can influence the formation and polarity of nearby small-scale magnetic features. These scenarios have measurable signatures, which we test for here using magnetograms from the Narrowband Filter Imager (NFI) on Hinode. We find no statistical tendency for newly-formed magnetic features to cluster around or away from network concentrations, nor do we find any statistical relationship between their polarities. We conclude that there is no shallow or "surface" dynamo on the spatial scales observable by Hinode/NFI. In light of these results, we offer a scenario in which the sub-surface field in a deep solar dynamo is stretched and distorted via turbulence, allowing the field to emerge at random locations on the photosphere.

Interaction between a central outflow and a surrounding wind is common in astrophysical sources powered by accretion. Understanding how the interaction might help to collimate the inner central outflow is of interest for assessing astrophysical jet formation paradigms. In this context, we studied the interaction between two nested supersonic plasma flows generated by focusing a long pulse high-energy laser beam onto a solid target. A nested geometry was created by shaping the energy distribution at the focal spot with a dedicated phase plate. Optical and X-ray diagnostics were used to study the interacting flows. Experimental results and numerical hydrodynamic simulations indeed show the formation of strongly collimated jets. Our work experimentally confirms the "shock-focused inertial confinement" mechanism proposed in previous theoretical astrophysics investigations.

This document presents (off-line) computing requrements and challenges for Cosmic Frontier science, covering the areas of data management, analysis, and simulations. We invite contributions to extend the range of covered topics and to enhance the current descriptions.

We explore Implicit Monte Carlo (IMC) and Discrete Diffusion Monte Carlo (DDMC) for radiation transport in high-velocity outflows with structured opacity. The IMC method is a stochastic computational technique for nonlinear radiation transport. IMC is partially implicit in time and may suffer in efficiency when tracking Monte Carlo particles through optically thick materials. DDMC accelerates IMC in diffusive domains. Abdikamalov extended IMC and DDMC to multigroup, velocity-dependent transport with the intent of modeling neutrino dynamics in core-collapse supernovae. Densmore has also formulated a multifrequency extension to the originally grey DDMC method. We rigorously formulate IMC and DDMC over a high-velocity Lagrangian grid for possible application to photon transport in the post-explosion phase of Type Ia supernovae. This formulation includes an analysis that yields an additional factor in the standard IMC-to-DDMC spatial interface condition. To our knowledge the new boundary condition is distinct from others presented in prior DDMC literature. The method is suitable for a variety of opacity distributions and may be applied to semi-relativistic radiation transport in simple fluids and geometries. Additionally, we test the code, called SuperNu, using an analytic solution having static material, as well as with a manufactured solution for moving material with structured opacities. Finally, we demonstrate with a simple source and 10 group logarithmic wavelength grid that IMC-DDMC performs better than pure IMC in terms of accuracy and speed when there are large disparities between the magnitudes of opacities in adjacent groups. We also present and test our implementation of the new boundary condition.

We propose a robust, quantitative method to compare the synthetic light curves of a Type Ia Supernova (SNIa) explosion model with a large set of observed SNeIa, and derive a figure of merit for the explosion model's agreement with observations. The synthetic light curves are fit with the data-driven model SALT2 which returns values for stretch, color, and magnitude at peak brightness, as well as a goodness-of-fit parameter. Each fit is performed multiple times with different choices of filter bands and epoch range in order to quantify the systematic uncertainty on the fitted parameters. We use a parametric population model for the distribution of observed SNIa parameters from large surveys, and extend it to represent red, dim, and bright outliers found in a low-redshift SNIa data set. We discuss the potential uncertainties of this population model and find it to be reliable given the current uncertainties on cosmological parameters. Using our population model, we assign each set of fitted parameters a likelihood of being observed in nature, and a figure of merit based on this likelihood. We define a second figure of merit based on the quality of the light curve fit, and combine the two measures into an overall figure of merit for each explosion model. We compute figures of merit for a variety of 1D, 2D and 3D explosion models and show that our evaluation method allows meaningful inferences across a wide range of light curve quality and fitted parameters.

The removal of magnetic flux from the quiet-sun photosphere is important for maintaining the statistical steady-state of the magnetic field there, for determining the magnetic flux budget of the Sun, and for estimating the rate of energy injected into the upper solar atmosphere. Magnetic feature death is a measurable proxy for the removal of detectable flux. We used the SWAMIS feature tracking code to understand how nearly 20000 detected magnetic features die in an hour-long sequence of Hinode/SOT/NFI magnetograms of a region of quiet Sun. Of the feature deaths that remove visible magnetic flux from the photosphere, the vast majority do so by a process that merely disperses the previously-detected flux so that it is too small and too weak to be detected. The behavior of the ensemble average of these dispersals is not consistent with a model of simple planar diffusion, suggesting that the dispersal is constrained by the evolving photospheric velocity field. We introduce the concept of the partial lifetime of magnetic features, and show that the partial lifetime due to Cancellation of magnetic flux, 22 h, is 3 times slower than previous measurements of the flux turnover time. This indicates that prior feature-based estimates of the flux replacement time may be too short, in contrast with the tendency for this quantity to decrease as resolution and instrumentation have improved. This suggests that dispersal of flux to smaller scales is more important for the replacement of magnetic fields in the quiet Sun than observed bipolar cancellation. We conclude that processes on spatial scales smaller than those visible to Hinode dominate the processes of flux emergence and cancellation, and therefore also the quantity of magnetic flux that threads the photosphere.

A common model of the explosion mechanism of Type Ia supernovae is based on a delayed detonation of a white dwarf. A variety of models differ primarily in the method by which the deflagration leads to a detonation. A common feature of the models, however, is that all of them involve the propagation of the detonation through a white dwarf that is either expanding or contracting, where the stellar internal velocity profile depends on both time and space. In this work, we investigate the effects of the pre-detonation stellar internal velocity profile and the post-detonation velocity of expansion on the production of alpha-particle nuclei, including Ni56, which are the primary nuclei produced by the detonation wave. We perform one-dimensional hydrodynamic simulations of the explosion phase of the white dwarf for center and off-center detonations with five different stellar velocity profiles at the onset of the detonation. We observe two distinct post-detonation expansion phases: rarefaction and bulk expansion. Almost all the burning to Ni56 occurs only in the rarefaction phase, and its expansion time scale is influenced by pre-existing flow structure in the star, in particular by the pre-detonation stellar velocity profile. We find that the mass fractions of the alpha-particle nuclei, including Ni56, are tight functions of the empirical physical parameter rho_up/v_down, where rho_up is the mass density immediately upstream of the detonation wave front and v_down is the velocity of the flow immediately downstream of the detonation wave front. We also find that v_down depends on the pre-detonation flow velocity. We conclude that the properties of the pre-existing flow, in particular the internal stellar velocity profile, influence the final isotopic composition of burned matter produced by the detonation.

Supersonic plasma outflows driven by multi-beam, high-energy lasers, such as Omega and NIF, have been and will be used as platforms for a variety of laboratory astrophysics experiments. Here we propose a new way of launching high density and high velocity, plasma jets using multiple intense laser beams in a hollow ring formation. We show that such jets provide a more flexible and versatile platform for future laboratory astrophysics experiments. Using high resolution hydrodynamic simulations, we demonstrate that the collimated jets can achieve much higher density, temperature and velocity when multiple laser beams are focused to form a hollow ring pattern at the target, instead of focused onto a single spot. We carried out simulations with different ring radii and studied their effects on the jet properties. Implications for laboratory collisionless shock experiments are discussed.

We study the magnetic properties of small-scale transients in coronal hole. We found all brightening events are associated with bipolar regions and caused by magnetic flux emergence followed by cancellation with the pre-existing and newly emerging magnetic flux. In the coronal hole, 19 of 22 events have a single stable polarity which does not change its position in time. In eleven cases this is the dominant polarity. The dominant flux of the coronal hole form the largest concentration of magnetic flux in terms of size while the opposite polarity is distributed in small concentrations. In the coronal hole the number of magnetic elements of the dominant polarity is four times higher than the non-dominant one. The supergranulation configuration appears to preserve its general shape during approximately nine hours of observations although the large concentrations in the network did evolve and were slightly displaced, and their strength either increased or decreased. The emission fluctuations seen in the X-ray bright points are associated with reoccurring magnetic cancellation in the footpoints. Unique observations of an X-ray jet reveal similar magnetic behaviour in the footpoints, i.e. cancellation of the opposite polarity magnetic flux. We found that the magnetic flux cancellation rate during the jet is much higher than in bright points. Not all magnetic cancellations result in an X-ray enhancement, suggesting that there is a threshold of the amount of magnetic flux involved in a cancellation above which brightening would occur at X-ray temperatures. Our study demonstrates that the magnetic flux in coronal holes is continuously recycled through magnetic reconnection which is responsible for the formation of numerous small-scale transient events. The open magnetic flux forming the coronal-hole phenomenon is largely involved in these transient features.

We describe the detonation mechanism comprising the "Pulsationally Assisted" Gravitationally Confined Detonation (GCD) model of Type Ia supernovae SNe Ia. This model is analogous to the previous GCD model reported in Jordan et al.(2008); however, the chosen initial conditions produce a substantively different detonation mechanism, resulting from a larger energy release during the deflagration phase. The resulting final kinetic energy and nickel-56 yields conform better to observational values than is the case for the "classical" GCD models. In the present class of models, the ignition of a deflagration phase leads to a rising, burning plume of ash. The ash breaks out of the surface of the white dwarf, flows laterally around the star, and converges on the collision region at the antipodal point from where it broke out. The amount of energy released during the deflagration phase is enough to cause the star to rapidly expand, so that when the ash reaches the antipodal point, the surface density is too low to initiate a detonation. Instead, as the ash flows into the collision region (while mixing with surface fuel), the star reaches its maximally expanded state and then contracts. The stellar contraction acts to increase the density of the star, including the density in the collision region. This both raises the temperature and density of the fuel-ash mixture in the collision region and ultimately leads to thermodynamic conditions that are necessary for the Zel'dovich gradient mechanism to produce a detonation. We demonstrate feasibility of this scenario with three 3-dimensional (3D), full star simulations of this model using the FLASH code. We characterized the simulations by the energy released during the deflagration phase, which ranged from 38% to 78% of the white dwarf's binding energy. We show that the necessary conditions for detonation are achieved in all three of the models.

The discovery of the accelerated expansion of the universe using Type Ia supernovae (SNe Ia) has stimulated a tremendous amount of interest in the use of SNe Type Ia events as standard cosmological candles, and as a probe of the fundamental physics of dark energy. Recent observations of SNe Ia have indicated a significant population difference depending on the host galaxy. These observational findings are consistent with SNe Ia Ni-56 production in star-forming spiral galaxies some 0.1 solar masses higher - and therefore more luminous than in elliptical galaxies. We present recent full-star, 3D simulations of Type Ia supernovae which may help explain the nature of this systematic variation in SNe Ia luminosities, as well as the nature of the Ia explosion mechanism. These insights may in turn eventually shed light on the mystery of dark energy itself.

Using a pulse-fit method, we investigate the spectral lags between the traditional gamma-ray band (50-400 keV) and the X-ray band (6-25 keV) for 8 GRBs with known redshifts (GRB 010921, GRB 020124, GRB 020127, GRB 021211, GRB 030528, GRB 040924, GRB 041006, GRB 050408) detected with the WXM and FREGATE instruments aboard the HETE-2 satellite. We find several relations for the individual GRB pulses between the spectral lag and other observables, such as the luminosity, pulse duration, and peak energy (Epeak). The obtained results are consistent with those for BATSE, indicating that the BATSE correlations are still valid at lower energies (6-25 keV). Furthermore, we find that the photon energy dependence for the spectral lags can reconcile the simple curvature effect model. We discuss the implication of these results from various points of view.

Several lines of evidence have suggested that the the galaxy cluster Cl~0024+17, an apparently relaxed system, is actually a collision of two clusters, the interaction occurring along our line of sight. In this paper we present a high-resolution $N$-body/hydrodynamics simulation of such a collision. We have created mock X-ray observations of our simulated system using MARX, a program that simulates the on-orbit performance of the Chandra X-ray Observatory. We analyze these simulated data to generate radial profiles of the surface brightness and temperature. At later times, $t = 2.0-3.0$ Gyr after the collision, the simulated surface brightness profiles are better fit by a superposition of two $\beta$-model profiles than a single profile, in agreement with the observations of Cl~0024+17. In general, due to projection effects, much of the post-collision density and temperature structure of the clusters is not seen in the observations. In particular, the observed temperatures from spectral fitting are much lower than the temperature of the hottest gas. We determine from our fitted profiles that if the system is modeled as a single cluster, the hydrostatic mass estimate is a factor $\sim$2-3 less than the actual mass, but if the system is modeled as two galaxy clusters in superposition, a hydrostatic mass estimation can be made which is accurate to within $\sim$10%. We examine some implications of these results for galaxy cluster X-ray surveys.

Several lines of evidence suggest that the galaxy cluster Cl0024+17, an apparently relaxed system, is actually a collision of two clusters, the interaction occurring along our line of sight. Recent lensing observations suggest the presence of a ring-like dark matter structure, which has been interpreted as the result of such a collision. In this paper we present $N$-body simulations of cluster collisions along the line of sight to investigate the detectability of such features. We use realistic dark matter density profiles as determined from cosmological simulations. Our simulations show a "shoulder" in the dark matter distribution after the collision, but no ring feature even when the initial particle velocity distribution is highly tangentially anisotropic ($\sigma_\theta/\sigma_r >> 1$). Only when the initial particle velocity distribution is circular do our simulations show such a feature. Even modestly anisotropic velocity distributions are inconsistent with the halo velocity distributions seen in cosmological simulations, and would require highly fine-tuned initial conditions. Our investigation leaves us without an explanation for the dark matter ring-like feature in Cl 0024+17 suggested by lensing observations.

We study the initiation of the detonation in the gravitationally confined detonation (GCD) model of Type Ia supernovae (SNe Ia). Initiation of the detonation occurs spontaneously in a region where the length scale of the temperature gradient extending from a flow (in which carbon burning is already occurring) into unburned fuel is commensurate to the range of critical length scales which have been derived from 1D simulations that resolve the initiation of a detonation. By increasing the maximum resolution in a truncated cone that encompasses this region, beginning somewhat before initiation of the detonation occurs, we successfully simulate in situ the first gradient-initiated detonation in a whole-star simulation. The detonation emerges when a compression wave overruns a pocket of fuel situated in a Kelvin-Helmholtz cusp at the leading edge of the inwardly directed jet of burning carbon. The compression wave pre-conditions the temperature in the fuel in such a way that the Zel'dovich gradient mechanism can operate and a detonation ensues. We explore the dependence of the length scale of the temperature gradient on spatial resolution and discuss the implications for the robustness of this detonation mechanism. We find that the time and the location at which initiation of the detonation occurs varies with resolution. In particular, initiation of a detonation had not yet occurred in our highest resolution simulation by the time we ended the simulation because of the computational demand it required. We suggest that the turbulent shear layer surrounding the inwardly directed jet provides the most favorable physical conditions, and therefore the most likely location, for initiation of a detonation in the GCD model.

Some explosion models for Type Ia supernovae (SN Ia), such as the gravitationally confined detonation (GCD) or the double detonation sub-Chandrasekhar (DDSC) models, rely on the spontaneous initiation of a detonation in the degenerate C/O material of a white dwarf. The length scales pertinent to the initiation of the detonation are notoriously unresolved in multi-dimensional stellar simulations, prompting the use of results of 1D simulations at higher resolution, such as the ones performed for this work, as guidelines for deciding whether or not conditions reached in the higher dimensional full star simulations successfully would lead to the onset of a detonation. Spontaneous initiation relies on the existence of a suitable gradient in self-ignition (induction) times of the fuel, which we set up with a spatially localized non-uniformity of temperature -- a hot spot. We determine the critical (smallest) sizes of such hot spots that still marginally result in a detonation in white dwarf matter by integrating the reactive Euler equations with the hydrodynamics code FLASH. We quantify the dependences of the critical sizes of such hot spots on composition, background temperature, peak temperature, geometry, and functional form of the temperature disturbance, many of which were hitherto largely unexplored in the literature. We discuss the implications of our results in the context of modeling of SNe Ia.

Aims: Taking advantage of the forthcoming Catalog of the HETE-2 mission, the aim of this paper is to evaluate the main properties of HETE-2 GRBs - the E_peak, the T_90 and the E_iso - in their source frames and to derive their unbiased distribution. Methods: We first construct a complete sample containing all the bursts localized by the WXM on-board HETE-2, which are selected with a uniform criterion and whose observed parameters can be constrained. We then derive the intrinsic E_peak, T_90 and E_iso distributions using their redshift when it is available, or their pseudo-redshift otherwise. We finally compute the number of GRB (N_Vmax) within the visibility volume (V_max) of each GRB, in order to derive a weight for each detected burst accounting both for the detection significance and the star formation history of the universe. Results: The unbiased distributions obtained clearly show the predominence of X-ray flashes (XRFs) in the global GRB population. We also derive the rate of local GRBs: R0^H2 > 11 Gpc-3 yr-1, which is intermediate between the local rate obtained by considering only the high-luminosity bursts (~1 Gpc-3 yr-1) and that obtained by including the low-luminosity bursts (>200 Gpc-3 yr-1).

We report a correlation based on a spectral simulation study of the prompt emission spectra of gamma-ray bursts (GRBs) detected by the Swift Burst Alert Telescope (BAT). The correlation is between the Epeak energy, which is the peak energy in the \nu F_\nu spectrum, and the photon index (\Gamma) derived from a simple power-law model. The Epeak - \Gamma relation, assuming the typical smoothly broken power-law spectrum of GRBs, is \log Epeak = 3.258 - 0.829\Gamma (1.3 < \Gamma < 2.3). We take into account not only a range of Epeak energies and fluences, but also distributions for both the low-energy photon index and the high-energy photon index in the smoothly broken power-law model. The distribution of burst durations in the BAT GRB sample is also included in the simulation. Our correlation is consistent with the index observed by BAT and Epeak measured by the BAT, and by other GRB instruments. Since about 85% of GRBs observed by the BAT are acceptably fit with the simple power-law model because of the relatively narrow energy range of the BAT, this relationship can be used to estimate Epeak when it is located within the BAT energy range.

We study the gravitationally confined detonation (GCD) model of Type Ia supernovae through the detonation phase and into homologous expansion. In the GCD model, a detonation is triggered by the surface flow due to single point, off-center flame ignition in carbon-oxygen white dwarfs. The simulations are unique in terms of the degree to which non-idealized physics is used to treat the reactive flow, including weak reaction rates and a time dependent treatment of material in nuclear statistical equilibrium (NSE). Careful attention is paid to accurately calculating the final composition of material which is burned to NSE and frozen out in the rapid expansion following the passage of a detonation wave over the high density core of the white dwarf; and an efficient method for nucleosynthesis post-processing is developed which obviates the need for costly network calculations along tracer particle thermodynamic trajectories. Observational diagnostics are presented for the explosion models, including abundance stratifications and integrated yields. We find that for all of the ignition conditions studied here, a self regulating process comprised of neutronization and stellar expansion results in final \iso{Ni}{56} masses of $\sim$1.1\msun. But, more energetic models result in larger total NSE and stable Fe peak yields. The total yield of intermediate mass elements is $\sim0.1$\msun and the explosion energies are all around 1.5$\times10^{51}$ ergs. The explosion models are briefly compared to the inferred properties of recent Type Ia supernova observations. The potential for surface detonation models to produce lower luminosity (lower \iso{Ni}{56} mass) supernovae is discussed.

GRB 041006 was detected by HETE-2 at 12:18:08 UT on 06 October 2004. This GRB displays a soft X-ray emission, a precursor before the onset of the main event, and also a soft X-ray tail after the end of the main peak. The light curves in four different energy bands display different features; At higher energy bands several peaks are seen in the light curve, while at lower energy bands a single broader bump dominates. It is expected that these different features are the result of a mixture of several components each of which has different energetics and variability. To reveal the nature of each component, we analysed the time resolved spectra and they are successfully resolved into several components. We also found that these components can be classified into two distinct classes; One is a component which has an exponential decay of $E_{p}$ with a characteristic timescale shorter than $\sim$ 30 sec, and its spectrum is well represented by a broken power law function, which is frequently observed in many prompt GRB emissions, so it should have an internal-shock origin. Another is a component whose $E_{p}$ is almost unchanged with characteristic timescale longer than $\sim$ 60 sec, and shows a very soft emission and slower variability. The spectrum of the soft component is characterized by either a broken power law or a black body spectrum. This component might originate from a relatively wider and lower velocity jet or a photosphere of the fireball. By assuming that the soft component is a thermal emission, the radiation radius is initially $4.4 \times 10^{6}$ km, which is a typical radius of a blue supergiant, and its expansion velocity is $2.4 \times 10^{5}$ km/s in the source frame.

We report the results of a series of three-dimensional (3-D) simulations of the deflagration phase of the gravitationally confined detonation mechanism for Type Ia supernovae. In this mechanism, ignition occurs at one or several off-center points, resulting in a burning bubble of hot ash that rises rapidly, breaks through the surface of the star, and collides at a point opposite breakout on the stellar surface. We find that detonation conditions are robustly reached in our 3-D simulations for a range of initial conditions and resolutions. Detonation conditions are achieved as the result of an inwardly-directed jet that is produced by the compression of unburnt surface material when the surface flow collides with itself. A high-velocity outwardly-directed jet is also produced. The initial conditions explored in this paper lead to conditions at detonation that can be expected to produce large amounts of $^{56}$Ni and small amounts of intermediate mass elements. These particular simulations are therefore relevant only to high luminosity Type Ia supernovae. Recent observations of Type Ia supernovae imply a compositional structure that is qualitatively consistent with that expected from these simulations.

We describe and discuss the spectral and temporal characteristics of the prompt emission and X-ray afterglow emission of X-ray flashes (XRFs) and X-ray-rich gamma-ray bursts (XRRs) detected and observed by Swift between December 2004 and September 2006. We compare these characteristics to a sample of conventional classical gamma-ray bursts (C-GRBs) observed during the same period. We confirm the correlation between Epeak_obs and fluence noted by others and find further evidence that XRFs, XRRs and C-GRBs form a continuum. We also confirm that our known redshift sample is consistent with the correlation between the peak energy in the GRB rest frame (Epeak_src) and the isotropic radiated energy (Eiso), so called the Epeak_src-Eiso relation. The spectral properties of X-ray afterglows of XRFs and C-GRBs are similar, but the temporal properties of XRFs and C-GRBs are quite different. We found that the light curves of C-GRB afterglows show a break to steeper indices (shallow-to-steep break) at much earlier times than do XRF afterglows. Moreover, the overall luminosity of XRF X-ray afterglows is systematically smaller by a factor of two or more compared to that of C-GRBs. These distinct differences between the X-ray afterglows of XRFs and C-GRBs may be the key to understanding not only the mysterious shallow-to-steep break in X-ray afterglow light curves, but also the unique nature of XRFs.

We develop an improved method for tracking the nuclear flame during the deflagration phase of a Type Ia supernova, and apply it to study the variation in outcomes expected from the gravitationally confined detonation (GCD) paradigm. A simplified 3-stage burning model and a non-static ash state are integrated with an artificially thickened advection-diffusion-reaction (ADR) flame front in order to provide an accurate but highly efficient representation of the energy release and electron capture in and after the unresolvable flame. We demonstrate that both our ADR and energy release methods do not generate significant acoustic noise, as has been a problem with previous ADR-based schemes. We proceed to model aspects of the deflagration, particularly the role of buoyancy of the hot ash, and find that our methods are reasonably well-behaved with respect to numerical resolution. We show that if a detonation occurs in material swept up by the material ejected by the first rising bubble but gravitationally confined to the white dwarf (WD) surface (the GCD paradigm), the density structure of the WD at detonation is systematically correlated with the distance of the deflagration ignition point from the center of the star. Coupled to a suitably stochastic ignition process, this correlation may provide a plausible explanation for the variety of nickel masses seen in Type Ia Supernovae.

Feature tracking and recognition are increasingly common tools for data analysis, but are typically implemented on an ad-hoc basis by individual research groups, limiting the usefulness of derived results when selection effects and algorithmic differences are not controlled. Specific results that are affected include the solar magnetic turnover time, the distributions of sizes, strengths, and lifetimes of magnetic features, and the physics of both small scale flux emergence and the small-scale dynamo. In this paper, we present the results of a detailed comparison between four tracking codes applied to a single set of data from SOHO/MDI, describe the interplay between desired tracking behavior and parameterization of tracking algorithms, and make recommendations for feature selection and tracking practice in future work.

A long X-ray flash was detected and localized by the instruments aboard the High Energy Transient Explorer II (HETE-2) at 00:03:30 UT on 2004 September 16. The position was reported to the GRB Coordinates Network (GCN) approximately 2 hours after the burst. This burst consists of two peaks separated by 200 s, with durations of 110 s and 60 s. We have analyzed the energy spectra of the 1st and 2nd peaks observed with the Wide Field X-Ray Monitor (WXM) and the French Gamma Telescope (FREGATE). We discuss the origin of the 2nd peak in terms of flux variabilities and timescales. We find that it is most likely part of the prompt emission, and is explained by the long-acting engine model. This feature is similar to some bright X-ray flares detected in the early afterglow phase of bursts observed by the Swift satellite.

A class of low-mass X-ray binary sources are known to be both X-ray burst sources and millisecond pulsars at the same time. A new source of this class was discovered by High Energy Transient Explorer 2 (HETE-2) on 14 June 2005 as a source of type-I X-ray bursts, which was named HETE J1900.1-2455. Five X-ray bursts from HETE J1900.1-2455 were observed during the summer of 2005. The time resolved spectral analysis of these bursts have revealed that their spectra are consistent with the blackbody radiation throughout the bursts. The bursts show the indication of radius expansion. The bolometric flux remains almost constant during the photospheric radius expansion while blackbody temperature dropped during the same period. Assuming that the flux reached to the Eddington limit on a standard 1.4 solar mass neutron star with a helium atmosphere, we estimate the distance to the source to be $\sim$ 4 kpc.

We analyze Swift gamma-ray bursts (GRBs) and X-ray afterglows for three GRBs with spectroscopic redshift determinations -- GRB 050401, XRF 050416a, and GRB 050525a. We find that the relation between spectral peak energy and isotropic energy of prompt emissions (the Amati relation) is consistent with that for the bursts observed in pre-Swift era. However, we find that the X-ray afterglow lightcurves, which extend up to 10-70 days, show no sign of the jet break that is expected in the standard framework of collimated outflows. We do so by showing that none of the X-ray afterglow lightcurves in our sample satisfies the relation between the spectral and temporal indices that is predicted for the phase after jet break. The jet break time can be predicted by inverting the tight empirical relation between the peak energy of the spectrum and the collimation-corrected energy of the prompt emission (the Ghirlanda relation). We find that there are no temporal breaks within the predicted time intervals in X-ray band. This requires either that the Ghirlanda relation has a larger scatter than previously thought, that the temporal break in X-rays is masked by some additional source of X-ray emission, or that it does not happen because of some unknown reason.

We develop and calibrate a realistic model flame for hydrodynamical simulations of deflagrations in white dwarf (Type Ia) supernovae. Our flame model builds on the advection-diffusion-reaction model of Khokhlov and includes electron screening and Coulomb corrections to the equation of state in a self-consistent way. We calibrate this model flame--its energetics and timescales for energy release and neutronization--with self-heating reaction network calculations that include both these Coulomb effects and up-to-date weak interactions. The burned material evolves post-flame due to both weak interactions and hydrodynamic changes in density and temperature. We develop a scheme to follow the evolution, including neutronization, of the NSE state subsequent to the passage of the flame front. As a result, our model flame is suitable for deflagration simulations over a wide range of initial central densities and can track the temperature and electron fraction of the burned material through the explosion and into the expansion of the ejecta.

We measure the large-scale real-space power spectrum P(k) using luminous red galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS) and use this measurement to sharpen constraints on cosmological parameters from the Wilkinson Microwave Anisotropy Probe (WMAP). We employ a matrix-based power spectrum estimation method using Pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 20 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.01h/Mpc < k < 0.2h/Mpc. Results from the LRG and main galaxy samples are consistent, with the former providing higher signal-to-noise. Our results are robust to omitting angular and radial density fluctuations and are consistent between different parts of the sky. They provide a striking confirmation of the predicted large-scale LCDM power spectrum. Combining only SDSS LRG and WMAP data places robust constraints on many cosmological parameters that complement prior analyses of multiple data sets. The LRGs provide independent cross-checks on Om and the baryon fraction in good agreement with WMAP. Within the context of flat LCDM models, our LRG measurements complement WMAP by sharpening the constraints on the matter density, the neutrino density and the tensor amplitude by about a factor of two, giving Omega_m=0.24+-0.02 (1 sigma), sum m_nu < 0.9 eV (95%) and r<0.3 (95%). Baryon oscillations are clearly detected and provide a robust measurement of the comoving distance to the median survey redshift z=0.35 independent of curvature and dark energy properties. Within the LCDM framework, our power spectrum measurement improves the evidence for spatial flatness, sharpening the curvature constraint Omega_tot=1.05+-0.05 from WMAP alone to Omega_tot=1.003+-0.010. Assuming Omega_tot=1, the equation of state parameter is constrained to w=-0.94+-0.09, indicating the potential for more ambitious future LRG measurements to provide precision tests of the nature of dark energy. All these constraints are essentially independent of scales k>0.1h/Mpc and associated nonlinear complications, yet agree well with more aggressive published analyses where nonlinear modeling is crucial.

Two new magnetic white dwarf accretion binaries with extremely low mass-transfer rates have been discovered in the course of the Sloan Digital Sky Survey. Measured magnetic fields are 42 MG and 57 MG, and one system orbits with a period of just 82 min. The new systems therefore significantly expand the range in properties exhibited by the small class. The measured accretion rates are very low, 0.6-5x10^{-13} Msun/yr, and multiple visits spanning more than a year confirm that this is not a short-lived characteristic. It is becoming increasingly clear that the low-mdot magnetic white dwarf binaries accrete by nearly complete magnetic capture of the stellar wind from the secondary star rather than by Roche lobe overflow. The accretion rates therefore provide some of the first realistic estimates of the total wind loss rates from M dwarfs. Although one or more of the eight systems known to date may be interrupted or possibly even extinct Polars, several lines of evidence suggest that most are pre-Polars whose evolution has not yet brought the secondaries into contact with their Roche surfaces. Considering the difficulties of identifying binaries over a wide range in field strength and accretion rate, it is quite possible that the space density of wind-accreting magnetic binaries exceeds that of the classical X-ray emitting, Roche-lobe overflow Polars.

In this work, we present a multi-wavelength study of XRF 040912, aimed at measuring its distance scale and the intrinsic burst properties. We performed a detailed spectral and temporal analysis of both the prompt and the afterglow emission and we estimated the distance scale of the likely host galaxy. We then used the currently available sample of XRFs with known distance to discuss the connection between XRFs and classical Gamma-ray Bursts (GRBs). We found that the prompt emission properties unambiguously identify this burst as an XRF, with an observed peak energy of E_p=17+/-13 keV and a burst fluence ratio S(2-30keV)/S(30-400keV)>1. A non-fading optical source with R~24 mag and with an apparently extended morphology is spatially consistent with the X-ray afterglow, likely the host galaxy. XRF 040912 is a very dark burst since no afterglow optical counterpart is detected down to R>25 mag (3 sigma limiting magnitude) at 13.6 hours after the burst. The host galaxy spectrum detected from 3800A to 10000A, shows a single emission line at 9552A. The lack of any other strong emission lines blue-ward of the detected one and the absence of the Ly alpha cut-off down to 3800A are consistent with the hypothesis of the [OII] line at redshift z=1.563+/-0.001. The intrinsic spectral properties rank this XRF among the soft GRBs in the E_peak-E_iso diagram. Similar results were obtained for most XRFs at known redshift. Only XRF 060218 and XRF 020903 represent a good example of instrinsic XRF(i-XRF) and are possibly associated with a different progenitor population. This scenario may calls for a new definition of XRFs.

GRB 051022 was detected at 13:07:58 on 22 October 2005 by HETE-2. The location of GRB 051022 was determined immediately by the flight localization system. This burst contains multiple pulses and has a rather long duration of about 190 seconds. The detections of candidate X-ray and radio afterglows were reported, whereas no optical afterglow was found. The optical spectroscopic observations of the host galaxy revealed the redshift z = 0.8. Using the data derived by HETE-2 observation of the prompt emission, we found the absorption N_H = 8.8 -2.9/+3.1 x 10^22 cm^-2 and the visual extinction A_V = 49 -16/+17 mag in the host galaxy. If this is the case, no detection of any optical transient would be quite reasonable. The absorption derived by the Swift XRT observations of the afterglow is fully consistent with those obtained from the early HETE-2 observation of the prompt emission. Our analysis implies an interpretation that the absorbing medium could be outside external shock at R ~ 10^16 cm, which may be a dusty molecular cloud.

We present Follow-Up Network for Gamma-Ray Bursts (FUN GRB) Collaboration observations of the optical afterglow of GRB 021211 made between 143 seconds and 102 days after the burst. Our unique data set includes the earliest filtered detections and color information for an afterglow in the pre-Swift era. We find that the afterglow is best described by (1) propagation through a wind-swept medium, (2) a cooling break that is blueward of the observed optical frequencies, and (3) a hard electron energy distribution. However, superimposed on this "standard model" behavior we find significant chromatic variations within the first few hours after the burst. We consider possible reasons for these variations, including the possibility that they are due to a dust echo. Finally, we constrain physical parameters that describe the afterglow and surrounding medium for a variety of scenarios and find that GRB 021211's afterglow is faint for a combination of 3-4 reasons: (1) a low fraction of energy in relativistic electrons, (2) a low density for the wind-swept medium, implying either a low mass-loss rate and/or a high wind velocity for the progenitor, (3) a wide opening/viewing angle for the jet, and possibly (4) moderate source frame extinction. The jet appears to be significantly far from equipartition and magnetically dominated. More extreme versions of this might explain the darkness of many afterglows in the Swift era.

Here we report the localizations and properties of four short-duration GRBs localized by the High Energy Transient Explorer 2 satellite (HETE-2): GRBs 010326B, 040802, 051211 and 060121, all of which were detected by the French Gamma Telescope (Fregate) and localized with the Wide-field X-ray Monitor (WXM) and/or Soft X-ray Camera (SXC) instruments. We discuss eight possible criteria for determining whether these GRBs are "short population bursts" (SPBs) or "long population bursts" (LPBs). These criteria are (1) duration, (2) pulse widths, (3) spectral hardness, (4) spectral lag, (5) energy Egamma radiated in gamma rays (or equivalently, the kinetic energy E_KE of the GRB jet), (6) existence of a long, soft bump following the burst, (7) location of the burst in the host galaxy, and (8) type of host galaxy. In particular, we have developed a likelihood method for determining the probability that a burst is an SPB or a LPB on the basis of its T90 duration alone. A striking feature of the resulting probability distribution is that the T90 duration at which a burst has an equal probability of being a SPB or a LPB is T90 = 5 s, not T90 = 2 s, as is often used. All four short-duration bursts discussed in detail in this paper have T90 durations in the Fregate 30-400 keV energy band of 1.90, 2.31, 4.25, and 1.97 sec, respectively, yielding probabilities P(S|T90) = 0.97, 0.91, 0.60, and 0.95 that these bursts are SPBs on the basis of their T90 durations alone. All four bursts also have spectral lags consistent with zero. These results provide strong evidence that all four GRBs are SPBs (abstract continues).

Thanks to the rapid position notice and response by HETE-2 and Swift, the X-ray afterglow emissions have been found for four recent short gamma-ray bursts (GRBs; GRB 050509b, GRB 050709, GRB 050724, and GRB 050813). The positions of three out of four short GRBs are coincident with galaxies with no current or recent star formation. This discovery tightens the case for a different origin for short and long GRBs. On the other hand, from the prompt emission point of view, a short GRB shows a harder spectrum comparing to that of the long duration GRBs according to the BATSE observations. We investigate the prompt emission properties of four short GRBs observed by Swift/BAT. We found that the hardness of all four BAT short GRBs is in between the BATSE range for short and long GRBs. We will discuss the spectral properties of short GRBs including the short GRB sample of Konus-Wind and HETE-2 to understand the hard nature of the BATSE short GRBs.

We present a theoretical study of double compact objects as potential short/hard gamma-ray burst (GRB) progenitors. An updated population synthesis code StarTrack is used to calculate properties of double neutron stars and black-hole neutron star binaries. We obtain their formation rates, estimate merger times and finally predict their most likely merger locations and afterglow properties for different types of host galaxies. Our results serve for a direct comparison with the recent HETE-II and SWIFT observations of several short bursts, for which afterglows and host galaxies were detected. We also discuss the possible constraints these observations put on the evolutionary models of double compact object formation. We emphasize that our double compact object models can successfully reproduce at the same time short GRBs within both young, star-forming galaxies (e.g., GRB 050709 and GRB 051221A), as well as within old, elliptical hosts (e.g., GRB 050724 and probably GRB 050509B)

Since the discovery of the first short-hard gamma-ray burst afterglows in 2005, the handful of observed events have been found to be embedded in nearby (z < 1), bright underlying galaxies. We present multiwavelength observations of the short-duration burst GRB 060121, which is the first observed to clearly outshine its host galaxy (by a factor >10^2). A photometric redshift for this event places the progenitor at a most probable redshift of z = 4.6, with a less probable scenario of z = 1.7. In either case, GRB 060121 could be the farthermost short-duration GRB detected to date and implies an isotropic-equivalent energy release in gamma-rays comparable to that seen in long-duration bursts. We discuss the implications of the released energy on the nature of the progenitor. These results suggest that GRB 060121 may belong to a family of energetic short-duration events, lying at z > 1 and whose optical afterglows would outshine their host galaxies, unlike the first short-duration GRBs observed in 2005. The possibility of GRB 060121 being an intermediate duration burst is also discussed.

We propose an improved version of the redshift indicator developed by Atteia (2003), which gets rid of the dependence on the burst duration and provides better estimates for high-redshift GRBs. We present the derivation and the definition of this redshift indicator, then its calibration with 17 GRBs with known redshifts detected by HETE-2 and 2 more detected by Konus-Wind. We also provide an estimation of the redshift for 59 bursts, and we finally discuss the redshift distribution of HETE-bursts and the possible other applications of this redshift indicator.

Gamma-Ray Bursts (GRBs) fall into two classes: short-hard and long-soft bursts. The latter are now known to have X-ray and optical afterglows, to occur at cosmological distances in star-forming galaxies, and to be associated with the explosion of massive stars. In contrast, the distance scale, the energy scale, and the progenitors of short bursts have remained a mystery. Here we report the discovery of a short-hard burst whose accurate localization has led to follow-up observations that have identified the X-ray afterglow and (for the first time) the optical afterglow of a short-hard burst. These, in turn, have led to identification of the host galaxy of the burst as a late-type galaxy at z=0.16 showing that at least some short-hard bursts occur at cosmological distances in the outskirts of galaxies, and are likely to be due to the merging of compact binaries.

In 2000, Lamb and Reichart predicted that gamma-ray bursts (GRBs) and their afterglows occur in sufficient numbers and at sufficient brightnesses at very high redshifts (z > 5) to eventually replace quasars as the preferred probe of element formation and reionization in the early universe and to be used to characterize the star-formation history of the early universe, perhaps back to when the first stars formed. Here we report the discovery of the afterglow of GRB 050904 and the identification of GRB 050904 as the first very high redshift GRB. We measure its redshift to be 6.39(+0.11,-0.12), which is consistent with the reported spectroscopic redshift (6.29 +/- 0.01). Furthermore, just redward of Ly-alpha the flux is suppressed by a factor of three on the first night, but returns to expected levels by the fourth night. We propose that this is due to absorption by molecular hydrogen that was excited to rovibrational states by the GRB's prompt emission, but was then overtaken by the jet. Now that very high redshift GRBs have been shown to exist, and at least in this case the afterglow was very bright, observing programs that are designed to capitalize on this science will likely drive a new era of study of the early universe, using GRBs as probes.

[Abridged] Very recently, relations between the peak energy of Gamma-Ray Burst burst spectra, the isotropic-equivalent energy of the burst, and the radiated energy of the burst have been found. In a way that is exactly analogous to the way in which the relation between the peak luminosity and the rate of decline of the light curve of Type Ia supernovae can be used to make Type Ia supernovae excellent standard candles for cosmology, so too, the relations between Epk, Eiso, and Egamma point toward a methodology for using GRBs as excellent standard candles for cosmology. In addition, GRBs occur over the broad redshift range from z=0.1 to at least z=4.5, and both they and their afterglows are easily detectable out to z > 8. Thus GRBs show great promise as cosmological ``yardsticks'' to measure the rate of expansion of the universe over time, and therefore the properties of dark energy (i.e., Omega_M, Omega_Lambda, w_0, and w_a). We describe a concept for a possible MIDEX-class mission dedicated to using GRBs to constrain the properties of dark energy that would obtain these quantities for > 800 bursts in the redshift range 0.1 >~ z <~ 10 during a 2-year mission. This burst sample would enable both Omega_M and w_0 to be determined to +/- 0.07 and +/- 0.06 (68% CL), respectively, and w_a to be significantly constrained.

The amount and nature of dark energy (DE) can be tightly constrained by measuring the spatial correlation features and evolution of a sample of ~ 100,000 galaxy clusters over the redshift range 0<z < 1.5. Such an X-ray survey will discover all collapsed structures with mass above 3.5e14 solar masss at redshifts z<2 (i.e. the full range where such objects are expected) in the high Galactic latitude sky. Above this mass threshold the tight correlations between X-ray observables and mass allow direct interpretation of the data. We describe the constraints on Dark Energy that can be inferred from such a survey, using powerful self-calibration techniques to relate the X-ray observables (luminosity and temperature) to the underlying mass.

HETE-II and BeppoSAX have produced a sample of GRBs and XRFs with known redshifts and $E_{pk}$. This sample provides four important empirical constraints on the nature of the source jets: Log $E_{iso}$ is approximately uniformly distributed over several orders of magnitude; the inferred prompt energy Log $E_{\gamma}$ is narrowly distributed; the Amati relation holds between $E_{iso}$ and $E_{pk}$; and the Ghirlanda relation holds between $E_{\gamma}$ and $E_{pk}$. We explore the implications of these constraints for GRB jet structure during the prompt emission phase. We infer the underlying angular profiles from the first two of the above constraints assuming all jets have the same profile and total energy, and show that such ``universal jet'' models cannot satisfy both constraints. We introduce a general and efficient method for calculating relativistic emission distributions and $E_{pk}$ distributions from jets with arbitrary (smooth) angular jet profiles. We also exhibit explicit analytical formulas for emission from top-hat jets (which are not smooth). We use these methods to exhibit $E_{pk}$ and $E_{iso}$ as a function of viewing angle, for several interesting families of GRB jet profiles. We use the same methods to calculate expected frequency distributions of $E_{iso}$ and $E_{\gamma}$ for the same families of models. We then proceed to explore the behavior of universal jet models under a range of profile shapes and parameters, to map the extent to which these models can conform to the above four empirical constraints.

HETE-2 has provided strong evidence that the properties of X-Ray Flashes (XRFs), X-ray-rich GRBs, and GRBs form a continuum, and therefore that these three kinds of bursts are the same phenomenon. A key feature found by HETE-2 is that the density of bursts is roughly constant per logarithmic interval in burst fluence S_E and observed spectral peak energy Ep_obs, and in isotropic-equivalent energy Eiso and rest frame spectral peak energy Epeak. In this paper, we explore a unified jet model of all three kinds of bursts, using population synthesis simulations of the bursts and detailed modeling of the instruments that detect them. We show that both a variable jet opening-angle model in which the emissivity is a constant independent of the angle relative to the jet axis and a universal jet model in which the emissivity is a power-law function of the angle relative to the jet axis can explain the observed properties of GRBs reasonably well. However, if one tries to account for the properties of all three kinds of bursts in a unified picture, the extra degree of freedom available in the variable jet opening-angle model enables it to explain the observations reasonably well while the power-law universal jet model cannot. The variable jet opening-angle model of XRFs, X-ray-rich GRBs, and GRBs implies that the energy Egamma radiated in gamma rays is ~ 100 times less than has been thought, and that most GRBs have very small jet opening angles (~ half a degree). It also implies that there are ~ 10^4 - 10^5 more bursts with very small jet opening angles for every burst that is observable. If this is the case, the rate of GRBs could be comparable to the rate of Type Ic core collapse supernovae.

We show analytically that GRB jets with an emissivity profile given by the Fisher distribution, epsilon(theta) = A * e^(B cos(theta)), have the unique property of producing equal numbers of bursts per logarithmic interval in Eiso, and therefore in most burst properties. Since this broad distribution of burst properties is a key feature found by HETE-2, a Fisher-shaped universal jet model can explain many of the observed properties of XRFs, X-ray-rich GRBs, and GRBs reasonably well, in contrast to a power-law universal model. For small viewing angles, the Fisher distribution can be approximated by a Gaussian, whose properties have been explored by Zhang et al. (2004). We also show that the Fisher universal jet model produces a broad distribution in the inferred radiated energy Egamma_inf, in contrast to the narrow distribution predicted by the uniform variable opening-angle jet model Lamb, Donaghy & Graziani (2005). Here we present Monte Carlo simulations of both a Fisher-shaped universal jet model and a Fisher-shaped variable opening-angle jet model.

We present a likelihood approach to modeling multi-dimensional GRB Epeak--fluence--redshift data that naturally incorporates instrument detection thresholds. The treatment of instrument thresholds is essential for analyzing evidence for GRB evolution. The method described here compares the data to a uniform jet model, in which the jet parameters are allowed to vary with redshift. Data from different experiments may be modeled jointly. In addition, BATSE data (for which no redshift information is available) may be incorporated by ascribing to each event a likelihood derived from the full model by integrating the probability density over the unknown redshift. The loss of redshift information is mitigated by the large number of available bursts. We discuss the implementation of the method, and validation of it using simulated data.

One third of all HETE-2--localized bursts are X-Ray Flashes (XRFs), a class of events first identified by Heise in which the fluence in the 2-30 keV energy band exceeds that in the 30-400 keV energy band. We summarize recent HETE-2 and other results on the properties of XRFs. These results show that the properties of XRFs, X-ray-rich gamma-ray bursts (GRBs), and GRBs form a continuum, and thus provide evidence that all three kinds of bursts are closely related phenomena. As the most extreme burst population, XRFs provide severe constraints on burst models and unique insights into the structure of GRB jets, the GRB rate, and the nature of Type Ib/Ic supernovae. We briefly mention a number of the physical models that have been proposed to explain XRFs. We then consider two fundamentally different classes of phenomenological jet models: universal jet models, in which it is posited that all GRBs jets are identical and that differences in the observed properties of the bursts are due entirely to differences in the viewing angle; and variable-opening angle jet models, in which it is posited that GRB jets have a distribution of jet opening angles and that differences in the observed properties of the bursts are due to differences in the emissivity and spectra of jets having different opening angles. We consider three shapes for the emissivity as a function of the viewing angle theta_v from the axis of the jet: power-law, top hat (or uniform), and Gaussian (or Fisher). We then discuss the effect of relativistic beaming on each of these models. We show that observations can distinguish between these various models.

We present the "C4 Cluster Catalog", a new sample of 748 clusters of galaxies identified in the spectroscopic sample of the Second Data Release (DR2) of the Sloan Digital Sky Survey (SDSS). The C4 cluster--finding algorithm identifies clusters as overdensities in a seven-dimensional position and color space, thus minimizing projection effects which plagued previous optical clusters selection. The present C4 catalog covers ~2600 square degrees of sky with groups containing 10 members to massive clusters having over 200 cluster members with redshifts. We provide cluster properties like sky location, mean redshift, galaxy membership, summed r--band optical luminosity (L_r), velocity dispersion, and measures of substructure. We use new mock galaxy catalogs to investigate the sensitivity to the various algorithm parameters, as well as to quantify purity and completeness. These mock catalogs indicate that the C4 catalog is ~90% complete and 95% pure above M_200 = 1x10^14 solar masses and within 0.03 <=z <= 0.12. The C4 algorithm finds 98% of X-ray identified clusters and 90% of Abell clusters within 0.03 <= z <= 0.12. We show that the L_r of a cluster is a more robust estimator of the halo mass (M_200) than the line-of-sight velocity dispersion or the richness of the cluster. L_r. The final SDSS data will provide ~2500 C4 clusters and will represent one of the largest and most homogeneous samples of local clusters.

GRB 020124 and GRB 030323 constitute half the sample of gamma-ray bursts with a measured redshift greater than 3. This paper presents the temporal and spectral properties of these two gamma-ray bursts detected and localized with HETE-2. While they have nearly identical redshifts (z=3.20 for GRB 020124, and z=3.37 for GRB 030323), these two GRBs span about an order of magnitude in fluence, thus sampling distinct regions of the GRB luminosity function. The properties of these two bursts are compared with those of the bulk of the GRB population detected by HETE-2. We also discuss the energetics of GRB 020124 and GRB 030323 and show that they are compatible with the Epeak - Eiso relation discovered by Amati et al. (2002). Finally, we compute the maximum redshifts at which these bursts could have been detected by HETE-2 and we address various issues connected with the detection and localization of high-z GRBs.

We present high resolution X-ray spectroscopy of two recent GRB afterglows observed with the Low Energy Transmission Gratings on Chandra. The afterglows to GRBs 030328 and 041006 are detected beginning 15.33 and 16.8 hours after each burst, respectively, and are observed to fade in time during each ~90 ksec observation. We fit for the continuum emission in each full data set and for the data sliced into half and quarter time sections. For both afterglows, the continuum emission is well described by an absorbed power-law model, and the model parameters describing the absorption and spectral slope do not appear to evolve in time. We perform a careful search for deviations from the model continua for the full and time-sliced data and find no evidence for significant (>~3sigma) narrow emission/absorption lines or edges. The lack of detections implies that line emission--if it is a general feature in GRB X-ray afterglows--occurs early (t~0.3 days in the source frame) and/or is short-lived (dt<~10 ksec). We also comment on synchrotron fireball models for the X-ray and optical data.

We report High Energy Transient Explorer 2 (HETE-2) Wide Field X-ray Monitor/French Gamma Telescope observations of XRF030723 along with observations of the XRF afterglow made using the 6.5m Magellan Clay telescope and the Chandra X-ray Observatory. The observed peak energy E_pk_obs of the nu F_nu burst spectrum is found to lie within (or below) the WXM 2-25 keV passband at 98.5% confidence, and no counts are detected above 30 keV. Our best fit value is E_pk_obs=8.4+3.5/-3.4 keV. The ratio of X-ray to Gamma-ray flux for the burst follows a correlation found for GRBs observed with HETE-2, and the duration of the burst is similar to that typical of long-duration GRBs. If we require that the burst isotropic equivalent energy E_iso and E_pk_rest satisfy the relation discovered by Amati et al. (2002), a redshift of z=0.38+0.36/-0.18 can be determined, in agreement with constraints determined from optical observations. We are able to fit the X-ray afterglow spectrum and to measure its temporal fade. Although the best-fit fade is shallower than the concurrent fade in the optical, the spectral similarity between the two bands indicates that the X-ray fade may actually trace the optical fade. If this is the case, the late time rebrightening observed in the optical cannot be due to a supernova bump. We interpret the prompt and afterglow X-ray emission as arising from a jetted GRB observed off-axis and possibly viewed through a complex circumburst medium due to a progenitor wind.

An exceptionally intense gamma-ray burst, GRB030329, was detected and localized by the instruments on board the High Energy Transient Explorer satellite (HETE) at 11:37:14 UT on 29 March 2003. The burst consisted of two \~10s pulses of roughly equal brightness and an X-ray tail lasting >100s. The energy fluence in the 30-400 keV energy band was 1.08e-4 erg/cm2, making GRB030329 one of the brightest GRBs ever detected. Communication of a 2 arcmin error box 73 minutes after the burst allowed the rapid detection of a counterpart in the optical, X-ray, radio and the ensuing discovery of a supernova with most unusual characteristics. Analyses of the burst lightcurves reveal the presence of a distinct, bright, soft X-ray component underlying the main GRB: the 2-10 keV fluence of this component is ~7e-6 erg/cm2. The main pulses of GRB030329 were preceded by two soft, faint, non-thermal bumps. We present details of the HETE observations of GRB030329.

We describe and discuss the global properties of 45 gamma-ray bursts (GRBs) observed by HETE-2 during the first three years of its mission, focusing on the properties of X-Ray Flashes (XRFs) and X-ray-rich GRBs (XRRs). We find that the numbers of XRFs, XRRs, and GRBs are comparable. We find that the durations and the sky distributions of XRFs and XRRs are similar to those of GRBs. We also find that the spectral properties of XRFs and XRRs are similar to those of GRBs, except that the values of the peak energy $E^{\rm obs}_{\rm peak}$ of the burst spectrum in $\nu F_\nu$, the peak energy flux $\Fp$, and the energy fluence $S_E$ of XRFs are much smaller -- and those of XRRs are smaller -- than those of GRBs. Finally, we find that the distributions of all three kinds of bursts form a continuum in the [$S_E$(2-30 keV),$S_E$(30-400) keV]-plane, the [$S_E$(2-400 keV), $E_{\rm peak}$]-plane, and the [$F_{\rm peak}$(50-300 keV), $E_{\rm peak}$]-plane. These results provide strong evidence that all three kinds of bursts arise from the same phenomenon.

The Sloan Digital Sky Survey has revealed an overdensity of luminous red giant stars ~ 3 degrees (40 projected kpc) to the northeast of M31, which we have called Andromeda NE. The line-of-sight distance to Andromeda NE is within approximately 50 kpc of M31; Andromeda NE is not a physically unrelated projection. Andromeda NE has a g-band absolute magnitude of ~ -11.6 and central surface brightness of ~ 29 mag/sq.arcsec, making it nearly two orders of magnitude more diffuse than any known Local Group dwarf galaxy at that luminosity. Based on its distance and morphology, Andromeda NE is likely undergoing tidal disruption. Andromeda NE's red giant branch color is unlike that of M31's present-day outer disk or the stellar stream reported by Ibata et al. (2001), arguing against a direct link between Andromeda NE and these structures. However, Andromeda NE has a red giant branch color similar to that of the G1 clump; it is possible that these structures are both material torn off of M31's disk in the distant past, or that these are both part of one ancient stellar stream.

We report the discovery of five very cool white dwarfs in the Sloan Digital Sky Survey (SDSS). Four are ultracool, exhibiting strong collision induced absorption (CIA) from molecular hydrogen and are similar in color to the three previously known coolest white dwarfs, SDSS J1337+00, LHS 3250 and LHS 1402. The fifth, an ultracool white dwarf candidate, shows milder CIA flux suppression and has a color and spectral shape similar to WD 0346+246. All five new white dwarfs are faint (g > 18.9) and have significant proper motions. One of the new ultracool white dwarfs, SDSS J0947, appears to be in a binary system with a slightly warmer (T_{eff} ~ 5000K) white dwarf companion.

We present a new mechanism for Type Ia supernova explosions in massive white dwarfs. The proposed scenario follows from relaxing the assumption of symmetry in the model and involves a detonation created in an unconfined environment. The explosion begins with an essentially central ignition of stellar material initiating a deflagration. This deflagration results in the formation of a buoyantly-driven bubble of hot material that reaches the stellar surface at supersonic speeds. The bubble breakout forms a strong pressure wave that laterally accelerates fuel-rich outer stellar layers. This material, confined by gravity to the white dwarf, races along the stellar surface and is focused at the location opposite to the point of the bubble breakout. These streams of nuclear fuel carry enough mass and energy to trigger a detonation just above the stellar surface. The flow conditions at that moment support a detonation that will incinerate the white dwarf and result in an energetic explosion. The stellar expansion following the deflagration redistributes stellar mass in a way that ensures production of intermediate mass and iron group elements consistent with observations. The ejecta will have a strongly layered structure with a mild amount of asymmetry following from the early deflagration phase. This asymmetry, combined with the amount of stellar expansion determined by details of the evolution (principally the energetics of deflagration, timing of detonation, and structure of the progenitor), can be expected to create a family of mildly diverse Type Ia supernova explosions.

We present the first high-resolution three-dimensional simulations of the deflagration phase of Type Ia supernovae that treat the entire massive white dwarf. We report the results of simulations in which ignition of the nuclear burning occurs slightly off-center. The subsequent evolution of the nuclear burning is surprisingly asymmetric with a growing bubble of hot ash rapidly rising to the stellar surface. Upon reaching the surface, the mass of burned material is $\approx 0.075 M_\sun$ and the kinetic energy is $4.3 \times 10^{49}$ ergs. The velocity of the top of the rising bubble approaches 8000 km s$^{-1}$. The amount of the asymmetry found in the model offers a natural explanation for the observed diversity in Type Ia supernovae. Our study strongly disfavors the classic central-ignition pure deflagration scenario by showing that the result is highly sensitive to details of the initial conditions.

We present the full spectroscopic white dwarf and hot subdwarf sample from the SDSS first data release, DR1. We find 2551 white dwarf stars of various types, 240 hot subdwarf stars, and an additional 144 objects we have identified as uncertain white dwarf stars. Of the white dwarf stars, 1888 are non-magnetic DA types and 171, non-magnetic DBs. The remaining (492) objects consist of all different types of white dwarf stars: DO, DQ, DC, DH, DZ, hybrid stars like DAB, etc., and those with non-degenerate companions. We fit the DA and DB spectra with a grid of models to determine the Teff and log(g) for each object. For all objects, we provide coordinates, proper motions, SDSS photometric magnitudes, and enough information to retrieve the spectrum/image from the SDSS public database. This catalog nearly doubles the known sample of spectroscopically-identified white dwarf stars. In the DR1 imaged area of the sky, we increase the known sample of white dwarf stars by a factor of 8.5. We also comment on several particularly interesting objects in this sample.

HETE-2 has provided strong evidence that the properties of X-Ray Flashes (XRFs) and GRBs form a continuum, and therefore that these two types of bursts are the same phenomenon. We show that both the structured jet and the uniform jet models can explain the observed properties of GRBs reasonably well. However, if one tries to account for the properties of both XRFs and GRBs in a unified picture, the uniform jet model works reasonably well while the structured jet model fails utterly. The uniform jet model of XRFs and GRBs implies that most GRBs have very small jet opening angles (~ half a degree). This suggests that magnetic fields play a crucial role in GRB jets. The model also implies that the energy radiated in gamma rays is ~100 times smaller than has been thought. Most importantly, the model implies that there are ~10^4 -10^5 more bursts with very small jet opening angles for every such burst we see. Thus the rate of GRBs could be comparable to the rate of Type Ic core collapse supernovae. Accurate, rapid localizations of many XRFs, leading to identification of their X-ray and optical afterglows and the determination of their redshifts, will be required in order to confirm or rule out these profound implications.

The X-ray-rich GRB030528 was detected by the HETE satellite and its localization was rapidly disseminated. However, early optical observations failed to detect a counterpart source. In a 2-epoch ToO observation with Chandra, we discovered a fading X-ray source likely counterpart to GRB030528. The source brightness was typical of X-ray afterglows observed at similar epochs. Other observers detected an IR source at a location consistent with the X-ray source. The X-ray spectrum is not consistent with a large absorbing column.

The X-ray-flash XRF030723 was detected by the HETE satellite and rapidly disseminated, allowing for an optical transient to be detected ~1 day after the burst. We discuss observations in the optical with Magellan, which confirmed the fade of the optical transient. In a 2-epoch ToO observation with Chandra, we discovered a fading X-ray source spatially coincident with the optical transient. We present spectral fits to the X-ray data. We also discuss the possibility that the source underwent a rebrightening in the X-rays, as was observed in the optical. We find that the significance of a possible rebrightening is very low (~1 sigma).

We review the significance determination for emission lines in the Chandra HETGS spectrum for GRB020813, and we report on a search for additional lines in high resolution Chandra spectra. No previously unclaimed features are found. We also discuss the significance of lines sets reportedly discovered using XMM data for GRB011211 and GRB030227. We find that these features are likely of modest, though not negligible, significance.

We present evidence for a ring of stars in the plane of the Milky Way, extending at least from l = 180 deg to l = 227 deg with turnoff magnitude $g \sim 19.5$; the ring could encircle the Galaxy. We infer that the low Galactic latitude structure is at a fairly constant distance of R = 18 +/- 2 kpc from the Galactic Center above the Galactic plane, and has R = 20 +/- 2 kpc in the region sampled below the Galactic plane. The evidence includes five hundred SDSS spectroscopic radial velocities of stars within 30 degrees of the plane. The velocity dispersion of the stars associated with this structure is found to be 27 km/s at (l,b) = (198, -27), 22 km/s at (l,b) = (225, 28), 30 km/s at (l,b) = (188, 24), and 30 km/s at (l,b) = (182, 27) degrees. The structure rotates in the same prograde direction as the Galactic disk stars, but with a circular velocity of 110+/-25 km/s. The narrow measured velocity dispersion is inconsistent with power law spheroid or thick disk populations. We compare the velocity dispersion in this structure with the velocity dispersion of stars in the Sagittarius dwarf galaxy tidal stream, for which we measure a velocity dispersion of 20 km/s at (l, b) = (165, -55) deg. We estimate a preliminary metallicity from the Ca II (K) line and color of the turnoff stars of [Fe/H] = -1.6 with a dispersion of 0.3 dex and note that the turnoff color is consistent with that of the spheroid population. We interpret our measurements as evidence for a tidally disrupted satellite of $2 \times 10^7$ to $5 \times 10^8 \rm M_\odot$ which rings the Galaxy.

After taking into account threshold effects, we find that the isotropic-equivalent energies E_iso and luminosities L_iso of gamma-ray bursts (GRBs) are correlated with redshift at the 5% and 0.9% signficance levels, respectively. Our results are based on 10 BeppoSAX GRBs and 11 HETE-2 GRBs with known redshifts. Our results suggest that the isotropic-equivalent energies and luminosities of GRBs increase with redshift. They strengthen earlier clues to this effect from analyses of the BATSE catalog of GRBs, using the variability of burst time histories as an estimator of burst luminosities (and therefore redshifts), and from an analysis of BeppoSAX bursts only. If the isotropic-equivalent energies and luminosities of GRBs really do increase with redshift, it suggests that GRB jets at high redshifts may be narrower and thus the cores of GRB progenitor stars at high redshifts may be rotating more rapidly. It also suggests that GRBs at very high redshifts may be more luminous -- and therefore easier to detect -- than has been thought, which would make GRBs a more powerful probe of cosmology and the early universe than has been thought.

HETE-2 has provided strong evidence that the properties of X-Ray Flashes (XRFs) and GRBs form a continuum, and therefore that these two types of bursts are the same phenomenon. We show that both the structured jet and the uniform jet models can explain the observed properties of GRBs reasonably well. However, if one tries to account for the properties of both XRFs and GRBs in a unified picture, the uniform jet model works reasonably well while the structured jet model does not. The uniform jet model of XRFs and GRBs implies that most GRBs have very small jet opening angles (~ half a degree). This suggests that magnetic fields play a crucial role in GRB jets. The model also implies that the energy radiated in gamma rays is ~100 times smaller than has been thought. Most importantly, the model implies that there are ~10^4 -10^5 more bursts with very small jet opening angles for every such burst we see. Thus the rate of GRBs could be comparable to the rate of Type Ic core collapse supernovae. Determination of the spectral parameters and redshifts of many more XRFs will be required in order to confirm or rule out the uniform jet model and its implications. HETE-2 is ideally suited to do this (it has localized 16 XRFs in ~2 years), whereas Swift is less so. The unique insights into the structure of GRBs jets, the rate of GRBs, and the nature of Type Ic supernovae that XRFs may provide therefore constitute a compelling scientific case for continuing HETE-2 during the Swift mission.

The HETE-2 mission has been highly productive. It has observed more than 250 GRBs so far. It is currently localizing 25 - 30 GRBs per year, and has localized 43 GRBs to date. Twenty-one of these localizations have led to the detection of X-ray, optical, or radio afterglows, and as of now, 11 of the bursts with afterglows have known redshifts. HETE-2 has confirmed the connection between GRBs and Type Ic supernovae, a singular achievement and certainly one of the scientific highlights of the mission so far. It has provided evidence that the isotropic-equivalent energies and luminosities of GRBs are correlated with redshift, implying that GRBs and their progenitors evolve strongly with redshift. Both of these results have profound implications for the nature of GRB progenitors and for the use of GRBs as a probe of cosmology and the early universe. HETE-2 has placed severe constraints on any X-ray or optical afterglow of a short GRB. It is also solving the mystery of "optically dark' GRBs, and revealing the nature of X-ray flashes.

The HETE-2 mission has been highly productive. It has observed more than 250 GRBs so far. It is currently localizing 25 - 30 GRBs per year, and has localized 43 GRBs to date. Twenty-one of these localizations have led to the detection of X-ray, optical, or radio afterglows, and as of now, 11 of the bursts with afterglows have redshift determinations. HETE-2 has also observed more than 45 bursts from soft gamma-ray repeaters, and more than 700 X-ray bursts. HETE-2 has confirmed the connection between GRBs and Type Ic supernovae, a singular achievement and certainly one of the scientific highlights of the mission so far. It has provided evidence that the isotropic-equivalent energies and luminosities of GRBs may be correlated with redshift; such a correlation would imply that GRBs and their progenitors evolve strongly with redshift. Both of these results have profound implications for the nature of GRB progenitors and for the use of GRBs as a probe of cosmology and the early universe. HETE-2 has placed severe constraints on any X-ray or optical afterglow of a short GRB. It has made it possible to explore the previously unknown behavior optical afterglows at very early times, and has opened up the era of high-resolution spectroscopy of GRB optical afterglows. It is also solving the mystery of "optically dark" GRBs, and revealing the nature of X-ray flashes (XRFs).

Gravitational lensing is a powerful tool for the study of the distribution of dark matter in the Universe. The cold-dark-matter model of the formation of large-scale structures predicts the existence of quasars gravitationally lensed by concentrations of dark matter so massive that the quasar images would be split by over 7 arcsec. Numerous searches for large-separation lensed quasars have, however, been unsuccessful. All of the roughly 70 lensed quasars known, including the first lensed quasar discovered, have smaller separations that can be explained in terms of galaxy-scale concentrations of baryonic matter. Although gravitationally lensed galaxies with large separations are known, quasars are more useful cosmological probes because of the simplicity of the resulting lens systems. Here we report the discovery of a lensed quasar, SDSS J1004+4112, which has a maximum separation between the components of 14.62 arcsec. Such a large separation means that the lensing object must be dominated by dark matter. Our results are fully consistent with theoretical expectations based on the cold-dark-matter model.

HETE-2 has provided new evidence that gamma-ray bursts may evolve with redshift. We investigate the consequences of this possibility for the unified jet model of XRFs and GRBs. We find that burst evolution with redshift can be naturally explained within the unified jet model, and the resulting model provides excellent agreement with existing HETE-2 and BeppoSAX data sets. In addition, this evolution model produces reasonable fits to the BATSE peak photon number flux distribution -- something that cannot be easily done without redshift evolution.

We have applied the proposed variability -> luminosity indicator to ten "X-Ray Flashes" (XRFs) observed by the Wide-Field Cameras on BeppoSAX for which BATSE survey data exists. Our results suggest that the variability -> luminosity indicator probably works for XRFs. Assuming this to be so, we find that the luminosity and redshift distributions of XRFs are consistent with those of long-duration gamma-ray bursts (GRBs), and therefore most XRFs are probably not very high redshift GRBs. The fact that XRFs and GRBs have similar luminosity and redshift distributions suggests that XRFs and long-duration GRBs are produced by a similar mechanism.

GRBs are expected to occur at redshifts far higher than the highest quasar redshifts so far detected. And unlike quasars, GRB afterglows may provide "clean" probes of the epoch of reionization because no complications (such as the presence of a strong Ly-alpha emission line) or "proximity effects" (such as the Stromgren sphere produced by ionizing photons from the quasar) are expected. Thus NIR and optical observations of GRB afterglows may provide unique information about the epoch of reionization. In particular, the flux at wavelengths shortward of Ly-alpha provides a direct measure of the density fluctuations of the IGM at the GRB redshift, while the flux at wavelengths longward of Ly-alpha provides an integrated measure of the number of ionizing photons produced by stars in the host galaxy of the GRB up until the burst occurs. A comparison of the sizes of the Stromgren spheres produced by stars in the host galaxies of GRBs and by quasars then provides an estimate of the relative contributions of star formation and quasars to reionization. We use detailed calculations of the expected shape of the GRB afterglow spectrum in the vicinity of Ly-alpha for GRBs at a variety of redshifts to illustrate these points.

The HETE-2 FREGATE and WXM instruments detected a short, hard GRB at 00:26:18.72 UT on 31 May 2002. A preliminary localization was reported as a GCN Position Notice 88 min after the burst, and a refined localization was disseminated 123 minutes later. An IPN localization of the burst was reported 18 hours after the GRB, and a refined IPN localization was disseminated ~5 days after the burst. The final IPN localization, disseminated on 25 July 2002, is a diamond-shaped region centered on RA=15h 15m 11.18s, Dec=-19o 24' 27.08" (J2000), and has an area of ~9 square arcminutes (99.7% confidence region). The prompt localization of the burst by HETE-2, coupled with the refinement of the localization by the IPN, made possible the most sensitive follow-up observations to date of a short, hard GRB at radio, optical, and X-ray wavelengths. The time history of GRB020531 at high (>30 keV) energies consists of a short, intense spike followed by a much less intense secondary peak, which is characteristic of many short, hard bursts. The duration of the burst increases with decreasing energy and the spectrum of the burst evolves from hard to soft, behaviors which are similar to those of long GRBs. This suggests that short, hard GRBs are closely related to long GRBs.

The availability of a few dozen GRB redshifts now allows studies of the intrinsic properties of these high energy transients. Amati et al. recently discovered a correlation between Epeak, the intrinsic peak energy of the $\nu f \nu$ spectrum, and Eiso, the isotropic equivalent energy radiated by the source. Lamb et al. have shown that HETE-2 data confirm and extend this correlation. We discuss here one of the consequences of this correlation: the existence of a 'spectral standard candle', which can be used to construct a simple redshift indicator for GRBs.

We report HETE-2 WXM/FREGATE observations of the X-ray flash, XRF 020903. This event was extremely soft: the ratio log(S_X/S_gamma) = 0.7, where S_X and S_gamma are the fluences in the 2-30 and 30-400 keV energy bands, is the most extreme value observed so far by HETE-2. In addition, the spectrum has an observed peak energy E^{obs}_{peak} < 5.0 keV (99.7 % probability upper limit) and no photons were detected above ~10 keV. The burst is shorter at higher energies, which is similar to the behavior of long GRBs. We consider the possibility that the burst lies at very high redshift and that the low value of E^{obs}_{peak} is due to the cosmological redshift, and show that this is very unlikely. We find that the properties of XRF 020903 are consistent with the relation between the fluences S(7-30 keV) and S(30-400 keV) found by Barraud et al. for GRBs and X-ray-rich GRBs, and are consistent with the extension by a decade of the hardness-intensity correlation (Mallozzi et al. 1995) found by the same authors. Assuming that XRF 020903 lies at a redshift z = 0.25 as implied by the host galaxy of the candidate optical and radio afterglows of this burst, we find that the properties of XRF 020903 are consistent with an extension by a factor ~300 of the relation between the isotropic-equivalent energy E_iso and the peak E_peak of the nu F_nu spectrum (in the source frame of the burst) found by Amati et al. for GRBs. The results presented in this paper therefore provide evidence that XRFs, X-ray-rich GRBs, and GRBs form a continuum and are a single phenomenon. The results also impose strong constraints on models of XRFs and X-ray-rich GRBs.

The Wide-field X-ray Monitor (WXM) is one of the scientific instruments carried on the High Energy Transient Explorer 2 (HETE-2) satellite launched on 2000 October 9. HETE-2 is an international mission consisting of a small satellite dedicated to provide broad-band observations and accurate localizations of gamma-ray bursts (GRBs). A unique feature of this mission is its capability to determine and transmit GRB coordinates in almost real-time through the burst alert network. The WXM consists of three elements: four identical Xe-filled one-dimensional position-sensitive proportional counters, two sets of one-dimensional coded apertures, and the main electronics. The WXM counters are sensitive to X-rays between 2 keV and 25 keV within a field-of-view of about 1.5 sr, with a total detector area of about 350 cm$^2$. The in-flight triggering and localization capability can produce a real-time GRB location of several to 30 arcmin accuracy, with a limiting sensitivity of $10^{-7}$ erg cm$^{-2}$. In this report, the details of the mechanical structure, electronics, on-board software, ground and in-flight calibration, and in-flight performance of the WXM are discussed.

Many classical nova ejecta are enriched in CNO and Ne. Rosner et al. recently suggested that the enrichment might originate in the resonant interaction between large-scale shear flows in the accreted H/He envelope and gravity waves at the interface between the envelope and the underlying C/O white dwarf. The shear flow amplifies the waves, which eventually form cusps and break. This wave breaking injects a spray of C/O into the superincumbent H/He. Using two-dimensional simulations, we formulate a quantitative expression for the amount of C/O per unit area that can be entrained into the H/He at saturation. The fraction of the envelope that is enriched depends on the horizontal distribution of shear velocity and the density contrast between the C/O white dwarf and the H/He layer but is roughly independent of the vertical shape of the shear profile. Using this parameterization for the mixed mass, we then perform several one-dimensional Lagrangian calculations of an accreting white dwarf envelope and consider two scenarios: that the wave breaking and mixing is driven by the convective flows; and that the mixing occurs prior to the onset of convection. In the absence of enrichment prior to ignition, the base of the convective zone, as calculated from mixing-length theory with the Ledoux instability criterion, does not reach the C/O interface. As a result, there is no additional mixing, and the runaway is slow. In contrast, the formation of a mixed layer during the accretion of H/He, prior to ignition, causes a more violent runaway. The envelope can be enriched by <25% of C/O by mass (consistent with that observed in some ejecta) for shear velocities, over the surface, with Mach numbers <0.4.

Using a sample of over 25000 spectroscopically confirmed quasars from the Sloan Digital Sky Survey, we show how quasar variability in the rest frame optical/UV regime depends upon rest frame time lag, luminosity, rest wavelength, redshift, the presence of radio and X-ray emission, and the presence of broad absorption line systems. The time dependence of variability (the structure function) is well-fit by a single power law on timescales from days to years. There is an anti-correlation of variability amplitude with rest wavelength, and quasars are systematically bluer when brighter at all redshifts. There is a strong anti-correlation of variability with quasar luminosity. There is also a significant positive correlation of variability amplitude with redshift, indicating evolution of the quasar population or the variability mechanism. We parameterize all of these relationships. Quasars with RASS X-ray detections are significantly more variable (at optical/UV wavelengths) than those without, and radio loud quasars are marginally more variable than their radio weak counterparts. We find no significant difference in the variability of quasars with and without broad absorption line troughs. Models involving multiple discrete events or gravitational microlensing are unlikely by themselves to account for the data. So-called accretion disk instability models are promising, but more quantitative predictions are needed.

HETE-2 has confirmed the connection between GRBs and Type Ic supernovae. Thus we now know that the progenitors of long GRBs are massive stars. HETE-2 has also provided strong evidence that the properties of X-Ray Flashes (XRFs) and GRBs form a continuum, and therefore that these two types of bursts are the same phenomenon. We show that both the structured jet and the uniform jet models can explain the observed properties of GRBs reasonably well. However, if one tries to account for the properties of both XRFs and GRBs in a unified picture, the uniform jet model works reasonably well while the structured jet model fails utterly. The uniform jet model of XRFs and GRBs implies that most GRBs have very small jet opening angles (~ half a degree). This suggests that magnetic fields play a crucial role in GRB jets. The model also implies that the energy radiated in gamma rays is ~100 times smaller than has been thought. Most importantly, the model implies that there are ~10^4 -10^5 more bursts with very small jet opening angles for every such burst we see. Thus the rate of GRBs could be comparable to the rate of Type Ic core collapse supernovae. Accurate, rapid localizations of many XRFs, leading to identification of their X-ray and optical afterglows and the determination of their redshifts, will be required in order to confirm or rule out these profound implications. HETE-2 is ideally suited to do this (it has localized 16 XRFs in ~2 years), whereas Swift is not. The unique insights into the structure of GRB jets, the rate of GRBs, and the nature of Type Ic supernovae that XRFs may provide therefore constitute a compelling scientific case for continuing HETE-2 during the Swift mission.

The Sloan Digital Sky Survey (SDSS) first data release provides a database of 106000 unique galaxies in the main galaxy sample with measured spectra. A sample of star-forming (SF) galaxies are identified from among the 3079 of these having 1.4 GHz luminosities from FIRST, by using optical spectral diagnostics. Using 1.4 GHz luminosities as a reference star formation rate (SFR) estimator insensitive to obscuration effects, the SFRs derived from the measured SDSS Halpha, [OII] and u-band luminosities, as well as far-infrared luminosities from IRAS, are compared. It is established that straightforward corrections for obscuration and aperture effects reliably bring the SDSS emission line and photometric SFR estimates into agreement with those at 1.4 GHz, although considerable scatter (~60%) remains in the relations. It thus appears feasible to perform detailed investigations of star formation for large and varied samples of SF galaxies through the available spectroscopic and photometric measurements from the SDSS. We provide herein exact prescriptions for determining the SFR for SDSS galaxies. The expected strong correlation between [OII] and Halpha line fluxes for SF galaxies is seen, but with a median line flux ratio F_[OII]/F_Halpha=0.23, about a factor of two smaller than that found in the sample of Kennicutt (1992). This correlation, used in deriving the [OII] SFRs, is consistent with the luminosity-dependent relation found by Jansen et al. (2001). The median obscuration for the SDSS SF systems is found to be A_Halpha=1.2 mag, while for the radio detected sample the median obscuration is notably higher, 1.6 mag, and with a broader distribution.

A new overdensity of A-colored stars in distant parts of the Milky Way's stellar halo, at a dereddened SDSS magnitude of g_0 = 20.3, is presented. Identification of associated variable RR Lyrae candidates supports the claim that these are blue horizontal branch stars. The inferred distance of these stars from the Galactic center is 90 kpc, assuming the absolute magnitude of these stars is M_g_0 = 0.7 and that the Sun is 8.5 kpc from the Galactic center. The new tidal debris is within 10 kpc of same plane as other confirmed tidal debris from the disruption of the Sagittarius dwarf galaxy, and could be associated with the trailing tidal arm. Distances to the Sagittarius stream estimated from M stars are about 13% smaller than our inferred distances. The tidal debris has a width of at least 10 degrees, and is traced for more than 20 degrees across the sky. The globular cluster NGC 2419 is located within the detected tidal debris, and may also have once been associated with the Sagittarius dwarf galaxy.

A bright, x-ray-rich GRB was detected by HETE-2 at 11:18:34.03 UT on 11 Dec 2002. The WXM localization was to 14' and relayed to the GCN 22 s after the start of the burst. The ground SXC localization was within 2' of R.A. 08h 09m 00s, Dec 06d 44' 20" (J2000). GRB021211 consists of a single, FRED-like pulse with t90s of 2.3 s (85-400 keV) and 8.5 s (2-10 keV). The peak photon number and photon energy fluxes in the 2-400 keV band, are 34.0 +/- 1.8 ph cm^-2 s^-1 and 1.68 +/- 0.11 x 10^-6 erg cm^-2 s^-1, respectively. The energy fluences in the 2-30 keV and 30-400 kev energy bands are S_X = 1.36 +/- 0.05 x 10^-6 erg cm^-2 and S_gamma 2.17 +/- 0.15 x 10^-6 erg cm^-2, respectively. Thus GRB021211 is an X-ray-rich GRB (S_X/S_gamma = 0.63 > 0.32). The spectrum is well-fit by a Band function (alpha = -0.805, beta = -2.37, E_peak = 46.8 keV). The prompt localization allowed the detection of an optical afterglow for what would otherwise have been an ``optically dark'' GRB. GRB 021211 demonstrates that some fraction of burst afterglows are ``optically dark'' because their optical afterglows at times > 1 hr after the burst are very faint, and thus have often escaped detection. GRB 021211 shows that such ``optically dim'' bursts can have very bright afterglows at times < 20 min after the burst.

The first full year of operation following the commissioning year of the Sloan Digital Sky Survey has revealed a wide variety of newly discovered cataclysmic variables. We show the SDSS spectra of forty-two cataclysmic variables observed in 2002, of which thirty-five are new classifications, four are known dwarf novae (CT Hya, RZ Leo, T Leo and BZ UMa), one is a known CV identified from a previous quasar survey (Aqr1) and two are known ROSAT or FIRST discovered CVs (RX J09445+0357, FIRST J102347.6+003841). The SDSS positions, colors and spectra of all forty-two systems are presented. In addition, the results of follow-up studies of several of these objects identify the orbital periods, velocity curves and polarization that provide the system geometry and accretion properties. While most of the SDSS discovered systems are faint (>18th mag) with low accretion rates (as implied from their spectral characteristics), there are also a few bright objects which may have escaped previous surveys due to changes in the mass transfer rate.

We report on the detection of an emission line near 1.3 keV, which we associate with blue-shifted hydrogen-like sulfur (S XVI), in a 76.8 ksec Chandra HETGS spectrum of the afterglow of GRB 020813. The line is detected at 3.3 sigma significance. We also find marginal evidence for a line possibly due to hydrogen-like silicon (Si XIV) with the same blue-shift. A line from Fe is not detected, though a very low significance Ni feature may be present. A thermal model fits the data adequately, but a reflection model may provide a better fit. There is marginal evidence that the equivalent width of the S XVI line decrease as the burst fades. We infer from these results that a supernova likely occurred >~ 2 months prior to the GRB. We find no discrete or variable spectral features in the Chandra HETGS spectrum of the GRB 021004 afterglow.

Critically understanding the `standard candle'-like behavior of Type Ia supernovae requires understanding their explosion mechanism. One family of models for Type Ia Supernovae begins with a deflagration in a Carbon-Oxygen white dwarf which greatly accelerates through wrinkling and flame instabilities. While the planar speed and behavior of astrophysically-relevant flames is increasingly well understood, more complex behavior, such as the flame's response to stretch and curvature, has not been extensively explored in the astrophysical literature; this behavior can greatly enhance or suppress instabilities and local flame-wrinkling, which in turn can increase or decrease the bulk burning rate. In this paper, we explore the effects of curvature on both nuclear flames and simpler model flames to understand the effect of curvature on the flame structure and speed.

An initial assessment is made of white dwarf and hot subdwarf stars observed in the Sloan Digital Sky Survey. In a small area of sky (190 square degrees), observed much like the full survey will be, 269 white dwarfs and 56 hot subdwarfs are identified spectroscopically where only 44 white dwarfs and 5 hot subdwarfs were known previously. Most are ordinary DA (hydrogen atmosphere) and DB (helium) types. In addition, in the full survey to date, a number of WDs have been found with uncommon spectral types. Among these are blue DQ stars displaying lines of atomic carbon; red DQ stars showing molecular bands of C_2 with a wide variety of strengths; DZ stars where Ca and occasionally Mg, Na, and/or Fe lines are detected; and magnetic WDs with a wide range of magnetic field strengths in DA, DB, DQ, and (probably) DZ spectral types. Photometry alone allows identification of stars hotter than 12000 K, and the density of these stars for 15<g<20 is found to be ~2.2 deg^{-2} at Galactic latitudes 29-62 deg. Spectra are obtained for roughly half of these hot stars. The spectra show that, for 15<g<17, 40% of hot stars are WDs and the fraction of WDs rises to ~90% at g=20. The remainder are hot sdB and sdO stars.

We present a spectral analysis of 35 GRBs detected with the HETE-2 gamma-ray detectors (the FREGATE instrument) in the energy range 7-400 keV. The GRB sample analyzed is made of GRBs localized with the Wide Field X-ray Monitor onboard HETE-2 or with the GRB Interplanetary Network. We derive the spectral parameters of the time-integrated spectra, and present the distribution of the low-energy photon index, alpha, and of the peak energy, e_peak . We then discuss the existence and nature of the recently discovered X-Ray Flashes and their relationship with classical GRBs.

The Sloan Digital Sky Survey (SDSS) has scanned the entire region containing the Draco dwarf spheroidal galaxy to 23rd magnitude in g*. We present a catalog of stars found in a 453 square arcminute, elliptical region centered on the Draco dwarf spheroidal galaxy. Objects in the catalog are matched with five previously published catalogs. The catalog contains SDSS photometry for 5634 individual objects, and also the photometry from matches to any of the other catalogs. A comparison of the photometry between catalogs allows us to identify 142 candidate variable objects. One hundred and twelve of the suspected variables have colors consistent with RR Lyrae variables.

We combine Sloan Digital Sky Survey spectra of 22,000 luminous, red, bulge-dominated galaxies to get high S/N average spectra in the rest-frame optical and ultraviolet (2600A to 7000A). The average spectra of these massive, quiescent galaxies are early-type with weak emission lines and with absorption lines indicating an apparent excess of alpha elements over solar abundance ratios. We make average spectra of subsamples selected by luminosity, environment and redshift. The average spectra are remarkable in their similarity. What variations do exist in the average spectra as a function of luminosity and environment are found to form a nearly one-parameter family in spectrum space. We present a high signal-to-noise ratio spectrum of the variation. We measure the properties of the variation with a modified version of the Lick index system and compare to model spectra from stellar population syntheses. The variation may be a combination of age and chemical abundance differences, but the conservative conclusion is that the quality of the data considerably exceeds the current state of the models.

The mass function of clusters of galaxies is determined from 400 deg^2 of early commissioning imaging data of the Sloan Digital Sky Survey; ~300 clusters in the redshift range z = 0.1 - 0.2 are used. Clusters are selected using two independent selection methods: a Matched Filter and a red-sequence color magnitude technique. The two methods yield consistent results. The cluster mass function is compared with large-scale cosmological simulations. We find a best-fit cluster normalization relation of sigma_8*omega_m^0.6 = 0.33 +- 0.03 (for 0.1 ~< omega_m ~< 0.4), or equivalently sigma_8 = (0.16/omega_m)^0.6. The amplitude of this relation is significantly lower than the previous canonical value, implying that either omega_m is lower than previously expected (omega_m = 0.16 if sigma_8 = 1) or sigma_8 is lower than expected (sigma_8 = 0.7 if omega_m = 0.3). The best-fit mass function parameters are omega_m = 0.19 (+0.08,-0.07) and sigma_8 = 0.9 (+0.3,-0.2). High values of omega_m (>= 0.4) and low sigma_8 (=< 0.6) are excluded at >~ 2 sigma.

We show that the clumpy structure of star-forming regions can naturally explain the fact that 50-70% of GRB afterglows are optically ``dark.'' We also show that dust echos from the GRB and its afterglow, produced by the clumpy structure of the star-forming region in which the GRB occurs, can lead to temporal variability and peaks in the NIR, optical, and UV lightcurves of GRB afterglows. We note that the detection of GRB ``orphan'' afterglows would provide strong evidence that the star-forming regions in which GRBs occur are clumpy.

We use the possible Cepheid-like luminosity estimator for the long-duration gamma-ray bursts (GRBs) developed by Reichart et al. (2000) to estimate the intrinsic luminosity, and thus the redshift, of 907 long-duration GRBs from the BATSE 4B catalog. We describe a method based on Bayesian inference which allows us to infer the intrinsic GRB burst rate as a function of redshift for bursts with estimated intrinsic luminosities and redshifts. We apply this method to the above sample of long-duration GRBs, and present some preliminary results.

We exploit the 14 gamma-ray bursts (GRBs) with known redshifts $z$ and the 7 GRBs for which there are constraints on $z$ to determine the GRB rate $R_{\rm GRB}(z)$, using a method based on Bayesian inference. We find that, despite the qualitative differences between the observed GRB rate and estimates of the SFR in the universe, current data are consistent with $R_{\rm GRB}(z)$ being proportional to the SFR.

We show that, if the long GRBs are produced by the collapse of massive stars, GRBs and their afterglows may provide a powerful probe of cosmology and the early universe.

We describe a new method of transient point source localization for coded-aperture X-ray detectors that we have applied to data from the HETE Wide-Field X-Ray Monitor (WXM). The method is based upon the calculation of the likelihood function and its interpretation as a probability density for the transient source location by an application of Bayes' Theorem. The method gives a point estimate of the source location by finding the maximum of this probability density, and credible regions for the source location by choosing suitable contours of constant probability density. We describe the application of this method to data from the WXM, and give examples of GRB localizations which illustrate the results that can be obtained using this method.

WXM gives GRB localizations in instrument coordinates. WXM localizations must be converted to celestial coordinates using spacecraft aspect information obtained by the optical cameras on HETE. We must therefore accurately determine the alignment of the WXM boresight with respect to that of the optical cameras, in order to accurately determine the celestial coordinates of WXM burst locations. We use a seven-parameter model that treats as free parameters the three Euler angles of a pure rotation, two horizontal shifts of the coded-aperture masks with respect to the detectors, and the heights of the masks above the two detectors. We determine the alignment by fitting the model to a set of 252 WXM localizations of Sco X-1 obtained between 23 April and 28 June 2001. We estimate the systematic error in WXM GRB locations by comparing the actual and the calculated locations of Sco X-1. We find that the systematic error corresponding to a 68.3% confidence region is 1.7$'$, and the systematic error corresponding to a 90% confidence region is 2.4$'$. We find that this astrometric solution also provides a satisfactory fit to an independent sample of SGR and XRB events. These results are consistent with the astrometric calibration and the systematic error in WXM localizations derived independently using the RIKEN localization method.

(Abridged) We present in this paper a detailed analysis of the effect of environment on the star-formation activity of galaxies within the EDR of the SDSS. We have used the Halpha emission line to derive the star-formation rate (SFR) for each galaxy within a volume-limited sample of 8598 galaxies with 0.05 < z < 0.095 and M(r)<= -20.45. We find that the SFR of galaxies is strongly correlated with the local (projected) galaxy density and thus we present here the density-SFR relation that is analogous to the density-morphology relation. The effect of density on the SFR of galaxies is seen in three ways. First, the overall distribution of SFRs is shifted to lower values in dense environments compared with the field population. Second, the effect is most noticeable for the strongly star-forming galaxies in the 75th percentile of the SFR distribution. Third, there is a ``break'' (or characteristic density) in the density-SFR relation at a local galaxy density of 1h-2 Mpc-2. To understand this break further, we have studied the SFR of galaxies as a function of clustercentric radius from 17 clusters and groups objectively selected from the SDSS EDR data. The distribution of SFRs of cluster galaxies begins to change, compared with the field population, at a clustercentric radius of 3-4 virial radii, which is consistent with the characteristic break in density that we observe in the density-SFR relation. Our tests suggest that the density-morphology relation alone is unlikely to explain the density-SFR relation we observe. Taken all together, these works demonstrate that the decrease in SFR of galaxies in dense environments is a universal phenomenon over a wide range in density (from 0.08 to 10h-2 Mpc-2) and redshift (out to z = 0.5).

The effective optical depth in the Ly-alpha forest region of 1061 low-resolution QSO spectra drawn from the SDSS database decreases with decreasing redshift over the range 2.5 < z < 4. Although the evolution is relatively smooth, tau_eff ~ (1+z)^{3.8 pm 0.2}, at z ~ 3.2 the effective optical depth decreases suddenly, by about ten percent with respect to this smoother evolution. It climbs back to the original smooth scaling again by z ~ 2.9. We describe two techniques, one of which is new, for quantifying this evolution which give consistent results. A variety of tests show that the feature is not likely to be a consequence of how the QSO sample was selected, nor the result of flux calibration or other systematic effects. Other authors have argued that, at this same epoch, the temperature of the IGM also shows a departure from an otherwise smooth decrease with time. These features in the evolution of the temperature and the optical depth are signatures of the reionization of He II.

Using photometry and spectroscopy of 144,609 galaxies from the Sloan Digital Sky Survey, we present bivariate distributions of pairs of seven galaxy properties: four optical colors, surface brightness, radial profile shape as measured by the Sersic index, and absolute magnitude. In addition, we present the dependence of local galaxy density (smoothed on 8 h^{-1} Mpc scales) on all of these properties. Several classic, well-known relations among galaxy properties are evident at extremely high signal-to-noise ratio: the color-color relations of galaxies, the color-magnitude relations, the magnitude-surface brightness relation, and the dependence of density on color and absolute magnitude. We show that most of the i-band luminosity density in the universe is in the absolute magnitude and surface brightness ranges used. Some of the relationships between parameters, in particular the color--magnitude relations, show stronger correlations for exponential galaxies and concentrated galaxies taken separately than for all galaxies taken together. We provide a simple set of fits of the dependence of galaxy properties on luminosity for these two sets of galaxies.

The {\it HETE-2} (hereafter \HETE) French Gamma Telescope (FREGATE) and the Wide-field X-ray Monitor (WXM) instruments detected a short ($t_{50} = 360$ msec in the FREGATE 85-300 keV energy band), hard gamma-ray burst (GRB) that occurred at 1578.72 SOD (00:26:18.72 UT) on 31 May 2002. The WXM flight localization software produced a valid location in spacecraft (relative) coordinates. However, since no on-board real-time star camera aspect was available, an absolute localization could not be disseminated. A preliminary localization was reported as a GCN Position Notice at 01:54:22 UT, 88 min after the burst. Further ground analysis produced a refined localization, which can be expressed as a 90% confidence rectangle that is 67 arcminutes in RA and 43 arcminutes in Dec (90% confidence region), centered at RA = +15$^{\rm h}$ 14$^{\rm m}$ 45$^{\rm s}$, Dec = -19$^\circ$ 21\arcmin 35\arcsec (J2000). An IPN localization of the burst was disseminated 18 hours after the GRB (Hurley et al. 2002b). A refined IPN localization was disseminated $\approx$ 5 days after the burst. This hexagonal-shaped localization error region is centered on RA = 15$^{\rm h}$ 15$^{\rm m}$ 03.57$^{\rm s}$, -19$^\circ$ 24\arcmin 51.00\arcsec (J2000), and has an area of $\approx$ 22 square arcminutes (99.7% confidence region). The prompt localization of this short, hard GRB by \HETE and the anti-Sun pointing of the \HETE instruments, coupled with the refinement of the localization by the IPN, has made possible rapid follow-up observations of the burst at radio, optical, and X-ray wavelengths.

We report the discovery of a highly luminous optical transient (OT), SDSS J124602.54+011318.8, associated with a galaxy at a redshift of 0.385. In this paper we consider the possibility that the OT may be a GRB afterglow. Three sets of images and two sets of spectra were obtained as part of the normal operations of the Sloan Digital Sky Survey (SDSS). In the first two image sets, observed two nights apart, the object appears as a point source at $r^{*}\approx 17$. The third image set, observed about 410 days later, shows an extended source which is more than 2.5 magnitudes fainter. The spectra were observed about 400 and 670 days after the first two image sets, and both show an apparently normal galaxy at a redshift of 0.385. Associating the OT with the galaxy, the absolute magnitude was $M_{r^*}=-24.8$, which is over 4 magnitudes brighter than the most luminous supernova ever measured. The spectral energy distributions of the galaxy-subtracted OT derived from the first two image sets are well-fit by single power-laws with indices of $\beta_{\nu}=-0.92$ and -1.29 respectively, similar to most GRB afterglows. Based upon the luminosity of the OT, non-detections in contemporaneous ROTSE-I images, and the change in spectral slope, the OT, if an afterglow, was likely discovered early during a ``plateau'' or slowly-fading phase. The discovery of a GRB afterglow at this stage of the SDSS is consistent with expectations, but only if the optical emission is much less strongly beamed than the gamma-rays. We emphasize that other explanations for the OT cannot be ruled out; a recent follow-up study by [galyam02] provides strong evidence that this source is in fact an unusual AGN.

Galaxy masses can be estimated by a variety of methods; each applicable in different circumstances, and each suffering from different systematic uncertainties. Confirmation of results obtained by one technique with analysis by another is particularly important. Recent SDSS weak lensing measurements of the projected-mass correlation function reveal a linear relation between galaxy luminosities and the depth of their dark matter halos (measured on 260 \hinv kpc scales). In this work we use an entirely independent dynamical method to confirm these results. We begin by assembling a sample of 618 relatively isolated host galaxies, surrounded by a total of 1225 substantially fainter satellites. We observe the mean dynamical effect of these hosts on the motions of their satellites by assembling velocity difference histograms. Dividing the sample by host properties, we find significant variations in satellite velocity dispersion with host luminosity. We quantify these variations using a simple dynamical model, measuring \mtsd a dynamical mass within 260 \hinv kpc. The appropriateness of this mass reconstruction is checked by conducting a similar analysis within an N-body simulation. Comparison between the dynamical and lensing mass-to-light scalings shows reasonable agreement, providing some quantitative confirmation for the lensing results.

We report on exploratory Chandra observations of the three highest redshift quasars known (z = 5.82, 5.99, and 6.28), all found in the Sloan Digital Sky Survey. These data, combined with a previous XMM-Newton observation of a z = 5.74 quasar, form a complete set of color-selected, z > 5.7 quasars. X-ray emission is detected from all of the quasars at levels that indicate that the X-ray to optical flux ratios of z ~ 6 optically selected quasars are similar to those of lower redshift quasars. The observations demonstrate that it will be feasible to obtain quality X-ray spectra of z ~ 6 quasars with current and future X-ray missions.

We present new VLT spectroscopic observations of the most distant quasar known, SDSS J1030+0524 at z=6.28 which was recently discovered by the Sloan Digital Sky Survey. We confirm the presence of a complete Gunn-Peterson trough caused by neutral hydrogen in the intergalactic medium. There is no detectable flux over the wavelength range from 8450 to 8710 A. We set an improved limit on the drop of the flux level blueward of the Ly alpha line: a factor of > 200. Below 8450 A the spectrum shows a rise in flux, with a large fraction (> 60 %) of the total emission produced by few narrow features of transmitted flux. We discuss the obvious proximity effect around this quasar, with the presence of transmitted flux with many absorption features in a region of about 23h^{-1} comoving Mpc. If assuming the surrounding medium were completely neutral, the size of this region would imply a quasar lifetime of \~1.3x10^7 years. We also present near-IR spectroscopy of both SDSS J1030+0524 and of SDSS J1306+05, the second most distant quasar known at redshift 6.0. We combine measurements of the CIV line and limits on the HeII emission with the NV line measurements from the optical spectra to derive line ratios, and by implication the abundances of these early quasar environments. The results are indistinguishable from those of lower redshift quasars and indicate little or no evolution in the abundances from z ~ 6 to z ~ 2. The line ratios suggest supersolar metallicities, implying that the first stars around the quasars must have formed at least a few hundreds of Myrs prior to the observation, i.e. at redshift higher than 8.

On September 21 at 18950.56 SOD (05:15:50.56) UT the FREGATE gamma-ray instrument on the High Energy Transient Explorer (HETE) detected a bright gamma-ray burst (GRB). The burst was also seen by the X-detector on the WXM X-ray instrument and was therefore well-localized in the X direction; however, the burst was outside the fully-coded field-of-view of the WXM Y-detector, and therefore information on the Y direction of the burst was limited. Cross-correlation of the HETE and Ulysses time histories yielded an Interplanetary Network (IPN) annulus that crosses the HETE error strip at a ~45 degree angle. The intersection of the HETE error strip and the IPN annulus produces a diamond-shaped error region for the location of the burst having an area of 310 square arcminutes. Based on the FREGATE and WXM light curves, the duration of the burst is characterized by a t90 = 18.4 s in the WXM 4 - 25 keV energy range, and 23.8 s and 21.8 s in the FREGATE 6 - 40 and 32 - 400 keV energy ranges, respectively. The fluence of the burst in these same energy ranges is 4.8 10^{-6}, 5.5 10^{-6}, and 11.4 10^{-6} erg cm^{-2}, respectively. Subsequent optical and radio observations by ground-based observers have identified the afterglow of GRB010921 and determined an apparent redshift of z = 0.450.

We present multi-instrument optical observations of the High Energy Transient Explorer (HETE-2)/Interplanetary Network (IPN) error box of GRB 010921. This event was the first gamma ray burst (GRB) localized by HETE-2 which has resulted in the detection of an optical afterglow. In this paper we report the earliest known observations of the GRB010921 field, taken with the 0.11-m Livermore Optical Transient Imaging System (LOTIS) telescope, and the earliest known detection of the GRB010921 optical afterglow, using the 0.5-m Sloan Digital Sky Survey Photometric Telescope (SDSS PT). Observations with the LOTIS telescope began during a routine sky patrol 52 minutes after the burst. Observations were made with the SDSS PT, the 0.6-m Super-LOTIS telescope, and the 1.34-m Tautenburg Schmidt telescope at 21.3, 21.8, and 37.5 hours after the GRB, respectively. In addition, the host galaxy was observed with the USNOFS 1.0-m telescope 56 days after the burst. We find that at later times (t > 1 day after the burst), the optical afterglow exhibited a power-law decline with a slope of $\alpha = 1.75 \pm 0.28$. However, our earliest observations show that this power-law decline can not have extended to early times (t < 0.035 day).

We identify new structures in the halo of the Milky Way Galaxy from positions, colors and magnitudes of five million stars detected in the Sloan Digital Sky Survey. Most of these stars are within 1.26 degrees of the celestial equator. We present color-magnitude diagrams (CMDs) for stars in two previously discovered, tidally disrupted structures. The CMDs and turnoff colors are consistent with those of the Sagittarius dwarf galaxy, as had been predicted. In one direction, we are even able to detect a clump of red stars, similar to that of the Sagittarius dwarf, from stars spread across 110 square degrees of sky. Focusing on stars with the colors of F turnoff objects, we identify at least five additional overdensities of stars. Four of these may be pieces of the same halo structure, which would cover a region of the sky at least 40 degrees in diameter, at a distance of 11 kpc from the Sun (18 kpc from the center of the Galaxy). The turnoff is significantly bluer than that of thick disk stars, and closer to the Galactic plane than a power-law spheroid. We suggest two models to explain this new structure. One possibility is that this new structure could be a new dwarf satellite of the Milky Way, hidden in the Galactic plane, and in the process of being tidally disrupted. The other possibility is that it could be part of a disk-like distribution of stars which is metal-poor, with a scale height of approximately 2 kpc and a scale length of approximately 10 kpc. The fifth overdensity, which is 20 kpc away, is some distance from the Sagittarius dwarf streamer orbit and is not associated with any known structure in the Galactic plane. It is likely that there are many smaller streams of stars in the Galactic halo.

We identify close pairs of galaxies from 278 deg^2 of Sloan Digital Sky Survey commissioning imaging data. The pairs are drawn from a sample of 330,041 galaxies with 18 < r^* < 20. We determine the angular correlation function of galaxy pairs, and find it to be stronger than the correlation function of single galaxies by a factor of 2.9 +/- 0.4. The two correlation functions have the same logarithmic slope of 0.77. We invert Limber's equation to estimate the three-dimensional correlation functions; we find clustering lengths of r_0= 4.2 +/- 0.4 h^{-1} Mpc for galaxies and 7.8 +/- 0.7 h^{-1} Mpc for galaxy pairs. These results agree well with the global richness dependence of the correlation functions of galaxy systems.

A sample of nearly 9000 early-type galaxies, in the redshift range 0.01<z<0.3, was selected from the Sloan Digital Sky Survey using morphological and spectral criteria. The sample was used to study how early-type galaxy observables, including luminosity L, effective radius R, surface brightness I, color, and velocity dispersion V, are correlated with one another. Measurement biases are understood with mock catalogs which reproduce all of the observed scaling relations. At any given redshift, the intrinsic distribution of luminosities, sizes and velocity dispersions in our sample are all approximately Gaussian. In the r* band L ~ V^3.91, L ~ R^1.58, R ~ I^(-0.75), and the Fundamental Plane relation is R ~ V^(1.49) I^(-0.75). These relations are approximately the same in the g*, i* and z* bands. At fixed luminosity, the mass-to-light ratio scales as M/L ~ L^0.14. The g*-r* color scales as V^0.25. Color also correlates with magnitude and size, but these correlations are entirely due to the L-V and R-V relations. Chemical evolution and star formation histories are investigated using co-added spectra of similar objects in our sample. Chemical abundances correlate primarily with velocity dispersion. At fixed V, the higher redshift population is bluer, is weaker in Mg2, and is stronger in Hbeta than the population nearby. In addition, the population at higher redshifts is slightly more luminous. These differences are consistent with that of a passively evolving population which formed the bulk of its stars about 9 Gyrs ago. The Fundamental Plane suggests that galaxies in dense regions are slightly different from those in less dense regions, but the co-added spectra and color--magnitude relations show no statistically significant dependence on environment.

The commissioning year of the Sloan Digital Sky Survey has demonstrated that many cataclysmic variables have been missed in previous surveys with brighter limits. We report the identification of 22 cataclysmic variables, of which 19 are new discoveries and 3 are known systems (SW UMa, BH Lyn and Vir4). A compendium of positions, colors and characteristics of these systems obtained from the SDSS photometry and spectroscopy is presented along with data obtained during follow-up studies with the Apache Point Observatory (APO) and Manastash Ridge Observatory (MRO) telescopes. We have determined orbital periods for 3 of the new systems: two show dwarf nova outbursts, and the third is a likely magnetic system with eclipses of its region of line emission. Based on these results, we expect the completed survey to locate at least 400 new CVs. Most of these will be faint systems with low accretion rates that will provide new constraints on binary evolution models.

We present moderate resolution Keck spectroscopy of quasars at z=5.82, 5.99 and 6.28, discovered by the Sloan Digital Sky Survey (SDSS). We find that the Ly Alpha absorption in the spectra of these quasars evolves strongly with redshift. To z~5.7, the Ly Alpha absorption evolves as expected from an extrapolation from lower redshifts. However, in the highest redshift object, SDSSp J103027.10+052455.0 (z=6.28), the average transmitted flux is 0.0038+-0.0026 times that of the continuum level over 8450 A < lambda < 8710 A (5.95<z(abs)<6.16), consistent with zero flux. Thus the flux level drops by a factor of >150, and is consistent with zero flux in the Ly Alpha forest region immediately blueward of the Ly Alpha emission line, compared with a drop by a factor of ~10 at z(abs)~5.3. A similar break is seen at Ly Beta; because of the decreased oscillator strength of this transition, this allows us to put a considerably stronger limit, tau(eff) > 20, on the optical depth to Ly Alpha absorption at z=6. This is a clear detection of a complete Gunn-Peterson trough, caused by neutral hydrogen in the intergalactic medium. Even a small neutral hydrogen fraction in the intergalactic medium would result in an undetectable flux in the Ly Alpha forest region. Therefore, the existence of the Gunn-Peterson trough by itself does not indicate that the quasar is observed prior to the reionization epoch. However, the fast evolution of the mean absorption in these high-redshift quasars suggests that the mean ionizing background along the line of sight to this quasar has declined significantly from z~5 to 6, and the universe is approaching the reionization epoch at z~6.

We show that, if the long GRBs are produced by the collapse of massive stars, GRBs and their afterglows may provide a powerful probe of cosmology and the early universe.

Wavelets are scaleable, oscillatory functions that deviate from zero only within a limited spatial regime and have average value zero. In addition to their use as source characterizers, wavelet functions are rapidly gaining currency within the source detection field. Wavelet-based source detection involves the correlation of scaled wavelet functions with binned, two-dimensional image data. If the chosen wavelet function exhibits the property of vanishing moments, significantly non-zero correlation coefficients will be observed only where there are high-order variations in the data; e.g., they will be observed in the vicinity of sources. In this paper, we describe the mission-independent, wavelet-based source detection algorithm WAVDETECT, part of the CIAO software package. Aspects of our algorithm include: (1) the computation of local, exposure-corrected normalized (i.e. flat-fielded) background maps; (2) the correction for exposure variations within the field-of-view; (3) its applicability within the low-counts regime, as it does not require a minimum number of background counts per pixel for the accurate computation of source detection thresholds; (4) the generation of a source list in a manner that does not depend upon a detailed knowledge of the point spread function (PSF) shape; and (5) error analysis. These features make our algorithm considerably more general than previous methods developed for the analysis of X-ray image data, especially in the low count regime. We demonstrate the algorithm's robustness by applying it to various images.

We show that, if the long GRBs are produced by the collapse of massive stars, GRBs and their afterglows may provide a powerful probe of cosmology and the early universe.

We present new measurements of scaling laws relating the luminosity of galaxies to the amplitude and shape of their dark matter halos. Early imaging and spectroscopic data from the Sloan Digital Sky Survey are used to make weak lensing measurements of the surface mass density contrast Delta Sigma_+ around classes of lens objects. This surface mass density contrast as a function of radius is a measure of the galaxy-mass correlation function (GMCF). Because spectroscopic redshifts are available for all lens objects, the mass and distance scales are well constrained. The GMCF measured around ~31,000 lenses is well fit by a power law of the form Delta Sigma_+ = (2.5+0.7-0.6) (R/1 Mpc)^{-0.8+-0.2} h M_sun pc^-2. We compare this GMCF to galaxy luminosity, type, and environment, and find that it varies strongly with all three. We quantify these variations by comparing the normalization of a fit to the inner 260 h^-1 kpc, M_260, to the galaxy luminosity. While M_260 is not strongly related to luminosity in bluest band u', there is a simple, linear relation between M_260 and luminosity in redder bands (g', r', i', and z'). We test the universality of these mass-to-light scalings by independently measuring them for spiral and elliptical galaxies,and for galaxies in a variety of environments. We find remarkable consistency in these determinations in the red bands, especially i' and z'. This consistency across a wide range of systems suggests that the measured scaling represents an excellent cosmic average, and that the integrated star formation history of galaxies is strongly related to the dark matter environments in which they form.

We exploit the 14 gamma-ray bursts (GRBs) with known redshifts z and the 7 GRBs for which there are constraints on z to determine the GRB rate R_{GRB}(z), using a method based on Bayesian inference. We find that, despite the qualitative differences between the observed GRB rate and estimates of the SFR in the universe, current data are consistent with R_{GRB}(z) being proportional to the SFR.

We show that the clumpy structure of star-forming regions can naturally explain the fact that 50-70% of GRB afterglows are optical``dark.'' We also show that dust echos from the GRB and its afterglow, produced by the clumpy structure of the star-forming region in which the GRB occurs, can lead to temporal variability and peaks in the NIR, optical, and UV lightcurves of GRB afterglows.

We present a Bayesian analysis of the energies and arrival times of the neutrinos from supernova SN 1987A detected by the Kamiokande II, IMB, and Baksan detectors, and find strong evidence for two components in the neutrino signal: a long time scale component from thermal Kelvin-Helmholtz cooling of the nascent neutron star, and a brief (~< 1 s), softer component similar to that expected from emission by accreting material in the delayed supernova scenario. In the context of this model, we show that the data constrain the electron antineutrino rest mass to be less than 5.7 eV with 95% probability. Our analysis takes advantage of significant advances that have occured in the years since the detections in both our understanding of the supernova mechanism and our ability to analyze sparse data. As a result there are substantial differences between our inferences and those found in earlier studies. We find that two-component models for the neutrino signal make the data >100 times more probable than single-component models. In addition, the radius and binding energy of the nascent neutron star implied by single-component models deviates significantly from the values predicted by current neutron star models, whereas those implied by models with an accretion component are in complete agreement with the predictions. As a result, two-component models are hundreds to thousands of times more probable than single-component models. The neutrino data thus provide the first direct observational evidence in favor of the delayed supernova scenario over the prompt scenario. (Abridged abstract)

The discovery of an optical counterpart to GRB010222 (detected by BeppoSAX; Piro 2001) was announced 4.4 hrs after the burst by Henden (2001a). The Sloan Digital Sky Survey's 0.5m photometric telescope (PT) and 2.5m survey telescope were used to observe the afterglow of GRB010222 starting 4.8 hours after the GRB. The 0.5m PT observed the afterglow in five, 300 sec g' band exposures over the course of half an hour, measuring a temporal decay rate in this short period of F_nu \propto t^{-1.0+/-0.5}. The 2.5m camera imaged the counterpart nearly simultaneously in five filters (u' g' r' i' z'), with r' = 18.74+/-0.02 at 12:10 UT. These multicolor observations, corrected for reddening and the afterglow's temporal decay, are well fit by the power-law F_nu \propto nu^{-0.90+/-0.03} with the exception of the u' band UV flux which is 20% below this slope. We examine possible interpretations of this spectral shape, including source extinction in a star forming region.

We present a possible Cepheid-like luminosity estimator for the long-duration gamma-ray bursts based on the variability of their light curves. We also present a preliminary application of this luminosity estimator to 907 long-duration bursts from the BATSE catalog.

We present a possible Cepheid-like luminosity estimator for the long-duration gamma-ray bursts based on the variability of their light curves.

We present a possible Cepheid-like luminosity estimator for the long gamma-ray bursts based on the variability of their light curves. To construct the luminosity estimator, we use CGRO/BATSE data for 13 bursts, Wind/KONUS data for 5 bursts, Ulysses/GRB data for 1 burst, and NEAR/XGRS data for 1 burst. Spectroscopic redshifts, peak fluxes, and high resolution light curves are available for 11 of these bursts; partial information is available for the remaining 9 bursts. We find that the isotropic-equivalent luminosities L of these bursts positively correlate with a rigorously-constructed measure V of the variability of their light curves. We fit a model to these data that accommodates both intrinsic scatter (statistical variance) and extrinsic scatter (sample variance). If one excludes GRB 980425 from the fit on the grounds that its association with SN 1998bw at a redshift of z = 0.0085 is not secure, the luminosity estimator spans approx. 2.5 orders of magnitude in L, and the slope of the correlation between L and V is positive with a probability of 1 - 1.4 x 10^-4 (3.8 sigma). Although GRB 980425 is excluded from this fit, its L and V values are consistent with the fitted model, which suggests that GRB 980425 may well be associated with SN 1998bw, and that GRB 980425 and the cosmological bursts may share a common physical origin. If one includes GRB 980425 in the fit, the luminosity estimator spans approx. 6.3 orders of magnitude in L, and the slope of the correlation is positive with a probability of 1 - 9.3 x 10^-7 (4.9 sigma). Independently of whether or not GRB 980425 should be included in the fit, its light curve is unique in that it is much less variable than the other approx. 17 light curves in our sample for which the signal-to-noise is reasonably good.

Early data taken during commissioning of the SDSS have resulted in the discovery of a very cool white dwarf. It appears to have stronger collision induced absorption from molecular hydrogen than any other known white dwarf, suggesting it has a cooler temperature than any other. While its distance is presently unknown, it has a surprisingly small proper motion, making it unlikely to be a halo star. An analysis of white dwarf cooling times suggests that this object may be a low-mass star with a helium core. The SDSS imaging and spectroscopy also recovered LHS 3250, the coolest previously known white dwarf, indicating that the SDSS will be an effective tool for identifying these extreme objects.

We discuss some of the computational challenges encountered in simulating the evolution of clusters of galaxies. Eulerian adaptive mesh refinement (AMR) techniques can successfully address these challenges but are currently being used by only a few groups. We describe our publicly available AMR code, FLASH, which uses an object-oriented framework to manage its AMR library, physics modules, and automated verification. We outline the development of the FLASH framework to include collisionless particles, permitting it to be used for cluster simulation.

The discoveries that GRBs have X-ray, optical and radio afterglows have connected the study of GRBs to the rest of astronomy, and revolutionized the field. In this review, I discuss the implications that the observation of these afterglows have for burst energies and luminosities, and for models of the bursts and their afterglows. I describe recent evidence linking the long, softer, smoother GRBs detected by BeppoSAX and core collapse supernovae. Finally, I summarize recent work showing that, if these GRBs are due to the collapse of massive stars, they may provide a powerful probe of the very high redshift universe.

The recent evidence that the light curves and the spectra of the afterglows of GRB 970228 and GRB 980326 appear to contain a SN component, in addition to a relativistic shock wave component, provide direct clues that at least the long GRBs detected by BeppoSAX are a result of the collapse of massive stars. We show that, if this is indeed the case, GRBs and their afterglows provide a powerful probe of the very high redshift (VHR) (z > 5) universe. We first establish that GRBs and their afterglows are both detectable out to VHRs. We then show that one expects GRBs to occur out to at least z ~ 10 and possibly z \~ 15 - 20, redshifts that are far larger than those expected for the most distant quasars. This implies that there are large numbers of GRBs with peak photon number fluxes below the detection thresholds of BATSE and HETE-2, and even below the detection threshold of Swift. The mere detection of VHR GRBs would give us our first information about the earliest generations of stars. We show that GRBs and their afterglows can be used as beacons to locate core collapse SNe at redshifts z >> 1, and to study the properties of these SNe. We describe the expected properties of the absorption-line systems and the Ly-alpha forest in the spectra of GRB afterglows, and discuss various strategies for determining the redshifts of VHR GRBs. We then show how the absorption-line systems and the Ly-alpha forest visible in the spectra of GRB afterglows can be used to trace the evolution of metallicity in the universe, and to probe the large-scale structure of the universe at VHRs. Finally, we show how measurement of the Ly-alpha break in the spectra of GRB afterglows can be used to constrain, or possibly measure, the epoch at which re-ionization of the universe occurred, using the Gunn-Peterson test.

There is increasingly strong evidence that gamma-ray bursts (GRBs) are associated with star-forming galaxies, and occur near or in the star-forming regions of these galaxies. These associations provide indirect evidence that at least the long GRBs detected by BeppoSAX are a result of the collapse of massive stars. The recent evidence that the light curves and the spectra of the afterglows of GRB 970228 and GRB 980326 appear to contain a supernova component, in addition to a relativistic shock wave component, provide more direct clues that this is the case. Here we establish that GRBs and their afterglows are both detectable out to very high redshifts (z > 5).

We show that, if many GRBs are indeed produced by the collapse of massive stars, GRBs and their afterglows provide a powerful probe of the very high redshift (z > 5) universe.

The discovery of high-amplitude brightness oscillations during type I X-ray bursts from six low-mass X-ray binaries has provided a powerful new tool to study the properties of matter at supranuclear densities, the effects of strong gravity, and the propagation of thermonuclear burning. There is substantial evidence that these brightness oscillations are produced by spin modulation of one or two localized hot spots confined to the stellar surface. It is therefore important to calculate the expected light curves produced by such hot spots under various physical assumptions, so that comparison with the observed light curves may most sensitively yield information about the underlying physical quantities. In this paper we make general relativistic calculations of the light curves and oscillation amplitudes produced by a rotating neutron star with one or two hot spots as a function of spot size, stellar compactness, rotational velocity at the stellar surface, spot location, orientation of the line of sight of the observer, and the angular dependence of the surface specific intensity. We find that stellar rotation and beaming of the emission tend to increase the observed oscillation amplitudes whereas greater compactness and larger spot size tend to decrease them. By applying these results to 4U 1636--536, we show that this source must have two emitting spots and place strong constraints on the neutron star's magnetic field geometry. We also show that the data on the phase lags between photons of different energies in the persistent pulsations in SAX J1808--58 can be fit well with a model in which the observed hard leads are due to Doppler beaming.

Many thermonuclear X-ray bursts exhibit brightness oscillations. The brightness oscillations are thought to be due to the combined effects of non-uniform nuclear burning and rotation of the neutron star. The waveforms of the oscillations contain information about the size and number of burning regions. They also contain substantial information about the mass and radius of the star, and hence about strong gravity and the equation of state of matter at supranuclear densities. We have written general relativistic ray-tracing codes that compute the waveforms and spectra of rotating hot spots as a function of photon energy. Using these codes, we survey the effect on the oscillation waveform and amplitude of parameters such as the compactness of the star, the spot size, the surface rotation velocity, and whether there are one or two spots. We also fit phase lag versus photon energy curves to data from the millisecond X-ray pulsar, SAX J1808--3658.

We describe a new hybrid N-body/hydrodynamical code based on the particle-mesh (PM) method and the piecewise-parabolic method (PPM) for use in solving problems related to the evolution of large-scale structure, galaxy clusters, and individual galaxies. The code, named COSMOS, possesses several new features which distinguish it from other PM-PPM codes. In particular, to solve the Poisson equation we have written a new multigrid solver which can determine the gravitational potential of isolated matter distributions and which properly takes into account the finite-volume discretization required by PPM. All components of the code are constructed to work with a nonuniform mesh, preserving second-order spatial differences. The PPM code uses vacuum boundary conditions for isolated problems, preventing inflows when appropriate. The PM code uses a second-order variable-timestep time integration scheme. Radiative cooling and cosmological expansion terms are included. COSMOS has been implemented for parallel computers using the Parallel Virtual Machine (PVM) library, and it features a modular design which simplifies the addition of new physics and the configuration of the code for different types of problems. We discuss the equations solved by COSMOS and describe the algorithms used, with emphasis on these features. We also discuss the results of tests we have performed to establish that COSMOS works and to determine its range of validity.

We better determine the spectral flux distribution of the supernova candidate associated with GRB 970228 by modeling the spectral flux distribution of the host galaxy of this burst, fitting this model to measurements of the host galaxy, and using the fitted model to better subtract out the contribution of the host galaxy to measurements of the afterglow of this burst.

We better determine the spectral flux distribution of the supernova candidate associated with GRB 970228 by modeling the spectral flux distribution of the host galaxy of this burst, fitting this model to measurements of the host galaxy, and using the fitted model to better subtract out the contribution of the host galaxy to measurements of the afterglow of this burst. Furthermore, we discuss why the non-detection of a SN1998bw-like component to the afterglow of GRB 990510 does not necessarily imply that a SN is not associated with this burst. Finally, we discuss how bursts can be used as beacons to locate SNe out to redshifts of z = 4 - 5.

I review some of the highlights of the Rome Workshop on Gamma-Ray Bursts, and discuss some of the questions these results pose about the nature and origin of gamma-ray bursts.

We use a rigorous method, based on Bayesian inference, for calculating the odds favoring the hypothesis that any particular class of astronomical transients produce gamma-ray bursts over the hypothesis that they do not. We then apply this method to a sample of 83 Type Ia supernovae and a sample of 20 Type Ib-Ic supernovae. We find overwhelming odds against the hypothesis that all Type Ia supernovae produce gamma-ray bursts, whether at low redshift ($10^{9}:1$) or high-redshift ($10^{12}:1$), and very large odds ($6000:1$) against the hypothesis that all Type Ib, Ib/c, and Ic supernovae produce observable gamma-ray bursts. We find large odds ($34:1$) against the hypothesis that a fraction of Type Ia supernovae produce observable gamma-ray bursts, and moderate odds ($6:1$) against the hypothesis that a fraction of Type Ib-Ic supernovae produce observable bursts.

We find that the local galactic extinction towards the field of gamma-ray burst GRB970228 is $A_V=1.09^{+0.10}_{-0.20}$, which implies a substantial dimming and change in the spectral slope of the intrinsic GRB970228 afterglow. We measure a color $(V_{606}-I_{814})_{ST} = -0.18^{+0.51}_{-0.61}$ for the extended source coincident with the afterglow. Taking into account our measurement of the extinction toward this field, this color implies that the extended source is most likely a galaxy undergoing star formation, in agreement with our earlier conclusion (\cite{CL98}). In a separate analysis, we find that the inferred intrinsic spectrum of the GRB 980329 afterglow is consistent with the predictions of the simplest relativistic fireball model. We also find that the intrinsic spectrum of the afterglow is extincted both by dust (source frame $A_V \ga 1$ mag), and that the shape of the extinction curve is typical of young star-forming regions like the Orion Nebula but is not typical of older star-forming or starburst regions. The $\approx$ 2 mag drop between the $R$ and the $I$ bands can be explained by the far-ultraviolet non-linear component of the extinction curve if $3 \la z \la 4$, and by the 2175 \AA\ bump if $z \approx 2$; other redshifts are not consistent with the observational data, given our general model.

We demonstrate that models of resonant cyclotron radiation transfer in a strong field (i.e. cyclotron scattering) can account for spectral lines seen at two epochs, denoted S1 and S2, in the Ginga data for GRB870303. Using a generalized version of the Monte Carlo code of Wang et al. (1988,1989b), we model line formation by injecting continuum photons into a static plane-parallel slab of electrons threaded by a strong neutron star magnetic field (~ 10^12 G) which may be oriented at an arbitrary angle relative to the slab normal. We examine two source geometries, which we denote "1-0" and "1-1," with the numbers representing the relative electron column densities above and below the continuum photon source plane. We compare azimuthally symmetric models, i.e. models in which the magnetic field is parallel to the slab normal, with models having more general magnetic field orientations. If the bursting source has a simple dipole field, these two model classes represent line formation at the magnetic pole, or elsewhere on the stellar surface. We find that the data of S1 and S2, considered individually, are consistent with both geometries, and with all magnetic field orientations, with the exception that the S1 data clearly favor line formation away from a polar cap in the 1-1 geometry, with the best-fit model placing the line-forming region at the magnetic equator. Within both geometries, fits to the combined (S1+S2) data marginally favor models which feature equatorial line formation, and in which the observer's orientation with respect to the slab changes between the two epochs. We interpret this change as being due to neutron star rotation, and we place limits on the rotation period.

The Ginga data for the gamma-ray burst GRB870303 exhibit low-energy dips in two temporally distinct spectra, denoted S1 and S2. S1, spanning 4 s, exhibits a single line candidate at ~ 20 keV, while S2, spanning 9 s, exhibits apparently harmonically spaced line candidates at ~ 20 and 40 keV. We evaluate the statistical evidence for these lines, using phenomenological continuum and line models which in their details are independent of the distance scale to gamma-ray bursts. We employ the methodologies based on both frequentist and Bayesian statistical inference that we develop in Freeman et al. (1999b). These methodologies utilize the information present in the data to select the simplest model that adequately describes the data from among a wide range of continuum and continuum-plus-line(s) models. This ensures that the chosen model does not include free parameters that the data deem unnecessary and that would act to reduce the frequentist significance and Bayesian odds of the continuum-plus-line(s) model. We calculate the significance of the continuum-plus-line(s) models using the Chi-Square Maximum Likelihood Ratio test. We describe a parametrization of the exponentiated Gaussian absorption line shape that makes the probability surface in parameter space better-behaved, allowing us to estimate analytically the Bayesian odds. The significance of the continuum-plus-line models requested by the S1 and S2 data are 3.6 x 10^-5 and 1.7 x 10^-4 respectively, with the odds favoring them being 114:1 and 7:1. We also apply our methodology to the combined (S1+S2) data. The significance of the continuum-plus-lines model requested by the combined data is 4.2 x 10^-8, with the odds favoring it being 40,300:1.

We report I-band observations of the GRB 980329 field made on March 29 with the 1.34-m Tautenberg Schmidt telescope, R-, J- and K-band observations made on April 1 with the APO 3.5-m telescope, R- and I-band observations made on April 3 with the Mayall 4-m telescope at KPNO, and J- and K-band observations made between April 6 - 8 with the Keck-I 10-m telescope. We show that these and other reported measurements are consistent with a power-law fading of the optical/near infrared source that is coincident with the variable radio source VLA J0702+3850. This firmly establishes that this source is the afterglow of GRB 980329.

We present a rigorous method, based on Bayesian inference, for calculating the odds favoring the hypothesis that any particular class of astronomical transients produce gamma-ray bursts over the hypothesis that they do not. We then apply this method to a sample of 83 Type Ia supernovae and a sample of 20 Type Ib-Ic supernovae. We find overwhelming odds against the hypothesis that all Type Ia supernovae produce gamma-ray bursts, whether at low redshift ($10^{9}:1$) or high-redshift ($10^{12}:1$), and very large odds ($6000:1$) against the hypothesis that all Type Ib, Ib/c, and Ic supernovae produce observable gamma-ray bursts. We find large odds ($34:1$) against the hypothesis that a fraction of Type Ia supernovae produce observable gamma-ray bursts, and moderate odds ($6:1$) against the hypothesis that a fraction of Type Ib-Ic supernovae produce observable bursts. We have also re-analyzed both a corrected version of the Wang & Wheeler sample of Type Ib-Ic SNe and our larger sample of 20 Type Ib-Ic SNe, using a generalization of their frequentist method. We find no significant evidence in either case of a correlation between Type Ib-Ic SNe and GRBs, consistent with the very strong evidence against such a correlation that we find from our Bayesian analysis.

We present an R-band and J-band photometry of an optical transient which is likely to be associated with the gamma-ray burst event GRB 971214. Our first measurement took place 13 hours after the gamma-ray event. The brightness decayed with a power-law exponent alpha = -1.20 +- 0.02, which is similar to those of GRB 970228 and GRB 970508 which had exponents of alpha = -1.10 +- 0.04 and alpha = -1.141 +- 0.014 respectively. The transient decayed monotonically during the first four days following the gamma-ray event in contrast with the optical transient associated with GRB 970508 which increased in brightness, peaking two days after the burst, before settling to a power-law decay.

We carefully analyze the WFPC2 and STIS images of GRB970228. We measure magnitudes for the GRB970228 point source component in the WFPC2 images of $V=26.20^{+0.14}_{-0.13}$, $I_c=23.94^{+0.10}_{-0.09}$ and $V=26.52^{+0.16}_{-0.18}$, $I_c=24.31^{+0.15}_{-0.11}$ on March 26 and April 7, respectively; and $R_c=27.09^{+0.14}_{-0.14}$ on September 4 in the STIS image. For the extended component, we measure magnitudes of $R_c=25.48^{+0.22}_{-0.20}$ in the combined WFPC2 images and $R_c=25.54^{+0.33}_{-0.22}$ in the STIS image, which are consistent with no variation. This value is fainter than previously reported (Galama et al. 98) and modifies the previously assumed magnitudes for the optical transient when it faded to a level where the extended source component contribution was not negligible, alleviating the discrepancy to a power-law temporal behavior. We also measure a color of $V_{606}-I_{814}=-0.18^{+0.51}_{-0.61}$ for the extended source component. Taking into account the extinction measured in this field (Castander & Lamb 1998), this color implies that the extended source is most likely a galaxy with ongoing star formation.

We determine the local galactic extinction towards the field of gamma-ray burst GRB970228 using a variety of methods. We develop a maximum likelihood method for measuring the extinction by comparing galaxy counts in the field of interest to those in a field of known extinction, and apply this method to the GRB970228 field. We also measure the extinction by comparing the observed stellar spectral energy distributions of stars in the GRB970228 field to the spectral energy distribution of library spectra of the same spectral type. Finally we estimate the extinction using the Balmer emission line ratios of a galaxy in the GRB970228 field, and the neutral hydrogen column density and amount of infrared dust emission toward this field. Combining the results of these methods, we find a best-fit galactic extinction in the optical of $A_V=1.19^{+0.10}_{-0.17}$, which implies a a substantial dimming and change of the spectral slope of the intrinsic GRB970228 afterglow.

BATSE, Ulysses, and TGRS and KONUS on WIND detected four gamma-ray events within 1.8 days during 1996 October 27-29, consistent with coming from the same location on the sky. We assess the evidence that these events may be due to a series of bursts from a single source by calculating the probability that such a clustering in position and in time of occurrence might happen by chance. The calculation of this probability is afflicted by the usual problem of a posteriori statistics. We introduce a clustering statistic, which is formed from the "minimum circle radius" (i.e. the radius of the smallest circle that just encloses the positions of all the events) and the minimum time lapse (i.e. the time elapsed between the first and last event). We also introduce a second clustering statistic, which is formed from the "cluster likelihood function" and the minimum time lapse. We show that the use of these statistics largely eliminates the "a posteriori" nature of the problem. The two statistics yield significances of the clustering of $3.3\times 10^{-4}$ and $3.1\times 10^{-5}$, respectively, if we interpret the four events as four bursts, whereas the clustering is not significant if we interpret the four events as only three bursts. However, in the latter case one of the bursts is the longest ever observed by BATSE.

The IRAS 100 micron image of the GRB 970228 field shows that the amount of galactic dust in this direction is substantial and varies on arcminute angular scales. From an analysis of the observed surface density of galaxies in the $2.6' \times 2.6'$ HST WFPC image of the GRB 970228 field, we find $A_V = 1.1 \pm 0.10$. From an analysis of the observed spectra of three stars in the GRB 970228 field, we find $A_V = 1.71^{+0.20}_{-0.40}$. This value may represent the best estimate of the extinction in the direction of GRB 970228, since these three stars lie only $2.7''$, $16''$, and $42''$ away from the optical transient. If instead we combine the two results, we obtain a conservative value $A_V = 1.3 \pm 0.2$. This value is significantly larger than the values $A_V = 0.4 - 0.8$ used in papers to date. The value of $A_V$ that we find implies that, if the extended source in the burst error circle is extragalactic and therefore lies beyond the dust in our own galaxy, its optical spectrum is very blue: its observed color $(V-I_c)_{obs} \approx 0.65^{+0.74}_{-0.94}$ is consistent only with a starburst galaxy, an irregular galaxy at $z > 1.5$, or a spiral galaxy at $z > 2$. On the other hand, its observed color and surface brightness $\mu_V \approx 24.5$ arcsec$^{-2}$ are similar to those expected for the reflected light from a dust cloud in our own galaxy, if the cloud lies in front of most of the dust in this direction.

We present a new method of transient point source deconvolution for coded-aperture X-Ray detectors. Our method is based upon the calculation of the likelihood function and its interpretation as a probability density for the transient source position by an application of Bayes' Theorem. The method obtains point estimates of source positions by finding the maximum of this probability density, and interval estimates of prescribed probability by choosing suitable contours of constant probability density. We give the results of simulations that we performed to test the method. We also derive approximate analytic expressions for the predicted performance of the method. These estimates underline the intuitively plausible properties of the method and provide a sound quantitative basis for the design of coded-aperture systems.

We calculate the properties of gamma-ray burst spectral lines formed by resonant cyclotron scattering in a radiation-driven outflow. Most previous models of line formation in gamma-ray bursts are appropriate at the polar cap of a neutron star located no further than several hundred parsecs away. However, the BATSE brightness and sky distributions indicate that, if bursters are galactic, they are located in a corona at distances greater than 100 kpc. At these distances the burst luminosity exceeds the Eddington luminosity and the plasma in a polar line-forming region is ejected along the field lines. The variation of the magnetic field strength and plasma velocity with altitude in such an outflow would seem to prevent the formation of narrow features like those observed by Ginga and other instruments. However, this is not the case because the majority of scatters occur close to the stellar surface, at altitudes below $z \approx r_h$, where $r_h \sim 10^{5}$ cm is the size of the photon source. Consequently, the interpretation of the observed features as cyclotron lines does not rule out burst sources in a galactic corona. The outflow model predicts both absorption-like and single-peaked emission-like features. The latter do not occur in models with static line-forming regions, but have been observed in gamma-ray bursts by the Konus and BATSE instruments.

We use a Monte Carlo radiative transfer code to examine the dependence of the properties of cyclotron resonant scattering lines on the spatial geometry and the optical depth of the line-forming region. We focus most of our attention on a line-forming region that is a plane-parallel slab. We also consider a cylindrical line-forming region. In both cases, the line-forming region contains an electron-proton plasma at the equilibrium Compton temperature, $T_c$, and is threaded with a uniform magnetic field with strength $\sim 10^{12}$ gauss. We consider geometries in which the photon source illuminates the line-forming region from below, and in which the photon source is embedded in the line-forming region. Our findings have implications for accretion-powered pulsars and gamma-ray bursters. In particular, the absence of pronounced shoulders on each side of the cyclotron first harmonic line in the spectra of accretion-powered pulsars suggests that the line-forming region is either illuminated from below and outside, as would be the case if it were plasma suspended in the magnetosphere of the neutron star, or it has a large ($\tau_{To} \gtrsim 10$) optical depth. Also, the ability of a slab line-forming region in which the magnetic field is parallel to the slab to produce narrow lines with large $W_E$ suggests that the lines observed in the X-ray spectra of some gamma-ray bursts might be able to be formed in plasma trapped at the magnetic equator of a neutron star.

We present results from a series of optical (g and r-band) and near-infrared (K'-band) observations of the region of the sky including the entire XTE and ROSAT error circles for the ``Bursting Pulsar'' GRO J1744-28. These data were taken with the Astrophysical Research Consortium's 3.5-m telescope at Apache Point Observatory and with the 2.2-m telescope at the European Southern Observatory. We see no new object, nor any significant brightening of any known object, in these error circles, with the exception of an object detected in our 8 February 1996 image. This object has already been proposed as a near-infrared counterpart to GRO J1744-28. While it is seen in only two of our ten 8 February frames, there is no evidence that this is an instrumental artifact, suggesting the possibility of near-infrared flares from GRO J1744-28, similar to those that have been reported from the Rapid Burster. The distance to the ``Bursting Pulsar'' must be more than 2 kpc, and we suggest that it is more than 7 kpc.

We investigate the viability of Galactic corona models of gamma-ray bursts by calculating the spatial distribution expected for a population of high-velocity neutron stars born in the Galactic disk and moving in a gravitational potential that includes the Galactic bulge, disk, and a dark matter halo. We consider models in which the bursts radiate isotropically and in which the radiation is beamed. We place constraints on the models by comparing the resulting brightness and angular distributions with the data in the BATSE 3B catalog. We find that, if the burst sources radiate isotropically, the Galactic corona model can reproduce the BATSE peak flux and angular distributions for neutron star kick velocities > 800 km s$^{-1}$, source turn-on ages > 20 Myrs, and BATSE sampling distances 130 kpc < d_{max} < 350 kpc. If the radiation is beamed, no turn-on age is required and agreement with the BATSE data can be found provided that the width of the beam is less than 20 deg.

We present the results of hydrodynamical simulations of nonaxisymmetric gas flow past a gravitating compact object in two dimensions. Calculations were performed with uniform flow as well as with transverse velocity and density gradients. We find that the flow is highly nonsteady, exhibiting the ``flip-flop'' behavior seen in previous studies in which accretion disks form with alternating directions of rotation. We investigate the periodicity of the flip-flop behavior, and study the effects of spatial resolution on the results. We find that the flip-flop motion creates accretion torques which, in some cases, may be large enough to explain the erratic spin behavior observed in some massive X-ray pulsars.

We argue that the bursting, transient X-ray source GRO J1744-28 is a binary consisting of a very low-mass ($M \approx 0.2 M_\sun$), highly evolved giant star that is transferring mass by Roche lobe overflow onto a high-mass ($M \approx 1.8 M_\sun$) neutron star. We explore a picture in which the bursts are due to thermonuclear flashes in matter that has accreted onto the neutron star. We attribute the unusually hard spectra of the bursts and the high burst rate compared to normal Type I X-ray burst sources -- indeed the existence of unstable nuclear burning itself -- to the unusually strong surface magnetic field, $ B_s \approx 1 \times 10^{13}$ G, of the neutron star.

The cyclotron line in the spectrum of the accretion-powered pulsar Her X-1 offers an opportunity to assess the ability of the BATSE Spectroscopy Detectors (SDs) to detect lines like those seen in some GRBs. Preliminary analysis of an initial SD pulsar mode observation of Her X-1 indicated a cyclotron line at an energy of approximately 44 keV, rather than at the expected energy of approximately 36 keV. Our analysis of four SD pulsar mode observations of Her X-1 made during high-states of its 35 day cycle confirms this result. We consider a number of phenomenological models for the continuum spectrum and the cyclotron line. This ensures that we use the simplest models that adequately describe the data, and that our results are robust. We find modest evidence (significance Q ~ 10^-4-10^-2) for a line at approximately 44 keV in the data of the first observation. Joint fits to the four observations provide stronger evidence (Q ~ 10^-7-10^-4) for the line. Such a shift in the cyclotron line energy of an accretion-powered pulsar is unprecedented.

We analyze the systematic errors in the positions of bursts in the BATSE 1B, 2B and 3B catalogs, using a likelihood approach. We use the BATSE data in conjunction with 196 single IPN arcs. We assume circular Gaussian errors, and that the total error is the sum in quadrature of the systematic error $\sigma_{\rm sys}$ and statistical error $\sigma_{\rm stat}$, as prescribed by the BATSE catalog. We find that the 3B burst positions are inconsistent with the value $\sigma_{\rm sys} = 1.6^\circ$ stated in the BATSE 3B catalog.

There is mounting evidence that, if gamma-ray bursters are Galactic in origin, they are located in a Galactic corona at distances greater than 100 kpc. This has created a need to explore new models of cyclotron line formation. In most previous calculations the line-forming region was modeled as a static slab of plasma, optically thin to continuum scattering, and threaded by a magnetic field of the order $10^{12}$ gauss oriented normal to the slab. Such a model is appropriate, for example, for the magnetic polar cap of a neutron star with a dipole field. However, if bursters lie at distances farther than several hundred parsecs, the burst luminosity exceeds the magnetic Eddington luminosity, and the plasma in a line-forming region at the magnetic polar cap would be ejected relativistically along the field lines. Mitrofanov and Tsygan have modeled the dynamics of such an outflow, and Miller {\it et al.} have calculated the properties of the cyclotron second and third harmonics, approximating them as due to cyclotron absorption. Here we describe Monte Carlo calculations of cyclotron resonant scattering at the first three harmonics in a relativistic outflow from the magnetic polar cap, and show that such scattering can produce narrow lines like those observed by Ginga.

We develop a likelihood methodology which can be used to search for evidence of burst repetition in the BATSE catalog, and to study the properties of the repetition signal. We use a simplified model of burst repetition in which a number $N_{\rm r}$ of sources which repeat a fixed number of times $N_{\rm rep}$ are superposed upon a number $N_{\rm nr}$ of non-repeating sources. The instrument exposure is explicitly taken into account. By computing the likelihood for the data, we construct a probability distribution in parameter space that may be used to infer the probability that a repetition signal is present, and to estimate the values of the repetition parameters. The likelihood function contains contributions from all the bursts, irrespective of the size of their positional errors --- the more uncertain a burst's position is, the less constraining is its contribution. Thus this approach makes maximal use of the data, and avoids the ambiguities of sample selection associated with data cuts on error circle size. We present the results of tests of the technique using synthetic data sets.

We investigate the viability of the Galactic halo model of $\gamma$-ray bursts by calculating the spatial distribution of neutron stars born with high velocities in the Galactic disk, and comparing the resulting brightness and angular distribution with the BATSE data. We find that the Galactic halo model can reproduce the BATSE peak flux and angular distribution data for neutron star kick velocities $\simgreat 800$~\kms, source turn-on ages $\simgreat 10$~Myrs, and sampling depths $100{\rm~kpc}\simless d_{max} \simless 400$~kpc.

We describe observational evidence and theoretical calculations which support the high velocity neutron star model of gamma-ray bursts. We estimate the energetic requirements in this model, and discuss possible energy sources. we also consider radiative processes involved in the bursts.

We investigate clustering in the angular distribution of the 260 $\gamma$-ray bursts in the publicly available BATSE catalogue, using a nearest neighbour analysis and the measures of burst brightness $B$ and short time scale variability $V$ which we introduced earlier. We find that while all 260 bursts are only modestly clustered (Q-value = $1.8 \times 10^{-2}$), the 202 bursts in this sample for which the statistical error in their locations is $<9^\circ$ are significantly clustered on an angular scale $\approx 5^\circ$ (Q-value = $2.5 \times 10^{-4}$, taking into acoount having chosen the cutoff in the statistical error). We also find a significant correlation between bright type I bursts and faint type I and type II bursts on an angular scale $\approx 5^\circ$ (Q-value = $4.0 \times 10^{-3}$). This angular scale is smaller than the typical (statistical plus systematic) error in burst locations of $6.8^\circ$, suggesting multiple recurrences from individual sources. We conclude that ``classical'' $\gamma$-ray burst sources repeat on a time scale of months, and that many faint type I and II bursts come from the sources of bright type I bursts.

We investigate the angular distribution of the $\gamma$-ray bursts in the publicly available BATSE catalogue, using the measures of burst brightness $B$ and short time scale ($\simless$ 0.3 s) variability $V$ which we introduced earlier. We show that the 54 type I ($\log V \le -0.8$) bursts lying in the middle brightness range 490 counts $\le B \le$ 1250 counts (corresponding to 1/3 of all type I bursts) exhibit a Galactic dipole moment of $\langle \cos \theta \rangle = 0.204 \pm 0.079$ {\it and} a deviation of the Galactic quadrupole moment from 1/3 of $\langle \sin^2 b \rangle - 1/3 = -0.104 \pm 0.041$. Using Monte Carlo simulations which include the BATSE sky exposure map and taking into account division of the type I bursts into three equal samples, we find that the probability by chance of an isotropic distribution of 54 bursts exhibiting values of $\langle \cos \theta \rangle$ {\it and} the negative of $\langle \sin^2 b \rangle - 1/3$ which equal or exceed the observed values is $6.6 \times 10^{-5}$. We conclude that $\gamma$-ray bursts are Galactic in origin.

We have analyzed the 241 bursts for which peak counts $\C$ exist in the publicly available Burst and Transient Source Experiment (BATSE) catalog. Introducing peak counts in 1024 ms as a measure of burst brightness $\B$ and the ratio of peak counts in 64 and 1024 ms as a measure of short timescale variability $\V$, we find a statistically significant correlation between the brightness and the short timescale variability of \g-ray bursts. The bursts which are smoother on short timescales are both faint and bright, while the bursts which are variable on short timescales are faint only, suggesting the existence of two distinct morphological classes of bursts.