Analyzing light curves of a complete sample of bright Swift long gamma-ray bursts (LGRBs) of which the peak photon fluxes constructed with the bin width of 1 second in the Swift 15-350 keV energy band exceed 2.6 photons cm$^{-2}$s$^{-1}$, we confirm that there do exist the third class in GRBs in addition to short and long GRBs. Being different from previous works based on the duration, fluence, etc. our classification method is based on two properties both quantified with light curve shapes of the prompt emission: the Absolute Deviation from the Constant Luminosity of their cumulative light curve $ADCL$, and the ratio of the mean counts to the maximum counts $\bar C/C_{\rm max}$. These are independent of the distance and the jet opening angle. A cluster analysis via the Gaussian mixture model detects three subclasses: one consisting of LGRBs with small $ADCL$ and large $\bar C/C_{\rm max}$ values referred to as Type I, one with large $ADCL$ and large $\bar C/C_{\rm max}$ referred to as Type II, and one with intermediate $ADCL$ and small $\bar C/C_{\rm max}$, which is composed of contaminating short GRBs with the extended emission. This result is reinforced by different temporal and spectral indices of their X-ray afterglows. The difference is prominent in the temporal index of the steep decay phase in particular: the indices for Type I LGRBs distribute between $-6$ and $-3$ while those for Type II LGRBs between $-3$ and $-2$. From these properties, we propose a possible scenario with different central engines: an accreting black hole and a magnetar.

Overwhelming diversity of long gamma-ray bursts (LGRBs), discovered after the launch of {\it Swift} satellite, is a major obstacle to LGRB studies. Recently, it is shown that the prompt emission of LGRBs is classified into three subclasses: Type I, Type II LGRBs populating separate fundamental planes in a 3D space defined by the peak luminosity, the duration, and the spectral peak energy, and outliers not belonging to either of the planes. Here we show that Type I LGRBs (LGRBs I) exhibit different shapes of light curves from Type II LGRBs (LGRBs II). Furthermore, we demonstrate that this classification has uncovered a new scaling law in the light curves of LGRBs II spanning 8 orders of magnitude from the prompt to late X-ray afterglow emission. The scaled light curve has four distinct phases. The first phase has a characteristic time scale while the subsequent three phases exhibit power law behaviors with different exponents. We discuss its possible interpretation in terms of the emission from an optically thick fireball propagating in the cricumstellar matter at relativistic speeds and argue that the observed four phases correspond to its hydrodynamical phases. Our classification scheme can pin down the intrinsic luminosities of LGRBs II through the scaling law from a sample of polymorphic GRBs. Further refinement of this scheme and scaling law will make a subclass of LGRBs a standard candle as reliable and accurate as Type Ia supernovae in the more distant universe than supernovae can reach.

We analyzed correlations among the rest frame spectral peak energy $E_{\rm p}$, the observed frame 64ms peak isotropic luminosity $L_{\rm p}$ and the isotropic energy $E_{\rm iso}$ for 13 Short Gamma Ray Burst (SGRB) candidates having the measured redshift $z$, $T_{\rm 90}^{\rm obs}/(1+z)<2$ sec and well determined spectral parameters. A SGRB candidate is regarded as a misguided SGRB if it is located in the 3-$\sigma_{\rm int}$ dispersion region from the best-fit function of the $E_{\rm p}$--$E_{\rm iso}$ correlation for Long GRBs (LGRBs) while the others are regarded as secure SGRBs possibly from compact star mergers. Using 8 secure SGRBs out of 13 SGRB candidates, we tested whether $E_{\rm p}$--$E_{\rm iso}$ and $E_{\rm p}$--$L_{\rm p}$ correlations exist for SGRBs. We found that $E_{\rm p}$--$E_{\rm iso}$ correlation for SGRBs($E_{\rm iso} =10^{51.42 \pm 0.15}{\rm erg s^{-1}} ({E_{\rm p}}/{\rm 774.5 keV})^{1.58 \pm 0.28}$) seems to exist with the correlation coefficeint $r=0.91$ and chance probability $p=1.5\times10^{-3}$. We found also that the $E_{\rm p}$--$L_{\rm p}$ correlation for SGRBs($L_{\rm p} = 10^{52.29 \pm 0.066}{\rm erg s^{-1}} ({E_{\rm p}}/{\rm 774.5 keV})^{1.59 \pm 0.11}$) is tighter than $E_{\rm p}$--$E_{\rm iso}$ correlation since $r=0.98$ and $p=1.5\times10^{-5}$. Both correlations for SGRBs are dimmer than those of LGRBs for the same $E_{\rm p}$ by factors $\sim$100 ($E_{\rm p}$--$E_{\rm iso}$) and $\sim$ 5($E_{\rm p}$--$L_{\rm p}$). Applying the tighter $E_{\rm p}$--$L_{\rm p}$ correlation for SGRBs to 71 bright BATSE SGRBs, we found that pseudo redshift $z$ ranges from 0.097 to 2.258 with the mean $<z>$ of 1.05. The redshifts of SGRBs apparently cluster at lower redshift than those of LGRBs ($<z>\sim 2.2 $), which supports the merger scenario of SGRBs.

