We introduce a simple model to explain the time-reversed and stretched residuals in gamma-ray burst (GRB) pulse light curves. In this model an impactor wave in an expanding GRB jet accelerates from subluminal to superluminal velocities, or decelerates from superluminal to subluminal velocities. The impactor wave interacts with the surrounding medium to produce Cherenkov and/or other collisional radiation when traveling faster than the speed of light in this medium, and other mechanisms (such as thermalized Compton or synchrotron shock radiation) when traveling slower than the speed of light. These transitions create both a time-forward and a time-reversed set of light curve features through the process of Relativistic Image Doubling (RID). The model can account for a variety of unexplained yet observed GRB pulse behaviors including the amount of stretching observed in time-reversed GRB pulse residuals and the relationship between stretching factor and pulse asymmetry. The model is applicable to all GRB classes since similar pulse behaviors are observed in long/intermediate GRBs, short GRBs, and x-ray flares. The free model parameters are the impactor's Lorentz factor when moving subluminally, its Lorentz factor when moving superluminally, and the speed of light in the impacted medium.

We demonstrate how pulse structures in Short gamma-ray bursts (SGRBs), coupled with observations of GRB/GW 170817A, constrain the geometries of dying HMXRB systems composed of merging neutron stars.

The Fermi GBM catalog provides a large database with many measured variables that can be used to explore and verify gamma-ray burst classification results. We have used Principal Component Analysis and statistical clustering techniques to look for clustering in a sample of 801 gamma-ray bursts described by sixteen classification variables. The analysis recovers what appears to be the Short class and two long-duration classes that differ from one another via peak flux, with negligible variations in fluence, duration and spectral hardness. Neither class has properties entirely consistent with the Intermediate GRB class. Spectral hardness has been a critical Intermediate class property. Rather than providing spectral hardness, Fermi GBM provides a range of fitting variables for four different spectral models; it is not intuitive how these variables can be used to support or disprove previous GRB classification results.

Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing.

We demonstrate that the `smoke' of limited instrumental sensitivity smears out structure in gamma-ray burst (GRB) pulse light curves, giving each a triple-peaked appearance at moderate signal-to-noise and a simple monotonic appearance at low signal-to-noise. We minimize this effect by studying six very bright GRB pulses (signal-to-noise generally $> 100$), discovering surprisingly that each exhibits $\textit{complex time-reversible}$ wavelike residual structures. These `mirrored' wavelike structures can have large amplitudes, occur on short timescales, begin/end long before/after the onset of the monotonic pulse component, and have pulse spectra that generally evolve hard to soft, re-hardening at the time of each structural peak. Among other insights, these observations help explain the existence of negative pulse spectral lags, and allow us to conclude that GRB pulses are less common, more complex, and have longer durations than previously thought. Because structured emission mechanisms that can operate forwards and backwards in time seem unlikely, we look to $\textit{kinematic}$ behaviors to explain the time-reversed light curve structures. We conclude that each GRB pulse involves a single impactor interacting with an independent medium. Either the material is distributed in a bilaterally symmetric fashion, the impactor is structured in a bilaterally symmetric fashion, or the impactor's motion is reversed such that it returns along its original path of motion. The wavelike structure of the time-reversible component suggests that radiation is being both produced and absorbed/deflected dramatically, repeatedly, and abruptly from the monotonic component.

GRB170817A, associated with the LIGO-Virgo GW170817 neutron-star merger event, lacks the short duration and hard spectrum of a Short gamma-ray burst (GRB) expected from long-standing classification models. Correctly identifying the class to which this burst belongs requires comparison with other GRBs detected by the Fermi GBM. The aim of our analysis is to classify Fermi GRBs and to test whether or not GRB170817A belongs - as suggested - to the Short GRB class. The Fermi GBM catalog provides a large database with many measured variables that can be used to explore gamma-ray burst classification. We use statistical techniques to look for clustering in a sample of 1298 gamma-ray bursts described by duration and spectral hardness. Classification of the detected bursts shows that GRB170817A most likely belongs to the Intermediate, rather than the Short GRB class. We discuss this result in light of theoretical neutron-star merger models and existing GRB classification schemes. It appears that GRB classification schemes may not yet be linked to appropriate theoretical models, and that theoretical models may not yet adequately account for known GRB class properties. We conclude that GRB170817A may not fit into a simple phenomenological classification scheme.

