Gamma-ray burst (GRB) data suggest that the jets from GRBs in the high redshift universe are more narrowly collimated than those at lower redshifts. This implies that we detect relatively fewer long GRB progenitor systems (i.e. massive stars) at high redshifts, because a greater fraction of GRBs have their jets pointed away from us. As a result, estimates of the star formation rate (from the GRB rate) at high redshifts may be diminished if this effect is not taken into account. In this paper, we estimate the star formation rate (SFR) using the observed GRB rate, accounting for an evolving jet opening angle. We find that the SFR in the early universe (z > 3) can be up to an order of magnitude higher than the canonical estimates, depending on the severity of beaming angle evolution and the fraction of stars that make long gamma-ray bursts. Additionally, we find an excess in the SFR at low redshifts, although this lessens when accounting for evolution of the beaming angle. Finally, under the assumption that GRBs do in fact trace canonical forms of the cosmic SFR, we constrain the resulting fraction of stars that must produce GRBs, again accounting for jet beaming-angle evolution. We find this assumption suggests a high fraction of stars in the early universe producing GRBs - a result that may, in fact, support our initial assertion that GRBs do not trace canonical estimates of the SFR.

Jet opening angles of long gamma-ray bursts (lGRBs) appear to evolve in cosmic time, with lGRBs at higher redshifts being on average more narrowly beamed than those at lower redshifts. We examine the nature of this anti-correlation in the context of collimation by the progenitor stellar envelope. First, we show that the data indicate a strong correlation between gamma-ray luminosity and jet opening angle, and suggest this is a natural selection effect - only the most luminous GRBs are able to successfully launch jets with large opening angles. Then, by considering progenitor properties expected to evolve through cosmic time, we show that denser stars lead to more collimated jets; we argue that the apparent anti-correlation between opening angle and redshift can be accounted for if lGRB massive star progenitors at high redshifts have higher average density compared to those at lower redshifts. This may be viable for an evolving IMF - under the assumption that average density scales directly with mass, this relationship is consistent with the form of the IMF mass evolution suggested in the literature. The jet angle-redshift anti-correlation may also be explained if the lGRB progenitor population is dominated by massive stars at high redshift, while lower redshift lGRBs allow for a greater diversity of progenitor systems (that may fail to collimate the jet as acutely). Overall, however, we find both the jet angle-redshift anti-correlation and jet angle-luminosity correlation are consistent with the conditions of jet launch through, and collimation by, the envelope of a massive star progenitor.

We examine the relationship between a number of long gamma-ray burst (lGRB) properties (isotropic emitted energy, luminosity, intrinsic duration, jet opening angle) and redshift. We find that even when accounting for conservative detector flux limits, there appears to be a significant correlation between isotropic equivalent energy and redshift, suggesting cosmological evolution of the lGRB progenitor. Analyzing a sub-sample of lGRBs with jet opening angle estimates, we find the beaming-corrected lGRB emitted energy does not correlate with redshift, but jet opening angle does. Additionally, we find a statistically significant anti-correlation between the intrinsic prompt duration and redshift, even when accounting for potential selection effects. We also find that - for a given redshift - isotropic energy is positively correlated with intrinsic prompt duration. None of these GRB properties appear to be correlated with galactic offset. From our selection-effect-corrected redshift distribution, we estimate a co-moving rate density for lGRBs, and compare this to the global cosmic star formation rate (SFR). We find the lGRB rate mildly exceeds the global star formation rate between a redshift of 3 and 5, and declines rapidly at redshifts above this (although we cannot constrain the lGRB rate above a redshift of about 6 due to sample incompleteness). We find the lGRB rate diverges significantly from the SFR at lower redshifts. We discuss both the correlations and lGRB rate density in terms of various lGRB progenitor models and their apparent preference for low-metallicity environments.

The current disparity in computational knowledge is a critical hindrance to the diversity and success of the field. Recommendations are outlined for policies and funding models to enable the growth and retention of a new generation of computational researchers that reflect the demographics of the undergraduate population in Astronomy and Physics.

