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.

We observed an X-ray afterglow of GRB 060904A with the Swift and Suzaku satellites. We found rapid spectral softening during both the prompt tail phase and the decline phase of an X-ray flare in the BAT and XRT data. The observed spectra were fit by power-law photon indices which rapidly changed from $\Gamma = 1.51^{+0.04}_{-0.03}$ to $\Gamma = 5.30^{+0.69}_{-0.59}$ within a few hundred seconds in the prompt tail. This is one of the steepest X-ray spectra ever observed, making it quite difficult to explain by simple electron acceleration and synchrotron radiation. Then, we applied an alternative spectral fitting using a broken power-law with exponential cutoff (BPEC) model. It is valid to consider the situation that the cutoff energy is equivalent to the synchrotron frequency of the maximum energy electrons in their energy distribution. Since the spectral cutoff appears in the soft X-ray band, we conclude the electron acceleration has been inefficient in the internal shocks of GRB 060904A. These cutoff spectra suddenly disappeared at the transition time from the prompt tail phase to the shallow decay one. After that, typical afterglow spectra with the photon indices of 2.0 are continuously and preciously monitored by both XRT and Suzaku/XIS up to 1 day since the burst trigger time. We could successfully trace the temporal history of two characteristic break energies (peak energy and cutoff energy) and they show the time dependence of $\propto t^{-3} \sim t^{-4}$ while the following afterglow spectra are quite stable. This fact indicates that the emitting material of prompt tail is due to completely different dynamics from the shallow decay component. Therefore we conclude the emission sites of two distinct phenomena obviously differ from each other.