We analyze pulse properties of Short gamma-ray bursts (GRBs) from a new catalog containing 434 pulses from 387 BATSE Time-Tagged Event (TTE) GRBs. Short GRB pulses exhibit correlated properties of duration, fluence, hardness, and amplitude, and they evolve hard-to-soft while undergoing similar triple-peaked light curves similar to those found in Long/Intermediate bursts. We classify pulse light curves using their temporal complexities, demonstrating that Short GRB pulses exhibit a range of complexities from smooth to highly variable. Most of the bright, hard, chaotic emission seen in complex pulses seems to represent a separate highly-variable emission component. Unlike Long/Intermediate bursts, as many as 90\% of Short GRBs are single-pulsed. However, emission in Short multi-pulsed bursts is coupled such that the first pulse's duration is a predictor of both the interpulse separation and subsequent pulse durations. These results strongly support the idea that external shocks produce the prompt emission seen in Short GRBs. The similarities between the triple-peaked structures and spectral evolution of Long, Short, and Intermediate GRBs then suggests that external shocks are responsible for the prompt emission observed in all GRB classes. In addition to these findings, we identify a new type of gamma-ray transient in which peak amplitudes occur at the end of the burst rather than at earlier times. Some of these "Crescendo" bursts are preceded by rapid-fire "Staccato" pulses, whereas the remaining are preceded by a variable episode that could be unresolved staccato pulses.

A characteristic energy of individual gamma-ray burst (GRB) spectra can in most cases be determined from the peak energy of the energy density spectra ($\nu{\cal F}_{\nu}$), called '$E_{\rm peak}$'. Distributions of $E_{\rm peak}$ have been compiled for time-resolved spectra from bright GRBs, and also time-averaged spectra and peak flux spectra for nearly every burst observed by CGRO-BATSE and Fermi-GBM. Even when determined by an instrument with a broad energy band, such as GBM (8 keV to 40 MeV), the distributions themselves peak at around 240 keV in the observer's frame, with a spread of roughly a decade in energy. $E_{\rm peak}$ can have considerable evolution (sometimes greater than one decade) within any given burst, as amply demonstrated by single pulses in GRB110721A and GRB130427A. Meanwhile, several luminosity or energy relations have been proposed to correlate with either the time-integrated or peak flux $E_{\rm peak}$. Thus, when discussing correlations with $E_{\rm peak}$, the question arises, "Which $E_{\rm peak}$?". A single burst may be characterized by any one of a number of values for $E_{\rm peak}$ that are associated with it. Using a single pulse simulation model with spectral evolution as a proxy for the type of spectral evolution observed in many bursts, we investigate how the time-averaged $E_{\rm peak}$ emerges from the spectral evolution within a single pulse, how this average naturally correlates with the peak flux derived $E_{\rm peak}$ in a burst and how the distribution in $E_{\rm peak}$ values from many bursts derives its surprisingly narrow width.

Isolated Swift gamma-ray burst (GRB) pulses, like their higher-energy BATSE counterparts, emit the bulk of their pulsed emission as a hard-to-soft component that can be fitted by the Norris et al. (2005) empirical pulse model. This signal is overlaid by a fainter, three-peaked signal that can be modeled by an empirical wave-like function (Hakkila and Preece, 2014): the two fits combine to reproduce GRB pulses with distinctive three-peaked shapes. The precursor peak appears on or before the pulse rise and is often the hardest component, the central peak is the brightest, and the decay peak converts exponentially decaying emission into a long, soft, power-law tail. Accounting for systematic instrumental differences, the general characteristics of the fitted pulses are remarkably similar. Isolated GRB pulses are dominated by hard-to-soft evolution; this is more pronounced for asymmetric pulses than for symmetric ones. Isolated GRB pulses can also exhibit intensity tracking behaviors that, when observed, are tied to the timing of the three peaks: pulses with the largest maximum hardnesses are hardest during the precursor, those with smaller maximum hardnesses are hardest during the central peak, and all pulses can re-harden during the central peak and/or during the decay peak. Since these behaviors are essentially seen in all isolated pulses, the distinction between "hard-to-soft" and "intensity-tracking" pulses really no longer applies. Additionally, the triple-peaked nature of isolated GRB pulses seems to indicate that energy is injected on three separate occasions during the pulse duration: theoretical pulse models need to account for this.