The combined detection of a gravitational-wave signal, kilonova, and short gamma-ray burst (sGRB) from GW170817 marked a scientific breakthrough in the field of multi-messenger astronomy. But even before GW170817, there have been a number of sGRBs with possible associated kilonova detections. In this work, we re-examine these "historical" sGRB afterglows with a combination of state-of-the-art afterglow and kilonova models. This allows us to include optical/near-infrared synchrotron emission produced by the sGRB as well as ultraviolet/optical/near-infrared emission powered by the radioactive decay of $r$-process elements (i.e., the kilonova). Fitting the lightcurves, we derive the velocity and the mass distribution as well as the composition of the ejected material. The posteriors on kilonova parameters obtained from the fit were turned into distributions for the peak magnitude of the kilonova emission in different bands and the time at which this peak occurs. From the sGRB with an associated kilonova, we found that the peak magnitude in H bands falls in the range [-16.2, -13.1] ($95\%$ of confidence) and occurs within $0.8-3.6\,\rm days$ after the sGRB prompt emission. In g band instead we obtain a peak magnitude in range [-16.8, -12.3] occurring within the first $18\,\rm hr$ after the sGRB prompt. From the luminosity distributions of GW170817/AT2017gfo, kilonova candidates GRB130603B, GRB050709 and GRB060614 (with the possible inclusion of GRB150101B) and the upper limits from all the other sGRBs not associated with any kilonova detection we obtain for the first time a kilonova luminosity function in different bands.

Our understanding of the engines and progenitors of gamma-ray bursts has expanded through the ages as a broader set of diagnostics has allowed us to test our understanding of these objects. Here we review the history of the growth in our understanding, focusing on 3 leading engines and 9 potential progenitors. The first gravitational wave detection of a short burst is the latest in a series of breakthrough observations shaping this understanding and we study the importance of multi-diagnostic, multi-messenger observations on these engines and their progenitors. Our understanding based on a detailed study of nearby bursts can be applied to make predictions for the trends expected as we begin to observe high redshift bursts and we discuss these trends.

GHz radio astronomy has played a fundamental role in the recent dazzling discovery of GW170817, a neutron star (NS)-NS merger observed in both gravitational waves (GWs) and light at all wavelengths. Here we show how the expected progress in sensitivity of ground-based GW detectors over the next decade calls for U.S.-based GHz radio arrays to be improved beyond current levels. We discuss specifically how several new scientific opportunities would emerge in multi-messenger time-domain astrophysics if a next generation GHz radio facility with sensitivity and resolution $10\times$ better than the current Jansky Very Large Array (VLA) were to work in tandem with ground-based GW detectors. These opportunities include probing the properties, structure, and size of relativistic jets and wide-angle ejecta from NS-NS mergers, as well as unraveling the physics of their progenitors via host galaxy studies.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

Under the assumption that a Gamma-ray Burst (GRB) is ultimately produced by a Blandford-Znajek (BZ) jet from a highly spinning black hole BH, we put limits on the magnetic field and BH mass needed to power observed long and short GRBs. For BHs in the range of $2-10 M_{\odot}$ (for long GRBs) and $0.5-4 M_{\odot}$ (for short GRBs), we find magnetic fields in the range of $5x10^{14} \lesssim B \lesssim 10^{17}$ G are needed to power the observed GRBs. Under the simple assumption of flux conservation, we estimate the magnetic fields of the progenitor systems for both long and short GRBs, finding that single massive star progenitors require fields $\sim 10^{6}$ G and NS merger systems require fields $\sim 10^{15}$ G. We also discuss the implications and consequences of high magnetic fields in GRB BH-disk systems, in terms of MRI field growth and magnetically arrested disks. Finally, we examine the conditions under which the progenitor systems can retain enough angular momentum to create BHs spinning rapidly enough to power BZ jets.

We extend our study of energetic radio loud and quiet gamma-ray bursts (GRBs), suggesting these GRBs potentially come from two separate progenitor systems. We expand the sample from our previous paper (Lloyd-Ronning & Fryer, 2017) and find our results are strengthened - radio quiet GRBs have significantly shorter intrinsic prompt duration, and are also less energetic on average. However, the tenuous correlation between isotropic energy and intrinsic duration in the radio dark sample remains tenuous and is slightly weakened by adding more bursts. Interestingly, we find an anti-correlation between the intrinsic duration and redshift in the radio bright sample but not the radio dark sample, further supporting that these two samples may come from separate progenitors. We also find that very high energy (0.1 - 100 GeV) extended emission is only present in the radio loud sample. There is no significant difference between the presence of X-ray/optical plateaus or the average jet opening angles between the two samples. We explore the interpretation of these results in the context of different progenitor models. The data are consistent with the radio loud GRBs coming from a Helium-merger system and the radio quiet GRBs coming from a collapsar system, but may also reflect other dichotomies in the inner engine such as a neutron star versus black hole core.

