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.

We present a new polynomial-free prolongation scheme for Adaptive Mesh Refinement (AMR) simulations of compressible and incompressible computational fluid dynamics. The new method is constructed using a multi-dimensional kernel-based Gaussian Process (GP) prolongation model. The formulation for this scheme was inspired by the GP methods introduced by A. Reyes et al. (A New Class of High-Order Methods for Fluid Dynamics Simulation using Gaussian Process Modeling, Journal of Scientific Computing, 76 (2017), 443-480; A variable high-order shock-capturing finite difference method with GP-WENO, Journal of Computational Physics, 381 (2019), 189-217). In this paper, we extend the previous GP interpolations and reconstructions to a new GP-based AMR prolongation method that delivers a high-order accurate prolongation of data from coarse to fine grids on AMR grid hierarchies. In compressible flow simulations special care is necessary to handle shocks and discontinuities in a stable manner. To meet this, we utilize the shock handling strategy using the GP-based smoothness indicators developed in the previous GP work by A. Reyes et al. We demonstrate the efficacy of the GP-AMR method in a series of testsuite problems using the AMReX library, in which the GP-AMR method has been implemented.

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.

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 gamma-ray burst (GRB) rate is essential for revealing the connection between GRBs, supernovae and stellar evolution. Additionally, the GRB rate at high redshift provides a strong probe of star formation history in the early universe. While hundreds of GRBs are observed by Swift, it remains difficult to determine the intrinsic GRB rate due to the complex trigger algorithm of Swift. Current studies of the GRB rate usually approximate the Swift trigger algorithm by a single detection threshold. However, unlike the previously flown GRB instruments, Swift has over 500 trigger criteria based on photon count rate and additional image threshold for localization. To investigate possible systematic biases and explore the intrinsic GRB properties, we develop a program that is capable of simulating all the rate trigger criteria and mimicking the image threshold. Our simulations show that adopting the complex trigger algorithm of Swift increases the detection rate of dim bursts. As a result, our simulations suggest bursts need to be dimmer than previously expected to avoid over-producing the number of detections and to match with Swift observations. Moreover, our results indicate that these dim bursts are more likely to be high redshift events than low-luminosity GRBs. This would imply an even higher cosmic GRB rate at large redshifts than previous expectations based on star-formation rate measurements, unless other factors, such as the luminosity evolution, are taken into account. The GRB rate from our best result gives a total number of 4571^{+829}_{-1584} GRBs per year that are beamed toward us in the whole universe. SPECIAL NOTE (2015.05.16): This new version incorporates an erratum. All the GRB rate normalizations ($R_{\rm GRB}(z=0)$) should be a factor of 2 smaller than previously reported. Please refer to the Appendix for more details. We sincerely apologize for the mistake.

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 apply the mathematical formalism of vector spherical harmonics decomposition to convective stellar velocity fields from multi-dimensional hydrodynamics simulations, and show that the resulting power spectra furnish a robust and stable statistical description of stellar convective turbulence. Analysis of the power spectra help identify key physical parameters of the convective process such as the dominant scale of the turbulent motions that influence the structure of massive evolved pre-supernova stars. We introduce the numerical method that can be used to calculate vector spherical harmonics power spectra from 2D and 3D convective shell simulation data. Using this method we study the properties of oxygen shell burning and convection for a 15 Msun star simulated by the hydrodynamics code FLASH in 2D and 3D. We discuss the importance of realistic initial conditions to achieving successful core-collapse supernova explosions in multi-dimensional simulations. We show that the calculated power spectra can be used to generate realizations of the velocity fields of pre-supernova convective shells. We find that the slope of the solenoidal mode power spectrum remains mostly constant throughout the evolution of convection in the oxygen shell in both 2D and 3D simulations. We also find that the characteristic radial scales of the convective elements are smaller in 3D than in 2D while the angular scales are larger in 3D.

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

Self-gravity computation by multipole expansion is a common approach in problems such as core-collapse and Type Ia supernovae, where single large condensations of mass must be treated. The standard formulation of multipole self-gravity suffers from two significant sources of error, which we correct in the formulation presented in this article. The first source of error is due to the numerical approximation that effectively places grid cell mass at the central point of the cell, then computes the gravitational potential at that point, resulting in a convergence failure of the multipole expansion. We describe a new scheme that avoids this problem by computing gravitational potential at cell faces. The second source of error is due to sub-optimal choice of location for the expansion center, which results in angular power at high multipole $l$ values in the gravitational field, requiring a high --- and expensive --- value of multipole cutoff \lmax. By introducing a global measure of angular power in the gravitational field, we show that the optimal coordinate for the expansion is the square-density-weighted mean location. We subject our new multipole self-gravity algorithm to two rigorous test problems: MacLaurin spheroids for which exact analytic solutions are known, and core-collapse supernovae. We show that key observables of the core-collapse simulations, particularly shock expansion, proto-neutron star motion, and momentum conservation, are extremely sensitive to the accuracy of the multipole gravity, and the accuracy of their computation is greatly improved by our reformulated solver.