Several large structures, including the Sloan Great Wall, the Huge Large Quasar Group, and a large gamma-ray burst cluster referred to as the Hercules-Corona Borealis Great Wall, appear to exceed the maximum structural size predicted by Universal inflationary models. The existence of very large structures such as these might necessitate cosmological model modifications. Gamma-ray bursts are the most luminous sources found in nature. They are associated with the stellar endpoints of massive stars and are found in and near distant galaxies. Since they are viable indicators of the dense part of the Universe containing normal matter, the spatial distribution of gamma-ray bursts can serve as tracers of Universal large-scale structure. An increased sample size of gamma-ray bursts with known redshift provides us with the opportunity to validate or invalidate the existence of the Hercules-Corona Borealis Great Wall. Nearest-neighbour tests are used to search the larger sample for evidence of clustering and a bootstrap point-radius method is used to estimate the angular cluster size. The potential influence of angular sampling biasing is studied to determine the viability of the results. The larger gamma-ray burst database further supports the existence of a statistically significant gamma-ray burst cluster at 1.6 < z < 2.1 with an estimated angular size of 2000-3000 Mpc. Although small number statistics limit our angular resolution and do not rule out the existence of adjacent and/or line-of-sight smaller structures, these structures must still clump together in order for us to see the large gamma-ray burst cluster detected here. This cluster provides support for the existence of very large-scale universal heterogeneities.

Swift's remarkable ability to quickly localize gamma-ray bursts has led to the accumulation of a sizable burst sample for which both angular locations and redshifts are measured. This sample has become large enough that it can potentially be used to probe angular anisotropies indicative of large-scale universal structure. In a previous work, a large clustering of gamma-ray bursts at redshift z about 2 was reported in the general direction of the constellations of Hercules and Corona Borealis. Since that report, a 42 per cent increase in the number of z about 2 gamma-ray bursts has been observed, warranting an updated analysis. Surprisingly, the cluster is more pronounced now than it was when it was first reported.

According to the cosmological principle, Universal large-scale structure is homogeneous and isotropic. The observable Universe, however, shows complex structures even on very large scales. The recent discoveries of structures significantly exceeding the transition scale of 370 Mpc pose a challenge to the cosmological principle. We report here the discovery of the largest regular formation in the observable Universe; a ring with a diameter of 1720 Mpc, displayed by 9 gamma ray bursts (GRBs), exceeding by a factor of five the transition scale to the homogeneous and isotropic distribution. The ring has a major diameter of $43^o$ and a minor diameter of $30^o$ at a distance of 2770 Mpc in the 0.78<z<0.86 redshift range, with a probability of $2\times 10^{-6}$ of being the result of a random fluctuation in the GRB count rate. Evidence suggests that this feature is the projection of a shell onto the plane of the sky. Voids and string-like formations are common outcomes of large-scale structure. However, these structures have maximum sizes of 150 Mpc, which are an order of magnitude smaller than the observed GRB ring diameter. Evidence in support of the shell interpretation requires that temporal information of the transient GRBs be included in the analysis. This ring-shaped feature is large enough to contradict the cosmological principle. The physical mechanism responsible for causing it is unknown.

The Fermi GBM Catalog has been recently published. Previous classification analyses of the BATSE, RHESSI, BeppoSAX, and Swift databases found three types of gamma-ray bursts. Now we analyzed the GBM catalog to classify the GRBs. PCA and Multiclustering analysis revealed three groups. Validation of these groups, in terms of the observed variables, shows that one of the groups coincides with the short GRBs. The other two groups split the long class into a bright and dim part, as defined by the peak flux. Additional analysis is needed to determine whether this splitting is only a mathematical byproduct of the analysis or has some real physical meaning.

A study of a set of well-isolated pulses in Long and Intermediate gamma ray burst light curves indicates that simple pulse models having smooth and monotonic pulse rise and decay regions are inadequate. Examining the residuals of fits of pulses to such models suggests the following patterns of departure from smooth pulses: three separate wavelike peaks found in the residuals of each pulse (the precursor peak, the central peak, and the decay peak) combine with the underlying Norris et al. (2005) pulse model to produce five distinct regions in the temporal evolution of each pulse. The Precursor Shelf occurs prior to or concurrent with the exponential Rapid Rise. The pulse reaches maximum intensity at the Peak Plateau, then undergoes a Rapid Decay. The decay gradually slows into an Extended Tail. Despite these distinct temporal segments, the pulses studied are almost universally characterized by hard to soft spectral evolution, arguing that the new pulse features reflect a single evolution, rather than being artifacts of pulse overlap. The fluctuations can give a single pulse the appearance of having up to three distinct localized peaks, leading to ambiguities in pulse-fitting if an incorrect pulse model is used. The approach demonstrates that complex GRBs may be composed of fewer pulses than indicated by the number of peaks. The large degree of similar spectro-temporal behavior within gamma-ray burst pulses indicates that a single process is responsible for producing pulses spanning a tremendous range of durations, luminosities, and spectral hardnesses, and the correlated characteristics of the wavelike peaks are related to the pulse asymmetry, suggesting kinematic origins that seem supportive of relativistic shocks.