We investigate the expected radio emission from the reverse shock of short GRBs, using the afterglow parameters derived from the observed short GRB light curves. In light of recent results suggesting that in some cases the radio afterglow is due to emission from the reverse shock, we examine the extent to which this component is detectable for short GRBs. In some GRBs, the standard synchrotron shock model predicts detectable radio emission from the reverse shock when none was seen. Because many physical parameters play a role in these estimates, our results highlight the need to more deeply explore the fundamental processes involved in GRB particle acceleration and emission. However, with more rapid follow-up, we can test our standard model of GRBs, which predicts an early, radio bright reverse shock in many cases.

We investigate the expected radio emission from short GRBs, guided by the observed distributions of their afterglow parameters. Our motivation is to explore the radio signatures of the forward and reverse shock in a sGRB jet, as well as the off-axis jet component and radio emission from the dynamical ejecta. We estimate the fraction of GRBs detectable in the radio that may be coincident with a GW signal from a neutron star (NS) merger. At the distance to which aLIGO at design specification is sensitive to GWs from NS mergers, we may detect up to 20 % of off-axis jet radio components of these events, and up to about 10 % of events with radio emission from the tidal ejecta shocking with the external medium. We discuss our results in light of the Next Generation Very Large Array, expected to make significant progress in detecting radio signals from sGRBs.

The simultaneous detection of gravitational and electromagnetic waves from a binary neutron star merger has both solidified the link between neutron star mergers and short-duration gamma-ray bursts (GRBs) and demonstrated the ability of astronomers to follow-up the gravitational wave detection to place constraints on the ejecta from these mergers as well as the nature of the GRB engine and its surroundings. As the sensitivity of aLIGO and VIRGO increases, it is likely that a growing number of such detections will occur in the next few years, leading to a sufficiently-large number of events to constrain the populations of these GRB events. While long-duration GRBs originate from massive stars and thus are located near their stellar nurseries, binary neutron stars may merge on much longer timescales, and thus may have had time to migrate appreciably. The strength and character of the electromagnetic afterglow emission of binary neutron star mergers is a sensitive function of the circum-merger environment. Though the explosion sites of short GRBs have been explored in the literature, the question has yet to be fully addressed in its cosmological context. We present cosmological simulations following the evolution of a galaxy cluster including star formation combined with binary population synthesis models to self-consistently track the locations and environmental gas densities of compact binary merger sites throughout the cosmic web. We present probability distributions for densities as a function of redshift and discuss model sensitivity to population synthesis model assumptions.

We offer an explanation for the anti-correlation between the minimum variability timescale ($MTS$) in the prompt emission light curve of gamma-ray bursts (GRBs) and the estimated bulk Lorentz factor of these GRBs, in the context of a magnetically arrested disk (MAD) model. In particular, we show that previously derived limits on the maximum available energy per baryon in a Blandford-Znajek jet leads to a relationship between the characteristic MAD timescale, $t_{MAD}$, in GRBs and the maximum bulk Lorentz factor: $t_{MAD} \propto \Gamma^{-6}$, somewhat steeper than (although within the error bars of) the fitted relationship found in the GRB data. Similarly, the MAD model also naturally accounts for the observed anti-correlation between $MTS$ and gamma-ray luminosity $L$ in the GRB data, and we estimate the accretion rates of the GRB disks (given these luminosities) in the context of this model. Both of these correlations ($MTS-\Gamma$ and $MTS-L$) are also observed in the AGN data, and we discuss the implications of our results in the context of both GRB and blazar systems.

We investigate the intrinsic properties of a sample of bright (E_iso > 10^52 erg) gamma-ray bursts, comparing those with and without radio afterglows. We find that the sample of bursts with no radio afterglows has a significantly shorter mean intrinsic duration of the prompt gamma-ray radiation, and the distribution of this duration is significantly different from those bursts with a radio afterglow. Although the sample with no radio afterglow has on average lower isotropic energy, the lack of radio afterglow does not appear to be a result of simply energetics of the burst, but a reflection of a separate physical phenomenon likely related to the circumburst density profile. We also find a weak correlation between the isotropic gamma-ray energy and intrinsic duration in the sample with no radio afterglow, but not in the sample which have observed radio afterglows. We give possible explanations for why there may exist a sample of GRBs with no radio afterglow depending on whether the radio emission comes from the forward or reverse shock, and why these bursts appear to have intrinsically shorter prompt emission durations. We discuss how our results may have implications for progenitor models of GRBs.