The long gamma-ray burst (GRB) rate is essential for revealing the connection between GRBs, supernovae and stellar evolution. Additionally, the GRB rate at high redshift provides a strong probe of star formation history in the early universe. While hundreds of GRBs are observed by Swift, it remains difficult to determine the intrinsic GRB rate due to the complex trigger algorithm of Swift. Current studies usually approximate the Swift trigger algorithm by a single detection threshold. However, unlike the previously flown GRB instruments, Swift has over 500 trigger criteria based on photon count rate and additional image threshold for localization. To investigate possible systematic biases and explore the intrinsic GRB properties, we developed a program that is capable of simulating all the rate trigger criteria and mimicking the image trigger threshold. We use this program to search for the intrinsic GRB rate. Our simulations show that adopting the complex trigger algorithm of Swift increases the detection rate of dim bursts. As a result, we find that either the GRB rate is much higher than previously expected at large redshift, or the luminosity evolution is non-negligible. We will discuss the best results of the GRB rate in our search, and their impact on the star-formation history.

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.

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.

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 2002 discovery of the "Amati Relation" of GRB spectra created the possibility that this and other correlations of GRB phenomenology might be used to make GRBs into standard candles. One recurring apparent difficulty with this program has been that some of the primary observational quantities to be fit as "data" - the isotropic-equivalent prompt energy $E_{iso}$ and the collimation-corrected "total" prompt energy energy $E_{\gamma}$ - depend for their construction on the very cosmological models that they are supposed to help constrain. This is the so-called "circularity problem" of standard candle GRBs. This paper is intended to point out that the circularity problem is not in fact a problem at all, except to the extent that it amounts to a self-inflicted wound. It arises essentially because of an unfortunate choice of data variables, such as $E_{iso}$, which are unnecessarily model-dependent. If, instead, the empirical correlations of GRB phenomenology which are formulated in source-variables are {\it mapped to the primitive observational variables} (such as fluence) and compared to the observations in that space, then all circularity disappears. I also indicate here a set of procedures for encoding high-dimensional empirical correlations in a "Gaussian Tube" smeared model that includes both the correlation and its intrinsic scatter, and how that source-variable model may easily be mapped to the space of primitive observables and fashioned into a likelihood. I discuss the projections of such Gaussian tubes into sub-spaces, which may be used to incorporate data from GRB events that may lack some element of the data (for example, GRBs without ascertained jet-break times). In this way, a large set of inhomogeneously observed GRBs may be assimilated into a single analysis, so long as each possesses at least two correlated data attributes.

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.

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

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.

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.

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.

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

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.

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

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

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.

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.

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

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

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.

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

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.

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.

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.

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.

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.

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.

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.

The traditional design of trigger algorithms for GRB experiments requires the specification of the background and burst samples in terms of acquisition times that are of fixed duration and of fixed elapsed time from each other. One such set of acquisition times is required for each characteristic timescale of GRB variation that one desires to detect. One then slides each set through the trigger data searching for samples that maximize the signal-to-noise of the background-subtracted burst sample. Here we describe a new triggering approach in which the times at which the background and burst samples are acquired are allowed to vary dynamically. Two background samples bracket a burst sample. The background and burst durations and elapsed time between them are allowed to be free parameters, which are maximized using the downhill simplex method. This produces great flexibility in the timescales that are available for detecting GRBs.

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

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.

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.

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 report the discovery of the optical and radio afterglow of GRB 010921, the first gamma-ray burst afterglow to be found from a localization by the High Energy Transient Explorer (HETE) satellite. We present optical spectroscopy of the host galaxy which we find to be a dusty and apparently normal star-forming galaxy at z = 0.451. The unusually steep optical spectral slope of the afterglow can be explained by heavy extinction, A_V > 0.5 mag, along the line of sight to the GRB. Dust with similar A_V for the the host galaxy as a whole appears to be required by the measurement of a Balmer decrement in the spectrum of the host galaxy. Thanks to the low redshift, continued observations of the afterglow will enable the strongest constraints, to date, on the existence of a possible underlying supernova.

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

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

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.

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

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 study the resonant divergences that occur in quantum scattering cross-sections in the presence of a strong external magnetic field. We demonstrate that all such divergences may be eliminated by introducing radiative corrections to the leading-order scattering amplitudes. These corrections impose a choice of basis states that must be used in scattering calculations: electron states must diagonalize the mass operator, while photon states must diagonalize the polarization operator. The radiative corrections introduce natural line-widths into the energy denominators of all propagators, as well as into the time-development exponentials of all scattering states corresponding to external lines. Since initial and final scattering states may now decay, it is logically necessary to compute scattering amplitudes for a finite time-lapse between the preparation of the initial state and the measurement of the final state. Strict energy conservation, which appeared in previous formulations of the theory, must thus be abandoned. We exhibit the generic formulae for the scattering cross-sections in two useful limits, corresponding to the cases where either the initial states or the final states are stable, and discuss the application of the general formula when neither of these limits applies.

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.