Engineering PapersSearch

Engineering topics

Hakkila, Jon

Publications and source records attributed to Hakkila, Jon.

At least 19 records

An Artificial Intelligence Classification Tool and Its Application to Gamma-Ray Bursts

Despite being the most energetic phenomenon in the known universe, the astrophysics of gamma-ray bursts (GRBs) has still proven difficult to understand. It has only been within the past five years that the GRB distance scale has been firmly established, on the basis of a few dozen bursts with x-ray, optical, and radio afterglows. The afterglows indicate source redshifts of z=1 to z=5, total energy outputs of roughly 10(exp 52) ergs, and energy confined to the far x-ray to near gamma-ray regime of the electromagnetic spectrum. The multi-wavelength afterglow observations have thus far provided more insight on the nature of the GRB mechanism than the GRB observations; far more papers have been written about the few observed gamma-ray burst afterglows in the past few years than about the thousands of detected gamma-ray bursts. One reason the GRB central engine is still so poorly understood is that GRBs have complex, overlapping characteristics that do not appear to be produced by one homogeneous process. At least two subclasses have been found on the basis of duration, spectral hardness, and fluence (time integrated flux); Class 1 bursts are softer, longer, and brighter than Class 2 bursts (with two second durations indicating a rough division). A third GRB subclass, overlapping the other two, has been identified using statistical clustering techniques; Class 3 bursts are intermediate between Class 1 and Class 2 bursts in brightness and duration, but are softer than Class 1 bursts. We are developing a tool to aid scientists in the study of GRB properties. In the process of developing this tool, we are building a large gamma-ray burst classification database. We are also scientifically analyzing some GRB data as we develop the tool. Tool development thus proceeds in tandem with the dataset for which it is being designed. The tool invokes a modified KDD (Knowledge Discovery in Databases) process, which is described as follows.

Hakkila, Jon

How Sample Completeness Affects Gamma-Ray Burst Classification

Unsupervised pattern-recognition algorithms support the existence of three gamma-ray burst classes: class 1 (long, large-fluence bursts of intermediate spectral hardness), class 2 (short, small-fluence, hard bursts), and class 3 (soft bursts of intermediate durations and fluences). The algorithms surprisingly assign larger membership to class 3 than to either of the other two classes. A known systematic bias has been previously used to explain the existence of class 3 in terms of class 1 ; this bias allows the fluences and durations of some bursts to be underestimated, as recently shown by Hakkila et al. We show that this bias primarily affects only the longest bursts and cannot explain the bulk of the class 3 properties. We resolve the question of class 3's 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 3. 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.

Hakkila, Jon

How Sample Completeness Affects Gamma-Ray Burst Classification

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. 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 reserves the duration distribution of faint bursts and correlates either 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.

Hakkila, Jon

Gamma-Ray Burst Class Properties

Guided by the supervised pattern recognition algorithm C4.5 developed by Quinlan in 1986, we examine the three gamma-ray burst classes identified by Mukherjee et al. in 1998. C4.5 provides strong statistical support for this classification. However, with C4.5 and our knowledge of the Burst and Transient Source Experiment (BATSE) instrument, we demonstrate that class 3 (intermediate fluence, intermediate duration, soft) does not have to be a distinct source population: statistical/systematic errors in measuring burst attributes combined with the well-known hardness/intensity correlation can cause low peak flux class 1 (high fluence, long, intermediate hardness) bursts to take on class 3 characteristics naturally. Based on our hypothesis that the third class is not a distinct one, we provide rules so that future events can be placed in either class 1 or class 2 (low fluence, short, hard). We find that the two classes are relatively distinct on the basis of Band's work in 1993 on spectral parameters alpha, beta, and E (sub peak) alone. Although this does not indicate a better basis for classification, it does suggest that different physical conditions exist for class 1 and class 2 bursts. In the process of studying burst class characteristics, we identify a new bias affecting burst fluence and duration measurements. Using a simple model of how burst duration can be underestimated, we show how this fluence duration bias can affect BATSE measurements and demonstrate the type of effect it can have on the BATSE fluence versus peak flux diagram.

Hakkila, Jon

Gamma-Ray Burst Subclasses: Separating Instrumental Biases from Potential Source Populations

Instrumental effects mar gamma-ray burst observations. These effects can alter whether or not bursts trigger, as well as changing measured burst properties. The Burst And Transient Source Experiment (BATSE) on NASA's Compton Gamma-Ray Observatory (CGRO) has provided the largest database of gamma-ray bursts observed to date. It also provides an excellent laboratory for studying instrumental biases. We examine trigger biases and burst property biases as they pertain to previously identified gamma-ray burst classes. We also study some burst class properties in the new age of burst afterglow studies.

Hakkila, Jon

