Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: it is generally agreed that low - frequency ( radio ) emission from relativistic jets in active galactic nuclei ( agn ) is synchrotron emission from nonthermal electrons . however , it is still debated whether the high - frequency ( x - ray & @xmath0-ray ) emission is due solely to compton scattering by the electrons responsible for the low - frequency synchrotron emission or due to hadronic processes resulting from protons being co - accelerated with the electrons ( e.g. , * ? ? ? if hadronic models are valid , agn would be an expected source of neutrinos , as the proton interactions responsible for the production of neutral pions which result in @xmath0-rays also produce charged pions , which result in secondary electrons , positrons , and neutrinos ( e.g. , * ? ? ? can we expect agn to be detectable neutrino sources at the level that starburst systems of similar @xmath0-ray luminosity are expected to be ? observations by the space - based @xmath0-ray telescope _ fermi _ and future observations by the hard x - ray space telescope nustar may help to further resolve this question .    building on the approach of previous models ( e.g. * * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) , we developed a semi - analytic model for cosmic ray interactions for the starburst galaxy m82 ( * ? ? ? * hereafter yegz ) . in developing the model , we emphasized the important role of galactic winds while limiting the number of free parameters and keeping our models observationally informed . now , we aim to determine whether the model is applicable to other starburst systems and explore benefits and limitations of such models due to the presence of agn . the starburst galaxy ngc 253 has been detected in @xmath0-rays by both space - based ( _ fermi _ ) and ground - based ( hess ) telescopes and possibly contains an active central supermassive black hole ( smbh ) @xcite . seyfert 2 galaxies ngc 1068 and ngc 4945 have also been detected with _ fermi_. though ngc 1068 ( @xmath1 mpc ; nasa / ipac extragalactic database ) is more distant than the nearby ngc 253 ( @xmath2 mpc ; * ? ? ? * ) and ngc 4945 ( @xmath3 mpc ; * ? ? ? * ) , all three have @xmath0-ray fluxes of the same order of magnitude @xcite .    in this paper , we determine whether our model is applicable to a combined starburst / agn environment . we select the starburst galaxy ngc 253 with a possible low - activity agn and the high - activity seyfert galaxy ngc 1068 . applying our model to ngc 253 and ngc 1068 , we seek the best - fit cosmic ray interaction model which accurately reproduces both radio and @xmath0-ray spectra within the confines of the observed properties of the interstellar medium ( ism ) and surrounding radiation field . additionally , we make a prediction of the neutrino fluxes from both galaxies . llclc distance & 3.9 mpc & 1 & 14.4 mpc & 7 + central molecular zone ( cmz ) radius & 150 pc & 2 & 200 pc & 8 + molecular gas mass & @xmath4 @xmath5 & 3 , 4 & @xmath6 @xmath5 & 9 + average ism density & 500 - 1400 @xmath7 & & 250 @xmath7 & + ir luminosity & @xmath8 @xmath9 & 5 & @xmath10 @xmath9 & 10 + radiation field energy density & 500 - 2000 ev  @xmath7 & & @xmath11 ev  @xmath7 & + sn explosion rate & 0.1 yr@xmath12 & 6 & 0.07 yr@xmath12 & 11 + sn explosion energy & 10@xmath13 ergs & & 10@xmath13 ergs & + sn energy transferred to cr & 4 - 20% & & 10% & + ratio of primary protons & 50 & & 50 & +  to electrons ( n@xmath14/n@xmath15 ) & & & & + slope of primary cr & 2.2/2.3 & & 2.0 & +  source function & & & & +    the next section describes the observed properties of ngc 253 and ngc 1068 . section 3 details how the population of energetic particles was computed . section 4 contains the results of our model for ngc 253 and ngc 1068 . in section 5 , we compare and contrast the results for ngc 253 and ngc 1068 and discuss predictions for the neutrino flux and the resulting implications . in section 6 , we present concluding remarks . ngc 253 is a giant , barred spiral galaxy with a central starburst region . tidal interactions are the cause of many starbursts in massive galaxies , and while there is no obvious companion or trigger for the starburst in ngc 253 , kinematic studies suggest ngc 253 may have been involved in a past merger @xcite . models indicate that the current high level of nuclear star formation has been ongoing for at least 20 - 30 myr @xcite . the center of ngc 253 is host to several compact objects which have been observed over a variety of wavelengths @xcite . the strongest non - thermal , compact radio source ( th2 ) has long been considered a possible agn candidate @xcite . additionally , x - ray emission from the nucleus may be powered by an obscured agn or a combination of supernova remnants and fluorescent line emission from molecular clouds @xcite . however , new observations with nustar and _ chandra _ are unable to rule out the possibility that these compact x - ray sources are ultraluminous x - ray sources ( ulxs ) and not a low - luminosity agn @xcite .    one of the other major features of ngc 253 is the galactic wind emanating from the nuclear region ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? this galactic wind has been detected in both h@xmath16 and x - rays . observations of the position and size of the outflow at both wavelengths match up almost exactly @xcite . optical observations give an outflow velocity on the order of a few 100 km  s@xmath12 ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? additionally , @xcite suggest that , like star formation , the wind is quasi - steady .    due to the obscured nature of the starburst nucleus of ngc 253 ( e.g. , * ? ? ? * ; * ? ? ? * ) , determining the properties of the nuclear region has been difficult . radio observations yield supernova rates ranging from 0.03 - 0.3 yr@xmath12 @xcite . in this work , we adopt a rate of 0.1 yr@xmath12 based on @xcite ( see table 1 for additional parameters ) . additionally , estimates of the molecular gas mass range from @xmath17 @xmath5 ( * ? ? ? * and references therein ) . we adopt a value of @xmath18 @xmath5 for the molecular gas mass for this paper , based on work by @xcite . we adopted an ionized gas mass for ngc 253 by scaling from the ratio of ionized to molecular gas in m82 , which gives a mass of @xmath19 @xmath5 . ngc 253 is well observed in the radio for both extended emission and the central starburst core . observations by @xcite and @xcite at 333 mhz and from 1 to 7 ghz show that the radio spectrum turns down only slightly in the core at low frequencies , unlike the radio spectrum for the core of m82 , which is completely absorbed at low frequencies . high frequency observations from @xcite at 18.5 and 22 ghz and from @xcite at 94 ghz require a free - free emission contribution to the flux in addition to non - thermal synchrotron emission . observed radio and @xmath0-ray spectra are presented with model spectra in section 4 . ngc 253 has been observed in @xmath0-rays at both gev and tev energies by _ fermi _ and hess . @xcite present four data points and two upper limits over 0.2 to 200 gev from analysis of 30 months of _ fermi _ data . @xcite present a combined analysis of five data points and one upper limit over 0.2 tev to 30 tev from hess data with the _ fermi _ data and fit a single power - law to both data sets , resulting in a spectral index of @xmath20 . @xcite present a more current reduction of the _ fermi _ data which we use for comparison with our models instead of the _ fermi _ data shown in @xcite . ngc 1068 is the closest , brightest seyfert galaxy . the total bolometric luminosity of the galaxy is @xmath21 @xmath9 ; half of which is accounted for by the circumnuclear disk ( cnd ) and agn @xcite . the galaxy supports a kiloparsec - scale , non - thermal radio jet @xcite and a compact radio nucleus @xcite . the complex chemical and kinematic properties of the nuclear region have been well - studied at multiple wavelengths @xcite . the circumnuclear starburst ring surrounding the 2.3 kpc stellar bar @xcite is illuminated by star formation in the past 10 - 40 myr @xcite . the stellar population inside the cnd , the inner @xmath22200 pc , developed in two bursts at 30 and 300 myr ago @xcite . @xcite give the star formation rate ( sfr ) for the inner galaxy as @xmath22 10 @xmath5 yr@xmath12 . this gives an upper bound on the supernova rate of @xmath220.07 yr@xmath12 for the nuclear region . we convert from sfr to supernova rate with by analysis similar to @xcite . both starburst regions are rich in molecular gas with masses of @xmath23 @xmath5 for the circumnuclear ring and @xmath24 @xmath5 for the cnd @xcite . co , hcn , and hco@xmath25 observations show that the gas is moderately dense ( @xmath26 @xmath7 ) with a kinetic temperature of @xmath27 k @xcite .      as the nearest and brightest seyfert galaxy , ngc 1068 has been extensively studied in the radio . on the kiloparsec - scale , radio observations have mapped the star - forming spiral arms and the radio jets . there are extensive observations specifically targeting the agn on the parsec - scale @xcite . however , as even observations of the inner hundred parsecs have a dominant contribution from the agn @xcite , specifically the radio jet and counter - jet , there are essentially no observations of the central molecular zone ( cmz ) alone . radio maps show little evidence of ongoing star - formation in the central @xmath22150 pc . @xmath0-ray data for ngc 1068 are more limited than for ngc 253 . while there are four data points and two upper limits at gev energies , there is only a single upper limit for tev energies from hess @xcite . the _ fermi _ data range from 0.2 gev to 300 gev and the hess upper limit is at 0.2 tev . diffusion is thought to be of minor importance in star - forming environments with high gas densities and galactic winds . as such , following the approach in yegz , a cosmic ray spectrum can be calculated from a single quantity including the source function ( dependent on supernova rate ) and the lifetime . because we only consider radiative and advective losses , our model is not applicable to spatially extended regions , such as the inner spiral arms in ngc 1068 , where cosmic ray diffusion is important . thus , we apply our model specifically to the cmzs of galaxies . supernovae are the assumed drivers of cosmic ray acceleration , and so we adopt a power - law source function ( @xmath28 ) such that @xmath29 where @xmath30 is the supernova rate , @xmath31 is the volume of the starburst region , @xmath32 is the fraction of the supernova energy transferred to cosmic rays , and @xmath33 is 10@xmath13 ergs , the typical energy from a supernova explosion . thus , the spectrum for primary protons is given by @xmath34 where @xmath35 is the minimum cosmic ray energy , here taken to be @xmath36 gev . @xmath37 is the cosmic ray lifetime composed of an energy loss timescale ( due to radiative and collisional losses ) and an energy independent timescale : @xmath38 where @xmath39 is the scale height of the starburst region and @xmath40 is the speed of the particles in the wind of the starburst region ( see also * ? ? ? energy losses include ionization , pion production , bremsstrahlung , synchrotron emission , and the inverse compton effect . pion production from proton - proton collisions is the dominant energy loss for cosmic ray protons above the 1.3 gev threshold . as the lifetimes of both charged and neutral pions are very short , pions quickly decay into secondary electrons and positrons , neutrinos , and @xmath0-rays . the primary proton source function can be used to calculate the source function for charged and neutral pions @xcite , and the pion source function can then be used to calculate the spectrum of secondary electrons and positrons @xcite . one of the main assumptions of the model is that cosmic rays sample the mean density of the interstellar medium . the density of the interstellar medium sampled by the cosmic rays is vital to the model as it affects pion production and thus both the @xmath0-ray and radio spectra . the majority of the ism is hot , low density gas in which pockets of cold , very dense molecular gas and warm ionized gas can be found . in @xcite , we assumed a simple two phase model of cold , molecular gas and hot , ionized gas of negligible density . we found that for a wide range of cosmic ray injection conditions , magnetic field properties , and wind speeds , the cosmic rays sample the mean density of the ism to within a factor of two . the effect of a third phase of dense , warm , ionized gas on cosmic ray sampling will be considered in future work .    from cosmic ray proton and primary and secondary electron / positron populations , we can calculate nonthermal radio synchrotron emission and @xmath0-ray emission by neutral pion decay , bremsstrahlung , and inverse compton ( see yegz for further details ) for comparison with observed data . primary variables include spectral index ( @xmath41 ) , cosmic ray acceleration efficiency ( @xmath32 ) , molecular gas mass ( @xmath42 ) , magnetic field strength ( @xmath43 ) , wind ( advection ) speed ( @xmath40 ) , ionized gas density ( @xmath44 ) , and interstellar radiation field energy density ( @xmath45 ) .      in yegz , we found that observations of the starburst core in m82 showed a complete turn down of the radio emission at low frequencies . we attributed this to severe free - free absorption in the core and adopted a model for the starburst core with a cylindrical wall of warm , ionized gas surrounding a lower density hot medium . because some observations at low frequencies did not have enough resolution to distinguish the core from the halo in m82 , we also included a component for radio emission from the halo . observations for ngc 253 do not show the same severe free - free absorption at low frequencies as in m82 . however , there is significant free - free emission at high frequencies , as is also the case in m82 . thus , we expect that at least some fraction of the radio synchrotron emission is still absorbed by ionized gas . to account for this , we keep the thin shell geometry ( with shell width @xmath46 and radius @xmath47 ) for the nucleus , but assume only a small fraction ( @xmath48 ) of the synchrotron emission is absorbed by the ionized gas . then , the radiative intensity of the absorbed synchrotron emission is given by @xmath49 where @xmath50 is the radius of the hot , diffuse gas ( @xmath51 ) , @xmath52 is the emission coefficient for the hot gas ( @xmath53 ) , and @xmath54 is the absorption coefficient for the warm , ionized gas ( @xmath55 ) . the majority of the radio spectrum is unabsorbed synchrotron emission from the hot , diffuse gas : @xmath56 additionally , we have some portion of the radio spectrum composed of free - free emission from the ionized gas clouds : @xmath57 where @xmath58 is the emission coefficient for the warm , ionized gas ( @xmath59 ) . thus , the total emergent intensity from the starburst region is @xmath60 for an absorption fraction of @xmath61 , this reduces to @xmath62 which is the expression for the radio emission that we used for modeling m82 ( yegz ) . note that unlike in yegz , we do not include a disk / halo component for radio emission in this paper as radio observations of the disk are separable from the halo ( e.g. , * ? ? ? * ; * ? ? ? we anticipate that inverse compton scattering is a vital process for producing @xmath0-rays in ngc 253 and ngc 1068 . for inverse compton scattering , the @xmath0-ray source function is given by @xcite @xmath63 with @xcite @xmath64,\ ] ] where @xmath65 is the energy of the resulting @xmath0-ray , @xmath66 is the energy of the incident photon , @xmath0 is the energy of the electron . here , the combined cosmic ray electron / positron spectrum , @xmath67 , is in units of @xmath7 . the function @xmath68 is part of the klein - nishina cross section and is given by @xcite @xmath69 where @xmath70 for the blackbody spectrum , @xmath71 , we use an isotropic , diluted , modified blackbody spectrum ( * ? ? ? * and references therein ) @xmath72 where @xmath73 is a spatial dilution factor ( given by the normalization @xmath74 ) and @xmath75 corresponds to @xmath76 hz . as with our previous model for m82 ( yegz ) , we test a range of parameters to find the best fits to both the radio and @xmath0-ray spectra for the starburst core of ngc 253 ( see table 2 ) . we choose to vary parameters such as magnetic field strength ( @xmath77 @xmath78 g ) and radiation field energy density ( @xmath79 ev  @xmath7 ) as they are not well - constrained by observations . we vary acceleration efficiency ( @xmath80 ) as it is equivalent to varying the supernova rate and test two spectral indices ( @xmath81 ) . as the observed @xmath0-ray spectrum for ngc 253 is steeper than that of m82 , we choose larger spectral indices and do not test for @xmath82 . we also test a variety of wind ( advection ) speeds ( @xmath83 km  s@xmath12 ) in order to study the no - wind case and wind speeds similar to those derived from optical and x - ray observations . ll magnetic field strength ( @xmath43 ) & 50 - 350 @xmath78 g + wind ( advection ) speed ( @xmath40 ) & 0 - 1000 km  s@xmath12 + ionized gas density ( @xmath44 ) & 50 - 500 @xmath7 + absorption fraction ( @xmath48 ) & 0.1 - 1.0 +    in regards to the molecular gas mass , there is a wide a range of observed values and there is not yet a clear consensus on which is most accurate . varying the gas mass would have a similar ( but not identical ) effect to varying the acceleration efficiency or supernova rate . in the case of the cosmic ray electrons , varying molecular gas mass affects the losses due to ionization and bremsstrahlung but not inverse comtpon or synchrotron emission . though we adopt the value of @xmath18 @xmath5 , there is likely at least a factor of two uncertainty in this value . as such , we also test models at @xmath84 @xmath5 . while we vary the molecular gas mass , we keep the ionized gas mass constant . however , we do vary the ionized gas density ( @xmath85 @xmath7 ) , and thus the filling fraction of ionized gas , as it has a significant impact on free - free absorption and emission in the radio spectrum .    to distinguish between models , we use @xmath86 tests to compare to the observed data . we can see from the contour plots of the @xmath86 values that the @xmath0-ray and radio data constrain the value for the magnetic field strength and wind speed in different ways ( see figure 1 ) . there is a degeneracy for the magnetic field strength and wind speed for the radio data such that as the magnetic field strength increases , an increase in wind speed allows for a similarly good fit . the reverse is true for the @xmath0-rays . in the radio spectrum , an increase in magnetic field strength results in an increase in synchrotron emission such that the number of electrons must be reduced ( see yegz for further explanation ) . however , in the @xmath0-ray spectrum , an increase in the magnetic field strength leads to a reduction of the inverse compton flux as fewer electrons are available as a result of the shorter energy loss lifetimes . thus , the wind speed must be reduced to allow for more electrons to be available to produce inverse compton emission .    because the results are essentially non - complementary , we were not able to find a joint solution ( a best - fit model for both the @xmath0-ray and radio spectra based on a common set of parameters ) for the radio and @xmath0-ray data assuming @xmath18 @xmath5 for the molecular mass . this is highlighted in figure 2 , where the combined reduced @xmath86 result is dominated by the radio data . for the radio spectrum , optimal wind speeds are essentially the no wind case ( @xmath40 ) which is not supported by observations of ngc 253 . for the @xmath0-ray spectrum , optimal wind speeds range from @xmath87 km  s@xmath12 , which is slightly higher than optical determinations of the wind speeds . additionally , while @xmath0-ray models allow for a wide range of magnetic field strengths , radio models constrain the magnetic field strength to @xmath88 @xmath78 g . these field strengths correspond to @xmath222000 to @xmath223000 ev  @xmath7 . cosmic ray energy densities are nearly an order of magnitude below this , @xmath22200 to @xmath22300 ev  @xmath7 for @xmath89 . while the magnetic field energy density lies at the high end of the range of the assumed radiation field energy densities , the cosmic ray energy densities are within a factor of two of the lower end of the range . taken at face value , these model results suggest that for the system the magnetic field and radiation field are in approximate equipartition while the cosmic rays are sub - equipartition . in addition to limiting the possible magnetic field strength and wind speed values , the @xmath0-ray data severely constrain the possible combinations of spectral index and acceleration efficiency . while the radio data is fit equally well by all combinations of acceleration efficiency and spectral index , the @xmath0-ray spectrum is easily overestimated because of the large radiation fields and high densities in starburst regions . as such , the majority of the combinations of acceleration efficiency and spectral index do not result in a minimized ( best - fit ) @xmath86 value for the @xmath0-ray spectrum . so , to keep the combination of pion decay and inverse compton emission in check , our model rules out acceleration efficiencies of @xmath90 and @xmath91 for both spectral indices . we also rule out the spectral index of @xmath92 for all acceleration efficiencies , as the best - fit models underestimate the @xmath0-ray spectrum at tev energies ( see figure 3 ) . thus , the only combination of acceleration efficiency and spectral index which both minimizes and fits the tev energy data is @xmath89 , @xmath93 . equivalently , a larger acceleration efficiency pared with a lower supernova rate would have the same result , as we constrain their product . the goodness of our @xmath86 fits also change slightly depending on the radiation field energy density . the best - fit @xmath0-ray spectrum occurs for an energy density of 2000 ev  @xmath7 for a spectral index of @xmath89 . the resulting @xmath0-ray spectrum is dominated by pion decay emission but still has a significant contribution from inverse compton emission at tev energies ( see figure 3 ) .    aside from the parameters discussed thus far , we also vary the fraction of synchrotron emission absorbed in ionized gas clouds ( see table 2 ) . observations show that at low frequencies , the radio spectrum for ngc 253 levels off , unlike the radio spectrum for m82 which turns over at low frequencies . in the case of m82 , to fit the radio spectrum properly , we required that all of the synchrotron emission from the hot , diffuse medium was absorbed in ionized gas clouds such that @xmath94 ( yegz ) . however , as the radio spectrum for ngc 253 only flattens , we varied the absorption fraction to find what absorption fraction was required to accurately fit the spectrum . we found that models with an absorption fraction of up to @xmath95 agreed with all compared radio data points ( see figure 4 ) , including models with no absorption fraction ( @xmath96 ) . the best - fit models were for an absorption fraction of @xmath97 , which was selected for use for the remainder of our models . in yegz , we found that while m82 was an electron calorimeter , it was , at best , only a partial ( @xmath2250% ) proton calorimeter . for ngc 253 , we find that it is an electron calorimeter and that it can be a proton calorimeter , depending upon the selected parameters . we based our conclusions on a comparison between energy - loss and advection timescales . in ngc 253 , we find that the average wind speeds for the best - fit @xmath0-ray models are on the order of @xmath87 km  s@xmath12 . the best - fit wind speed decreases depending on the assumed magnetic field strength . comparing timescales , we find that ngc 253 is a @xmath2250% proton calorimeter . however , the best - fits for the radio spectrum were for the no wind case . in this case , ngc 253 would be a complete proton calorimeter ( cf . while the @xmath0-ray data were readily fit with our model , the radio spectrum presented much more of a challenge . as previously stated , we were not able to find a joint solution for the @xmath0-ray and radio models . while our best - fit for the @xmath0-ray spectrum had a reduced @xmath86 value of 2.0 , the best - fit for the radio spectrum had a reduced @xmath86 value of 50 . in addition to the preferred wind speed being unphysical , the best - fit model does not match up well with the @xcite data . although we did not explicitly vary molecular gas mass in our models , we did originally test a lower gas mass . models tested with a lower gas mass of @xmath98 @xmath5 resulted in significantly better radio fits ( see figures 2 and 4 ) . our best - fit model with the lower mass has a reduced @xmath86 value of 12 . this reduction in @xmath86 with the lower mass is due to higher bremsstrahlung losses for higher masses . thus , for the lower mass , more electrons are available to produce synchrotron emission at the low frequencies instead of bremsstrahlung . the discrepancy between fits for the different gas masses is , in large part , due to the limitations of the one - zone model . while there are many computational advantages to assuming a simple single zone , there are also disadvantages . in yegz , we noted that by assuming a single ionized gas density in our model , we were unable to fit both the high - frequency radio data and the low - frequency turn down . in ngc 253 , our one zone model has reached the calorimeter limit for the higher gas mass . as such , while we are producing the same amount of secondary electrons and positrons at both masses , there are fewer electrons / positrons available to produce inverse compton and synchrotron emission due to an increase in bremsstrahlung . the cmz dust temperature is the key to modeling the @xmath0-ray observations for ngc 1068 as it determines the radiation field that inverse compton emission depends on . observations by @xcite show a blackbody spectrum with temperatures in the range of @xmath99 for the inner cnd . when assuming a radiation field from dust with @xmath100 k , the photon number is significantly decreased such that our models produce negligible inverse compton @xmath0-ray emission . however , this dust temperature is attributed to the dusty torus of the agn nucleus and likely does not dominate the larger , surrounding cmz . as such , we assume that the dust temperature in the cmz is on par with the temperatures of molecular gas in the region , @xmath22100 k , and we use this to determine the radiation field spectrum .    as with ngc 253 , we intended to test a variety of different sets of parameters with which to model ngc 1068 . however , we found ngc 1068 significantly harder to model than ngc 253 . the upper bound on the supernova rate produces a @xmath0-ray spectrum that is lower by a factor of only a few ( see figure 5 ) . however , a lower bound on the supernova rate results in a @xmath0-ray spectrum that is nearly two order of magnitudes lower than the observed data . because we were underestimating the @xmath0-ray emission , we selected parameters to maximize the inverse compton emission ( a magnetic field strength of @xmath101 @xmath78 g and a radiation field energy density of @xmath102 ev  @xmath7 ) and the pion decay emission and bremsstrahlung ( a wind speed of @xmath103 km  s@xmath12 ) . even selecting parameters to augment the @xmath0-ray emission , without invoking an extra source of cosmic rays , we were not able to produce a model which agrees with the _ fermi _ observations to better than a factor of a few .    further complicating matters is the radio spectrum for ngc 1068 . while the galaxy has been extensively observed in the radio spectrum , the presence of a radio jet greatly overshadows any emission not originating from the agn or its jets . ultimately , we chose to compare our radio models against a few different radio observations . first , we plot our radio models against radio observations of the agn core ( s1 ) found in @xcite . though the observations of the agn are on the parsec - scale , they do highlight the fact that they are of a fundamentally different nature than the emission from star - forming regions . we also plot radio observations of the region labeled s2 @xcite , which could be the counter - jet as seen through the star - forming disk . additionally , we use the 18  cm radio map in @xcite to compare with our models . while the outermost contour includes the spiral arms , the second contour should give us an upper limit on the radio continuum in the cmz excluding the component from the agn and the jets . taking the brightness level ( of 74 mjy / beam ) and multiplying by the area of the cmz , we obtain an upper limit of @xmath22 13 mjy ( see red squared in figure 5 ) . while nearly all of the energy of the cosmic ray electrons goes into inverse compton emission , some small fraction of the energy is lost due to bremsstrahlung and synchrotron emission . assuming a moderate starburst magnetic field strength of 200  @xmath78 g results in a significant amount of synchrotron emission . as the upper limit ( red square ) that we plot is nearly an order of magnitude below the model , we must be significantly overestimating the synchrotron emission in the cmz . thus , while assuming the upper limit on the supernova rate results in a @xmath0-ray spectrum which comes within a factor of a few of the observed data , the radio spectrum disagrees significantly with observations ( see figure 5 ) . previously , we developed and tested a model for cosmic ray interactions in starburst nuclei ( yegz ) . in testing our model against other nuclei and combined starburst / agn environments , we have applied our updated model to two giant , barred spiral galaxies ; specifically , we have modeled both ngc 253 , with an unambiguous starburst nucleus , and ngc 1068 , the archetypal seyfert 2 galaxy . the @xmath0-ray spectrum for ngc 253 is well fit with a combination of emission from neutral pion decay , bremsstrahlung , and inverse compton , requiring both hadronic and leptonic emission mechanisms .    while we find a joint solution for the @xmath0-ray and radio spectra for the lower molecular gas mass for the ngc 253 cmz , we were unable to do the same at a higher gas mass ( @xmath18 @xmath5 ) . @xcite find best - fit models for both the @xmath0-ray and radio spectra for similarly large densities ( @xmath104 @xmath7 ) . however , they do not include wind speed as a free parameter and fix their advection / diffusion timescale at 10 myr ( which is equivalent to a wind speed of @xmath22 5 km  s@xmath12 for a half - thickness of 50 pc ) . this would allow for a good fit to the radio spectrum at higher masses and their lower radiation field would also allow for a similarly good fit for the @xmath0-ray spectrum . other models include those of @xcite who find joint solutions when assuming parameters closer to our lower gas mass and wind speeds between 300 and 600 km  s@xmath12 . we believe that the data support our choice of gas masses in the @xmath4 @xmath5 range as well as wind speeds of @xmath105 km  s@xmath12 ( see * ? ? ? * ; * ? ? ? * ) and that our models therefore contain astrophysically reasonable parameters .    though located nearly four times further away , ngc 1068 has a @xmath0-ray flux comparable to that of ngc 253 . in comparing our model with observations , we find that the inner starburst region in ngc 1068 does not supply a sufficient cosmic ray population to produce the observed @xmath0-ray emission by itself . this agrees with analysis of the supernova rate , total gas mass , and @xmath0-ray luminosity by @xcite for ngc 1068 in comparison to other starburst galaxies . while a scenario which boosts the cosmic ray population in the cmz , such as a diffusion of cosmic rays from the circumnuclear ring into the cmz , is possible , it is most likely that the high - energy emission results from interactions in the jet . an agn jet origin of the @xmath0-ray emission in ngc 1068 agrees with calculations of @xcite who propose a leptonic external inverse compton ( eic ) model for the high - energy emission . while proton - proton interactions dominate as a cosmic ray energy loss mechanism in the ism , proton-@xmath0 interactions will be dominant in the jet due to the lack of molecular gas necessary for proton - proton interactions . proton - proton interactions result in both charged and neutral pions . while neutral pions decay quickly into @xmath0-rays , charged pions decay into muons and neutrinos followed by electrons and positrons and more neutrinos . thus , galaxies whose @xmath0-ray spectra are dominated by or are significantly contributed to by pion decay will also be sources of neutrinos . as proton-@xmath0 interactions still result in pion production , agn could also be a source for neutrinos as starbursts are , if @xmath0-ray emission in the jet is hadronic .    while a detection of neutrinos from agn would help to determine the origin of @xmath0-ray emission , such a detection will also provide vital clues for starburst galaxies . observations of neutrinos provide direct information about the nature of the @xmath0-ray spectrum and the cosmic ray proton population . better constraining the contribution of neutral pion decay to the @xmath0-ray spectrum , also better constrains the contribution of bremsstrahlung and inverse compton emission . in the case of inverse compton , this would be another diagnostic for determining the radiation field , which can be uncertain in such environments .    to make a neutrino flux prediction , we start with the source function for muon neutrinos which results from charged pion decay ( @xmath106 and @xmath107 ) depends on the source function for pions : @xmath108 where @xmath109 is the source functions for pions ( see yegz for details ) and @xmath110 @xcite . we found that the @xmath0-ray spectra for ngc 253 had a significant contribution from pion decay . while this results in a significant neutrino flux ( see figure 6 ) , the steep spectrum of ngc 253 means that m82 is a more likely candidate for eventual detection of the energetic neutrinos for which icecube has good sensitivity @xcite . though all of our models for ngc 1068 underestimated the @xmath0-ray emission , we find that the corresponding neutrino flux from the starburst region is comparable to that of ngc 253 at tev energies ( see figure 6 ) . however , if the high - energy emission in the jet is hadronic , ngc 1068 may still be a detectable source of neutrinos , particularly if the flat @xmath0-ray spectrum seen in figure 5 continues to higher energies . at energies below 1 gev , models predict a significant difference in the spectra for inverse compton and pion decay @xmath0-ray emission @xcite . thus , in the future , hard x - ray observations , such as those by nustar , should be able to distinguish between leptonic and hadronic emission mechanisms for @xmath0-ray production in such galaxies and determine whether these galaxies are likely to be detectable in neutrinos . in this paper , we explore the differences between starburst galactic nuclei and agn and their impact on cosmic ray interactions . as a typical nuclear starburst environment , ngc 253 can be fit with our model that includes inverse compton emission . we find that for an assumed gas mass of @xmath111 @xmath5 , the @xmath0-ray and radio spectra can not uniquely determine the magnetic field strength or wind speed . additionally , the @xmath0-ray spectrum also severely limits the possible combinations of spectral index ( to @xmath89 ) and acceleration efficiency ( to @xmath93 ) due to the powerful combination of effects due to inverse compton and pion decay emission . past models for ngc 253 have assumed a fixed wind ( advection ) speed which presupposes calorimetry ( e.g. * ? ? ? based on our results , a wind is a critical component for determining the degree to which a starburst system is a calorimeter . while the best - fits to the @xmath0-ray spectrum result in a @xmath2250% proton calorimeter , fits to the radio spectrum result in a complete proton calorimeter . additionally , our fits result in systems below equipartition ( in regards to the cosmic ray energy densities ) with strong magnetic fields ( see also * ? ? ? the lack of a solution which fits the @xmath0-ray and radio spectra together with a single set of parameters and the problems in fitting the radio spectrum highlight the limitations of this one - zone model . we , however , can obtain a good fit to @xmath0-rays and the radio by assuming a lower gas mass . understanding of the ism conditions is vital to applying these models to starbursts . as such , while the model works well for a late phase starburst such as m82 , it does not as easily fit an early phase starburst such as ngc 253 .    as an agn environment containing a possible starburst region , ngc 1068 is not fit with our self - consistent starburst model . while we find a few models that come close to predicting the observed @xmath0-ray spectrum for ngc 1068 , all of our models underestimate the high - energy emission and we are unable to find any models of the inner starburst region that agree with the observed data . thus , we find that the cosmic ray population from the starburst region alone is not enough to accurately fit the observed @xmath0-ray emission and even those models in rough agreement with high - energy observations result in an overestimation of the radio spectrum . we conclude that the observed @xmath0-ray emission likely occurs in the agn or its jets .    additionally , we find that starburst galaxies such as ngc 253 are potential sources for high energy neutrino detection . however , agn with high @xmath112 such as ngc 1068 will not be luminous neutrino sources unless their high - energy emission is of hadronic origin . this work was supported in part by nsf ast-0907837 , nsf phy-0821899 ( to the center for magnetic self - organization in laboratory and astrophysical plasmas ) , and nsf phy-0969061 ( to the icecube collaboration ) . we thank alberto bolatto for conversations concerning ngc 253 . we also thank francis halzen for his help and support . abbasi , r. , abdou , y. , abu - zayyad , t. , et al . 2011 , , 732 , 18 abramowski , a. , acero , f. , aharonian , f. , et al . 2012 , , 757 , 158 ackermann , m. , ajello , m. , allafort , a. , et al . 2012 , , 755 , 164 anchordoqui , l.  a. , hooper , d. , sarkar , s. , & taylor , a.  m.  2008 , aph , 29 , 1 blumenthal , g.  r. , & gould , r.  j.  1970 , rev . , 42 , 273 bock , j.  j. , neugebauer , g. , matthews , k. , et al . 2000 , , 120 , 2904 boettcher , e. , zweibel , e.  g. , yoast - hull , t.  m. , & gallagher , j.  s. , iii  2013 , , 779 , 12 carilli , c.  l.  1996 , , 305 , 402 condon , j.  j.  1992 , , 30 , 575 davidge , t.  j.  2010 , , 725 , 1342 davies , r.  i. , sugai , h. , & ward , m.  j.  1998 , , 300 , 388 domingo - santamara , e. , & torres , d.  f.  2005 , , 444 , 403 engelbracht , c.  w. , rieke , m.  j. , rieke , g.  h. , kelly , d.  m. , & achtermann , j.  m.  1998 , , 505 , 639 exposito , j. , gratadour , d. , clnet , y. , & rouan , d.  2011 , , 533 , 63 fabbiano , g. , & trinchieri , g.  1984 , , 286 , 491 forbes , d.  a. , polehampton , e. , stevens , i.  r. , brodie , j.  p. , & ward , m.  j.  2000 , , 312 , 689 gallimore , j.  f. , baum , s.  a. , & odea , c.  p.  1996a , , 464 , 198 gallimore , j.  f. , baum , s.  a. , odea , c.  p. , & pedlar , a.  1996b , , 458 , 136 gallimore , j.  f. , baum , s.  a. , & odea , c.  p.  2004 , , 613 , 794 hailey - dunsheath , s. , nikola , t. , stacey , g.  j. , et al . 2008 , , 689 , l109 hailey - dunsheath , s. , sturm , e. , fischer , j. , et al .  2012 , , 755 , 57 halzen , f.  2006 , epjc , 46 , 669 heesen , v. , beck , r. , krause , m. , & dettmar , r .- j . 2009a , , 494 , 563 heesen , v. , krause , m. , beck , r. , & dettmar , r .- j .  2009b , , 506 , 1123 hnig , s.  f. , prieto , m.  a. , & beckert , t.  2008 , , 485 , 33 kamenetzky , j. , glenn , j. , maloney , p.  r. , et al . 2011 , , 731 , 83 karachentsev , i.  d. , grebel , e.  k. , sharina , m.  e. , et al . 2003 , , 404 , 93 karachentsev , i.  d. , tully , r.  b. , dolphin , a. , et al . 2007 , , 133 , 504 kelner , s.  r. , aharonian , f.  a. , & bugayov , v.  v.  2006 , , 74 , 034018 kornei , k.  a. , & mccrady , n.  2009 , , 697 , 1180 krips , m. , martn , s. , eckart , a. , et al . 2011 , , 736 , 37 lacki , b.  c. , & beck , r.  2013 , , 430 , 3171 lacki , b.  c. , thompson , t.  a. , & quataert , e.  2010 , , 717 , 1 lacki , b.  c. , thompson , t.  a. , quataert , e. , loeb , a. , & waxman , e.  2011 , , 734 , 107 lehmer , b.  d. , wik , d.  r. , hornschemeier , a.  e. , et al . 2013 , , 771 , 134 lenain , j .- p . , ricci , c. , trler , m. , dorner , d. , & walter , r.  2010 , , 524 , 72 lenc , e. , & tingay , s.  j.  2006 , , 132 , 1333 mitsuishi , i. , yamasaki , n.  y. , & takei , y.  2011 , , 742 , l31 mller - snchez , f. , gonzlez - martn , o. , fernndez - ontiveros , j.  a. , acosta - pulido , j.  a. , & prieto , m.  a.  2010 , , 716 , 1166 muxlow , t.  w.  b. , pedlar , a. , holloway , a.  j. , gallimore , j.  f. , & antonucci , r.  r.  j.  1996 , , 278 , 854 paglione , t.  a.  d. , & abrahams , r.  d.  2012 , , 755 , 106 paglione , t.  a.  d. , marscher , a.  p. , jackson , j.  m. , & bertsch , d.  l.  1996 , , 460 , 295 peng , r. , zhou , s. , whiteoak , j.  b. , lo , k.  y. , & sutton , e.  c.  1996 , , 470 , 821 persic , m. , rephaeli , y. , & arieli , y.  2008 , , 486 , 143 ricci , r. , prandoni , i. , gruppioni , c. , sault , r.  j. , & de zotti , g.  2006 , , 445 , 465 rybicki , g.  b. , & lightman , a.  p.  1979 , radiative processes in astrophysics ( new york : wiley ) sakamoto , k. , mao , r .- q . , matsushita , s. , et al . 2011 , , 735 , 19 schinnerer , e. , eckart , a. , tacconi , l.  j. , genzel , r. , & downes , d.  2000 , , 533 , 850 schlickeiser , r.  2002 , cosmic ray astrophysics ( new york : springer ) schlickeiser , r. , vainio , r. , bttcher , m. , et al . 2002 , , 393 , 69 schulz , h. , & wegner , g.  1992 , , 266 , 167 scoville , n.  z. , matthews , k. , carico , d.  p. , & sanders , d.  b.  1988 , , 327 , l61 spinoglio , l. , pereira - santaella , m. , busquet , g. , et al . 2012 , , 758 , 108 stecker , f.  w.  1979 , , 228 , 919 strickland , d.  k. , heckman , t.  m. , colbert , e.  j.  m. , hoopes , c.  g. , & weaver , k.  a.  2004 , , 151 , 193 strickland , d.  k. , heckman , t.  m. , weaver , k.  a. , dahlem , m.  2000 , , 120 , 2965 storchi - bergmann , t. , riffel , r.  a. , riffel , r. , et al . 2012 , , 755 , 87 telesco , c.  m. , & harper , d.  a.  1980 , , 235 , 392 thompson , t.  a. , quataert , e. , & waxman , e.  2007 , , 654 , 219 thronson , h.  a. , jr . , hereld , m. , majewski , s. , et al . 1989 , , 343 , 158 torres , d.  f.  2004 , , 617 , 966 ulrich , m .- h .  1978 , , 219 , 424 ulvestad , j.  s. , & antonucci , r.  r.  j.  1997 , , 488 , 621 van buren , d. , & greenhouse , m.  a.  1994 , , 431 , 640 waller , w.  h. , kleinmann , s.  g. , ricker , g.  r.  1988 , , 95 , 1057 wei , a. , kovcs , a. , gsten , r. , et al . 2008 , , 490 , 77 westmoquette , m.  s. , smith , l.  j. , & gallagher , j.  s. , iii  2011 , , 414 , 3719 williams , p.  k.  g. , & bower , g.  c.  2010 , , 710 , 1462 wilson , a.  s. , & ulvestad , j.  s.  1983 , , 275 , 8 yoast - hull , t.  m. , everett , j.  e. , gallagher , j.  s. , iii , & zweibel , e.  g.  2013 , , 768 , 53 ( yegz )
Abstract:

Answer: the galaxies m82 , ngc 253 , ngc 1068 , and ngc 4945 have been detected in @xmath0-rays by _ fermi_. previously , we developed and tested a model for cosmic ray interactions in the starburst galaxy m82 . now , we aim to explore the differences between starburst and active galactic nuclei ( agn ) environments by applying our self - consistent model to the starburst galaxy ngc 253 and the seyfert galaxy ngc 1068 . assuming constant cosmic - ray acceleration efficiency by supernova remnants with milky - way parameters , we calculate the cosmic - ray proton and primary and secondary electron / positron populations , predict the radio and @xmath0-ray spectra , and compare with published measurements . we find that our models easily fits the observed @xmath0-ray spectrum for ngc 253 while constraining the cosmic ray source spectral index and acceleration efficiency . however , we encountered difficultly modeling the observed radio data and constraining the speed of the galactic wind and the magnetic field strength , unless the gas mass is less than currently preferred values . additionally , our starburst model consistently underestimates the observed @xmath0-ray flux and overestimates the radio flux for ngc 1068 ; these issues would be resolved if the agn is the primary source of @xmath0-rays . we discuss the implications of these results and make predictions for the neutrino fluxes for both galaxies .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: as very few gamma - ray burst ( grb ) sources have astronomical counterparts at other wavebands , empirical studies of grbs have been largely restricted to the analysis of their gamma ray properties : bulk properties such as fluence and spectral hardness , and evolution of these properties within a burst event ( fishman & meegan 1995 ) . while bursts exhibit a vast range of complex temporal behaviors , their bulk properties appear simpler and amenable to straightforward statistical analyses . studies fall into two categories : examination of whether grb bulk properties comprise a homogeneous population or are divided into distinct classes ; and search for relationships between bulk properties . both types of study may lead to astrophysical insight , just as the distinction between main sequence stars and red giants and the measurement of a luminosity - mass relation along the main sequence assisted the development of stellar astrophysics early in the century . the most widely accepted taxonomy of grbs is the division between short - hard and long - soft bursts proposed by dezalay et al . ( 1992 ) and kouveliotou et al . ( 1993 , henceforth k93 ) . k93 noticed a bimodality in the burst duration variable @xmath1 ( time within which 90% of the flux arrived ) , suggesting the presence of two distinct types of bursts separated at @xmath2 sec . the short bursts have systematically harder gamma - ray spectra than longer bursts . the two groups seemed indistinguishable in most other bulk properties , although the larger group of long - soft bursts may have a subclass with a different fluence distribution ( i.e. , different @xmath3 ; katz & canel 1996 ) and the groups may have different galactic latitude distributions ( belli 1997 ) . other researchers point to small groups of bursts with distinctive properties such as the soft - gamma repeaters ( norris et al . 1991 ) , two possible classes with differing short - timescale variability ( lamb , graziani & smith 1993 ) , fast - rise exponential - decay bursts ( bhat et al . 1994 ) , and two types of bursts with different ratios of total fluence and @xmath4300 kev fluence ( pendleton et al . 1997 ) .    a variety of relationships between burst properties have also been reported . norris et al . ( 1995 ) find an anti - correlation between @xmath1 ( calculated after wavelet thresholding ) and peak intensity , consistent with a cosmological time dilation . however , a positive correlation between @xmath1 and total fluence is also seen which does not agree with the simplest cosmological interpretation ( lee & petrosian 1997 ) . additional reported relationships include : @xmath1 correlated with peak heights ( lestrade 1994 ) , peak energy correlated with peak flux ( mallozzi et al . 1995 ) , and peak duration anticorrelated with gamma - ray energy ( fenimore et al . 1995 ) . most of these studies suffer from a failure to treat all of the bulk property variables in an unbiased and quantitative way . astronomers typically examine univariate or bivariate distributions , sometimes constructing composite variables ( such as hardness ratios ) with pre - determined relationships to include one or two additional variables . but it is quite possible that the complex astrophysics producing grbs will not manifest themselves in simple bivariate plots , just as the division between short - hard and long - soft bursts is not evident in spectral variables alone ( pendleton et al . grb catalogs , like most multiwavelength astronomical catalogs , are multivariate databases and should be treated with multivariate statistical methods that can objectively and effectively uncover structure involving many variables ( feigelson & babu 1997 ) . two previous studies take a fully multivariate approach to understanding grb bulk properties . baumgart ( 1994 ) constructs a neural network taxonomy of 99 grbs from the pvo satellite using 26 variables representing both bulk burst properties and detailed temporal characteristics ( e.g. number of peaks , fractal dimension , wavelet transform crossings ) and finds two or three distinct grb classes . bagoly et al . ( 1997 ) perform principal components and factor analyses of nine bulk property variables using 625 grbs from the batse 3b catalog . they find that the relationships in the database are determined principally by only three variables : an appropriately weighted fluence , a weighted burst duration , and ( to a lesser extent ) flux in the highest energy bin . we note , however , that it can be dangerous to look for correlations prior to classifying ( or establishing the homogeneity of ) the population . while the anticorrelation between hardness ratio and burst duration seen in full samples ( k93 ) may be the manifestation of a single astrophysical process , it may alternatively reflect differences between distinct processes . the latter possibility is suggested by a reported hardness - duration positive correlation within the long - soft class of bursts ( dezalay et al . 1996 ; horack & hakkila 1997 ) . most multivariate analyses thus begin with a study of homogeneity and classification , and then investigate the variance - covariance structure ( i.e. correlations ) within each class . this paper describes a multivariate analysis of grbs from the third batse catalog ( meegan et al . 1996 ) . after defining the sample (  2 ) , we start with a simple statistical description of the variables and their bivariate relationships for the entire dataset (  3 ) . we then seek distinct types of clusters in two ways . first , a standard nonparametric agglomerative hierarchical clustering analysis is performed (  4 ) which reveals three distinct classes . the statistical significance of the third cluster is validated , under gaussian assumptions , with manova tests . second , a parametric maximum likelihood model - based clustering procedure is adopted which reveals the same three groups and indicates strong evidence for the presence for three rather than two groups (  5 ) . the variance - covariance structure of each group is then examined (  6 ) . results are synthesized in the discussion (  7 ) .    throughout the paper , we discuss our mathematical techniques to help the reader understand the complexities of multivariate analysis . from the vast literature in this subject , we recommend the following monographs for interested readers : johnson & wichern ( 1992 ) and jobson ( 1992 ) for overviews of applied multivariate analysis ; hartigan ( 1975 ) , jain & dubes ( 1988 ) and kaufman & rousseeuw ( 1990 ) for multivariate clustering algorithms ; murtagh & heck ( 1987 ) and , more briefly , babu & feigelson ( 1996 ) and feigelson & babu ( 1997 ) , for multivariate methodology in astronomy . our sample is drawn from the third catalog of the _ burst and transient source experiment _ ( batse ) on board the compton gamma ray observatory . this 3b catalog has 1122 grbs detected by batse between 1991 april 19 and 1994 september 19 . the catalog is presented and fully described by meegan et al . our database was extracted from the on - line database www.batse.msfc.nasa.gov/data/grb/catalog in may 1996 , which provides many properties of each burst . there are roughly eleven variables of potential astrophysical interest : two measures of location in galactic coordinates , @xmath5 and @xmath6 ; two measures of burst durations , the times within which 50% ( @xmath7 ) and 90% ( @xmath1 ) of the flux arrives ; three peak fluxes @xmath8 , @xmath9 and @xmath10 measured in 64 ms , 256 ms and 1024 ms bins respectively ; and four time - integrated fluences @xmath11-@xmath12 in the 20 - 50 kev , 50 - 100 kev , 100 - 300 kev and @xmath13 kev spectral channels respectively . researchers commonly consider three composite variables : the total fluence , @xmath14 , and two measures of spectral hardness derived from the ratios of channel fluences , @xmath15 and @xmath16 . due to the limitations of available multivariate statistical techniques , we ignore other variables of potential relevance including the heteroscedastic measurement errors of each quantity ( i.e. , errors that differ from point to point ) and truncation values associated with batse triggering operations .    of the 1122 listed bursts , 807 have data on all the variables described above . ten bursts listed with zero fluences were eliminated . our sample thus has 797 batse grbs . for some analyses , we also used a subset of 644 bursts with ` debiased ' durations , @xmath17 . here the durations are modified to account for the effect that brighter bursts will have signal above the noise for longer periods than fainter bursts with the same time profiles ( j. norris , private communication ) . statistical analyses in   3 , 4 and 6 were conducted within the statistical analysis system _ sas / stat _ , a very large and widely distributed commercial statistical software package ( sas institute inc . _ sas / stat _ procedures _ cluster _ , _ glm _ , _ princomp _ and _ varclus _ were used . the analysis in  5 was performed with the _ mclust _ software ( banfield & raftery 1993 ; fraley 1998 ) , which is interfaced to the splus statistical package ( splus version 3.4 ; mathsoft , inc . 1996 ) and its extensions . further information is provided at http://stat.washington.edu/fraley/software.html and http://lib.stat.cmu.edu/general/mclust . for multivariate data visualization , we used the xgobi ( swayne , cook & buja 1991 ) program , available from the _ statlib _ software archive at http://lib.stat.cmu.edu/general/xgobi . hypertext links to a variety of public domain software for multivariate analysis , classification and visualization are available at the penn state _ statcodes _ web site , http://www.astro.psu.edu/statcodes . we are faced with a multivariate database of 797 objects and 15 variables ( 11 variables from the catalog , 3 composite variables , and @xmath18 ) . two initial problems are frequently faced in analyses of multivariate databases . first , variables with incompatible units and ranges must be compared . units can be removed by normalization ( e.g. replacing @xmath11 by @xmath19 ) , by standardization ( e.g. replacing @xmath11 by @xmath20 where @xmath21 is the sample standard deviation ) , or by taking logarithms . second , the dimensionality of the problem should be reduced , as many of the variables are closely interrelated either by construction or by astronomical circumstance . although there are no mathematical rules regulating reduction of dimensionality , it can usefully be guided by a correlation matrix showing bivariate relationships and a principal components analysis showing multivariate relationships that are mainly responsible for structured variance in the data . scientific reasoning can also be used to eliminate consideration of variables . we conducted a preliminary examination of data representations , correlation matrices and bivariate plots , and principal components analyses to facilitate choice of variables . when no mathematical preference arose , we selected variables most commonly used by previous researchers to facilitate comparison of results . our choices were as follows . we use log variables , rather than normalized or standardized variables . we kept information on burst fluence and spectra through @xmath22 and hardness ratios rather than through the original fluences @xmath11-@xmath12 . we initially eliminated @xmath8 and @xmath10 from consideration , and later eliminated @xmath9 when we found that its main contribution to the clustering process was noise . we chose to remove the location variables @xmath23 , already established by other researchers to be random for the entire sample , but use them later to test for isotropy of subsamples . the debiased @xmath24 is used only in special tests . our analysis was thus performed in six or fewer dimensions using log @xmath7 , log @xmath1 , log @xmath22 , log @xmath9 , log @xmath25 and log @xmath26 . tables 1 and 2 give basic statistics for these six variables : means , standard deviations , and bivariate values of pearson s linear correlation coefficient @xmath27 . for @xmath28 and assuming bivariate normal populations , any @xmath29 implies that a correlation between the two variables exists at a two - tailed significance level @xmath30 ( beyer 1968 , pp . 389 and 283 ) . but from an astrophysical perspective , we might consider any relationship with @xmath31 to be of little interest . figure 1 shows the bivariate scatter plots . the correlation structure of the entire sample ( table 2 ) shows that the two measures of duration and the two measures of spectral hardness have correlation near unity , indicating that they are nearly redundant . @xmath22 and @xmath9 are quite dissimilar : @xmath22 shows a strong correlation with burst duration ( e.g. lee & petrosian 1997 ) and no relation to hardness , while @xmath9 shows no relation to duration but is mildly correlated with hardness ( mallozzi et al . burst duration is anticorrelated with hardness ( k93 ; fenimore et al . the cosmological anticorrelation between duration and peak flux reported by norris et al . ( 1995 ) is statistically significant , but accounts for only a few percent of the variance between these variables . the correlation matrix based on the debiased @xmath18 values yields very similar results . however , the scatter plots ( figure 1 ) show a more complex story . first , many plots show inhomogeneous distributions inconsistent with the unimodal multinormal ( i.e. multivariate gaussian ) population assumed by pearson s @xmath27 . the distributions often seem bimodal with asymmetrical non - gaussian shapes . one outlier burst is also seen in several projections . we therefore consider the hypothesis that the sample consists of two or more distinct classes , and proceed to find the ` clusters ' using well - established methods . agglomerative hierarchical clustering is a procedure based on the successive merging of proximate pairs of clusters of objects . it produces a clustering tree or dendrogram starting with n clusters of 1 member ( or a coarse partition based on prior knowledge ) and ending with one cluster of n members . unfortunately , there are many possible ways to proceed ; mathematics provides little guidance among the choices and no probabilistic evaluation of the results without the imposition of additional assumptions . the scientist must make four decisions to fully define the clustering procedure :    1 .   creating _ unit - free variables _ is essential for meaningful treatment of objects in multivariate space (  3 ) . a favorite choice by statisticians is standardization , where each variable is normalized by the standard deviation of the sample . astronomers more commonly make logarithmic transformations or construct ratios of variables sharing the same units . we follow the tradition of grb researchers by measuring spectral hardness with ratios of fluences having the same units , and making logarithmic transformations of all variables . the _ metric _ defines the meaning of proximity between two objects or clusters . common choices are the simple euclidean distance between unit - free variables and the squares of euclidean distances . we chose the former option for most of the analyses in this section . 3 .   several _ merging procedures _ can be used . one might begin by merging the clusters with the nearest neighbors . this is called single linkage clustering and is most familiar in astronomy where it is frequently called the friends - of - friends algorithm . it tends to produce long stringy clusters , and is equivalent to a well - known divisive clustering procedure known as pruning the minimal spanning tree . complete linkage proceeds by maximizing the distance between clusters , and leads to evenly bifurcating dendrograms . for most of our analysis , we choose average linkage where the distance between two clusters is the average of the distances between pairs of observations , where each member of the pair comes from a different cluster . this is a compromise between single and complete linkage and tends to give compact clusters . specifically , the distance between clusters k and l is given by ( e.g. sas institute inc . 1989 , pp . 529ff ; johnson & wichern 1992 , pp . @xmath32 where the bar indicates an unweighted mean , @xmath33-th cluster . another popular choice is ward s minimum variance criterion where the distance between the two clusters is the anova ( analysis of variance ) sum of squares between two clusters added up over all variables ( ward 1963 ) , @xmath34 if the sample is generated by a mixture of multinormal ( i.e. multidimensional gaussian ) distributions where each distribution has covariance matrix of the form @xmath35 , this method joins clusters to maximize the likelihood at each level of the hierarchy , and so is a special case of the model - based clustering methodology to be discussed in  5 . 4 . as the procedure gives a hierarchy from @xmath36 clusters with 1 object down to one cluster with @xmath36 objects , the user must choose _ how many clusters to report _ as scientifically important clusters . this choice can be assisted by examination of two statistics . the squared correlation coefficient , @xmath37 , states the fraction of the total variance accounted for by a partition into @xmath38 clusters , @xmath39 the squared semi - partial correlation coefficient , @xmath40 measures the difference in the variance between the resulting cluster and the immediate parent clusters normalized by the total sample variance , @xmath41 @xmath37 thus tells how much of the scatter is explained by a given level of clustering , and @xmath40 tells how much improvement is achieved between levels . we emphasize again that there is no mathematically ` best ' choice , although extensive experience with problems in many fields has led to a preference for certain combinations ( e.g. standardized variables and ward s minimum variance criterion ) . we conducted extensive experiments with different choices . the last several levels of the clustering tree for the 797 grbs using the six unit - free variables shown in table 1 , average linkage and a euclidean metric are shown in figure 2 ( left panel ) with details in table 3a . the action taken at each level is indicated in column 2 of table 3 , which may refer to a level higher in the tree which ( for brevity ) is not shown here . two types of mergers are seen : the incorporation of ` twigs ' of one or a few grbs into a large preexisting ` trunk ' ( levels 1 , 3 , 4 , 5 and 7 ) ; and the union of two substantial branches into a single larger trunk ( levels 2 , 6 and 8) . the first type has little effect on the variance of the sample with @xmath42% . the single grb brought into the main trunk at level 1 is the distant outlier seen in several panels of figure 1 . the level 2 merger of clusters with 190 and 606 members is clearly the most important structure , accounting for roughly 53% of the variance of the entire sample . this is the bifurcation of the sample into two classes easily seen in figure 1 and noted by k93 and others . the principal finding that is not immediately obvious from figure 1 is the structure indicated at level 6 . the main trunk of 599 bursts ( plus a few twigs to be merged later ) is divided into groups of 93 and 506 bursts . this division accounts for 10% of the total variance of the sample , indicated in both the @xmath37 and @xmath40 values . we found that the twigs in the tree structure disappear if the peak flux @xmath9 variable is omitted and the analysis is made in 5-dimensional space ( figure 2 , right panel , and table 3b ) . here the largest cluster of 593 members is formed by the union of clusters with 107 and 486 bursts , again accounting for 10% of the sample variance . it is possible that @xmath9 is a nuisance variable irrelevant to the basic astrophysics of grbs , producing noisy ` twigs ' seen in table 3a and figure 2 ( left panel ) . we tested many variants of hierarchical clustering . we replaced average linkage hypothesis with complete linkage , single linkage and ward s minimum variance criterion . the ward s criterion computation , for example , gave three clusters with 468 , 145 and 184 bursts . we clustered using nonparametric density estimation based on the 100 nearest neighbors , and clustered using the principal components rather than the observed variables . various methods were tried with both the observed @xmath1 values and debiased @xmath18 values , with little effect on the results . all methods showed two strong clusters and the outlier but , in some cases , the third cluster appeared only weakly .    to proceed further , we choose a single clustering structure for detailed study : the 5-dimensional average linkage analysis ( table 3b ) with three clusters  class i with 486 bursts , class ii with 203 bursts , class iii with 107 bursts . class iv consisting of the single outlier is ignored because of independent evidence that its properties are due to data of poor quality (  7 ) . the membership of these clusters is given in table 4 , and four projections of the clusters onto two - dimensional scatter plots are shown in figure 3 . these are frames from the ` grand tour ' movie of the 5-dimensional dataset provided by the xgobi software where each cluster is ` brushed ' with a different symbol . note that , in general , there is no reason why classification structure should be most evident in projections parallel to the variable axes shown in figure 1 . it is more important that the clusters show cohesion in many projections of the dataset . the grand tour of the 797 grbs shows that classes i , ii and the outlier are very distinct in most projections . class iii often lies between classes i and ii ( e.g. top panels of figure 3 ) , but in other projections is offset from the line between classes i and ii ( e.g. bottom panels of figure 3 ) . it also appears elongated along some projections , while the larger classes i and ii appear roughly hyperspherical . this analyses described here provide considerable evidence for three major clusters and an outlier . but as some nonparametric clustering procedures did not find a strong third cluster , there is some some worry that class iii is simply a group of bursts with properties intermediate between the classes i and ii . while nonparametric hierarchical clustering methods can not address this question , it can be investigated with parametric methods . mathematically well - founded methods for evaluating the statistical significance of a proposed multivariate classification scheme are available under the assumption that the population is a multinormal mixture ; that is , the objects of each class are drawn from multivariate gaussians . all relationships between the variables must thus be linear ( as in table 2 ) , although the relationships may differ between clusters . there is no requirement of sphericity , so that clusters may have shapes akin to pancakes or cigars with arbitrary orientations in multidimensional space . the separate existence of each of the postulated subpopulations can be tested using multivariate analysis of variance ( manova ) . the results of our manova calculations are summarized in table 5 . the columns give the values of the three manova statistics , followed by details for the wilks @xmath0 : the corresponding value of the @xmath58 statistic , the numerator and denominator degrees of freedom for that @xmath58 value ; and the resulting @xmath60 value . details for pillai s and hotelling - lawley s traces are omitted , but give similar results in all cases . the first row tests the null hypothesis that the classes i , ii and iii have the same mean , the second row tests the equality of classes i and ii , and so forth . the @xmath58 values are very high in all cases , indicating that the clusters are different with extremely high statistical significance ( @xmath61 ) . as the tails of the distribution are poorly determined unless the sample size is extremely large . ] this is a clear _ quantitative _ demonstration that at least two clusters exist among grbs ( for a univariate test , see ashman , bird & zepf 1994 ) , which was _ qualitatively _ reported by delazay et al . ( 1992 ) and k93 . the other rows in table 5 test the hypotheses that each proposed class has the same mean as each other class . one problem with these manova tests is that they are conditional on the classification that has been found using the clustering algorithm . because the clustering algorithm is constructed to find groups that are different from one another , tests such as these tend to be biased towards finding structure , perhaps where none exists . although the manova results seem to indicate very strong evidence of structure , they can not be taken as definitive for this reason . tests arising from model - based clustering can overcome this problem , as discussed in the following section . in the previous section , we conducted a hierarchical clustering analysis without making assumptions regarding the shapes of the clusters , but needed the parametric assumption of normality to estimate the statistical significance of the resulting classification . it is reasonable to conduct the entire analysis , both clustering and validation , within a model - based framework . we report here an analysis of this type again assuming that the grb population consists of a mixture of multivariate gaussian classes . early development of this model for clustering is discussed in mclachlan & basford ( 1988 ) ; we use more recent developments here . first , an initial classification for each possible number of clusters is found via agglomerative hierarchical clustering ( murtagh & raftery 1984 , banfield & raftery 1993 , fraley 1998 ) . next , the em ( expectation - maximization ) algorithm is used to refine partitions obtained from hierarchical clustering ( celeux & govaert 1995 , dasgupta & raftery 1998 ) . finally , the bayesian information criterion is used to select the ` best ' partition among those associated with different numbers of clusters ( dasgupta & raftery 1998 ) .    in the model considered here , the @xmath43-dimensional observations @xmath62 are drawn from @xmath38 multinormal groups , each of which is characterized by a vector of parameters @xmath63 for @xmath64 . our goals are : to determine the number of grb types , @xmath38 ; to determine the cluster assignment of each burst ; and to estimate the mean @xmath65 and covariance matrix @xmath66 for each cluster . following fraley ( 1998 ) , the density of an observation @xmath67 from the @xmath68th subpopulation is expressed as follows : @xmath69 where @xmath70 means multivariate normal . we estimate the parameters using the principle of maximum likelihood . in the hierarchical clustering phase , we use the classification likelihood @xmath71 where @xmath72 represents the observations and @xmath73 is the cluster assignment : @xmath74 when @xmath75 comes from the @xmath68-th group . equivalently , @xmath76 where @xmath77 is the set of indices corresponding to observations belonging to the @xmath68-th group . the method used here for maximizing the likelihood ( fraley 1998 ) and implemented in the _ mclust _ code involves parameterization of the @xmath78 matrices in terms of their eigenvectors and eigenvalues ( analogous to a principal components analysis ) , and iterative relocation of the clusters using the em algorithm . the em ( estimation - maximization ) algorithm ( dempster , laird & rubin 1977 ) , one of the most successful methods in modern statistics , is a procedure for iteratively maximizing likelihoods in a wide variety of circumstances . for example , the lucy - richardson algorithm in astronomical image restoration is the em algorithm . in the present application , we apply em to the mixture likelihood @xmath79 where @xmath80 are mixing probabilities associated with each group . for a given number @xmath38 of components in the mixture , we use em to estimate the conditional probability that observation @xmath62 belongs to the @xmath68-th group for each @xmath45 and selected @xmath68 via maximum likelihood . although the computational procedure has some limitations ( e.g. convergence of the em iterations is not guaranteed ; clusters can not be extremely small ) , it is generally efficient and effective for gaussian clustering problems when started from reasonable partitions such as those produced by hierarchical agglomeration . we use the bayes factor to assess the evidence for a given number of clusters against a different number of clusters . the bayes factor , defined in the context of bayesian statistics , is the posterior odds for one model against the other when the prior odds are equal to one ( i.e. , when one does not favor one model over the other _ a priori _ ) . kass & raftery ( 1995 ) review the use of bayes factors in adjudicating between competing scientific hypotheses on the basis of data . the bayes factor for a model @xmath81 against a competing model @xmath82 ( say , for three @xmath83 two classes of grbs ) is defined as @xmath84 where @xmath85 ) for @xmath86 is obtained by integrating the likelihood times the prior density over the parameters of the model . it can be viewed as a likelihood ratio , but it differs from the usual frequentist ratio that underlies the likelihood ratio test in that the latter is obtained by maximizing ( rather than integrating ) the likelihood over the model parameters .    twice the logarithm of the bayes factor can be approximated by the bayesian information criterion or bic ( schwarz 1978 ) , @xmath87 where @xmath88 is the likelihood and @xmath89 is the number of parameters for one mixture model , and similarly for @xmath90 and @xmath91 . the bic measures the balance between the improvement in the likelihood and the number of model parameters needed to achieve that likelihood . while the absolute value of the bic is not informative , differences between the bic values for two competing models provide estimates of the evidence in the data for one model against another . conventionally , bic differences @xmath92 represent weak evidence , differences between 2 and 6 represent positive evidence , @xmath93 strong evidence , and @xmath94 very strong evidence ( jeffreys 1961 , appendix b ; kass & raftery 1995 ) . the use of the bic in choosing clusters in a mixture or clustering model is discussed by roeder & wasserman ( 1997 ) and dasgupta & raftery ( 1998 ) . bayes factors and bic have the advantage that they can be used to assess the evidence _ for _ a null hypothesis , unlike standard significance tests which can only reject a null hypothesis . they can also easily be used to compare non - nested models , again unlike standard significance tests which require competing models to be nested .      to reduce the dimensionality of the problem and the complexity of the calculation , we eliminated the highly redundant @xmath7 and @xmath26 variables ( see figure 1 ) and considered only the 3 variables @xmath1 , @xmath22 and @xmath25 for the sample of 797 batse grbs . the _ mclust _ model - based clustering procedure described above was run for trials of @xmath95 groups . the resulting values of bic(@xmath38 ) are plotted in figure 4 . the maximum bic is achieved for three classes . most importantly , the bic value for @xmath96 is @xmath97 units above that for @xmath98 . this corresponds to strong evidence indeed for the presence of three groups rather than two . this result strongly confirms the analysis in  4 indicating the existence of three clusters , and this time is free of the problem that the manova tests are conditional on the estimated partition . the result here takes account of the fact that the partition is not known in advance . we have also calculated the bic for @xmath99 with various constraints on the covariance matrix * @xmath100 * such as hypersphericity and uniformly shaped ellipsoids . spherical clusters give poor fits . uniform ellipsoids give good fits with 4 and 8 clusters . but in all cases , the maximum likelihood assuming 2 clusters is much lower than the likelihood of @xmath101 clusters . the cluster assignment vector @xmath102 for the @xmath96 model with unconstrained * @xmath100 * is given in table 4 . over 85% of the assignments are the same as those obtained from the nonparametric hierarchical clustering procedure in  4 , so that we note only differences between the two clustering results with @xmath103 and @xmath104 markings . all but one of the 96 assignment differences move bursts from classes i and ii into class iii . the close agreement between the cluster assignments in the two methods reinforces confidences in the conclusions from both of them . we can now examine the properties of grbs within each cluster with reasonable confidence that the populations are distinct from each other but internally homogeneous . these properties become inputs to astrophysical theories seeking to explain grb bulk properties . table 6a lists the means and standard deviations of the principal variables for each cluster based on both the nonparametric and model - based clustering procedures . the two methods give very similar results . the three types are well - separated in the burst duration variables : cluster i bursts have the longest durations around @xmath105-@xmath106 seconds , cluster ii bursts have the shortest durations below 1 second , and cluster iii bursts have intermediate durations around 2 seconds . this is shown clearly in figure 5 , which projects each class onto the univariate @xmath1 axis . cluster iii bursts are also intermediate in their fluences , although their fluence distribution overlaps that of the fainter class ii bursts . the hardness ratios of all three clusters overlap considerably , but class iii bursts have the softest spectra and class i bursts have intermediate spectra . we can thus classify the types in the three principal dimensions duration / fluence / spectrum ( table 6c ) : class i is long / bright / intermediate , class ii is short / faint / hard , and class iii is intermediate / intermediate / soft . a major constraint for the astrophysical interpretation of grbs has been the remarkable isotropy of their spatial distribution in the celestial sphere . it is possible that , while the bulk of grbs are isotropic and have an inferred extragalactic origin , some class of grbs have significant anisotropy which would reflect a galactic origin ( see lamb 1995 ) . we apply four statistical tests for isotropy discussed by briggs ( 1993 ) and applied by briggs et al . ( 1996 ) to various subsamples of the 3rd batse catalog of grbs . the statistics are : @xmath107cos@xmath108 , where @xmath109 is the angle between a burst and the galactic center ; @xmath107sin@xmath110 , where @xmath6 is the galactic latitude ; rayleigh - watson @xmath111 ; and bingham @xmath112 . @xmath107cos@xmath108 tests the dipole moment around the galactic center , @xmath107sin@xmath110 tests the quadrapole moment with respect to the galactic plane , @xmath111 tests the dipole moment around any point in the celestial sphere , and @xmath112 tests the quadrapole moment around any plane or two poles . the expected values for the four statistics assuming random isotropic distribution on the sphere are 0 , 0 , 3 and 5 respectively . the asymptotic distributions of these statistics are known . table 6b shows the results of this analysis for clusters i - iii , kindly calculated for us by michael briggs . no deviations from isotropy are found . the @xmath107cos@xmath108 and @xmath107sin@xmath113 values lie within one standard deviation of the expected value for a random distribution . the @xmath111 and @xmath112 values must be larger than the expected value to indicate anisotropy . the only such case , class ii with @xmath1147.32 , has a deviation with very low significance ( @xmath115 ) . we thus do not confirm belli s ( 1997 ) report of significant differences in spatial distributions of burst classes i and iii , although we did not specifically test the galactic latitude distribution .    in principle , the relative populations of the three classes may be an important constraint on astrophysical theory . we find that class i contains more than half of the bursts with the remainder divided between class ii and class iii ( table 6c ) . but we do not believe our analysis gives a precise census for two reasons . first , the exact assignments of individual bursts to clusters depend on the detailed assumptions of the clustering algorithms . for example , class ii is larger than class iii in the 5-dimensional nonparametric procedure but is smaller in the 3-dimensional model - based procedure . second , the numbers of weaker bursts in classes ii and iii are strongly dependent on the details of the batse instrument s burst triggering process which produces a complicated truncation bias for fainter bursts .    we look for structure within each of the clusters by computing correlation coefficients similar to those in  3 for the entire sample . results are given in table 7 . here we see a systematic difference between the two clustering methodologies : nonparametric average linkage clustering tends to give stronger correlations between the variables than the model - based clustering . for example , in the nonparametric analysis we find significant positive correlations between total fluence and hardness in classes i and ii , and a correlation between duration and fluence in class ii . however , we attach more credence to the model - based results for this purpose than the average linkage results because the former method is specifically designed to provide optimal estimates of the within - group covariances given the clustering model . the model - based results do not give strong evidence for any non - zero correlations between variables , suggesting that the partition into three clusters explains all of the correlation between variables in the full dataset . we thus find , using clustering and validation methods with different mathematical underpinnings , that three classes of grbs are present in our large subset of the third batse catalog . most of the structure can be found using three fundamental burst properties , duration / fluence / spectrum . the class properties and relation to previous research can be briefly summarized as follows :    class i : :    these long / bright / intermediate bursts correspond to the well - known    populous long - soft class of k93 and others . within this group , we do    not confirm a hardness - duration correlation reported by dezalay et al . ( 1996 ) and horack & hakkila ( 1997 ) . class ii : :    this short / faint / hard group corresponds to the short - hard burst type    of k93 and others . fluence - duration and fluence - hardness correlations    may tentatively be present within the class . note that while the mean    location of this type is consistent in the two clustering schemes , its    size and population ( e.g. 1/2 or 1/4 that of class i ) differs between    clustering algorithms . class iii : :    the discovery of this group with intermediate / intermediate / soft    properties is the principal result of this study . the group is easily    distinguished in the projections of figure 3 , but can also be    discerned in some panels of figure 1 . for example , it lies between    classes i and ii in the @xmath116 ,    @xmath117 and @xmath118 scatter    plots . in the univariate @xmath1 distribution shown in    figure 5 , class iii accounts for most , but not all , of the bursts in    the small peak around @xmath119 sec between the major    short and long duration peaks . it is possible that our class iii is    related to the class of no - high - energy ( nhe ) burst and peaks discussed    by pendleton et al . these bursts have unusually weak    @xmath12 emission , soft @xmath120 kev spectra , and    low @xmath22 . however , the nhe class does not appear to    exhibit a clear duration segregation from other bursts as we find for    class iii . class iii does not appear to be the third cluster found by    baumgart ( 1994 , see his table 3 ) , but the high dimensionality of his    analysis prevents a simple comparison with our low dimensionality    study . outlier : :    batse trigger event 2757 , burst 3b 940114 , is the outlier in the    nonparametric analysis of  4 and is clearly visible in many    projections in figures 1 and 3 . it has an exceedingly soft hardness ratio and short    burst duration . but examination of the original batse database shows    that the @xmath121 fluxes are very weak with large    measurement uncertainties . the published 3b catalog gives only an    upper limit to its total fluence and no estimate of its hardness ratio    ( meegan et al . the unusual properties of this burst are thus    illusory and are due to its very weak fluence . the multivariate analysis described here is not comprehensive and may not have uncovered all of the structure in the third batse catalog of bulk grb properties . our reduction of dimensionality may have been too severe , omitting , for example , the potentially important @xmath12 as a distinct variable ( pendleton et al . 1997 ; bagoly et al . 1997 ) . many methodological options were not exercised . for example , it would be valuable to repeatedly apply the @xmath68-means partitioning algorithm to the database under the assumption that three clusters are present ( see murtagh 1992 for an astronomical application of this method ) , check for skewness or kurtosis in the clusters , and undertake an oblique decision tree analysis to give analytical formulation to hyperplanes separating the clusters ( see white 1997 ) . codes for these and many other multivariate techniques are publicly available through the web metasite _ statcodes _ at www.astro.psu.edu/statcodes .    however , the efforts described here are far more capable of finding and quantifying clustering in the database than most previous analyses (  1 ) . previous studies have been based on qualitative rather than quantitative procedures for identifying structures , and provide no statistical validation of their claims . it is thus not surprising that we uncovered structure missed by previous researchers . in particular , our confidence in the presence of a third cluster , class iii , is strong . two completely independent mathematical procedures (  4 and  5 ) found very similar structure , each validated with high statistical confidence . it is possible that the clustering reported here is indeed present in the database , but does not have an astrophysical origin . the complex triggering mechanism of the batse instrument mechanism and biases in bulk property values at low signal - to - noise ratio are two problems that probably affect the multivariate structure . we have investigated one manifestation of the latter effect using @xmath18 and found no effect on our results . instrumental biases generally affect the number of bursts found in some regions of the multivariate hyperspace ( thereby biasing log @xmath36 - log @xmath122 distributions ) and may alter the location of clusters , but are unlikely to cause the appearance of clustering that is not present in the underlying population . nonetheless , since the batse instrument identifies bursts on three separate timescales , it is possible that the third cluster here is related to a selection effect associated with the batse triggering mechanisms . we conclude that the third batse catalog shows three statistically significant types of bursts ( duration / fluence / spectrum ) : class i grbs are long / bright / intermediate , class ii grbs are short / faint / hard , and class iii grbs are intermediate / intermediate / soft . unless the separation of class iii from the other types is due to some subtle batse instrumental effect , these types are likely to be real and their existence should be considered an important input into astrophysical theories for grbs . for example , the three types may reflect different types of external environments and internal shocks in relativistic fireball models ( mszros & rees 1993 ; panaitescu & mszros 1998 ) . note that statistical anlaysis is unable to determine whether burst types represent fundamentally different astrophysical processes or distinct conditions within a single astrophysical model . our results can be confirmed and extended in two fashions . first , the analysis described here can be validated with several hundred more bursts collected by batse since the september 1994 cutoff in the database used here . second , following baumgart ( 1994 ) , the dimensionality of the problem can be enlarged to include detailed characteristics of the burst temporal behaviors . burst smoothness @xmath83 peakiness , characteristic wavelet scales , spectral evolution , and other parameters can be included . with this enlarged database , one can perform both an unsupervised exploratory cluster analysis similar to that described here , and manova - type analyses that assume the existence of the three groups to determine whether the clusters have distinctive temporal properties . this work was supported by nasa grant nagw-2120 and nsf grant dms-9626189 funding astrostatistical research at penn state , and onr grants n-00014 - 96 - 1 - 0192 and n-00014 - 96 - 1 - 0330 funding statistical research at washington . sm would like to thank the departments of astronomy & astrophysics and statistics at penn state for hospitality while this work was underway . we are very grateful to jay norris ( nasa - gsfc ) for providing debiased duration values and michael briggs ( nasa - msfc ) for performing isotropy calculations . we also thank drs . briggs , norris and peter mszros ( penn state ) for sharing many valuable insights into grb astronomy . pendleton , g. n. , paciesas , w. s. , briggs , m. s. , preece , r. d. , mallozzi , r. s. , meegan , c. a. , horack , j. m. , fishman , g. j. , band , d. l. , matteson , j. l. , skelton , r. t. , hakkila , j. , ford , l. a. , kouveliotou , c. & koshut , t. m. 1997 , apj , 489 , 175                  lrc & & + & mean & s.d@xmath123 + log @xmath7 ( sec ) & 0.55 & 0.92 + log @xmath1 ( sec ) & 0.96 & 0.92 + log @xmath22 ( erg @xmath124 ) & -5.61 & 0.76 + log @xmath9 ( photons s@xmath125 @xmath126 & 0.16 & 0.45 + log @xmath25 & 0.25 & 0.33 + log @xmath26 & 0.48 & 0.30 +    lrrrrrr & & & & & & + & log @xmath7 & log @xmath1 & log @xmath22 & log @xmath9 & log @xmath25 & log @xmath26 + log @xmath7 & 1.00 & & & & & + log @xmath1 & 0.97 & 1.00 & & & & + log @xmath22 & 0.63 & 0.66 & 1.00 & & & + log @xmath9 & -0.01 & 0.04 & 0.59 & 1.00 & & + log @xmath25 & -0.36 & -0.36 & 0.02 & 0.24 & 1.00 & + log @xmath26 & -0.35 & -0.35 & -0.00 & 0.19 & 0.96 & 1.00 +    rcrrr & & & & + level & merger & members & @xmath40 & @xmath37 +   + 8 & 10 + 15 & 506 & 0.08 & 0.65 + 7 & 14 + 137 & 93 & 0.00 & 0.65 + 6 & 8 + 7 & 599 & 0.10 & 0.55 + 5 & 9 + 266 & 188 & 0.00 & 0.55 + 4 & 5 + 26 & 190 & 0.00 & 0.55 + 3 & 6 + 12 & 606 & 0.01 & 0.54 + 2 & 3 + 4 & 796 & 0.53 & 0.01 + 1 & 2 + 616 & 797 & 0.00 & 0.00 + & & & & +   + 6 & 15 + 21 & 107 & 0.01 & 0.70 + 5 & 10 + 8 & 486 & 0.01 & 0.69 + 4 & 7 + 20 & 203 & 0.01 & 0.68 + 3 & 6 + 5 & 593 & 0.10 & 0.58 + 2 & 3 + 4 & 796 & 0.58 & 0.00 + 1 & 2 + 616 & 797 & 0.01 & 0.00 +       lcrrr & & & & + variable & method@xmath127 & + & & & & +   + log @xmath7 & np & 1.13 @xmath128 0.44 & -0.80 @xmath128 0.41 & 0.33 @xmath128 0.27 + log @xmath1 & np & 1.55 @xmath128 0.40 & -0.42 @xmath128 0.44 & 0.71 @xmath128 0.32 + & mb & 1.22 @xmath128 0.39 & -0.91 @xmath128 0.35 & 0.29 @xmath128 0.41 + log @xmath22 & np & -5.21 @xmath128 0.59 & -6.37 @xmath128 0.57 & -6.11 @xmath128 0.37 + & mb & -5.13 @xmath128 0.58 & -6.46 @xmath128 0.54 & -5.93 @xmath128 0.47 + log @xmath25 & np & 0.19 @xmath128 0.27 & 0.51 @xmath128 0.27 & 0.08 @xmath128 0.40 + & mb & 0.21 @xmath128 0.26 & 0.52 @xmath128 0.28 & 0.16 @xmath128 0.40 + log @xmath26 & np & 0.43 @xmath128 0.23 & 0.70 @xmath128 0.26 & 0.35 @xmath128 0.39 + & & & & +   + @xmath107cos @xmath109@xmath4 & np & 0.015 & -0.041 & 0.010 + @xmath107sin@xmath129 & np & -0.012 & -0.025 & 0.028 + rayleigh - watson & np & 0.39 & 1.28 & 1.80 + bingham & np & 2.02 & 7.32 & 1.95 + & & & & +   + number & np & 486 & 203 & 107 + & mb & 426 & 170 & 201 + duration & & long & short & intermediate + fluence & & bright & faint & intermediate + spectrum & & intermediate & hard & soft +   +    lrrrrr & & & & & + & log @xmath7 & log @xmath1 & log @xmath22 & log @xmath25 & log @xmath26 +   +   + log @xmath7 & 1.00 & & & & + log @xmath1 & 0.88 & 1.00 & & & + log @xmath22 & 0.10 & 0.22 & 1.00 & & + log @xmath25 & -0.11 & -0.08 & 0.39 & 1.00 & + log @xmath26 & -0.11 & -0.08 & 0.38 & 0.97 & 1.00 + & & & & & +   +   + log @xmath7 & 1.00 & & & & + log @xmath1 & n / a & n / a & & & + log @xmath22 & n / a & 0.01 & 1.00 & & + log @xmath25 & n / a & -0.01 & 0.06 & 1.00 & + log @xmath26 & n / a & n / a & n / a & n / a & n / a + & & & & & +   +   + log @xmath7 & 1.00 & & & & + log @xmath1 & n / a & n / a & & & + log @xmath22 & n / a & -0.03 & 1.00 & & + log @xmath25 & n / a & -0.05 & 0.02 & 1.00 & + log @xmath26 & n / a & n / a & n / a & n / a & n / a + & & & & & +   +   + log @xmath7 & 1.00 & & & & + log @xmath1 & 0.86 & 1.00 & & & + log @xmath22 & 0.02 & 0.06 & 1.00 & & + log @xmath25 & -0.24 & -0.34 & -0.16 & 1.00 & + log @xmath26 & -0.22 & -0.32 & -0.22 & 0.95 & 1.00 + & & & & & +   +   + log @xmath7 & 1.00 & & & & + log @xmath1 & n / a & n / a & & & + log @xmath22 & n / a & 0.03 & 1.00 & & + log @xmath25 & n / a & -0.01 & 0.07 & 1.00 & + log @xmath26 & n / a & n / a & n / a & n / a & n / a +
Abstract:

Answer: a multivariate analysis of gamma - ray burst ( grb ) bulk properties is presented to discriminate between distinct classes of grbs . several variables representing burst duration , fluence and spectral hardness are considered . two multivariate clustering procedures are used on a sample of 797 bursts from the third batse catalog : a nonparametric average linkage hierarchical agglomerative clustering procedure validated with wilks @xmath0 and other manova tests ; and a parametric maximum likelihood model - based clustering procedure assuming multinormal populations calculated with the em algorithm and validated with the bayesian information criterion . the two methods yield very similar results . the batse grb population consists of three classes with the following duration / fluence / spectrum bulk properties : class i with long / bright / intermediate bursts , class ii with short / hard / faint bursts , and class iii with intermediate / intermediate / soft bursts . one outlier with poor data is also present . classes i and ii correspond to those reported by kouveliotou et al . ( 1993 ) , but class iii is clearly defined here for the first time .