Context. Research over the past three decades has revolutionized cosmology while supporting the standard cosmological model. However, the cosmological principle of Universal homogeneity and isotropy has always been in question, since structures as large as the survey size have always been found each time the survey size has increased. Until 2013, the largest known structure in our Universe was the Sloan Great Wall, which is more than 400 Mpc long located approximately one billion light years away. Aims. Gamma-ray bursts are the most energetic explosions in the Universe. As they are associated with the stellar endpoints of massive stars and are found in and near distant galaxies, they are viable indicators of the dense part of the Universe containing normal matter. The spatial distribution of gamma-ray bursts can thus help expose the large scale structure of the Universe. Methods. As of July 2012, 283 GRB redshifts have been measured. Subdividing this sample into nine radial parts, each containing 31 GRBs, indicates that the GRB sample having 1.6 < z < 2.1 differs significantly from the others in that 14 of the 31 GRBs are concentrated in roughly 1/8 of the sky. A two-dimensional Kolmogorov-Smirnov test, a nearest-neighbour test, and a Bootstrap Point-Radius Method explore the significance of this clustering. Results. All tests used indicate that there is a statistically significant clustering of the GRB sample at 1.6 < z < 2.1. Furthermore, this angular excess cannot be entirely attributed to known selection biases, making its existence due to chance unlikely. Conclusions. This huge structure lies ten times farther away than the Sloan Great Wall, at a distance of approximately ten billion light years. The size of the structure defined by these GRBs is about 2000-3000 Mpc, or more than six times the size of the Sloan Great Wall.

Research over the past three decades has revolutionized the field of cosmology while supporting the standard cosmological model. However, the cosmological principle of Universal homogeneity and isotropy has always been in question, since structures as large as the survey size have always been found as the survey size has increased. Until now, the largest known structure in our Universe is the Sloan Great Wall (SGW), which is more than 400 Mpc long and located approximately one billion light-years away. Here we report the discovery of a structure at least six times larger than the Sloan Great Wall that is suggested by the distribution of gamma-ray bursts (GRBs). Gamma-ray bursts are the most energetic explosions in the Universe. They are associated with the stellar endpoints of massive stars and are found in and near distant galaxies. Therefore, they are very good indicators of the dense part of the Universe containing normal matter. As of July 2012, 283 GRB redshifts have been measured. If one subdivides this GRB sample into nine radial parts and compares the sky distributions of these subsamples (each containing 31 GRBs), one can observe that the fourth subsample (1.6 < z < 2.1) differs significantly from the others in that many of the GRBs are concentrated in the same angular area of the sky. Using the two-dimensional Kolmogorov-Smirnov test, the significance of this observation is found to be less than 0.05 per cent. Fourteen out of the 31 Gamma-Ray Bursts in this redshift band are concentrated in approximately 1/8 of the sky. The binomial probability to find such a deviation is p=0.0000055. This huge structure lies ten times farther away than the Sloan Great Wall, at a distance of approximately ten billion light-years. The size of the structure defined by these GRBs is about 2000-3000 Mpc, or more than six times the size of the largest known object (SGW) in the Universe.

Curvature effects in Gamma-ray bursts (GRBs) have long been a source of considerable interest. In a collimated relativistic GRB jet, photons that are off-axis relative to the observer arrive at later times than on-axis photons and are also expected to be spectrally softer. In this work, we invoke a relatively simple kinematic two-shell collision model for a uniform jet profile and compare its predictions to GRB prompt-emission data for observations that have been attributed to curvature effects such as the peak-flux--peak-frequency relation, i.e., the relation between the $\nu$F$_\nu$ flux and the spectral peak, E$_{pk}$ in the decay phase of a GRB pulse, and spectral lags. In addition, we explore the behavior of pulse widths with energy. We present the case of the single-pulse Fermi GRB 110920, as a test for the predictions of the model against observations.

