Studying the GRBs' gamma-ray spectra may reveal some physical information of Gamma-ray Bursts. The Fermi satellite observed more than two thousand GRBs. The FERMIGBRST catalog contains GRB parameters (peak energy, spectral indices, intensity) estimated for both the total emission (fluence), and the emission during the interval of the peak flux. We found a relationship with linear discriminant analysis between the spectral categories and the model independent physical data. We compared the Swift and Fermi spectral types. We found a connection between the Fermi fluence spectra and the Swift spectra but the result of the peak flux spectra can be disputable. We found that those GRBs which were observed by both Swift and Fermi can similarly discriminate as the complete Fermi sample. We concluded that the common observation probably did not mean any trace of selection effects in the spectral behavior of GRBs.

The afterglow of a gamma ray burst (GRB) can give us valuable insight into the properties of its host galaxy. To correctly interpret the spectra of the afterglow we need to have a good understanding of the foreground interstellar medium (ISM) in our own Galaxy. The common practice to correct for the foreground is to use neutral hydrogen (HI) data from the Leiden/Argentina/Bonn (LAB) survey. However, the poor spatial resolution of the single dish data may have a significant effect on the derived column densities. To investigate this, we present new high-resolution HI observations with the Australia Telescope Compact Array (ATCA) towards 4 GRBs. We combine the interferometric ATCA data with single dish data from the Galactic All Sky Survey (GASS) and derive new Galactic HI column densities towards the GRBs. We use these new foreground column densities to fit the Swift XRT X-ray spectra and calculate new intrinsic hydrogen column density values for the GRB host galaxies. We find that the new ATCA data shows higher Galactic HI column densities compared to the previous single dish data, which results in lower intrinsic column densities for the hosts. We investigate the line of sight optical depth near the GRBs and find that it may not be negligible towards one of the GRBs, which indicates that the intrinsic hydrogen column density of its host galaxy may be even lower. In addition, we compare our results to column densities derived from far-infrared data and find a reasonable agreement with the HI data.

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

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA).

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

Statistical studies of gamma-ray burst (GRB) spectra may result in important information on the physics of GRBs. The Fermi GBM catalog contains GRB parameters (peak energy, spectral indices, intensity) estimated fitting the gamma-ray SED of the total emission (fluence, flnc), and during the time of the peak flux pflx. Using contingency tables we studied the relationship of the models best fitting pflx and flnc time intervals. Our analysis revealed an ordering of the spectra into a power law - Comptonized - smoothly broken power law - Band series. This result was further supported by a correspondence analysis (CA) of the pflx and flnc spectra categorical variables. We performed a linear discriminant analysis (LDA) to find a relationship between categorical (spectral) and model independent physical data. LDA resulted in highly significant physical differences among the spectral types, that is more pronounced in the case of the pflx spectra, than for the flnc spectra. We interpreted this difference as caused by the temporal variation of the spectrum during the outburst. This spectral variability is confirmed by the differences in the low energy spectral index and peak energy, between the pflx and flnc spectra. We found that the synchrotron radiation is significant in GBM spectra. The mean low energy spectral index is close to the canonical value of {\alpha} = -2/3 during the peak flux. However, $\alpha$ is ~ -0.9 for the spectra of the fluences. We interpret this difference as showing that the effect of cooling is important only for the fluence spectra.

A significant number of the parameters of a gamma-ray burst (GRB) and its host galaxy are calculated from the afterglow. There are various methods obtaining extinction values for the necessary correction for galactic foreground. These are: galaxy counts, from HI 21 cm surveys, from spectroscopic measurements and colors of nearby Galactic stars, or using extinction maps calculated from infrared surveys towards the GRB. We demonstrate that AKARI Far-Infrared Surveyor sky surface brightness maps are useful uncovering the fine structure of the galactic foreground of GRBs. Galactic cirrus structures of a number of GRBs are calculated with a 2 arcminute resolution, and the results are compared to that of other methods.

The Fermi collaboration identified a possible electromagnetic counterpart of the gravitational wave event of September 14, 2015. Our goal is to provide an unsupervised data analysis algorithm to identify similar events in Fermi's Gamma-ray Burst Monitor CTTE data stream. We are looking for signals that are typically weak. Therefore, they can only be found by a careful analysis of count rates of all detectors and energy channels simultaneously. Our Automatized Detector Weight Optimization (ADWO) method consists of a search for the signal, and a test of its significance. We developed ADWO, a virtual detector analysis tool for multi-channel multi-detector signals, and performed successful searches for short transients in the data-streams. We have identified GRB150522B, as well as possible electromagnetic candidates of the transients GW150914 and LVT151012. ADWO is an independently developed, unsupervised data analysis tool that only relies on the raw data of the Fermi satellite. It can therefore provide a strong, independent test to any electromagnetic signal accompanying future gravitational wave observations.

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