We estimate the distance modulus to long gamma-ray bursts (LGRBs) using the Type I Fundamental Plane, a correlation between the spectral peak energy $E_{\rm p}$, the peak luminosity $L_{\rm p}$, and the luminosity time $T_{\rm L}$ ($\equiv E_{\rm iso}/L_{\rm p}$ where $E_{\rm iso}$ is isotropic energy) for small Absolute Deviation from Constant Luminosity(ADCL). The Type I Fundamental Plane of LGRBs is calibrated using 8 LGRBs with redshift $z<1.4$. To avoid any assumption on the cosmological model, we use the distance modulus of 557 Type Ia supernovae (SNeIa) from the Union 2 sample. This calibrated Type I Fundamental Plane is used to measure the distance moduli to 9 high-redshift LGRBs with the mean error $\bar \sigma_{\mu}=0.31$, which is comparable with that of SNe Ia $\bar \sigma_{\mu}=0.26$ where $\mu$ stands for the distance modulus. The Type I Fundamental Plane is so tight that our distance moduli have very small uncertainties. From those distance moduli, we obtained the constraint $\Omega_{\rm M}=0.22\pm0.04$ for flat $\Lambda$CDM universe. Adding 9 LGRBs distance moduli ($z>1.4$) to 557 SNeIa distance moduli ($z<1.4$) significantly improves the constraint for non-flat $\Lambda$CDM universe from ($\Omega_{\rm M}, \Omega_{\rm \Lambda}$)=($0.29\pm0.10$, $0.76\pm0.13$) for SNeIa only to ($\Omega_{\rm M}, \Omega_{\rm \Lambda}$)=($0.23\pm0.06$, $0.68\pm0.08$) for SNeIa and 9 LGRBs.

We argue a new classification scheme of long gamma-ray bursts (LGRBs) using the morphology of the cumulative light curve of the prompt emission. We parametrize the morphology by the absolute deviation from their constant luminosity ($ADCL$) and derive the value for 36 LGRBs which have spectropic redshifts, spectral parameters determined by the Band model, 1-second peak fluxes, fluences, and 64-msec resolution light curves whose peak counts are 10 times larger than background fluctuations. Then we devide the sample according to the value of ADCL into two groups ($ADCL < 0.17$ and $ADCL > 0.17$) and, for each group, derive the spectral peak energy $E_{\rm p}$ - peak luminosity $L_{\rm p}$ correlation and the Fundamental Plane of LGRBs, which is a correlation between the spectral peak energy $E_{\rm p}$, the luminosity time $T_{\rm L}$ ($\equiv E_{\rm iso}/L_{\rm p}$ where $E_{\rm iso}$ is isotropic energy) and the peak luminosity $L_{\rm p}$. We find that both of the correlations for both groups are statistically more significant compared with ones derived from all samples. The Fundamental Planes with small and large ADCL are given by $L_{\rm p}=10^{52.53\pm 0.01}(E_{\rm p}/10^{2.71}{\rm keV})^{1.84\pm 0.03} (T_{\rm L}/10^{0.86}{\rm sec})^{0.29\pm0.08}$ with $\chi^2_{\nu}=10.93/14$ and $L_{\rm p}=10^{52.98\pm0.08}(E_{\rm p}/10^{2.71}{\rm keV})^{1.82\pm 0.09} (T_{\rm L}/10^{0.86}{\rm sec})^{0.85\pm 0.27}$ with $\chi^2_{\nu}=7.58/8$, respectively. This fact implies the existence of subclasses of LGRBs characterized by the value of $ADCL$. Also there is a hint for the existence of the intermediate-$ADCL$ class which deviates from both fundamental planes. Both relations are so tight that our result provides a new accurate distance measurement scheme up to the high redshift universe.