A preliminary study of a set of well-isolated pulses in GRB light curves indicates that simple pulse models, with smooth and monotonic pulse rise and decay regions, are inadequate. Examining the residuals of fits of pulses to such models suggests the following patterns of departure from the smooth pulse model of Norris et al. (2005): A Precursor Shelf occurs prior to or concurrent with the exponential Rapid Rise. The pulse reaches maximum intensity at the Peak Plateau, then undergoes a Rapid Decay. The decay changes into an Extended Tail. Pulses are almost universally characterized by hard-to-soft evolution, arguing that the new pulse features reflect a single physical phenomenon, rather than artifacts of pulse overlap.

We demonstrate that distinguishable gamma-ray burst pulses exhibit similar behaviors as evidenced by correlations among the observable pulse properties of duration, peak luminosity, fluence, spectral hardness, energy-dependent lag, and asymmetry. Long and Short burst pulses exhibit these behaviors, suggesting that a similar process is responsible for producing all GRB pulses. That these properties correlate in the observer's frame indicates that intrinsic correlations are strong enough to not be diluted into insignificance by the dispersion in distances and redshift. We show how all correlated pulse characteristics can be explained by hard-to-soft pulse evolution, and we demonstrate that "intensity tracking" pulses not having these properties are not single pulses; they instead appear to be composed of two or more overlapping hard-to-soft pulses. In order to better understand pulse characteristics, we recognize that hard-to-soft evolution provides a more accurate definition of a pulse than its intensity variation. This realization, coupled with the observation that pulses begin near-simultaneously across a wide range of energies, leads us to conclude that the observed pulse emission represents the energy decay resulting from an initial injection, and that one simple and as yet unspecified physical mechanism is likely to be responsible for all gamma-ray burst pulses regardless of the environment in which they form and, if GRBs originate from different progenitors, then of the progenitors that supply them with energy.

We test the hypothesis that prompt gamma-ray burst pulse emission starts simultaneously at all energies (the Pulse Start Conjecture). Our analysis, using a sample of BATSE bursts observed with four channel, 64-ms data and performed using a pulse fit model, generally supports this hypothesis for the Long GRB class, although a few discrepant pulses belong to bursts observed during times characterized by low signal-to-noise, hidden pulses, and/or significant pulse overlap. The typical uncertainty in making this statement is < 0.4 s for pulses in Long GRBs (and < 0.2 s for 40% of the pulses) and perhaps < 0.1 s for pulses in Short GRBs. When considered along with the Epk decline found in GRB pulse evolution, this result implies that energy is injected at the beginning of each and every GRB pulse, and the subsequent spectral evolution, including the pulse peak intensity, represents radiated energy losses from this initial injection.

Continued study of the BATSE catalog verifies previously-identified correlations between pulse lag and pulse duration and corresponding anti-correlations between both properties and pulse peak flux for a large sample of Long GRB pulses; the study also finds correlations between pulse peak lags, pulse asymmetry, and pulse hardness. These correlations apparently can be used to delineate Long GRBs from Short ones. Correlated pulse properties represent constraints that can be used to guide theoretical modeling, whereas bulk prompt emission properties appear to be constructed by combining and smearing out pulse characteristics in ways that potentially lose valuable information.

Correlations among pulse properties in the prompt emission of long GRBs can potentially be used as cosmological distance indicators to estimate redshifts of GRBs to which these pulses belong. We demonstrate application of this technique to a sample of GRBs for which redshifts are not known. We also study the scatter of predicted redshifts of pulses found within individual bursts. We explore the characteristics of this scatter in hopes of identifying systematic corrections and/or pulse subsets that can be used to increase the technique's reliability.

We derive a new peak lag vs. peak luminosity relation in gamma-ray burst (GRB) pulses. We demonstrate conclusively that GRB spectral lags are pulse rather than burst properties and show how the lag vs. luminosity relation determined from CCF measurements of burst properties is essentially just a rough measure of this newly derived relation for individual pulses. We further show that most GRB pulses have correlated properties: short-lag pulses have shorter durations, are more luminous, and are harder within a burst than long-lag pulses. We also uncover a new pulse duration vs. pulse peak luminosity relation, and indicate that long-lag pulses often precede short-lag pulses. Although most pulse behaviors are supportive of internal shocks (including long-lag pulses), we identify some pulse shapes that could result from external shocks.