We present a model for the temporal variability of long gamma-ray bursts during the prompt phase (the highly variable first 100 seconds or so), in the context of a magnetically arrested disk (MAD) around a black hole. In this state, sufficient magnetic flux is held on to the black hole such that it stalls the accretion near the inner region of the disk. The system transitions in and out of the MAD state, which we relate to the variable luminosity of the GRB during the prompt phase, with a characteristic timescale defined by the free fall time in the region over which the accretion is arrested. We present simple analytic estimates of the relevant energetics and timescales, and compare them to gamma-ray burst observations. In particular, we show how this model can reproduce the characteristic one second time scale that emerges from various analyses of the prompt emission light curve. We also discuss how our model can accommodate the potentially physically important correlation between a burst quiescent time and the duration of its subsequent pulse (Ramirez-Ruiz & Merloni 2001).

We use the Los Alamos VPIC code to investigate particle acceleration in relativistic, unmagnetized, collisionless electron-ion plasmas. We run our simulations both with a realistic proton-to-electron mass ratio m_p/m_e = 1836, as well as commonly employed mass ratios of m_p/m_e =100 and 25, and show that results differ among the different cases. In particular, for the physically accurate mass ratio, electron acceleration occurs efficiently in a narrow region of a few hundred inertial lengths near the flow front, producing a power law dN/dgamma ~ gamma^(-p) with p ~ -2 developing over a few decades in energy, while acceleration is weak in the region far downstream. We find 20%, 10%, and 0.2% of the total energy given to the electrons for mass ratios of 25, 100, and 1836 respectively at a time of 2500 (w_p)^-1. Our simulations also show significant magnetic field generation just ahead of and behind the the flow front, with about 1% of the total energy going into the magnetic field for a mass ratio of 25 and 100, and 0.1% for a mass ratio of and 1836. In addition, lower mass ratios show significant fields much further downstream than in the realistic mass ratio case. Our results suggest the region and energetic extent of particle acceleration is directly related to the presence of magnetic field generation. Our work sheds light on the understanding of particle acceleration and emission in gamma-ray bursts, among other relativistic astrophysical outflows, but also underscores the necessity of optimizing numerical and physical parameters, as well as comparing among PIC codes before firm conclusions are drawn from these types of simulations.

Magnetic reconnection is a leading mechanism for dissipating magnetic energy and accelerating nonthermal particles in Poynting-flux dominated flows. In this letter, we investigate nonthermal particle acceleration during magnetic reconnection in a magnetically-dominated ion-electron plasma using fully kinetic simulations. For an ion-electron plasma with the total magnetization $\sigma_0=B^2/(4\pi n(m_i+m_e)c^2)$, the magnetization for each species is $\sigma_i \sim \sigma_0$ and $\sigma_e \sim (m_i/m_e) \sigma_0$, respectively. We have studied the magnetically dominated regime by varying $\sigma_{e} = 10^3 - 10^5$ with initial ion and electron temperatures $T_i = T_e = 5 - 20 m_ec^2$ and mass ratio $m_i/m_e = 1 - 1836$. The results demonstrate that reconnection quickly establishes power-law energy distributions for both electrons and ions within several ($2-3$) light-crossing times. For the cases with periodic boundary conditions, the power-law index is $1<s<2$ for both electrons and ions. The hard spectra limit the power-law energies for electrons and ions to be $\gamma_{be} \sim \sigma_e$ and $\gamma_{bi} \sim \sigma_i$, respectively. The main acceleration mechanism is a Fermi-like acceleration through the drift motions of charged particles. When comparing the spectra for electrons and ions in momentum space, the spectral indices $s_p$ are identical as predicted in Fermi acceleration. We also find that the bulk flow can carry a significant amount of energy during the simulations. We discuss the implication of this study in the context of Poynting-flux dominated jets and pulsar winds especially the applications for explaining the nonthermal high-energy emissions.