Short and long bursts were identified by the BATSE team in the early 90s. A decade ago there were some suggestions about the intermediate duration type of bursts. We are going to summarize recent analyses of the duration distributions of the Beppo-Sax and Swift data. Our conclusion is all the three satellites (CGRO, Swift, Beppo-Sax) can see the third type of the GRBs. The properties of the group members are very similar in the different data sets.

We studied the unbiased optical brightness distribution which was calculated from the survival analysis of host galaxies and its relationship with the Swift GRB data of the host galaxies observed by the Keck telescopes. Based on the sample obtained from merging the Swift GRB table and the Keck optical data we also studied the dependence of this distribution on the data of the GRBs. Finally, we compared the HGs distribution with standard galaxies distribution which is in the DEEP2 galaxies catalog.

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

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

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

The statistical tests - done by the authors - are surveyed, which verify the null-hypothesis of the intrinsic randomness in the angular distribution of gamma-ray bursts collected at BATSE Catalog. The tests use the counts-in-cells method, an analysis of spherical harmonics, a test based on the two-point correlation function and a method based on multiscale methods. The tests suggest that the intermediate subclass of gamma-ray bursts are distributed anisotropically.

One of the most important task of the Gamma-Ray Burst field is the classification of the bursts. Many researches have proven the existence of the third kind (intermediate duration) of GRBs in the BATSE data. Recent works have analyzed BeppoSax and Swift observations and can also identify three types of GRBs in the data sets. However, the class memberships are probabilistic we have enough observed redshifts to calculate the redshift and spatial distribution of the intermediate GRBs. They are significantly farther than the short bursts and seems to be closer than the long ones.

Gamma-ray bursts can be divided into three groups: short-, intermediate and long-duration bursts. While the progenitors of the short and long ones are relatively known, the progenitor objects of the intermediate-duration bursts are generally unknown, however, recent statistical studies suggest, that they should be related to the long-duration bursts. In this work we investigate whether there is a differ- ence between the global parameters of the X-ray flashes and intermediate-duration group of gamma-ray bursts. The statistical tests do not show any significant discrepancy regarding most of the parameters, except the BAT photon index, which is only a consequence of the defintion of the X-ray flashes.

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

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

We present a method for determining the background of Fermi GBM GRBs using the satellite positional information and a physical model. Since the polynomial fitting method typically used for GRBs is generally only indicative of the background over relatively short timescales, this method is particularly useful in the cases of long GRBs or those which have Autonomous Repoint Request (ARR) and a background with much variability on short timescales. We give a Direction Dependent Background Fitting (DDBF) method for separating the motion effects from the real data and calculate the duration (T90 and T50, as well as confidence intervals) of the nine example bursts, from which two resulted an ARR. We also summarize the features of our method and compare it qualitatively with the official GBM Catalogue. Our background filtering method uses a model based on the physical information of the satellite position. Therefore, it has many advantages compared to previous methods. It can fit long background intervals, remove all the features caused by the rocking behaviour of the satellite, and search for long emissions or not-triggered events. Furthermore, many part of the fitting have now been automatised, and the method have been shown to work for both Sky Survey mode and ARR mode data. Future work will provide a burst catalogue with DDBF.

We investigated the main prompt and afterglow emission parameters of gamma-ray bursts detected by the Burst Alert Telescope (BAT) and X-Ray Telescope installed on the Swift satellite. Our aim was to look for differences or connections between the different types of gamma-ray bursts, so we compared the BAT fluences, 1-sec peak photon fluxes, photon indices, XRT early fluxes, initial temporal decay and spectral indices. We found that there might be a connection between the XRT initial decay index and XRT early flux/BAT photon index. Using statistical tools we also determined that beside the duration and hardness ratios, the means of the \gamma- and X-ray--fluences and the \gamma-ray photon index differ significantly between the three types of bursts.

We studied the distribution of the GRBs, observed by the Fermi satellite, in the multidimensional parameter space consisting of the duration, Fluence, Peak flux and Peak energy (if it was available). About 10% of the Fermi bursts was observed also by the Swift satellite. We did not find significant differences between the Peak flux and Peak energy of GRBs observed and not observed also by the Swift satellite. In contrast, those GRBs detected also by the Swift satellite had significantly greater Fluence and duration. We did a similar study for the GRBs detected by the Swift satellite. About 30% percent of these bursts was also measured by the Fermi satellite. We found a significant difference in the Fluence, Peak flux and Photon index but none in duration. These differences may be accounted for the different construction and observing strategy of the Fermi and Swift satellites.

The Fermi Gamma-ray Burst Monitor (GBM) detects gamma-rays in the energy range 8 keV - 40 MeV. We developed a new background fitting process of these data, based on the motion of the satellite. Here we summarize this method, called Direction Dependent Background Fitting (DDBF), regarding the GBM triggered catalog. We also give some preliminary results and compare the duration parameters with the 2-years Fermi Burst Catalog.