A simple scheme delineates Long GRBs with similar time histories using the Internal Luminosity Function power-law index and the spectral lag. Several generalizations are made about time history morphologies: 1) GRBs with long spectral lags contain fewer pulses that are broader than those found in bursts with short spectral lags, 2) short-lag bursts with small ILF power-law indices have many narrow pulses and are highly variable, while long-lag bursts with small ILF power-law indices are characterized by broad, smooth pulses and have low variability. GRB time history morphologies primarily identify intrinsic rather than extrinsic characteristics based on correlations with gamma-ray luminosity, afterglow luminosity, and numbers of pulses. These characteristics result because internal relativistic effects due to bulk Lorentz factor are larger than cosmological effects, and because the numbers and shapes of pulses indicate different efficiencies and forms of GRB energy release. Single-pulsed GRBs are characterized either by large ILF power indices (indicating a range of jet opening angles and Lorentz factors, with a FRED pulse shape), or they have long lags (large jet opening angles with low Lorentz factors, with either a FRED pulse shape or an unpeaked, smooth pulse shape). They also have lower-luminosity afterglows than multi-pulsed GRBs. GRBs with simple time histories are often associated with Type Ibc supernovae. This suggests that some single-pulsed GRBs contain single, beamed blast waves that are similar to and have characteristics that overlap with those of many supernovae. Such a connection may not exist between multi-pulsed GRBs and supernovae.

Currently, the best available probe of the early phase of gamma-ray burst (GRB) jet attributes is the prompt gamma-ray emission, in which several intrinsic and extrinsic variables determine GRB pulse evolution. Bright, usually complex bursts have many narrow pulses that are difficult to model due to overlap. However, the relatively simple, long spectral lag, wide-pulse bursts may have simpler physics and are easier to model. In this work we analyze the temporal and spectral behavior of wide pulses in 24 long-lag bursts, using a pulse model with two shape parameters -- width and asymmetry -- and the Band spectral model with three shape parameters. We find that pulses in long-lag bursts are distinguished both temporally and spectrally from those in bright bursts: the pulses in long spectral lag bursts are few in number, and ~ 100 times wider (10s of seconds), have systematically lower peaks in nu*F(nu), harder low-energy spectra and softer high-energy spectra. We find that these five pulse descriptors are essentially uncorrelated for our long-lag sample, suggesting that at least five parameters are needed to model burst temporal and spectral behavior. However, pulse width is strongly correlated with spectral lag; hence these two parameters may be viewed as mutual surrogates. We infer that accurate formulations for estimating GRB luminosity and total energy will depend on several gamma-ray attributes, at least for long-lag bursts. The prevalence of long-lag bursts near the BATSE trigger threshold, their predominantly low nF(n) spectral peaks, and relatively steep upper power-law spectral indices indicate that Swift will detect many such bursts.

We have identified two quiescent GRBs (bursts having two or more widely-separated distinct emission episodes) in which the post-quiescent emission exhibits distinctly different characteristics than the pre-quiescent emission. In these two cases (BATSE GRBs 960530 and 980125), the second emission episode has a longer lag, a smoother morphology, and softer spectral evolution than the first episode. Although the pre-quiescent emission satisfies the standard internal shock paradigm, we demonstrate that the post-quiescent emission is more consistent with external shocks. We infer that some observed soft, faint, long-lag GRBs are external shocks in which the internal shock signature is not observed. We further note that the peak luminosity ratio between quiescent episodes is not in agreement with the ratio predicated by the lag vs. peak luminosity relationship. We briefly discuss these observations in terms of current collapsar jet models.

Unsupervised pattern recognition algorithms support the existence of three gamma-ray burst classes; Class I (long, large fluence bursts of intermediate spectral hardness), Class II (short, small fluence, hard bursts), and Class III (soft bursts of intermediate durations and fluences). The algorithms surprisingly assign larger membership to Class III than to either of the other two classes. A known systematic bias has been previously used to explain the existence of Class III in terms of Class I; this bias allows the fluences and durations of some bursts to be underestimated (Hakkila et al., ApJ 538, 165, 2000). We show that this bias primarily affects only the longest bursts and cannot explain the bulk of the Class III properties. We resolve the question of Class III existence by demonstrating how samples obtained using standard trigger mechanisms fail to preserve the duration characteristics of small peak flux bursts. Sample incompleteness is thus primarily responsible for the existence of Class III. In order to avoid this incompleteness, we show how a new dual timescale peak flux can be defined in terms of peak flux and fluence. The dual timescale peak flux preserves the duration distribution of faint bursts and correlates better with spectral hardness (and presumably redshift) than either peak flux or fluence. The techniques presented here are generic and have applicability to the studies of other transient events. The results also indicate that pattern recognition algorithms are sensitive to sample completeness; this can influence the study of large astronomical databases such as those found in a Virtual Observatory.