Gamma-ray bursts (GRBs) are the ultimate cosmic lighthouses, capable of illuminating the universe at its earliest epochs. Could such events probe the properties of the first stars at z $\sim$ 20, the end of the cosmic Dark Ages? Previous studies of Population III GRBs only considered explosions in the diffuse relic H II regions of their progenitors, or bursts that are far more more energetic than those observed to date. But the processes that produce GRBs at the highest redshifts likely reset their local environments, creating much more complicated structures than those in which relativistic jets have been modeled so far. These structures can greatly affect the luminosity of the afterglow, and hence the redshift at which it can be detected. We have now simulated Population III GRB afterglows in H II regions, winds, and dense shells ejected by the star during the processes that produce the burst. Our model, which has been used in previous work, has been extended to include contributions from reverse shocks, inverse Compton cooling and the effects of sphericity and beaming in the blast wave, and is valid in a variety of circumjet density profiles. We find that GRBs with E$_{\mathrm{iso},\gamma} =$ 10$^{ 51}$ - 10$^{53}$ erg will be visible at z $\gtrsim$ 20 to the next generation of near infrared and radio observatories. In many cases, the environment of the burst, and hence progenitor type, can be inferred from the afterglow light curve. Although some Population III GRBs are visible to Swift and the Very Large Array now, the optimal strategy for their detection will be future missions like EXIST and JANUS, which have large survey areas and onboard X-ray and infrared telescopes that can track their near infrared flux from the moment of the burst, thereby identifying its redshift.

It is now generally accepted that long-duration gamma ray bursts (GRBs) are due to the collapse of massive rotating stars. The precise collapse process itself, however, is not yet fully understood. Strong winds, outbursts, and intense ionizing UV radiation from single stars or strongly interacting binaries are expected to destroy the molecular cloud cores that give birth to them and create highly complex circumburst environments for the explosion. Such environments might imprint features on GRB light curves that uniquely identify the nature of the progenitor and its collapse. We have performed numerical simulations of realistic environments for a variety of long-duration GRB progenitors with ZEUS-MP, and have developed an analytical method for calculating GRB light curves in these profiles. Though a full, three-dimensional, relativistic magnetohydrodynamical computational model is required to precisely describe the light curve from a GRB in complex environments, our method can provide a qualitative understanding of these phenomena. We find that, in the context of the standard afterglow model, massive shells around GRBs produce strong signatures in their light curves, and that this can distinguish them from those occurring in uniform media or steady winds. These features can constrain the mass of the shell and the properties of the wind before and after the ejection. Moreover, the interaction of the GRB with the circumburst shell is seen to produce features that are consistent with observed X-ray flares that are often attributed to delayed energy injection by the central engine. Our algorithm for computing light curves is also applicable to GRBs in a variety of environments such as those in high-redshift cosmological halos or protogalaxies, both of which will soon be targets of future surveys such as JANUS or Lobster.

Using measured X-ray luminosities to 17 Gamma-Ray Bursts (GRBs) during the afterglow phase and accounting for radiative losses, we calculate the kinetic energy of these bursts and investigate its relation to other GRB properties. We then use the observed radiated energy during the prompt phase to determine the radiative efficiency of these bursts, and explore how the efficiency relates to other GRB observables. We find that the kinetic energy in the afterglow phase is directly correlated with the radiated energy, total energy as well as possibly the jet opening angle and spectral peak energy. More importantly, we find the intriguing fact that the efficiency is correlated with the radiated energy, and mildly with the total energy, jet opening angle and spectral peak energy. XRF020903 also seems to follow the trends we find for our GRB sample. We discuss the implications of these results for the GRB radiation and jet models.

An observed correlation Ep \propto (E_{iso})^{1/2} extending from Gamma-ray Bursts (GRBs) to X-ray flashes (XRFs) poses problems both for a power-law universal jet model where the energy per solid angle decreases as the inverse square of the angle respect to the jet axis, and for a conical jet model with a uniform energy density within the jet beam and a sharp energy cut-off at the jet edge. Here we show that the current GRB-XRF prompt emission/afterglow data can be understood in terms of a picture in which the GRB-XRF jets are quasi-universal and structured, with a Gaussian-like or similar structure, i.e., one where the jet has a characteristic angle, with a mild variation of energy inside and a rapid (e.g. exponential) decrease of energy outside of it. A Monte Carlo simulation shows that the current data is compatible with such a quasi-universal Gaussian jet with a typical opening angle of 5.7(+3.4)(-2.1) degrees, and with a standard jet energy of about log (E_j/1 erg)= 51.1+-0.3. According to this model, the true-to-observed number ratio of the whole GRB-XRF population is about 14 with the current instrumental sensitivity.