The Swift satellite made a real break through with measuring simultaneously the gamma X-ray and optical data of GRBs, effectively. Although, the satellite measures the gamma, X-ray and optical properties almost in the same time a significant fractions of GRBs remain undetected in the optical domain. In a large number of cases only an upper bound is obtained. Survival analysis is a tool for studying samples where a part of the cases has only an upper (lower) limit. The obtained survival function may depend on some other variables. The Cox regression is a way to study these dependencies. We studied the dependence of the optical brightness (obtained by the UVOT) on the gamma and X-ray properties, measured by the BAT and XRT on board of the Swift satellite. We showed that the gamma peak flux has the greatest impact on the afterglow's optical brightness while the gamma photon index and the X-ray flux do not. This effect probably originates in the energetics of the jet launched from the central engine of the GRB which triggers the afterglow.

Since the early 1990s Gamma Ray Bursts have been accepted to be of extra-galactic origin due to the isotropic distribution observed by BATSE and the redshifts observed via absorption line spectroscopy. Nevertheless, upon further examination at least one case turned out to be of galactic origin. This particular event presented a Fast Rise, Exponential Decay (FRED) structure which leads us to believe that other FRED sources might also be Galactic. This study was set out to estimate the most probable degree of contamination by galactic sources that certain samples of FREDs have. In order to quantify the degree of anisotropy the average dipolar and quadripolar moments of each sample of GRBs with respect to the galactic plane were calculated. This was then compared to the probability distribution of simulated samples comprised of a combination of isotropically generated sources and galactic sources. We observe that the dipolar and quadripolar moments of the selected subsamples of FREDs are found more than two standard deviations outside those of random isotropically generated samples.The most probable degree of contamination by galactic sources for the FRED GRBs of the Swift catalog detected until February 2011 that do not have a known redshift is about 21 out of 77 sources which is roughly equal to 27%. Furthermore we observe, that by removing from this sample those bursts that may have any type of indirect redshift indicator and multiple peaks gives the most probable contamination increases up to 34% (17 out of 49 sources). It is probable that a high degree of contamination by galactic sources occurs among the single peak FREDs observed by Swift.

In the case of Gamma-Ray Bursts with measured redshift, we can calculate the k-correction to get the fluence and energy that were actually produced in the comoving system of the GRB. To achieve this we have to use well-fitted parameters of a GRB spectrum, available in the GCN database. The output of the calculations is the comoving isotropic energy E_iso, but this is not the endpoint: this data can be useful for estimating the {\Omega}M parameter of the Universe and for making a GRB Hubble diagram using Amati's relation.

The Fermi satellite has a particular motion during its flight which enables it to catch the gamma-ray bursts mostly well. The side-effect of this favourable feature is that the lightcurves of the GBM detectors are stressed by rapidly and extremely varying background. Before this data is processed, it needs to be separated from the background. The commonly used methods were useless for most cases of Fermi, so we developed a new technique based on the motion and orientation of the satellite. The background-free lightcurve can be used to perform statistical surveys, hence we showed the efficiency of our background-filtering method presenting a statistical analysis known from the literature.

Swift satellite measurements contributed substantially to the gamma-ray burst (GRB) redshift observations through fast slewing to the source of the GRBs. Still, a large number of bursts are without redshift. We study the celestial distribution of bursts with various methods and compare them to a random catalog using Monte-Carlo simulations. We find an anisotropy in the distribution of the intermediate class of bursts and find that the short and long population are distributed isotropically.

Gamma-ray bursts are the most luminous events in the Universe. Going beyond the short-long classification scheme we work in the context of three burst populations with the third group of intermediate duration and softest spectrum. We are looking for physical properties which discriminate the intermediate duration bursts from the other two classes. We use maximum likelihood fits to establish group memberships in the duration-hardness plane. To confirm these results we also use k-means and hierarchical clustering. We use Monte-Carlo simulations to test the significance of the existence of the intermediate group and we find it with 99.8% probability. The intermediate duration population has a significantly lower peak-flux (with 99.94% significance). Also, long bursts with measured redshift have higher peak-fluxes (with 98.6% significance) than long bursts without measured redshifts. As the third group is the softest, we argue that we have {related} them with X-ray flashes among the gamma-ray bursts. We give a new, probabilistic definition for this class of events.

With the early afterglow localizations of gamma-ray burst positions made by Swift, the clear delimitation of the prompt phase and the afterglow is not so obvious any more. It is important to see weather the two phases have the same origin or they stem from different parts of the progenitor system. We will combine the two kinds of gamma-ray burst data from the Swift-XRT instrument (windowed timing and photon counting modes) and from BAT. A thorough desription of the applied procedure is given. We apply various binning techniques to the different data: Bayes blocks, exponential binning and signal-to-noise type of binning. We present a handful of flux curves and some possible applications.