Gamma-ray bursts provide what is probably one of the messiest of all astrophysical data sets. Burst class properties are indistinct, as overlapping characteristics of individual bursts are convolved with effects of instrumental and sampling biases. Despite these complexities, data mining techniques have allowed new insights to be made about gamma-ray burst data. We demonstrate how data mining techniques have simultaneously allowed us to learn about gamma-ray burst detectors and data collection, cosmological effects in burst data, and properties of burst subclasses. We discuss the exciting future of this field, and the web-based tool we are developing (with support from the NASA AISR Program). We invite others to join us in AI-guided gamma-ray burst classification (http://grb.mnsu.edu/grb/).

We describe the design of a suite of software tools to allow users to query Gamma-Ray Burst (GRB) data and perform data mining expeditions. We call this suite of tools a shed (SHell for Expeditions using Datamining). Our schedule is to have a completed prototype (funded via the NASA AISRP) by February, 2002. Meanwhile, interested users will find a partially functioning tool shed at http://grb.mankato.msus.edu.

We use ESX, a product of Information Acumen Corporation, to perform unsupervised learning on a data set containing 797 gamma-ray bursts taken from the BATSE 3B catalog. Assuming all attributes to be distributed logNormally, Mukherjee et al. (1998) analyzed these same data using a statistical cluster analysis. Utilizing the logarithmic values for T90 duration, total fluence, and hardness ratio HR321 their results showed the instances formed three classes. Class I contained long/bright/intermediate bursts, class II consisted of short/faint/hard bursts and class III was represented by intermediate/intermediate/soft bursts. When ESX was presented with these data and restricted to forming a small number of classes, the two classes found by previous standard techniques were determined. However, when ESX was allowed to form more than two classes, four classes were created. One of the four classes contained a majority of short bursts, a second class consisted of mostly intermediate bursts, and the final two classes were subsets of the Class I (long) bursts determined by Mukherjee et al. We hypothesize that systematic biases may be responsible for this variation.

The three gamma-ray burst (GRB) classes identified by statistical clustering analysis (Mukherjee et al. 1998) are examined using the pattern recognition algorithm C4.5 (Quinlan 1986). Although the statistical existence of Class 3 (intermediate duration, intermediate fluence, soft) is supported, the properties of this class do not need to arise from a distinct source population. Class 3 properties can easily be produced from Class 1 (long, high fluence, intermediate hardness) by a combination of measurement error, hardness/intensity correlation, and a newly-identified BATSE bias (the fluence duration bias). Class 2 (short, low fluence, hard) does not appear to be related to Class 1.

The fluence duration bias causes fluences and durations of faint gamma-ray bursts to be systematically underestimated relative to their peak fluxes. Using Monte Carlo analysis, we demonstrate how this effect explains characteristics of structure of the fluence vs. 1024 ms peak flux diagram. Evidence of this bias exists in the BATSE fluence duration database, and provides a partial explanation for the existence of burst class properties.

A possible relationship between the presence of a radio afterglow and gamma-ray burst spectral hardness is discussed. The correlation is marginally significant; the spectral hardness of the bursts with radio afterglows apparently results from a combination of the break energy Ebreak and the high-energy spectral index beta. If valid, this relationship would indicate that the afterglow does carry information pertaining to the GRB central engine.

We introduce a definition of gamma-ray burst (GRB) duty cycle that describes the GRB's efficiency as an emitter; it is the GRB's average flux relative to the peak flux. This GRB duty cycle is easily described in terms of measured BATSE parameters; it is essentially fluence divided by the quantity peak flux times duration. Since fluence and duration are two of the three defining characteristics of the GRB classes identified by statistical clustering techniques (the other is spectral hardness), duty cycle is a potentially valuable probe for studying properties of these classes.

The Burst and Transient Source Experiment (BATSE) on the Compton Gamma Ray Observatory (CGRO) has triggered on 1637 cosmic gamma-ray bursts between 1991 April 19 and 1996 August 29. These events constitute the Fourth BATSE burst catalog. The current version (4Br) has been revised from the version first circulated on CD-ROM in September 1997 (4B) to include improved locations for a subset of bursts that have been reprocssed using additional data. A significant difference from previous BATSE catalogs is the inclusion of bursts from periods when the trigger energy range differed from the nominal 50-300 keV. We present tables of the burst occurrence times, locations, peak fluxes, fluences, and durations. In general, results from previous BATSE catalogs are confirmed here with greater statistical significance.

Angular sky exposure is presented for a number of published BATSE gamma-ray burst catalogs. A new algorithm was required due to telemetry gaps resulting from onboard tape recorder failures; the new algorithm improves the 1B Catalog exposure calculation. The most influential effects limiting BATSE's exposure are (1) deadtime due to triggering, (2) sky blockage by the Earth, and (3) trigger disabling when the spacecraft is in the SAA and over other specific Earth locations. Exposure has improved during the CGRO mission as a result of decreased Solar flares and magnetospheric particle events.

The Two-Point Angular Correlation Function is a standard analysis tool used to study angular anisotropies. Since BATSE's sky exposure (the angular sampling of gamma-ray bursts) is anisotropic, the TPACF should at some point identify anisotropies in BATSE burst catalogs due to sky exposure. The effects of BATSE sky exposure are thus explored here for BATSE 3B and 4B catalogs. Sky-exposure effects are found to be small.

Revised upper limits on gamma-ray burst repetition rates are found using the BATSE 3B and 4B catalogs. A statistical repetition model is assumed in which sources burst at a mean rate but in which BATSE observes bursts randomly from each source.

Artificial intelligence (AI) classifiers can be used to classify unknowns, refine existing classification parameters, and identify/screen out ineffectual parameters. We present an AI methodology for classifying new gamma-ray bursts, along with some preliminary results.

If gamma-ray bursts originate in the Galaxy, at some level there should be a galactic pattern in their distribution on the sky. We test published galactic models by comparing their dipole and quadrupole moments with the moments of the BATSE 3B catalog. While many models have moments that are too large, several models are in acceptable or good agreement with the data.

We tighten previous upper limits on gamma ray burst repetition by analyzing the angular power spectrum of the BATSE 3B catalog of 1122 bursts. At 95% confidence, we find that no more than 2% of all observed bursts can be labeled as repeaters, even if no sources are observed to repeat more than once. If a fraction f of all observed bursts can be labeled as repeaters that are observed to burst v times each, then all models with (v-1)f>0.05 are ruled out at 99% confidence, as compared to the best previous 99% limit (v-1)f>0.27. At 95% confidence, our new limit is (v-1)f>0.02. Thus even a cluster of 6 events from a single source would have caused excess power above that present in the 3B catalog. We conclude that the current BATSE data are consistent with no repetition of classical gamma ray bursts, and that any repeater model is severely constrained by the near perfect isotropy of their angular distribution.

We use dipole and quadrupole statistics to test the large-scale isotropy of the first 1005 gamma-ray bursts observed by the Burst and Transient Source Experiment (BATSE). In addition to the entire sample of 1005 gamma-ray bursts, many subsets are examined. We use a variety of dipole and quadrupole statistics to search for Galactic and other predicted anisotropies and for anisotropies in a coordinate-system independent manner. We find the gamma-ray burst locations to be consistent with isotropy, e.g., for the total sample the observed Galactic dipole moment <cos theta> differs from the value predicted for isotropy by 0.9 sigma and the observed Galactic quadrupole moment <sin^2 b - 1/3> by 0.3 sigma. We estimate for various models the anisotropies that could have been detected. If one-half of the locations were within 86 degrees of the Galactic center, or within 28 degrees of the Galactic plane, the ensuing dipole or quadrupole moment would have typically been detected at the 99% confidence level. We compare the observations with the dipole and quadrupole moments of various Galactic models. Several Galactic gamma-ray bursts models have moments within 2 sigma of the observations; most of the Galactic models proposed to date are no longer in acceptable agreement with the data. Although a spherical Dark Matter Halo distribution could be consistent with the data, the required core radius is larger than the core radius of the Dark Matter Halo used to explain the Galaxy's rotation curve. Gamma-ray bursts are much more isotropic than any observed Galactic population, strongly favoring but not requiring an origin at cosmological distances.