Two separate topics are discussed. (1) We describe the classifications of the long GRB early afterglow lightcurves within the framework of the fireball shock model, focusing on the interplay between the reverse and forward shock emission components. We will also provide evidence that the central engine of at least two bursts are entrained with strong magnetic fields, and discuss the implications of this result for our understanding of the GRB phenomenon; (2) We argue that the current gamma-ray burst (GRB) and X-ray flash (XRF) data are consistent with a picture that all GRB-XRF jets are structured and quasi-universal, with a typical Gaussian-like jet structure.

The idea that GRBs originate from uniform jets has been used to explain numerous observations of breaks in the GRB afterglow lightcurves. We explore the possibility that GRBs instead originate from a structured jet that may be quasi-universal, where the variation in the observed properties of GRBs is due to the variation in the observer viewing angle. We test how various models reproduce the jet data of Bloom, Frail, & Kulkarni (2003), which show a negative correlation between the isotropic energy output and the inferred jet opening angle (in a uniform jet configuration). We find, consistent with previous studies, that a power-law structure for the jet energy as a function of angle gives a good description. However, a Gaussian jet structure can also reproduce the data well, particularly if the parameters of the Gaussian are allowed some scatter. We place limits on the scatter of the parameters in both the Gaussian and power-law models needed to reproduce the data, and discuss how future observations will better distinguish between these models for the GRB jet structure. In particular, the Gaussian model predicts a turnover at small opening angles and in some cases a sharp cutoff at large angles, the former of which may already have been observed. We also discuss the predictions each model makes for the observed luminosity function of GRBs and compare these predictions with the existing data.

The variety of gamma-ray burst phenomenology could be largely attributable to differences in the opening angle of an isotropic outflow or to a standard type of event viewed from different orientations. Motivated by this currently popular idea, we study the effects of varying the energy per unit solid angle in the unsteady expelled outflow by an increase either in the bulk Lorentz factor or in the baryon loading. We apply these models to interpret the observed correlations between variability, luminosity and spectral peak energy and find that while the latter scenario fails to provide a good description, bursts produced by collisions between similar mass shells but with increasingly large Lorentz factors are both more variable and have larger peak spectral energies. We present detailed internal shock calculations confirming this interpretation and discuss the roles various timescales, radii and the optical thickness of the wind play in determining this wide range of behaviors. Finally, we discuss the variety of scenarios in which large variations of the source expansion velocity are naturally expected.

The variable activity of a gamma-ray burst (GRB) source is thought to be correlated with its absolute peak luminosity - a relation that, if confirmed, can be used to derive an independent estimate of the redshift of a GRB. We find that bursts with highly variable light curves have greater nu-F_nu spectral peak energies, when we transform these energies to the cosmological rest frame using the redshift estimates derived either from optical spectral features or from the luminosity-variability distance indicator itself. This positive correlation between peak energy and variability spans ~ 2 orders of magnitude and appears to accommodate GRB 980425, lending credibility to the association of this burst with SN 1998bw. The existence of such a correlation not only provides an interesting clue to the nature of this luminosity indicator but potentially reinforces the validity of the redshift estimates derived from this method. It also implies that the rest frame GRB peak energy is correlated with the intrinsic luminosity of the burst, as has been suggested in the past as an explanation of the observed hardness-intensity correlation in GRBs.

Using 220 Gamma-Ray Burst (GRB) redshifts and luminosities derived from the luminosity-variability relationship of Fenimore & Ramirez-Ruiz (2000), we show that there exists a significant correlation between the GRB luminosity and redshift. In particular, we find that the evolution of the average luminosity can be parameterized as L ~ (1+z)^(1.4 +- 0.5), where z is the burst redshift. We discuss the possible reasons behind this evolution and compare it to other known sources that exhibit similar behavior. In addition, we use non-parametric statistical techniques to independently estimate the distributions of the luminosity and redshift of bursts, accounting for the evolution (in contrast to previous studies which have assumed that the luminosity function is independent of redshift). We present these distributions and discuss their implications. Most significantly, we find a co-moving rate density of GRBs that continues to increase to (1+z) ~ 10. From this estimate of the GRB rate density, we then use the population synthesis codes of Fryer et al. (1999) to estimate the star formation rate at high redshifts, based on different progenitor models of GRBs. We find that no matter what the progenitor or population synthesis model, the star formation rate increases or remains constant to very high redshifts (z ~ 10).