We collect and reanalyze about 200 GRB data of prompt-emission with known redshift observed until the end of 2009, and select 101 GRBs which were well observed to have good spectral parameters to determine the spectral peak energy ($E_p$), 1-second peak luminosity ($L_p$) and isotropic energy ($E_{\rm iso}$). Using our newly-constructed database with 101 GRBs, we first revise the $E_p$--$L_p$ and $E_p$--$E_{\rm iso}$ correlations. The correlation coefficients of the revised correlations are 0.889 for 99 degree of freedom for the $E_p$--$L_p$ correlation and 0.867 for 96 degree of freedom for the $E_p$--$E_{\rm iso}$ correlation. These values correspond to the chance probability of $2.18 \times 10^{-35}$ and $4.27 \times 10^{-31}$, respectively. It is a very important issue whether these tight correlations are intrinsic property of GRBs or caused by some selection effect of observations. In this paper, we examine how the truncation of the detector sensitivity affects the correlations, and we conclude they are surely intrinsic properties of GRBs. Next we investigate origins of the dispersion of the correlations by studying their brightness and redshift dependence. Here the brightness (flux or fluence) dependence would be regarded as an estimator of the bias due to the detector threshold. We find a weak fluence-dependence in the $E_p$--$E_{\rm iso}$ correlations and a redshift dependence in the $E_p$--$L_p$ correlation both with 2 $\sigma$ statistical level. These two effects may contribute to the dispersion of the correlations which is larger than the statistical uncertainty. We discuss a possible reason of these dependence and give a future prospect to improve the correlations.

The accuracy and reliability of gamma-ray bursts (GRBs) as distance indicators are strongly restricted by their systematic errors which are larger than statistical errors. These systematic errors might come from either intrinsic variations of GRBs, or systematic errors in observations. In this paper, we consider the possible origins of systematic errors in the following observables, (i) the spectral peak energies (Ep) estimated by Cut-off power law (CPL) function, (ii) the peak luminosities (Lp) estimated by 1 second in observer time. Removing or correcting them, we reveal the true intrinsic variation of the Ep-TL-Lp relation of GRBs. Here TL is the third parameter of GRBs defined as TL ~ Eiso / Lp. Not only the time resolution of Lp is converted from observer time to GRB rest frame time, the time resolution with the largest likelihood is sought for. After removing obvious origin of systematic errors in observation mentioned above, there seems to be still remain some outliers. For this reason, we take account another origin of the systematic error as below, (iii) the contamination of short GRBs or other populations. To estimate the best fit parameters of the Ep-TL-Lp relations from data including outliers, we develop a new method which combine robust regression and an outlier identification technique. Using our new method for 18 GRBs with {\sigma}Ep/Ep < 0.1, we detect 6 outliers and find the Ep-TL-Lp relation become the tightest around 3 second.

In the \swift era, two-component jet models were introduced to explain the complex temporal profiles and the diversity of early afterglows. In this paper, we concentrate on the two-component jet model; first component is the conventional afterglow and second is the emission due to the late internal dissipation such as the late-prompt emission. We suggest herein that the two-component jet model can be probed by the existence of two optical peaks for orphan GRB afterglows. Each peak is caused by its respective jet as its relativistic beaming cone widens to encompass the off-axis line of sight. Typically, the first peak appears at $10^4-10^5$ s and the second at $10^5-10^6$ s. Furthermore, we expect to observe a single, bright X-ray peak at the same time as the first optical peak. Because orphan afterglows do not have prompt emission, it is necessary to monitor the entire sky every $10^4$ s in the X-ray regime. We can test the model with orphan afterglows through the X-ray all-sky survey collaboration and by using ground-based optical telescopes.

We reconsider correlations among the spectral peak energy ($E_p$), 1-second peak luminosity ($L_p$) and isotropic energy (\Eiso), using the database constructed by \citet{yonetoku10} which consists of 109 Gamma-Ray Bursts (GRBs) whose redshifts are known and $E_p$, $L_p$ and \Eiso are well determined. We divide the events into two groups by their data quality. One (gold data set) consists of GRBs with peak energies determined by the Band model with four free parameters. On the other hand, GRBs in the other group (bronze data set) have relatively poor energy spectra so that their peak energies were determined by the Band model with fixed spectral index (i.e. three free parameters) or by the Cut-off power law (CPL) model with three free parameters. Using only the gold data set we found the intrinsic dispersion in $\log L_p$ ($=\sigma_{\rm int}$) is 0.13 and 0.22 for \tsutsui correlation ($T_L \equiv E_{\rm iso}/L_p$) and \yonetoku correlation, respectively. We also find that GRBs in the bronze data set have systematically larger $E_p$ than expected by the correlations constructed with the gold data set. This means that the intrinsic dispersion of correlations among $E_p$, $L_p$, and \Eiso of GRBs depends on the quality of data set. At present, using \tsutsui correlation with gold data set, we would be able to determine the luminosity distance with $\sim 16%$ error, which might be useful to determine the nature of the dark energy at high redshift $z > 3$.