Gamma-ray bursts are usually classified through their high-energy emission into short-duration and long-duration bursts, which presumably reflect two different types of progenitors. However, it has been shown on statistical grounds that a third, intermediate population is needed in this classification scheme, although an extensive study of the properties of this class has so far not been done. The large amount of follow-up studies generated during the Swift era allows us to have a suficient sample to attempt a study of this third population through the properties of their prompt emission and their afterglows. Our study is focused on a sample of GRBs observed by Swift during its first four years of operation. The sample contains those bursts with measured redshift since this allows us to derive intrinsic properties. Intermediate bursts are less energetic and have dimmer afterglows than long GRBs, especially when considering the X-ray light curves, which are on average one order of magnitude fainter than long bursts. There is a less significant trend in the redshift distribution that places intermediate bursts closer than long bursts. Except for this, intermediate bursts show similar properties to long bursts. In particular, they follow the Epeak vs. Eiso correlation and have, on average, positive spectral lags with a distribution similar to that of long bursts. Like long GRBs, they normally have an associated supernova, although some intermediate bursts have shown no supernova component. This study shows that intermediate bursts are different from short bursts and, in spite of sharing many properties with long bursts, there are some differences between them as well. We suggest that the physical difference between intermediate and long bursts could be that for the first the ejecta are thin shells while for the latter they are thick shells.

Recently the anisotropy of the short gamma-ray bursts detected by BATSE was announced (Vavrek et al. 2008). The impact of this discovery on cosmology is discussed. It is shown that the anisotropy found may cause the breakdown of the cosmological principle.

Earlier classification analyses found three types of gamma-ray bursts (short, long and intermediate in duration) in the BATSE sample. Recent works have shown that these three groups are also present in the RHESSI and the BeppoSAX databases. The duration distribution analysis of the bursts observed by the Swift satellite also favors the three-component model. In this paper, we extend the analysis of the Swift data with spectral information. We show, using the spectral hardness and the duration simultaneously, that the maximum likelihood method favors the three-component against the two-component model. The likelihood also shows that a fourth component is not needed.

The Gamma Ray Inverse Problem is discussed. Four methods of spectral deconvolution are presented here and applied to the BATSE's MER data type. We compare these to the Band spectra.

As GRBs trace the high-z Universe, there is a non-negligible probability of a lensing effect being imprinted on the lightcurves of the bursts. We propose to search for lensed candidates with a cross-correlation method, by looking at bursts days to years apart coming from the same part of the sky. We look for similarities and hypothesize a Singular Isothermal Sphere (SIS) model for the lens. A lensed pair would enable us to constrain the mass of the lensing object. Our search did not reveal any gravitationally lensed events.

The power spectra of 457 short BATSE bursts were analyzed, focusing on the 64 ms lightcurves' tails in the low energy bands. Using MC simulations, 22 GRBs were identified with unusually high harmonic power above 0.03 Hz. The sky distribution of these bursts shows an extraordinarily strong dipole moment with a 99.994% significance.

GRB 080319B is one of the brightest and most extensively sampled bursts. It has good coverage at many wavelengths. Here we present the optical observations of the Palomar 60 inch telescope, which spans a long time interval after the burst. We augment the optical dataset with freely available Swift BAT and XRT observations reduced by us. We also compare our conclusions with the published parameters from the rich literature about this burst.

Two classes of gamma-ray bursts were identified in the BATSE catalogs characterized by their durations. There were also some indications for the existence of a third type of gamma-ray bursts. Swift satellite detectors have different spectral sensitivity than pre-Swift ones for GRBs. Therefore in this paper we analyze the bursts' duration distribution and also the duration-hardness bivariate distribution, published in The First BAT Catalog. Similarly to the BATSE data, to explain the BAT GRBs' duration distribution three components are needed. Although, the relative frequencies of the groups are different than they were in the BATSE GRB sample, the difference in the instrument spectral sensitivities can explain this bias in a natural way. This means theoretical models may have to explain three different type of gamma-ray bursts.

Two classes of gamma-ray bursts have been identified in the BATSE catalogs characterized by durations shorter and longer than about 2 seconds. There are, however, some indications for the existence of a third type of burst. Swift satellite detectors have different spectral sensitivity than pre-Swift ones for gamma-ray bursts. Therefore it is worth to reanalyze the durations and their distribution and also the classification of GRBs. Using The First BAT Catalog the maximum likelihood estimation was used to analyzed the duration distribution of GRBs. The three log-normal fit is significantly (99.54 % probability) better than the two for the duration distribution. Monte-Carlo simulations also confirm this probability (99.2 %).