The Fourth BATSE Gamma-Ray Burst Catalog

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 1997 September (4B) to include improved locations for a subset of bursts that have been reprocessed 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.

Paciesas, William S.

Gamma-Ray Burst Class Properties

Guided by the Supervised pattern recognition algorithm C4.5, we examine the three gamma-ray burst classes identified by Mukherjee et al. C4.5 provides strong statistical support for this classification. However, with C4.5 and our knowledge of the BATSE instrument, we demonstrate that Class 3 (intermediate fluence, intermediate duration, soft) does not have to be a distinct source population: statistical/systematic errors in measuring burst attributes combined with the well-known hardness/intensity correlation can cause low peak flux Class I (high fluence, long, intermediate hardness) bursts to take on Class 3 characteristics naturally. Based on our hypothesis that the third class is not a distinct one, we provide rules so that future events can be placed in either Class I or Class 2 (low fluence, short, hard). Using classified bursts from the BATSE 4B Catalog, we plot log(N>P) vs. log(P) curves and study spectral features of each class. We find that the two classes are relatively distinct on the basis of spectral parameters, alpha, Beta, and E(sub peak) alone. Although this does not indicate a better basis for classification, it does suggest that different physical conditions exist for Class I and Class 2 bursts. In the process of studying burst class characteristics, we identify a new bias that affects measurement of burst fluences and durations. Using a simple model of how burst duration can be underestimated, we generally characterize how this fluence duration bias affects BATSE measurements, and demonstrate the type of effect it can have on the BATSE fluence vs. peak flux diagram.

Hakkila, Jon

The Fourth BATSE Gamma-Ray Burst Catalog

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

Paciesas W. S.

Internal Luminosity Distribution of Bright Gamma-Ray Bursts and its Relation to Duration and Spectral Hardness

We present first results from a comprehensive investigation into the distribution of luminosity within the 50 brightest cosmic gamma-ray bursts detected by the Burst and Transient Source Experiment (BATSE). The internal luminosity function psi(L) is defined such that the quantity psi(L)dL represents the fraction of total emission time during which the burst possesses a luminosity between L and L + dL. For these brightest bursts, the psi(L) functions are quasi-power-law-like and decrease in amplitude with increasing luminosity. Through investigation of both individual psi(L) distributions and data from the ensemble of bursts, we demonstrate a high probability for correlation between the shape of the internal luminosity function as measured by the average logarithmic slope and the burst duration as measured by the T(sub 90) parameter and, with lower significance, between the shape of psi(L) and the burst photon-fluence spectral index. We furthermore demonstrate a correlation between burst hardness ratio and duration in these brightest bursts which is opposite to that of the entire gamma-ray burst ensemble.

Horack, John M.

The Third BATSE Gamma-Ray Burst Catalog

The Burst and Transient Source Experiment (BATSE) on the Compton Gamma Ray Observatory (CGRO) has triggered on 1122 cosmic gamma-ray bursts between 1991 April 19 and 1994 September 19. These events constitute the Third BATSE (3B) burst catalog. This catalog includes the events previously reported in the 2B catalog, which covered the time interval 1991 April 19 to 1993 March 9. 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. The angular distribution is consistent with isotropy. The mean galactic dipole and quadrupole moments are within 0.6 a and 0.3 a, respectively, of the values expected for isotropy. The intensity distribution is not consistent with a homogeneous distribution of burst sources, with V/V(sub max) = 0.33 +/- 0.01. The duration distribution (T(sub 90)) exhibits bimodality, with peaks at approx. 0.5 and approx. 30 s. There is no compelling evidence for burst repetition, but only weak limits can be placed on the repetition rate.

Meegan, Charles A.

Improved Limits on Gamma-Ray Burst Repetition from BATSE

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 upsilon times each, then all models with (upsilon - 1)f greater than or equal to 0.05 are ruled out at 99% confidence, as compared to the best previous 99% limit (upsilon - 1)f greater than or equal to 0.27. At 95% confidence, our new limit is (upsilon - 1)f greater than or equal to 0.02. Thus, even a cluster of six 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.

Tegmark, Max

Improved Limits on Gamma-Ray Burst Repetition from BATSE

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 nu times each, then all models with (nu - 1)f greater than or equal to 0.05 are ruled out at 99% confidence, as compared to the best previous 99% limit (nu - 1)f greater than or equal to 0.27. At 95% confidence, our new limit is (nu - 1)f greater than or equal to 0.02. Thus, even a cluster of six 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 dis- tribution.

Tegmark, Max

The Intensity Distribution for Gamma-Ray Bursts Observed with BATSE