We consider two empirical relations using data only from the prompt emission of Gamma-Ray Bursts (GRBs), peak energy ($E_p$) - peak luminosity ($L_p$) relation (so called Yonetoku relation) and $E_p$-isotropic energy ($E_{\rm iso}$) relation (so called Amati relation). We first suggest the independence of the two relations although they have been considered similar and dependent. From this viewpoint, we compare constraints on cosmological parameters, $\Omega_m$ and $\Omega_{\Lambda}$, from the Yonetoku and Amati relations calibrated by low-redshift GRBs with $z < 1.8$. We found that they are different in 1-$\sigma$ level, although they are still consistent in 2-$\sigma$ level. This and the fact that both Amati and Yonetoku relations have systematic errors larger than statistical errors suggest the existence of a hidden parameter of GRBs. We introduce the luminosity time $T_L$ defined by $T_L\equiv E_{\rm iso}/L_p$ as a hidden parameter to obtain a generalized Yonetoku relation as $(L_p/{10^{52} \rm{erg s^{-1}}}) = 10^{-3.88\pm0.09}(E_p/{\rm{keV}})^{1.84\pm0.04} (T_L/{\rm{s}})^{-0.34\pm0.04}$. The new relation has much smaller systematic error, 30%, and can be regarded as "Fundamental plane" of GRBs. We show a possible radiation model for this new relation. Finally we apply the new relation for high-redshift GRBs with $1.8 < z < 5.6$ to obtain $(\Omega_m,\Omega_{\Lambda}) = (0.16^{+0.04}_{-0.06},1.20^{+0.03}_{-0.09})$, which is consistent with the concordance cosmological model within 2-$\sigma$ level.

We extend the Hubble diagram up to $z = 5.6$ using 63 gamma-ray bursts (GRBs) via peak energy-peak luminosity relation (so called Yonetoku relation), and obtain constraints on cosmological parameters including dynamical dark energy parametrized by $P/\rho\equiv w(z) = w_0 + w_a \cdot z/(1+z)$. It is found that the current GRB data are consistent with the concordance model, ($\Omega_m = 0.28, \Omega_{\Lambda} = 0.72, w_0 = -1, w_a = 0$), within two sigma level. Although constraints from GRBs themselves are not so strong, they can improve the conventional constraints from SNeIa because GRBs have much higher redshifts. Further we estimate the constraints on the dark-energy parameters expected by future observations with GLAST (Gamma-ray Large Area Space Telescope) and \swift by Monte-Carlo simulation. Constraints would improve substantially with another 150 GRBs.

We calibrated the peak energy-peak luminosity relation of GRBs (so called Yonetoku relation) using 33 events with the redshift $z < 1.62$ without assuming any cosmological models. The luminosity distances to GRBs are estimated from those of large amount of Type Ia supernovae with $z<1.755$. This calibrated Yonetoku relation can be used as a new cosmic distance ladder toward higher redshifts. We determined the luminosity distances of 30 GRBs in $1.8 < z < 5.6$ using the calibrated relation and plotted the likelihood contour in $(\Omega_m,\Omega_\Lambda)$ plane. We obtained $(\Omega_m, \Omega_{\Lambda})= (0.37^{+0.14}_{-0.11}, 0.63^{+0.11}_{-0.14})$ for a flat universe. Since our method is free from the circularity problem, we can say that our universe in $1.8 < z < 5.6$ is compatible with the so called concordance cosmological model derived for $z < 1.8$. This suggests that the time variation of the dark energy is small or zero up to $z\sim 6$.

We compared redshifts $z_Y$ from Yonetoku relation and $z_{lag}$ from the lag-luminosity relation for 565 BASTE GRBs and were surprised to find that the correlation is very low. Assuming that the luminosity is a function of both $z_Y$ and the intrinsic spectral lag $\tau_{lag}$, we found a new redshift dependent lag-luminosity relation as $L=7.5\times 10^{50}{\rm erg/s}(1+z)^{2.53}\tau_{lag}^{-0.282}$ with the correlation coefficient of 0.77 and the chance probability of $7.9\times 10^{-75}$. To check the validity of this method, we examined the other luminosity indicator, Amati relation, using $z_Y$ and the observed fluence and found the correlation coefficient of 0.92 and the chance probability of $5.2\times 10^{-106}$. Although the spectral lag is computed from two channels of BATSE, our new lag-luminosity relation suggests that a possible lag-luminosity relation in the \swift era should also depend on redshift.