We discuss the prompt emission of Gamma-Ray Bursts in different spectral energy bands. First, we suggest that a three-part synchrotron emission model is a good description of the ~20 keV - 1 MeV gamma-ray emission of GRBs. We show that this model provides excellent fits to the data and naturally explains the observed global correlations between spectral parameters. In particular, we show there exists a negative correlation between between the peak of the nu-Fnu spectrum, Ep, and the low energy photon index \alpha for bursts with -2/3 < \alpha < 0, and suggest that this correlation is due to the mechanism responsible for producing \alpha's above the value of -2/3 - namely, a decreasing mean pitch angle of the electrons. We then discuss the physical origin of the increasing number of GRBs that are observed to peak in the X-ray energy band (~5-40 keV). Although either a cosmological (i.e. high redshift) or intrinsic interpretation for the low values of Ep is viable at this point, the data appear to suggest that intrinsic effects are playing the dominant role. Finally, we briefly comment on the prompt GRB optical emission (~ eV) and very high energy emission (>10 MeV), and how these spectral bands may be used to place additional constraints on the physics of gamma-ray bursts.

In this paper, we explore time resolved Gamma-Ray Burst (GRB) spectra in the context of the synchrotron emission model presented in Lloyd and Petrosian (2000; LP00). First, we show that our model - which involves three distinct emission regimes - can provide excellent fits to the time resolved spectra of GRBs, and we present these results for a few bursts. We then describe how the phenomenological Band spectrum (Band et al., 1993) can be interpreted in the context of our models based on the value of the low energy photon index $\alpha$. We discuss the types of correlations one would expect to observe among the Band parameters if these models are correct. We then compare these predictions to the existing data, combining a sample of 2,026 time resolved spectra (from approximately 80 bursts). We show that the correlations found in the data are consistent with the models, and discuss the constraints they place on the emission physics. In particular, we find a 4-sigma negative correlation between the peak of the nu-Fnu spectrum, Ep, and the low energy photon index $\alpha$ for bursts with $-2/3 < \alpha < 0$, in contrast to what is predicted by the instrumental effect discussed in LP00. We suggest that this correlation is simply due to the mechanism responsible for producing $\alpha$'s above the value of -2/3 - namely, a decreasing mean pitch angle of the electrons. We also show that Ep is correlated with the photon flux, and interpret this as a result of changing magnetic field or characteristic electron energy between emission episodes. Finally, we discuss the implications our results have on particle acceleration in GRBs, and prospects for further testing these models with the anticipated data from HETE-2, Swift and GLAST.

We investigate synchrotron emission models as the source of gamma-ray burst spectra. We show that including the possibility for synchrotron self-absorption, a ``smooth cutoff'' to the electron energy distribution, and an anisotropic distribution for the electron pitch angles produces a whole range of low energy spectral behavior. In addition, we show that the procedure of spectral fitting to GRB data over a finite bandwidth can introduce a spurious correlation between spectral parameters - in particular, the value of the peak of the nu F_nu spectrum, E_p, and the low energy photon spectral index alpha (the lower E_p is, the lower (softer) the fitted value of alpha will be). From this correlation and knowledge of the E_p distribution, we show how to derive the expected distribution of alpha. We show that optically thin synchrotron models with an isotropic electron pitch angle distribution can explain the distribution of alpha below alpha=-2/3. This agreement is achieved if we relax the unrealistic assumption of the presence of a sharp low energy cutoff in the spectrum of accelerated electrons, and allow for a more gradual break. We show that this low energy portion of the electron spectrum can be at most flat. We also show that optically thin synchrotron models with an anisotropic electron pitch angle distribution can explain all bursts with -2/3 < alpha <= 0$. The very few bursts with low energy spectral indices that fall above alpha=0 may be due the presence of a the synchrotron self-absorption frequency entering the lower end of the BATSE window. Our results also predict a particular relationship between alpha and E_p during the temporal evolution of a GRB. We give examples of spectral evolution in GRBs and discuss how the behavior are consistent with the above models.

We test the models of synchrotron emission presented in Part I of this series (Lloyd & Petrosian, these proceedings) against the distributions and evolution of GRB spectral parameters (particularly the low energy index, $\alpha$). With knowledge of the $E_{p}$ distribution and the correlation between $\alpha$ and $E_{p}$ presented in Part I, we show how to derive the expected distribution of $\alpha$ from fits to optically thin synchrotron spectra, and compare this with the observed distribution. We show that there is no difficulty explaining bursts below the ``line of death'', $\alpha < -2/3$, and that these bursts indicate that the spectrum of accelerated electrons must flatten or decline at low energies. Bursts with low energy spectral indices that fall above this limit are explained by the synchrotron self-absorption frequency entering the lower end of the BATSE window. Finally, we discuss a variety of spectral evolution behavior seen in GRBs and explain this behavior in the context of synchrotron emission from internal shocks.