The measured redshifts of gamma-ray bursts (GRBs), which were first detected by the Swift satellite, seem to be bigger on average than the redshifts of GRBs detected by other satellites. We analyzed the redshift distribution of GRBs triggered and observed by different satellites (Swift, HETE2, BeppoSax, Ulyssses). After considering the possible biases {significant difference was found at the p=95.70% level in the redshift distributions of GRBs measured by HETE and the Swift.

After the discovery of the anisotropy in the sky-distribution of intermediate gamma-ray bursts recently also the distribution of the short gamma-ray bursts is proven to be anisotropic. The impact of these behaviors on the validity of the cosmological principle is shortly discussed.

In a systematic search of the OTs at GRBs the Swift satellite determined only an upper limit of the apparent brightness in a significant fraction of cases. Combining these upper limits with the really measured OT brightness we obtained a sample well suited to survival analysis. Performing a Kaplan-Meier product limit estimation we obtained an unbiased cumulative distribution of the V visual brightness. The lg(N(V)) logarithmic cumulative distribution can be well fitted with a linear function of V in the form of lg(N(V))= 0.234 V + const. We studied the dependence of V on the gamma ray properties of the bursts. We tested the dependence on the fluence, T90 duration and peak flux. We found a dependence on the peak flux on the 99.7% significance level.

Using the discriminant analysis of the multivariate statistical analysis we compared the distribution of the physical quantities of the optically dark and bright GRBs, detected by the BAT and XRT on board of the Swift Satellite. We found that the GRBs having detected optical transients (OT) have systematically higher peak fluxes and lower HI column densities than those without OT.

We presented the results of several statistical tests of the randomness in the angular sky-distribution of gamma-ray bursts in BATSE Catalog. Thirteen different tests were presented based on Voronoi tesselation, Minimal spanning tree and Multifractal spectrum for five classes (short1, short2, intermediate, long1, long2) of gamma-ray bursts, separately. The long1 and long2 classes are distributed randomly. The intermediate subclass, in accordance with the earlier results of the authors, is distributed non-randomly. Concerning the short subclass earlier statistical tests also suggested some departure from the random distribution, but not on a high enough confidence level. The new tests presented in this article suggest also non-randomness here.

A sample of 197 long BATSE GRBs is studied statistically. In the sample 11 variables, describing for any burst the time behavior of the spectra and other quantities, are collected. The application of the factor analysis on this sample shows that five factors describe the sample satisfactorily. Both the pseudo-redshifts coming from the variability and the Amati-relation in its original form are disfavored.

In this paper we make an attempt to combine the two kinds of data from the Swift-XRT instrument (windowed timing and photon counting modes) and the from BAT. A thorough desription of the applied procedure will be given. We apply various binning techniques to the different data: Bayes blocks, exponential binning and signal-to-noise type of binning. We present a handful of lightcurves and some possible applications.

The measured redshifts of gamma-ray bursts (GRBs), which were first detected by the Swift satellite, seem to be bigger on average than the redshifts of GRBs detected by other satellites. We analyzed the redshift distribution of GRBs triggered and observed by different satellites (Swift[1], HETE2[2], BeppoSax, Ulyssses). After considering the possible biases significant difference was found at the p = 95.70% level in the redshift distributions of GRBs measured by HETE and the Swift.

We study statistically 197 long gamma-ray bursts, detected and measured in detail by the BATSE instrument of the Compton Gamma-Ray Observatory. In the sample 10 variables, describing for any burst the time behavior of the spectra and other quantities, are collected. The factor analysis method is used to find the latent random variables describing the temporal and spectral properties of GRBs. The application of this particular method to this sample indicates that five factors and the $\REpk$ spectral variable (the ratio of peak energies in the spectrum) describe the sample satisfactorily. Both the pseudo-redshifts inferred from the variability, and the Amati-relation in its original form, are disfavored.

Two classes of gamma-ray bursts have been identified in the BATSE catalogs characterized by durations shorter and longer than about 2 seconds. There are, however, some indications for the existence of a third class. Swift satellite detectors have different spectral sensitivity than pre-Swift ones for gamma-ray bursts. Therefore, it is worth to reanalyze the durations and their distribution. We analyze, the maximum likelihood estimation, the bursts duration distribution, published in The First BAT Catalog, whether it contains two, three or more groups. The three log-normal fit is significantly (99.54% probability) better than the two for the duration distribution. Monte-Carlo simulations also confirm this probability (99.2%). Similarly, in previous results we found that the fourth component is not needed. The relative frequencies of the distribution of the groups are 7% short 35% intermediate and 58% long. Although the relative frequencies of the groups are different than in the BATSE GRB sample, the difference in the instrument spectral sensitivities can explain this bias on a natural way. This means theoretical models may be needed to explain three different type of gamma-ray bursts.

We studied the complete randomness of the angular distribution of gamma-ray bursts (GRBs) detected by BATSE. Since GRBs seem to be a mixture of objects of different physical nature we divided the BATSE sample into 5 subsamples (short1, short2, intermediate, long1, long2) based on their durations and peak fluxes and studied the angular distributions separately. We used three methods, Voronoi tesselation, minimal spanning tree and multifractal spectra to search for non-randomness in the subsamples. To investigate the eventual non-randomness in the subsamples we defined 13 test-variables (9 from the Voronoi tesselation, 3 from the minimal spanning tree and one from the multifractal spectrum). Assuming that the point patterns obtained from the BATSE subsamples are fully random we made Monte Carlo simulations taking into account the BATSE's sky-exposure function. The MC simulations enabled us to test the null hypothesis i.e. that the angular distributions are fully random. We tested the randomness by binomial test and introducing squared Euclidean distances in the parameter space of the test-variables. We concluded that the short1, short2 groups deviate significantly (99.90%, 99.98%) from the fully randomness in the distribution of the squared Euclidean distances but it is not the case at the long samples. At the intermediate group the squared Euclidean distances also give significant deviation (98.51%).

Principal component analysis is a statistical method, which lowers the number of important variables in a data set. The use of this method for the bursts' spectra and afterglows is discussed in this paper. The analysis indicates that three principal components are enough among the eight ones to describe the variablity of the data. The correlation between spectral index alpha and the redshift suggests that the thermal emission component becomes more dominant at larger redshifts.

The low energy spectra of some gamma-ray bursts' show excess components beside the power-law dependence. The consequences of such a feature allows to estimate the gamma photometric redshift of the long gamma-ray bursts in the BATSE Catalog. There is good correlation between the measured optical and the estimated gamma photometric redshifts. The estimated redshift values for the long bright gamma-ray bursts are up to z=4, while for the the faint long bursts - which should be up to z=20 - the redshifts cannot be determined unambiguously with this method. The redshift distribution of all the gamma-ray bursts with known optical redshift agrees quite well with the BATSE based gamma photometric redshift distribution.

We provided separate bivariate log-normal distribution fits to the BATSE short and long burst samples using the durations and fluences. We show that these fits present an evidence for a power-law dependence between the fluence and the duration, with a statistically significant different index for the long and short groups. We argue that the effect is probably real, and the two subgroups are different physical phenomena. This may provide a potentially useful constraint for models of long and short bursts.

The origin of dark bursts - i.e. that have no observed afterglows in X-ray, optical/NIR and radio ranges - is unclear yet. Different possibilities - instrumental biases, very high redshifts, extinction in the host galaxies - are discussed and shown to be important. On the other hand, the dark bursts should not form a new subgroup of long gamma-ray bursts themselves.

Gamma-ray bursts can be divided into three groups ("short", "intermediate", "long") with respect to their durations. The third type of gamma-ray bursts - as known - has the intermediate duration. We show that the intermediate group is the softest one. An anticorrelation between the hardness and the duration is found for this subclass in contrast to the short and long groups.

We studied the statistical properties of the luminosity distance and internal extinction data of type Ia supernovae in the lists published by Tonry et al. (2003) and Barris et al. (2004). After selecting the luminosity distance in an empty Universe as a reference level we divided the sample into low $z<0.25$ and high $z \ge 0.25$ parts. We further divided these subsamples by the median of the internal extinction. Performing sign tests using the standardized residuals between the estimated logarithmic luminosity distances and those of an empty universe, on the four subsamples separately, we recognized that the residuals were distributed symmetrically in the low redshift region, independently from the internal extinction. On the contrary, the low extinction part of the data of $z \ge 0.25$ clearly showed an excess of the points with respect to an empty Universe which was not the case in the high extinction region. This diversity pointed to an interrelation between the estimated luminosity distance and internal extinction. To characterize quantitatively this interrelation we introduced a hidden variable making use of the technics of factor analysis. After subtracting that part of the residual which was explained by the hiddenmaking use of the technics of factor analysis. After subtracting that part of the residual which was explained by the hidden variable we obtained luminosity distances which were already free from interrelation with internal extinction. Fitting the corrected luminosity distances with cosmological models we concluded that the SN Ia data alone did not exclude the possibility of the $\Lambda=0$ solution.

The redshift distribution of gamma-ray bursts collected in the BATSE Catalog is compared with the star formation rate. We aim to clarify the accordance between them. We also study the case of comoving number density of bursts monotonously increasing up to redshift 6-20. A method independent of the models of the gamma-ray bursts is used. The short and the long subgroups are studied separately. The redshift distribution of the long bursts may be proportional to the star formation rate. For the short bursts this can also happen, but the proportionality is less evident. For the long bursts the monotonously increasing scenario is also less probable but still can occur. For the short bursts this alternative seems to be excluded.

Until 6 October 2005 sixteen redshifts have been measured of long gamma-ray bursts discovered by the Swift satellite. Further 45 redshifts have been measured of the long gamma-ray bursts discovered by other satellites. Here we perform five statistical tests comparing the redshift distributions of these two samples assuming - as the null hypothesis - identical distribution for the two samples. Three tests (Student's $t$-test, Mann-Whitney test, Kolmogorov-Smirnov test) reject the null hypothesis on the significance levels between 97.19 and 98.55%. Two different comparisons of the medians show extreme $(99.78-99.99994)$% significance levels of rejection. This means that the redshifts of the Swift sample and the redshifts of the non-Swift sample are distributed differently - in the Swift sample the redshifts are on average larger. This statistical result suggests that the long GRBs should on average be at the higher redshifts of the Swift sample.

Gamma-ray bursts can be divided into three groups ("short", "intermediate", "long") with respect to their durations. This classification is somewhat imprecise, since the subgroup of intermediate duration has an admixture of both short and long bursts. In this paper a physically more reasonable definition of the intermediate group is presented, using also the hardnesses of the bursts. It is shown again that the existence of the three groups is real, no further groups are needed. The intermediate group is the softest one. From this new definition it follows that 11% of all bursts belong to this group. An anticorrelation between the hardness and the duration is found for this subclass in contrast to the short and long groups. Despite this difference it is not clear yet whether this group represents a physically different phenomenon.

The gamma-ray burst's lightcurves have been analyzed using a special wavelet transformation. The applied wavelet base is based on a typical Fast Rise-Exponential Decay (FRED) pulse. The shape of the wavelet coefficients' total distribution is determined on the observational frequency grid. Our analysis indicates that the pulses in the long bursts' high energy channel lightcurves are more FRED-like than the lower ones, independently from the actual physical time-scale.

We present an analysis indicating that there is a correlation between the fluences and the durations of gamma-ray bursts, and provide arguments that this reflects a correlation between the total emitted energies and the intrinsic durations. For the short (long) bursts the total emitted energies are roughly proportional to the first (second) power of the intrinsic duration. This difference in the energy-duration relationship is statistically significant, and may provide an interesting constraint on models aiming to explain the short and long gamma-ray bursts.

It is shown that in the duration-hardness plane the GRBs of the third intermediate subgroup are well defined. Their durations are intermediate (i.e. roughly between 2 and 10 seconds), but their hardnesses are the lowest. They are even softer than the long bursts.

We analyze the spectral lags of a sample of bright gamma-ray burst pulses observed by CGRO BATSE and compare these with the results of high-resolution spectroscopical investigations. We find that pulses with hard spectra have the largest lags, and that there is a similar, but weaker correlation between hardness-intensity correlation index, eta, and lag. We also find that the lags differ considerably between pulses within a burst. Furthermore, the peak energy mainly decreases with increasing lag. Assuming a lag-luminosity relation as suggested by Norris et al., there will thus be a positive luminosity--peak-energy correlation. We also find that the hardness ratio, of the total flux in two channels, only weakly correlates with the spectral evolution parameters. These results are consistent with those found in the analytical and numerical analysis in Paper I. Finally, we find that for these bursts, dominated by a single pulse, there is a correlation between the observed energy-flux, F, and the inverse of the lag, t_lag: F propto t_lag^{-1}. We interpret this flux-lag relation found as a consequence of the lag-luminosity relation and that these bursts have to be relatively narrowly distributed in z. However, they still have to, mainly, lie beyond z ~ 0.01, since they do not coincide with the local super-cluster of galaxies. We discuss the observed correlations within the collapsar model, in which the collimation of the outflow varies. Both the thermal photospheric emission as well as non-thermal, optically-thin synchrotron emission should be important.

We argue that the distributions of both the intrinsic fluence and the intrinsic duration of the gamma-ray emission in gamma-ray bursts from the BATSE sample are well represented by log-normal distributions, in which the intrinsic dispersion is much larger than the cosmological time dilatation and redshift effects. We perform separate bivariate log-normal distribution fits to the BATSE short and long burst samples. The bivariate log-normal behaviour results in an ellipsoidal distribution, whose major axis determines an overall statistical relation between the fluence and the duration. We show that this fit provides evidence for a power-law dependence between the fluence and the duration, with a statistically significant different index for the long and short groups. We discuss possible biases, which might affect this result, and argue that the effect is probably real. This may provide a potentially useful constraint for models of long and short bursts.

It is known that the soft tail of the gamma-ray bursts' spectra show excesses from the exact power-law dependence. In this article we show that this departure can be detected in the peak flux ratios of different BATSE DISCSC energy channels. This effect allows to estimate the redshift of the bright long gamma-ray bursts in the BATSE Catalog. A verification of these redshifts is obtained for the 8 GRB which have both BATSE DISCSC data and measured optical spectroscopic redshifts. There is good correlation between the measured and esti redshifts, and the average error is $\Delta z \approx 0.33$. The method is similar to the photometric redshift estimation of galaxies in the optical range, hence it can be called as "gamma photometric redshift estimation". The estimated redshifts for the long bright gamma-ray bursts are up to $z \simeq 4$. For the the faint long bursts - which should be up to $z \simeq 20$ - the redshifts cannot be determined unambiguously with this method.

We argue that the distributions of both the intrinsic fluence and the intrinsic duration of the gamma-ray emission in gamma-ray bursts from the BATSE sample are well represented by log-normal distributions, in which the intrinsic dispersion is much larger than the cosmological time dilatation and redshift effects. We perform separate bivariate log-normal distribution fits to the BATSE short and long burst samples. The bivariate log-normal behavior results in an ellipsoidal distribution, whose major axis determines an overall statistical relation between the fluence and the duration. We show that this fit provides evidence for a power-law dependence between the fluence and the duration, with a statistically significant different index for the long and short groups. We discuss possible biases which might affect this result, and argue that the effect is probably real. This may provide a potentially useful constraint for models of long and short bursts.

In the article a test is developed, which allows to test the null-hypothesis of the intrinsic randomness in the angular distribution of gamma-ray bursts collected at the Current BATSE Catalog. The method is a modified version of the well-known counts-in-cells test, and fully eliminates the non-uniform sky-exposure function of BATSE instrument. Applying this method to the case of all gamma-ray bursts no intrinsic non-randomness is found. The test also did not find intrinsic non-randomnesses for the short and long gamma-ray bursts, respectively. On the other hand, using the method to the new intermediate subclass of gamma-ray bursts, the null-hypothesis of the intrinsic randomness for 181 intermediate gamma-ray bursts is rejected on the 96.4% confidence level. Taking 92 dimmer bursts from this subclass itself, we obtain the surprising result: This "dim" subclass of the intermediate subclass has an intrinsic non-randomness on the 99.3% confidence level. On the other hand, the 89 "bright" bursts show no intrinsic non-randomness.

The counts-in-cells and the two-point angular correlation function method are used to test the randomnesses in the angular distributions of both the all gamma-ray bursts collected at BATSE Catalog, and also their three subclasses ("short", "intermediate", "long"). The methods elimate the non-zero sky-exposure function of BATSE instrument. Both tests suggest intrinsic non-randomnesses for the intermediate subclass; for the remaining three cases only the correlation function method. The confidence levels are between 95% and 99.9%. Separating the GRBs into two parts ("dim half" and "bright half", respectively) we obtain the result that the "dim" half shows a non-randomness on the 99.3% confidence level from the counts-in-cells test.

This articles is concerned primarily with the intrinsic anisotropy in the angular distribution of 2281 gamma-ray bursts (GRBs) collected in Current BATSE Gamma-Ray Burst Catalog until the end of year 1998, and, second, with intrinsic anisotropies of three subclasses ("short", "intermediate", "long") of GRBs. Testing based on spherical harmonics of each class, in equatorial coordinates, is presented. Because the sky exposure function of BATSE instrument is not dependent on the right ascension $\alpha$, any non-zero spherical harmonic proportional either to $P_n^m(\sin \delta) \sin m \alpha$ or to $P_n^m(\sin \delta) \cos m\alpha$ with $m \neq 0$ ($\delta$ is the declination), immediately indicates an intrinsic non-zero term. It is a somewhat surprising result that the "intermediate" subclass shows an intrinsic anisotropy at the $97\%$ significance level caused by the high non-zero $P_3^1(\sin \delta) \sin m \alpha$ harmonic. The remaining two subclasses, and the full sample of GRBs, remain isotropic.

We have carried out a principal component analysis for 625 gamma-ray bursts in the BATSE 3B catalog for which non-zero values exist for the nine measured variables. This shows that only two out of the three basic quantities of duration, peak flux and fluence are independent, even if this relation is strongly affected by instrumental effects, and these two account for 91.6% of the total information content. The next most important variable is the fluence in the fourth energy channel (at energies above 320 keV). This has a larger variance and is less correlated with the fluences in the remaining three channels than the latter correlate among themselves. Thus a separate consideration of the fourth channel, and increased attention on the related hardness ratio $H43$ appears useful for future studies. The analysis gives the weights for the individual measurements needed to define a single duration, peak flux and fluence. It also shows that, in logarithmic variables, the hardness ratio $H32$ is significantly correlated with peak flux, while $H43$ is significantly anticorrelated with peak flux. The principal component analysis provides a potentially useful tool for estimating the improvement in information content to be achieved by considering alternative variables or performing various corrections on available measurements