The intensity distributions of gamma-ray bursts observed by BATSE from 19 April 1991 to 19 September 1994 are presented. For this data set, (V/V(sub max)) is 0.329 +/- 0.011, which is 15.5 sigma away from the value of 0.5 expected for a homogeneous distribution. Standard cosmological model parameters are obtained by fitting the differentially binned peak flux distribution expressed in units of photons cm(exp -2) s(exp -1) in the energy range 50-300 keV. The value of z calculated for a peak flux of 1 photon cm(exp -2) s(exp -1) is 0.8 +/- 0.33. The procedures used to produce the peak flux data and C(sub p)/C(sub lim) data are presented. The differences between the two representations of burst intensity are emphasized so that researchers can determine which type of data is most appropriate for their studies. The sky sensitivity correction as a function of intensity for the peak flux data is also described.

Pendleton, Geoffrey N.

BATSE Observations of the Large-Scale Isotropy of Gamma-Ray Bursts

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(exp 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 deg of the Galactic center, or within 28 deg 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.

Briggs, Michael S.

Constraints on galactic distributions of gamma-ray burst sources from BATSE observations

The paradigm that gamma-ray bursts originate from Galactic sources is studied in detail using the angular and intensity distributions observed by the Burst and Transient Source Experiment (BATSE) on NASA's Compton Gamma Ray Observatory (CGRO). Monte Carlo models of gamma-ray burst spatial distributions and luminosity functions are used to simulate bursts, which are then folded through mathematical models of BATSE selection effects. The observed and computed angular intensity distributions are analyzed using modifications of standard statistical homogeneity and isotropy studies. Analysis of the BATSE angular and intensity distributions greatly constrains the origins and luminosities of burst sources. In particular, it appears that no single population of sources confined to a Galactic disk, halo, or localized spiral arm satisfactorily explains BATSE observations and that effects of the burst luminosity function are secondary when considering such models. One family of models that still satisfies BATSE observations comprises sources located in an extended spherical Galactic corona. Coronal models are limited to small ranges of burst luminosity and core radius, and the allowed parameter space for such models shrinks with each new burst BATSE observes. Multiple-population models of bursts are found to work only if (1) the primary population accounts for the general isotropy and inhomogeneity seen in the BATSE observations and (2) secondary populations either have characteristics similar to the primary population or contain numbers that are small relative to the primary population.

Hakkila, Jon

Preliminary model constraints on the spatial and luminosity distributions of GRBs observed by BATSE

The GRB angular/intensity distributions observed by BATSE greatly constrain models of the spatial and luminosity source distributions. Single populations of Galactic Disk and Galactic Halo sources appear to be untenable, and only a limited subset of models made up of sources in a symmetric Galactic Corona satisfy the observations. Comments are made on other classes of models.

Hakkila, Jon

BATSE analysis techniques for probing the GRB spatial and luminosity distributions

The Burst And Transient Source Experiment (BATSE) has measured homogeneity and isotropy parameters from an increasingly large sample of observed gamma-ray bursts (GRBs), while also maintaining a summary of the way in which the sky has been sampled. Measurement of both of these are necessary for any study of the BATSE data statistically, as they take into account the most serious observational selection effects known in the study of GRBs: beam-smearing and inhomogeneous, anisotropic sky sampling. Knowledge of these effects is important to analysis of GRB angular and intensity distributions. In addition to determining that the bursts are local, it is hoped that analysis of such distributions will allow boundaries to be placed on the true GRB spatial distribution and luminosity function. The technique for studying GRB spatial and luminosity distributions is direct. Results of BATSE analyses are compared to Monte Carlo models parameterized by a variety of spatial and luminosity characteristics.

Hakkila, Jon

BATSE software for the analysis of the gamma ray burst spatial distribution

The Burst and Transient Source Experiment (BATSE) on the Gamma Ray Observatory (GRO) is designed to study astronomical gamma ray sources and to provide better positional, spectral, and time resolution about these objects than has previously been possible from one experiment. The procedure to be used in the analysis of the gamma ray burst spatial distribution is presented. Data is input from BATSE via the Gamma Ray Burst Catalog (listing individual burst positions, flux values, and associated errors) and the Sky Sensitivity Map (which summarizes observational selection effects in table format). A FORTRAN program generates Monte Carlo burst catalogs, which are models to be compared to the actual distribution. The Monte Carlo models are then filtered through the Sky Sensitivity Map so that they suffer from the same selection effects as the actual catalog data. Additionally, each burst position is converted into a probability distribution to mimic BATSE positional sensitivity. The Burst Catalog, Monte Carlo burst catalog, and Sky Sensitivity Map are then passed onto an IDL program that compares the catalogs for statistical significance. The Sky Sensitivity Map is used to estimate how often each sky area is observed above the minimum flux level in question. Each burst found in this sky area is then weighted according to the frequency with which this sky area is observed. The catalogs are then compared via tests of homogeneity (based on their radial distributions) and isotropy (based upon their angular distributions). The results of the statistical comparisons along with graphs and charts of the summaries, are output from the IDL program for study.

Hakkila, Jon