We investigate synchrotron emission models as the source of gamma ray burst spectra. We show that allowing for synchrotron self absorption and a ``smooth cutoff'' to the electron energy distribution produces a wide range of low energy spectral behavior. We show that there exists a correlation between the value of the peak of the $\nu F_{\nu}$ spectrum, $E_{p}$, and the low energy spectral index $\alpha$ as determined by spectral fits over a finite bandwidth. Finally, we discuss the implications of synchrotron emission from internal shocks for GRB spectral evolution.

We present results of correlation studies, examining the association between the peak of the nu F_nu spectrum of gamma ray bursts, E_p, with the burst's energy fluence and photon peak flux. We discuss methods to account for data truncation in E_p and fluence or flux when performing the correlation analyses. However, because bursts near the detector threshold are not usually able to provide reliable spectral parameters, we focus on results for the brightest bursts in which we can better understand the selection effects relevant to E_p and burst strength. We find that there is a strong correlation between total fluence and E_p. We discuss these results in terms of both cosmological and intrinsic effects. In particular, we show that for realistic distributions of the burst parameters, cosmological expansion alone cannot account for the correlation between E_p and total fluence; the observed correlation is likely a result of an intrinsic relation between the burst rest-frame peak energy and the total radiated energy. We investigate this latter scenario in the context of synchrotron radiation from external and internal shock models of GRBs. We find that the internal shock model is consistent with our interpretation of the correlation, while the external shock model cannot easily explain this intrinsic relation between peak energy and burst radiated energy.

There have been several suggestions of the existence of cosmological redshift signatures in the temporal and spectral characteristics of gamma-ray bursts. However, recent discoveries of afterglows and redshift measurements indicate the presence of broad ``luminosity functions'' which may overwhelm such weaker cosmological signatures. The primary goal of this paper is to determine if the intrinsic and cosmological dispersions can be separated. We have expanded the search for cosmological signatures to several other temporal and spectral features, have determined correlations which could arise from the cosmological redshift of the sources and have carried out tests to determine if the observed correlations can be due to cosmology. We find that the intrinsic dispersions are the dominant factors and that detection of cosmological signatures must await accumulation of a much larger number of identification with galaxies and measurements of redshifts.

We investigate the cosmological evolution of GRBs, using the total gamma ray fluence as a measure of the burst strength. This involves an understanding of the distributions of the spectral parameters of GRBs as well as the total fluence distribution - both of which are subject to detector selection effects. We present new non-parametric statistical techniques to account for these effects, and use these methods to estimate the true distribution of the peak of the nu F_nu spectrum, E_p, from the raw distribution. The distributions are obtained from four channel data and therefore are rough estimates. Here, we emphasize the methods and present qualitative results. Given its spectral parameters, we then calculate the total fluence for each burst, and compute its cumulative and differential distributions. We use these distributions to estimate the cosmological rate evolution of GRBs, for three cosmological models. Our two main conclusions are the following: 1) Given our estimates of the spectral parameters, we find that there may exist a significant population of high E_p bursts that are not detected by BATSE, 2) We find a GRB co-moving rate density quite different from that of other extragalactic objects; in particular, it is different from the recently determined star formation rate.

We argue that recent observations of afterglows and some theoretical considerations indicate that the gamma-ray fluence may be a better measure of strength of a burst than its peak flux. We discuss how the distribution of the fluence, or any physical quantity related to the peak counts (or flux), can be obtained with proper correction for threshhold effects. Using this method, we compare the cumulative and differential distributions of flux and fluence. We find some differences between these distributions, but more remarkable are the similarities between the distributions. The other striking feature is how different these distributions are from those of other extragalactic objects. Using the fluence distributions, we derive the expected comoving density rate evolution for different assumed total luminosity; we compare this with the GRB rate expected from the recently determined star formation rate.

The distributions of spectral characteristics and their correlations with fluence, peak flux, or duration, are essential in understanding the nature of GRBs. However, the selection effects involved in detecting GRBs can distort these distributions. Here, we discuss how to deal with selection effects involving the peak energy E_p of the GRB nu F_nu spectrum, which suffers from both an upper and lower threshold. We describe a new method to account for this double-sided truncation, and show that the true distribution of E_p is significantly different from the observed distribution.