Observationally, it has been reported that the densest stellar system in the Universe does not exceed a maximum stellar surface density, $\Sigma^{\max}_{*}$ = $3\times10^5$M$_{\odot}$pc$^{-2}$, throughout a wide physical scale ranging from star cluster to galaxy. This suggests there exists a fundamental physics which regulates the star formation and stellar density. However, factors that determine this maximum limit are not clear. In this study, we show that $\Sigma^{\max}_{*}$ of galaxies is not a constant as previous work reported, but actually depends on the stellar mass. We select galaxy sample from the Sloan Digital Sky Survey Data Release 12 at $z=0.01-0.5$. In contrast to a constant maximum predicted by theoretical models, $\Sigma^{\max}_{*}$ strongly depends on stellar mass especially for less massive galaxies with $\sim10^{10}$M$_{\odot}$. We also found that a majority of high-$\Sigma_{*}$ galaxies show red colours and low star-formation rates. These galaxies probably reach the $\Sigma^{\max}_{*}$ as a consequence of the galaxy evolution from blue star forming to red quiescent by quenching star formation. One possible explanation of the stellar-mass dependency of $\Sigma^{\max}_{*}$ is a mass dependent efficiency of stellar feedback. The stellar feedback could be relatively more efficient in a shallower gravitational potential, which terminates star formation quickly before the stellar system reaches a high stellar density.

Fast radio bursts (FRBs) are mysterious radio bursts with a time scale of approximately milliseconds. Two populations of FRB, namely repeating and non-repeating FRBs, are observationally identified. However, the differences between these two and their origins are still cloaked in mystery. Here we show the time-integrated luminosity-duration ($L_{\nu}$-$w_{\rm int,rest}$) relations and luminosity functions (LFs) of repeating and non-repeating FRBs in the FRB Catalogue project. These two populations are obviously separated in the $L_{\nu}$-$w_{\rm int,rest}$ plane with distinct LFs, i.e., repeating FRBs have relatively fainter $L_{\nu}$ and longer $w_{\rm int,rest}$ with a much lower LF. In contrast with non-repeating FRBs, repeating FRBs do not show any clear correlation between $L_{\nu}$ and $w_{\rm int,rest}$. These results suggest essentially different physical origins of the two. The faint ends of the LFs of repeating and non-repeating FRBs are higher than volumetric occurrence rates of neutron-star mergers and accretion-induced collapse (AIC) of white dwarfs, and are consistent with those of soft gamma-ray repeaters (SGRs), type Ia supernovae, magnetars, and white-dwarf mergers. This indicates two possibilities: either (i) faint non-repeating FRBs originate in neutron-star mergers or AIC and are actually repeating during the lifetime of the progenitor, or (ii) faint non-repeating FRBs originate in any of SGRs, type Ia supernovae, magnetars, and white-dwarf mergers. The bright ends of LFs of repeating and non-repeating FRBs are lower than any candidates of progenitors, suggesting that bright FRBs are produced from a very small fraction of the progenitors regardless of the repetition. Otherwise, they might originate in unknown progenitors.

Investigating the Gunn-Peterson trough of high redshift quasars (QSOs) is a powerful way to reveal the cosmic reionization. As one of such attempts, we perform a series of analyses to examine the absorption lines observed with one of the highest redshift QSOs, PSO J006.1240+39.2219, which we previously discovered at z = 6.62. Using the Subaru telescope, we obtained medium-resolution spectrum with a total exposure time of 7.5 hours. We calculate the Ly$\alpha$ transmission in different redshift bins to determine the near zone radius and the optical depth at 5.6$<$z$<$6.5. We find a sudden change in the Ly$\alpha$ transmission at 5.75$<$z$<$5.86, which is consistent with the result from the literature. The near zone radius of the QSO is 5.79$\pm$0.09 $p$Mpc, within the scatter of the near zone radii of other QSOs measured in previous studies. We also analyze the dark gap distribution to probe the neutral hydrogen fractions beyond the saturation limit of the Gunn-Peterson trough. We extend the measurement of the dark gaps to 5.7$<$z$<$6.3. We find that the gap widths increase with increasing redshifts, suggesting more neutral Universe at higher redshifts. However, these measurements strongly depend on the continuum modeling. As a continuum model-free attempt, we also perform the dark-pixel counting analysis, to find the upper limit of $\langle x_{\rm H I} \rangle \sim$0.6 (0.8) at $z<$5.8 ($z>$5.8). All three analyses based on this QSO show increasingly neutral hydrogen towards higher redshifts, adding precious measurements up to z$\sim$6.5.

Based on the cosmic shear data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), Kilbinger et al. (2013) obtained a constraint on the amplitude of matter fluctuations of $\sigma_8({\Omega_\mathrm{m}}/0.27)^{0.6}=0.79\pm0.03$ from the two-point correlation function (2PCF). This is $\approx3\sigma$ lower than the value $0.89\pm0.01$ derived from Planck data on cosmic microwave background (CMB) anisotropies. On the other hand, based on the same CFHTLenS data, but using the power spectrum, and performing a different analysis, Liu et al. (2015) obtained the higher value of $\sigma_8({\Omega_\mathrm{m}}/0.27)^{0.64}=0.87^{+0.05}_{-0.06}$. We here investigate the origin of this difference, by performing a fair side-by-side comparison of the 2PCF and power spectrum analyses on CFHTLenS data. We find that these two statistics indeed deliver different results, even when applied to the same data in an otherwise identical procedure. We identify excess power in the data on small scales ($\ell>5,000$) driving the larger values inferred from the power spectrum. We speculate on the possible origin of this excess small-scale power. More generally, our results highlight the utility of analysing the 2PCF and the power spectrum in tandem, to discover (and to help control) systematic errors.

We present the analysis of a new near-infrared (NIR) spectrum of a recently discovered $z=6.621$ quasar PSO J006+39 in an attempt to explore the early growth of supermassive black holes (SMBHs). This NIR (rest-frame ultraviolet, UV) spectrum shows blue continuum slope and rich metal emission lines in addition to Ly$\alpha$ line. We utilize the MgII line width and the rest frame luminosity $L_\text{3000\AA}$ to find the mass of SMBH ($M_\text{BH}$) to be $\sim 10^8 M_\odot$, making this one of the lowest mass quasars at high redshift. The power-law slope index ($\alpha_\lambda$) of the continuum emission is $-2.94\pm0.03$, significantly bluer than the slope of $\alpha_\lambda=-7/3$ predicted from standard thin disc models. We fit the spectral energy distribution (SED) using a model which can fit local SMBHs, which includes warm and hot Comptonisation powered by the accretion flow as well as an outer standard disc. The result shows that the very blue slope is probably produced by a small radial ($\sim230$ gravitational radius, $R_\text{g}$) extent of the standard accretion disc. All plausible SED models require that the source is super-Eddington ($L_\text{bol}/L_\text{Edd} \gtrsim 9$), so the apparently small disc may simply be the inner funnel of a puffed up flow, and clearly the SMBH in this quasar is in a rapid growth phase. We also utilize the rest-frame UV emission lines to probe the chemical abundance in the broad line region (BLR) of this quasar. We find that this quasar has super solar metallicity through photoionization model calculations.

In a recent study, using the distribution of galaxies in the north galactic pole of SDSS DR7 region enclosed in a 500$\mpch$ box, we carried out our ELUCID simulation (Wang et al. 2016, ELUCID III). Here we {\it light} the dark matter halos and subhalos in the reconstructed region in the simulation with galaxies in the SDSS observations using a novel {\it neighborhood} abundance matching method. Before we make use of thus established galaxy-subhalo connections in the ELUCID simulation to evaluate galaxy formation models, we set out to explore the reliability of such a link. For this purpose, we focus on the following a few aspects of galaxies: (1) the central-subhalo luminosity and mass relations; (2) the satellite fraction of galaxies; (3) the conditional luminosity function (CLF) and conditional stellar mass function (CSMF) of galaxies; and (4) the cross correlation functions between galaxies and the dark matter particles, most of which are measured separately for all, red and blue galaxy populations. We find that our neighborhood abundance matching method accurately reproduces the central-subhalo relations, satellite fraction, the CLFs and CSMFs and the biases of galaxies. These features ensure that thus established galaxy-subhalo connections will be very useful in constraining galaxy formation processes. And we provide some suggestions on the three levels of using the galaxy-subhalo pairs for galaxy formation constraints. The galaxy-subhalo links and the subhalo merger trees in the SDSS DR7 region extracted from our ELUCID simulation are available upon request.

Accurate reconstruction of the spatial distributions of the Point Spread Function (PSF) is crucial for high precision cosmic shear measurements. Nevertheless, current methods are not good at recovering the PSF fluctuations of high spatial frequencies. In general, the residual PSF fluctuations are spatially correlated, therefore can significantly contaminate the correlation functions of the weak lensing signals. We propose a method to correct for this contamination statistically, without any assumptions on the PSF and galaxy morphologies or their spatial distribution. We demonstrate our idea with the data from the W2 field of CFHTLenS.

In this work, we present a novel centroiding method based on Fourier space Phase Fitting(FPF) for Point Spread Function(PSF) reconstruction. We generate two sets of simulations to test our method. The first set is generated by GalSim with elliptical Moffat profile and strong anisotropy which shifts the center of the PSF. The second set of simulation is drawn from CFHT i band stellar imaging data. We find non-negligible anisotropy from CFHT stellar images, which leads to $\sim$0.08 scatter in unit of pixels using polynomial fitting method Vakili and Hogg (2016). And we apply FPF method to estimate the centroid in real space, this scatter reduces to $\sim$0.04 in SNR=200 CFHT like sample. In low SNR (50 and 100) CFHT like samples, the background noise dominates the shifting of the centroid, therefore the scatter estimated from different methods are similar. We compare polynomial fitting and FPF using GalSim simulation with optical anisotropy. We find that in all SNR$\sim$50, 100 and 200) samples, FPF performs better than polynomial fitting by a factor of $\sim$3. In general, we suggest that in real observations there are anisotropy which shift the centroid, and FPF method is a better way to accurately locate it.

As the second paper of a series on studying galaxy-galaxy lensing signals using the Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we present our measurement and modelling of the lensing signals around groups of galaxies. We divide the groups into four halo mass bins, and measure the signals around four different halo-center tracers: brightest central galaxy (BCG), luminosity-weighted center, number-weighted center and X-ray peak position. For X-ray and SDSS DR7 cross identified groups, we further split the groups into low and high X-ray emission subsamples, both of which are assigned with two halo-center tracers, BCGs and X-ray peak positions. The galaxy-galaxy lensing signals show that BCGs, among the four candidates, are the best halo-center tracers. We model the lensing signals using a combination of four contributions: off-centered NFW host halo profile, sub-halo contribution, stellar contribution, and projected 2-halo term. We sample the posterior of 5 parameters i.e., halo mass, concentration, off-centering distance, sub halo mass, and fraction of subhalos via a MCMC package using the galaxy-galaxy lensing signals. After taking into account the sampling effects (e.g. Eddington bias), we found the best fit halo masses obtained from lensing signals are quite consistent with those obtained in the group catalog based on an abundance matching method, except in the lowest mass bin. Subject headings: (cosmology:) gravitational lensing, galaxies: clusters: general

Halo bias is the one of the key ingredients of the halo models. It was shown at a given redshift to be only dependent, to the first order, on the halo mass. In this study, four types of cosmic web environments: clusters, filaments, sheets and voids are defined within a state of the art high resolution $N$-body simulation. Within those environments, we use both halo-dark matter cross-correlation and halo-halo auto correlation functions to probe the clustering properties of halos. The nature of the halo bias differs strongly among the four different cosmic web environments we describe. With respect to the overall population, halos in clusters have significantly lower biases in the {$10^{11.0}\sim 10^{13.5}h^{-1}\rm M_\odot$} mass range. In other environments however, halos show extremely enhanced biases up to a factor 10 in voids for halos of mass {$\sim 10^{12.0}h^{-1}\rm M_\odot$}. Such a strong cosmic web environment dependence in the halo bias may play an important role in future cosmological and galaxy formation studies. Within this cosmic web framework, the age dependency of halo bias is found to be only significant in clusters and filaments for relatively small halos $\la 10^{12.5}\msunh$.

Reconstruction of the point spread function (PSF) is a critical process in weak lensing measurement. We develop a real-data based and galaxy-oriented pipeline to compare the performances of various PSF reconstruction schemes. Making use of a large amount of the CFHTLenS data, the performances of three classes of interpolating schemes - polynomial, Kriging, and Shepard - are evaluated. We find that polynomial interpolations with optimal orders and domains perform the best. We quantify the effect of the residual PSF reconstruction error on shear recovery in terms of the multiplicative and additive biases, and their spatial correlations using the shear measurement method of Zhang et al. (2015). We find that the impact of PSF reconstruction uncertainty on the shear-shear correlation can be significantly reduced by cross correlating the shear estimators from different exposures. It takes only 0.2 stars (SNR > 100) per square arcmin on each exposure to reach the best performance of PSF interpolation, a requirement that is generally satisfied in the CFHTlenS data.

Strange quark matter (SQM) may be the true ground state of hadronic matter, indicating that the observed pulsars may actually be strange stars, but not neutron stars. According to this SQM hypothesis, the existence of a hydrostatically stable sequence of strange quark matter stars has been predicted, ranging from 1 --- 2 solar mass strange stars, to smaller strange dwarfs and even strange planets. While gravitational wave (GW) astronomy is expected to open a new window to the universe, it will shed light on the searching for SQM stars. Here we show that due to their extreme compactness, strange planets can spiral very close to their host strange stars, without being tidally disrupted. Like inspiraling neutron stars or black holes, these systems would serve as a new kind of sources for GW bursts, producing strong gravitational waves at the final stage. The events occurring in our local Universe can be detected by the upcoming gravitational wave detectors, such as Advanced LIGO and the Einstein Telescope. This effect provides a unique probe to SQM objects and is hopefully a powerful tool for testing the SQM hypothesis.

We present NuSTAR spectral and timing studies of the Supergiant Fast X-ray Transient (SFXT) IGR J17544-2619. The spectrum is well-described by a ~1 keV blackbody and a hard continuum component, as expected from an accreting X-ray pulsar. We detect a cyclotron line at 17 keV, confirming that the compact object in IGR J17544-2619 is indeed a neutron star. This is the first measurement of the magnetic field in a SFXT. The inferred magnetic field strength, B = (1.45 +/- 0.03) * 10^12 G * (1+z) is typical of neutron stars in X-ray binaries, and rules out a magnetar nature for the compact object. We do not find any significant pulsations in the source on time scales of 1-2000 s.

We present NuSTAR observations of the bright Seyfert 2 galaxy NGC 2110 obtained in 2012, when the source was at the highest flux level ever observed, and in 2013, when the source was at a more typical flux level. We include archival observations from other X-ray satellites, namely XMM-Newton, Suzaku, BeppoSAX, Chandra and Swift. Simultaneous NuSTAR and Swift broad band spectra (in the 3-80 keV range) indicate a cutoff energy $E_{\rm c}>210$ keV, with no detectable contribution from Compton reflection. NGC 2110 is one of the very few sources where no evidence for distant Compton thick scattering is found and, by using temporal information collected over more than a decade, we investigate variations of the iron K$\alpha$ line on time scales of years. The Fe K$\alpha$ line is likely the sum of two components: one constant (originating from distant Compton-thick material) and the other one variable and linearly correlated with the source flux (possibly arising from Compton-thin material much closer to the black hole).

In this paper, we set the new limits on the equation of state parameter (EoS) of dark energy with the observations of cosmic microwave background radiation (CMB) from Planck satellite, the type Ia supernovae from Pan-STARRS and the baryon acoustic oscillation (BAO). We consider two parametrization forms of EoS: a constant $w$ and time evolving $w(a)=w_0+w_a(1-a)$. The results show that with a constant EoS, $w=-1.141\pm{0.075}$ ($68\%~C.L.$), which is consistent with $\Lambda$CDM at about $2\sigma$ confidence level. For a time evolving $w(a)$ model, we get $w_0=-1.09^{+0.16}_{-0.18}$ ($1\sigma~C.L.$), $w_a=-0.34^{+0.87}_{-0.51}$ ($1\sigma~C.L.$), and in this case $\Lambda$CDM can be comparable with our observational data at $1\sigma$ confidence level. In order to do the parametrization independent analysis, additionally we adopt the so called principal component analysis (PCA) method, in which we divide redshift range into several bins and assume $w$ as a constant in each redshift bin (bin-w). In such bin-w scenario, we find that for most of the bins cosmological constant can be comparable with the data, however, there exists few bins which give $w$ deviating from $\Lambda$CDM at more than $2\sigma$ confidence level, which shows a weak hint for the time evolving behavior of dark energy. To further confirm this hint, we need more data with higher precision.

The GeV$-$TeV $\gamma-$ray line signal is the smoking gun signature of the dark matter annihilation or decay. The detection of such a signal is one of the main targets of some space-based telescopes, including Fermi-LAT and the upcoming CALET, DAMPE and Gamma-400. An important feature of the dark-matter-annihilation originated $\gamma-$ray line photons is their concentration at the center of the Galaxy. So far no reliable $\gamma-$ray line has been detected by Fermi-LAT and the upper limits on the cross section of annihilation into $\gamma-$rays have been reported. We use these upper limits to estimate the "maximal" number of $\gamma-$ray line photons detectable for Fermi-LAT, DAMPE and Gamma-400 and then investigate the spatial distribution of these photons. We show that usually the center of the distribution region will be offset from the Galactic centre (Sgr A$^{\star}$) due to the limited statistics. Such a result is almost independent of the dark matter distribution models and renders the reconstruction of the dark matter distribution with the $\gamma-$ray line signal very challenging for the foreseeable space-based detectors.

The ultracompact minihalos could be formed during the earlier epoch of the universe. The dark matter annihilation within them is very strong due to the steep density profile, $\rho \sim r^{-2.25}$. The high energy electrons and positrons from the dark matter annihilation can inverse Compton scatter (ICS) with the background photons, such as CMB photons, to acquire higher energy. On the other hand, the synchrotron radiation can also be produced when they meet the magnetic field. In this paper, we study the signals from the UCMHs due to the dark matter annihilation for the radio, X-ray and $\gamma$-ray band. We found that for the radio emission the UCMHs can provide one kind of source for the radio excess observed by ARCADE 2. But the X-ray signals due to the ICS effect or the $\gamma$-ray signals mainly due to the prompt emission from dark matter would exceed the present observations, such as Fermi, COMPTEL and CHANDRA. We found that the strongest limits on the fraction of UCMHs come from the X-ray observations and the constraints from the radio data are the weakest.

We present a multi-wavelength analysis of the very fast X-ray transient MAXI J0158-744, which was detected by MAXI/GSC on 2011 November 11. The subsequent exponential decline of the X-ray flux was followed with Swift observations, all of which revealed spectra with low temperatures (~100eV) indicating that MAXI J0158-744 is a new Supersoft Source (SSS). The Swift X-ray spectra near maximum show features around 0.8 keV that we interpret as possible absorption from OVIII, and emission from O, Fe, and Ne lines. We obtained SAAO and ESO optical spectra of the counterpart early in the outburst and several weeks later. The early spectrum is dominated by strong Balmer and HeI emission, together with weaker HeII emission. The later spectrum reveals absorption features that indicate a B1/2IIIe spectral type, and all spectral features are at velocities consistent with the Small Magellanic Cloud. At this distance, it is a luminous SSS (>10^37 erg/s) but whose brief peak luminosity of >10^39 erg/s in the 2-4 keV band makes it the brightest SSS yet seen at "hard" X-rays. We propose that MAXI J0158-744 is a Be-WD binary, and the first example to possibly enter ULX territory. The brief hard X-ray flash could possibly be a result of the interaction of the ejected nova shell with the B star wind in which the white dwarf (WD) is embedded. This makes MAXI J0158-744 only the third Be/WD system in the Magellanic Clouds, but it is by far the most luminous. The properties of MAXI J0158-744 give weight to previous suggestions that SSS in nearby galaxies are associated with early-type stellar systems.

We present a short Chandra observation that confirms a previous unidentified extended X-ray source, G308.3-1.4, as a new supernova remnant (SNR) in the Milky Way. Apart from identifying its SNR nature, a bright X-ray point source has also been discovered at the geometrical center. Its X-ray spectral properties are similar to those of a particular class of neutron star known as central compact objects (CCOs). On the other hand, the optical properties of this counterpart suggests it to be a late-type star. Together with the interesting ~ 1.4 hours X-ray periodicity found by Chandra, this system can possibly provide the first direct evidence of a compact binary survived in a supernova explosion.

Observational growth rate data had been derived from observations of redshift distortions in galaxy redshift surveys. Here we use the growth rate data to place constraints on the dark energy model parameters. By performing a joint analysis with the Type Ia supernova, baryon acoustic oscillation and cosmic microwave background data, it is found that the growth rate data are useful for improving the constraints. The joint constraints show that the $\Lambda$CDM model is still in good agreement with current observations, although a time-variant dark energy still cannot be ruled out. It is argued that the growth rate data are helpful for understanding the dark energy. With more accurate data available in the future, we will have a powerful tool for constraining the cosmological and dark energy parameters.

We investigate in detail some popular cosmological models in light of the latest observational data, including the Union2.1 supernovae compilation, the baryon acoustic oscillation measurements from the WiggleZ Dark Energy Survey, the cosmic microwave background information from the WMAP 7-year observations along with the observational Hubble parameter data. Based on the model selection statistics such as the Akaike and the Bayesian information criterias, we compare different models to assess the worth of them. We do not assume a flat universe in the fitting. Our results show that the concordance $\Lambda$CDM model remains the best one to explain the data, while the DGP model is clearly disfavored by the data. Among these models, those whose parameters can reduce themselves to the $\Lambda$CDM model provide good fits to the data. These results indicate that for the current data, there is no obvious evidence supporting any more complex models over the simplest $\Lambda$CDM model.

A new method is proposed for directly measuring the expansion rate of the universe through very precise measurement of the fluence of extremely stable sources. The method is based on the definition of the luminosity distance and its change along the time due to the cosmic expansion. It is argued that galaxies may be chosen as the targets of the observation to perform the measurement. We show that, by simultaneously increasing the observation time and physically adding the fluences from different galaxies, the requirement on the relative precision of the detector for an observation of 1 second on a single galaxy can be relaxed to $10^{-5}$. Benefiting from the abundance of galaxies in the universe, the method may be quite promising.

The large scatters of luminosity relations of gamma-ray bursts (GRBs) have been one of the most important reasons that prevent the extensive applications of GRBs in cosmology. Many efforts have been made to seek tight luminosity relations. With the latest sample of 116 GRBs with measured redshift and spectral parameters, we investigate 6 two-dimensional (2D) correlations and 14 derived three-dimensional (3D) correlations of GRBs to explore the possibility of decreasing the intrinsic scatters of the luminosity relations of GRBs. We find the 3D correlation of $E_{\mathrm{peak}}$--$\tau_{\mathrm{RT}}$--$L$ to be evidently tighter (at the $2 \sigma$ confidence level) than its corresponding 2D correlations, i.e., the $E_{\mathrm{peak}}$--$L$ and $\tau_{\mathrm{RT}}$--$L$ correlations. In addition, the coefficients before the logarithms of $E_{\mathrm{peak}}$ and $\tau_{\mathrm{RT}}$ in the $E_{\mathrm{peak}}$--$\tau_{\mathrm{RT}}$--$L$ correlation are almost exact opposites of each other. Inputting this situation as a prior reduces the relation to $L \propto (E_{\mathrm{peak}}^{'} / \tau_{\mathrm{RT}}^{'})^{0.842 \pm 0.064}$, where $E_{\mathrm{peak}}^{'}$ and $\tau_{\mathrm{RT}}^{'}$ denote the peak energy and minimum rise time in the GRB rest frame. We discuss how our findings can be interpreted/understood in the framework of the definition of the luminosity (energy released in units of time). Our argument about the connection between the luminosity relations of GRBs and the definition of the luminosity provides a clear direction for exploring tighter luminosity relations of GRBs in the future.

A new parameterization for the dark energy equation of state(EoS) is proposed and some of its cosmological consequences are also investigated. This new parameterization is the modification of Efstathiou' dark energy EoS parameterization. $w (z)$ is a well behaved function for $z\gg1$ and has same behavior in $z$ at low redshifts with Efstathiou' parameterization. In this parameterization there are two free parameter $w_0$ and $w_a$. We discuss the constraints on this model's parameters from current observational data. The best fit values of the cosmological parameters with $1\sigma$ confidence-level regions are: $\Omega_m=0.2735^{+0.0171}_{-0.0163}$, $w_0=-1.0537^{+0.1432}_{-0.1511}$ and $w_a=0.2738^{+0.8018}_{-0.8288}$.

ROSAT all-sky survey (RASS) data have provided another window to search for supernova remnants (SNRs). In reexamining this data archive, a list of unidentified extended X-ray objects have been suggested as promising SNR candidate. However, most of these targets have not yet been fully explored by the state-of-art X-ray observatories. For selecting a pilot target for a long-term identification campaign, we have observed the brightest candidate, G308.3-1.4, with Chandra X-ray observatory. An incomplete shell-like X-ray structure which well-correlated with the radio shell emission at 843 MHz has been revealed. The X-ray spectrum suggests the presence of a shock-heated plasma. All these evidences confirm G308.3-1.4 as a SNR. The brightest X-ray point source detected in this field-of-view is also the one locates closest to the geometrical center of G308.3-1.4, which has a soft spectrum. The intriguing temporal variability and the identification of optical/infrared counterpart rule out the possibility of an isolated neutron star. On the other hand, the spectral energy distribution from Ks band to R band suggests a late-type star. Together with a putative periodicity of \sim1.4 hrs, the interesting excesses in V, B bands and H-alpha suggest this source as a promising candidate of a compact binary survived in a supernova explosion (SN).

It has been proposed that ultracompact minihalos (UCMHs) might be formed in earlier epoch. If dark matter consists of Weakly Interacting Massive Particles (WIMPs), UCMHs can be treated as the {\gamma}-ray sources due to dark matter annihilation within them. In this paper, we investigate the contributions of UCMHs formed during three phase transi- tions (i.e., electroweak symmetry breaking, QCD confinement and e+ e- annihilation) to the extragalactic {\gamma}-ray background. Moreover, we use the Fermi-LAT observation data of the extragalactic {\gamma}-ray background to get the constraints on the current abundance of UCMHs produced during these phase transitions. We also compare these results with those obtained from Cosmic Microwave Background (CMB) observations and find that the constraints from the Fermi-LAT are more stringent than those from CMB

A new kind of dark matter structures, ultracompact minihalos (UCMHs) was proposed recently. They would be formed during the radiation dominated epoch if the large density perturbations are existent. Moreover, if the dark matter is made up of weakly interacting massive particles, the UCMHs can have effect on cosmological evolution because of the high density and dark matter annihilation within them. In this paper, one new parameter is introduced to consider the contributions of UCMHs due to the dark matter annihilation to the evolution of cosmology, and we use the current and future CMB observations to obtain the constraint on the new parameter and then the abundance of UCMHs. The final results are applicable for a wider range of dark matter parameters

Using the optical data from the Wide component of the CFHT Legacy Survey, and new ultraviolet data from GALEX, we study the colours and specific star formation rates (SSFR) of ~100 galaxy clusters at 0.16<z<0.36, over areas extending out to radii of r~7Mpc. We use a multicolour, statistical background subtraction method to study the galaxy population at this radius; thus our results pertain to those galaxies which constitute an excess over the average field density. We find that the average SSFR, and its distribution, of the star-forming galaxies (with SFR>0.7 M_sun/yr at z~0.2 and SFR>1.2 M_sun/yr at z~0.3) have no measurable dependence on the cluster-centric radius, and are consistent with the field values. However, the fraction of galaxies with SFR above these thresholds, and the fraction of optically blue galaxies, are lower for the overdense galaxy population in the cluster outskirts compared with the average field value, at all stellar masses M*>10^{9.8} M_sun and at all radii out to at least 7Mpc. Most interestingly, the fraction of blue galaxies that are forming stars at a rate below our UV detection limit is much higher in all radial bins around our cluster sample, compared with the general field value. This is most noticeable for massive galaxies M*>10^{10.7} M_sun; while almost all blue field galaxies of this mass have detectable star formation, this is true for less than 20% of the blue cluster galaxies, even at 7Mpc from the cluster centre. Our results support a scenario where galaxies are pre-processed in locally overdense regions, in a way that reduces their SFR below our UV detection limit, but not to zero.

We present results of deep radio continuum observations of two galactic globular clusters 47 Tucanae (47 Tuc) and Omega Centauri (Omega Cen) with Australia Telescope Compact Array (ATCA). No statistically significant evidence for radio emission was found from the central region for the two clusters. However, both clusters show a 2.5-sigma detection near the center that may be confirmed by future deeper radio observations. The 3-sigma upper limits of the radio observations is 20 and 40 microJy/beam for Omega Cen and 47 Tuc, respectively. By using the fundamental plane of accreting black holes which describes the relationship between radio luminosity, X-ray luminosity and black hole mass, we constrain the mass of a possible intermediate-mass black hole (IMBH) in the globualar clusters. We also compare our results with other globular clusters and discuss the existence of IMBHs in globular clusters.

Globular clusters (GCs) are emerging as a new class of gamma-ray emitters, thanks to the data obtained from the Fermi Gamma-ray Space Telescope. By now, eight GCs are known to emit gamma-rays at energies >100~MeV. Based on the stellar encounter rate of the GCs, we identify potential gamma-ray emitting GCs out of all known GCs that have not been studied in details before. In this paper, we report the discovery of a number of new gamma-ray GCs: Liller 1, NGC 6624, and NGC 6752, and evidence for gamma-ray emission from M80, NGC 6139, and NGC 6541, in which gamma-rays were found within the GC tidal radius. With one of the highest metallicity among all GCs in the Milky Way, the gamma-ray luminosity of Liller 1 is found to be the highest of all known gamma-ray GCs. In addition, we confirm a previous report of significant gamma-ray emitting region next to NGC 6441. We briefly discuss the observed offset of gamma-rays from some GC cores. The increasing number of known gamma-ray GCs at distances out to ~10 kpc is important for us to understand the gamma-ray emitting mechanism and provides an alternative probe to the underlying millisecond pulsar populations of the GCs.

We analyse the stellar and hot gas content of 18 nearby, low-mass galaxy clusters, detected in redshift space and selected to have a dynamical mass 3E14<M/Msun<6E14, as measured from the 2dF Galaxy Redshift Survey. We combine X-ray measurements from both Chandra and XMM with ground-based near-infrared observations from CTIO, AAT and CFHT to compare the mass in hot gas and stars to the dynamical mass and state of the clusters. Only 13 of the clusters are detected in X-ray emission, and for these systems we find that a range of 7-20 per cent of their baryonic mass, and <3 per cent of their dynamical mass, is detected in starlight, similar to what is observed in more massive clusters. In contrast, the five undetected clusters are underluminous in X-ray emission, by up to a factor 10, given their stellar mass. Although the velocity distribution of cluster members in these systems is indistinguishable from a Gaussian, all show subtle signs of being unrelaxed: either they lack a central, dominant galaxy, or the bright galaxy distribution is less concentrated and/or more elongated than the rest of the sample. Thus we conclude that low-mass clusters and groups selected from the velocity distribution of their galaxies exhibit a dichotomy in their hot gas properties. Either they are detected in X-ray, in which case they generally lie on the usual scaling relations, or they are completely undetected in X-ray emission. The non-detections may be partly related to the apparently young dynamical state of the clusters, but it remains a distinct possibility that some of these systems are exceptionally devoid of hot emitting gas as the result of its expulsion or rarefaction.

In this paper we generalize the flux ratio method Bailey et al. (2009) to the case of two luminosity indicators and search the optimal luminosity-flux ratio relations on a set of spectra whose phases are around not only the date of bright light but also other time. With these relations, a new method is proposed to constrain the host galaxy extinction of SN Ia and its distance. It is first applied to the low redshift supernovas and then to the high redshift ones. The results of the low redshift supernovas indicate that the flux ratio method can indeed give well constraint on the host galaxy extinction parameter E(B-V), but weaker constraints on R_{V}. The high redshift supernova spectra are processed by the same method as the low redshift ones besides some differences due to their high redshift. Among 16 high redshift supernovas, 15 are fitted very well except 03D1gt. Based on these distances, Hubble diagram is drew and the contents of the Universe are analyzed. It supports an acceleration behavior in the late Universe. Therefore, the flux ratio method can give constraints on the host galaxy extinction and supernova distance independently. We believe, through further studies, it may provide a precise tool to probe the acceleration of the Universe than before.

Gamma-ray bursts (GRBs) are the most luminous astrophysical events observed so far. They are conventionally classified into long and short ones depending on their time duration, $T_{90}$. Because of the advantage their high redshifts offer, many efforts have been made to apply GRBs to cosmology. The key to this is to find correlations between some measurable properties of GRBs and the energy or the luminosity of GRBs. These correlations are usually referred to as luminosity relations and are helpful in understanding the GRBs themselves. In this paper, we explored such correlations in the X-ray emission of GRBs. The X-ray emission of GRBs observed by \emph{Swift} has the exponential functional form in the prompt phase and relaxes to a power-law decay at time $T_p$. We have assumed a linear relation between $\log L_{X, p}$ (with $L_{X, p}$ being the X-ray luminosity at $T_p$) and $\log [T_p/(1+z)]$, but there is some evidence for curvature in the data and the true relationship between $L_{X, p}$ and $T_p/(1+z)$ may be a broken power law. The limited GRB sample used in our analysis is still not sufficient for us to conclude whether the break is real or just an illusion caused by outliers. We considered both cases in our analysis and discussed the implications of the luminosity relation, especially on the time duration of GRBs and their classification.

We present photometric, spectroscopic and weak lensing analysis of the large-scale structure and dynamics of the most X-ray luminous galaxy cluster known, RX J1347-1145, at z=0.451. We spectroscopically confirmed 47 new members with LDSS3 on the Magellan telescope. Together with previously known members, we measure a new velocity dispersion of 1163(+/-97) km/s. The mass inferred from our velocity dispersion is M200 = 1.16^{+0.32}_{-0.27}x10^{15} solar mass, with r200=1.85Mpc, under the assumption of a singular isothermal sphere. We also present a weak lensing analysis using deep CFHT data on this cluster, and find a deprojected mass of 1.47^{+0.46}_{-0.43}x10^{15} solar mass within r200, in excellent agreement with our dynamical estimate. Thus, our new dynamical mass estimate is consistent with that from weak lensing and X-ray studies in the literature, resolving a previously claimed discrepancy. We photometrically detect and spectroscopically confirm another massive cluster with sigma=780(+/-100) km/s and M200=3.4^{+1.4}_{-1.1}x10^{14} solar mass ~7Mpc south-west of RX J1347-1145, which we refer to as RXJ1347-SW. Our spectroscopic survey reveals a possible excess of galaxies in velocity space in the region between RX J1347-1145 and RXJ1347-SW; comparing with simulations, this excess appears consistent with that expected from a large filamentary structure traced by galaxies connecting these two clusters.

We study the optimal use of third order statistics in the analysis of weak lensing by large-scale structure. These higher order statistics have long been advocated as a powerful tool to break measured degeneracies between cosmological parameters. Using ray-tracing simulations, incorporating important survey features such as a realistic depth-dependent redshift distribution, we find that a joint two- and three-point correlation function analysis is a much stronger probe of cosmology than the skewness statistic. We compare different observing strategies, showing that for a limited survey time there is an optimal depth for the measurement of third-order statistics, which balances statistical noise and cosmic variance against signal amplitude. We find that the chosen CFHTLS observing strategy was optimal and forecast that a joint two- and three-point analysis of the completed CFHTLS-Wide will constrain the amplitude of the matter power spectrum $\sigma_8$ to 10% and the matter density parameter $\Omega_m$ to 17%, a factor of ~2.5 improvement on the two-point analysis alone. Our error analysis includes all non-Gaussian terms, finding that the coupling between cosmic variance and shot noise is a non-negligible contribution which should be included in any future analytical error calculations.

The large scatters of luminosity relations of gamma-ray bursts (GRBs) have been one of the most important reasons that prevent the extensive applications of GRBs in cosmology. In this paper, we extend the two-dimensional (2D) luminosity relations with $\tau_{\mathrm{lag}}$, $V$, $E_{\mathrm{peak}}$, and $\tau_{\mathrm{RT}}$ as the luminosity indicators to three dimensions (3D) using the same set of luminosity indicators to explore the possibility of decreasing the intrinsic scatters. We find that, for the 3D luminosity relations between the luminosity and an energy scale ($E_{\mathrm{peak}}$) and a time scale ($\tau_{\mathrm{lag}}$ or $\tau_{\mathrm{RT}}$), their intrinsic scatters are considerably smaller than those of corresponding 2D luminosity relations. Enlightened by the result and the definition of the luminosity (energy released in units of time), we discussed possible reasons behind, which may give us helpful suggestions on seeking more precise luminosity relations for GRBs in the future.

A long plateau phase and an amazing brightness have been observed in the Xray afterglow of GRB 060729. This peculiar light curve is likely due to long-term energy injection in external shock. Here we present a detailed numerical study on the energy injection process of magnetic dipole radiation from a strongly magnetized millisecond pulsar and model the multi-band afterglow observations. It is found that this model can successfully explain the long plateaus in the observed X-ray and optical afterglow light curves. The sharp break following the plateaus should be due to the rapid decline of the emission power of the central pulsar. At an even late time (~5*10^6s), an obvious jet break appears, which implies a relatively large half opening angle of theta~0.3 for the GRB ejecta. Due to the energy injection, the Lorentz factor of the outflow is still larger than two 10^7s post the GRB trigger, making the X-ray afterglow of this burst detectable by Chandra even 642 days after the burst.

We study a Chandra X-ray Observatory ACIS-S observation of the Galactic globular clusterM12. With a 26 ks exposure time, we detect 6 X-ray sources inside the half-mass radius (2.16 arcminutes) of which two are inside the core radius (0.72 arcminutes) of the cluster. If we assume these sources are all associated with globular cluster M12, the luminosity L_{X} among these sources between 0.3-7.0 keV varies roughly from 10^{30} to 10^{32} ergs s^{-1}. For identification, we also analyzed the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) and Wide Field and Planetary Camera 2 (WFPC2) data and identified the optical counterparts to five X-ray sources inside the HST ACS field of view. According to the X-ray and optical features, we found 2-5 candidate active binaries (ABs) or cataclysmic variables (CVs) and 0-3 background galaxies within the HST ACS field of view. Based on the assumption that the number of X-ray sources scales with the encounter rate and the mass of the globular cluster, we expect 2 X-ray source inside M12, and the expectation is consistent with our observational results. Therefore, the existence of identified X-ray sources (possible CVs or ABs) in M12 suggests the primordial origin of X-ray sources in globular clusters which is in agreement with previous studies.

In the fireball model, it is more physically realistic that gamma-ray burst (GRB) ejecta have a range of bulk Lorentz factors (assuming $M\propto \Gamma^{-s}$). The low Lorentz factor part of the ejecta will catch up with the high Lorentz factor part when the latter is decelerated by the surrounding medium to a comparable Lorentz factor. Such a process will develop a long-lasting weak reverse shock until the whole ejecta are shocked. Meanwhile, the forward shocked materials are gradually supplied with energy from the ejecta that are catching-up, and thus the temporal decay of the forward shock emission will be slower than that without an energy supply. However, the reverse shock may be strong. Here, we extend the standard reverse-forward shock model to the case of radially nonuniform ejecta. We show that this process can be classified into two cases: the thick shell case and the thin shell case. In the thin shell case, the reverse shock is weak and the temporal scaling law of the afterglow is the same as that in Sari & Mesz (2000). However, in the thick shell case, the reverse shock is strong and thus its emission dominates the afterglow in the high energy band. Our results also show slower decaying behavior of the afterglow due to the energy supply by low Lorentz factor materials, which may help the understanding of the plateau observed in the early optical and X-ray afterglows.

We calculate the diffuse high energy (TeV - PeV) neutrino emission from hyperflares of Soft-Gamma Repeaters (SGRs), like the hyperflare risen from \astrobj{SGR 1806-20} on December 27 of 2004, within the framework of the fireball model. The fireball model for gamma-ray bursts (GRBs) can explain well the main features of this hyperflare and the subsequent multi-frequency terglow emission. The expected rate, $\sim 20-100$ Gpc$^{-3}$day$^{-1}$, of such hyperflares is well in excess of the GRBs rate. Our result shows that the contribution to the diffuse TeV-PeV neutrino background from such hyperflares is less than 10% of the contribution from GRBs. We also discuss the high energy cosmic rays (CRs) from these sources.

Extensive multi-band afterglow data are available for GRB 980703. Especially, its radio afterglow was very bright and was monitored until more than 1000 days after the trigger time. Additionally, there is no obvious special feature, i.e., no rebrightenings, no plateau, and no special steep decay or slow decay in the multi-band afterglow light curves. All these conditions make GRB 980703 a precious sample in gamma-ray burst research. Here we use the observational data of GRB 980703 to test the standard fireball model in depth. It is found that the model can give a satisfactory explanation to the multi-band and overall afterglow light curves. The beaming angle of GRB 980703 is derived as ~0.23 radian, and the circum-burst medium density is ~ 27 cm-3. The total isotropic equivalent kinetic energy of the ejecta is ~ 3.8E52 ergs. A rest-frame extinction of Av ~ 2.5 mag in the host galaxy is also derived.

The global star formation rate has decreased significantly since z ~ 1, for reasons that are not well understood. Red-sequence galaxies, dominating in galaxy clusters, represent the population that have had their star formation shut off, and may therefore be the key to this problem. In this work, we select 127 rich galaxy clusters at 0.17<z<0.36, from 119 square degrees of the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) optical imaging data, and construct the r'-band red-sequence luminosity functions (LFs). We show that the faint end of the LF is very sensitive to how red-sequence galaxies are selected, and an optimal way to minimise the contamination from the blue cloud is to mirror galaxies on the redder side of the colour-magnitude relation (CMR). The LFs of our sample have a significant inflexion centred at Mr' ~- 18.5, suggesting a mixture of two populations. Combining our survey with low redshift samples constructed from the Sloan Digital Sky Survey, we show that there is no strong evolution of the faint end of the LF (or the red-sequence dwarf-to-giant ratio) over the redshift range 0.2 < z < 0.4, but from z ~ 0.2 to z ~ 0 the relative number of red-sequence dwarf galaxies has increased by a factor of ~3, implying a significant build-up of the faint end of the cluster red-sequence over the last 2.5 Gyr.

We quantify the Fisher information content of the cosmic shear survey two-point function as a function of noise and resolution. The two point information of dark matter saturates at the trans-linear scale. We investigate the impact of non-linear non-Gaussianity on the information content for lensing, which probes the same dark matter. To do so we heavily utilize N-body simulations in order to probe accurately the non-linear regime. While we find that even in a perfect survey, there is no clear saturation scale, we observe that non-linear growth induced non-Gaussianity could lead to a factor of ~4 reduction for the common Dark Energy figure of merit. This effect is however mitigated by realistic levels of shot noise and we find that for future surveys, the effect is closer to a factor of 1.5. To do so, we develop a new scheme to compute the relevant covariant matrix. It leads us to claim an unbiased estimator with an order of magnitude improvement in accuracy with only twice more simulations than previously used. Finally, we evaluate the error on the errors using bootstrap methods.

Wide area Large-Scale Structure (LSS) surveys are planning to map a substantial fraction of the visible universe to quantify dark energy through Baryon Acoustic Oscillations (BAO). At increasing redshift, for example that probed by proposed 21-cm intensity mapping surveys, gravitational lensing potentially limits the fidelity (Hui et al., 2007) because it distorts the apparent matter distribution. In this paper we show that these distortions can be reconstructed, and actually used to map the distribution of intervening dark matter. The lensing information for sources at z=1-3 allows accurate reconstruction of the gravitational potential on large scales, l <~ 100, which is well matched for Integrated Sachs-Wolfe (ISW) effect measurements of dark energy and its sound speed, and a strong constraint for modified gravity models of dark energy. We built an optimal quadratic lensing estimator for non-Gaussian sources, which is necessary for LSS. The phenomenon of "information saturation" (Rimes & Hamilton, 2005) saturates reconstruction at mildly non-linear scales, where the linear source power spectrum Delta^2 ~ 0.2-0.5. We find that steeper power spectra, saturate more quickly. We compute the effective number densities of independent lensing sources for LSS lensing, and find that they increase rapidly with redshifts. For LSS/21-cm sources at z ~ 2-4, the lensing reconstruction is limited by cosmic variance at l<~ 100.

Motivated by the fact that both SNe Ia and GRBs seem to prefer a dark energy EOS greater than -1 at redshifts $z \gtrsim 0.5$, we perform a careful investigation on this situation. We find that the deviation of dark energy from the cosmological constant at redshifts $z \gtrsim 0.5$ is large enough that we should pay close attention to it with future observational data. Such a deviation may arise from some biasing systematic errors in the handling of SNe Ia and/or GRBs or more interestingly from the nature of the dark energy itself.

In this paper, we revisit the weak gravitational lensing of CMB. The key point is that the deflection angle of CMB photon varies with the redshift. From this we obtain a set of modified power spectra. Compared with the ones calculated before, the temperature, the polarization E mode and the cross power spectra correction are of the order $0.01 %$, while the polarization B mode correction is about $0.4 %$ (to $l=2000$).

We investigate the spectrum of vector modes today which is generated at second order by density perturbations. The vector mode background that is generated by structure formation is small but in principle it contributes to the integrated Sachs-Wolfe effect, to redshift-space distortions and to weak lensing. We recover, clarify and extend previous results, and explain carefully why no vorticity is generated in the fluid at second order. The amplitude of the induced vector mode in the metric is around 1% that of the first-order scalars on small scales. We also calculate the power spectrum and the energy density of the vector part of the shear at second order.

The usual inflationary scenario predicts a Gaussian random primordial density fluctuation, different Fourier modes of which do not correlate with each other. In this paper we propose a correlation between these different modes. A simple case is that these different Fourier modes correlate with each other following a Gaussian function. For such a primordial density fluctuation we calculate the CMB angular power spectrum and find that its amplitude decreases but the decrease is different for different $l$. This feature can be used to constrain the the correlation strength from the real data.

Long-term precise timing of Galactic millisecond pulsars holds great promise for measuring the long-period (months-to-years) astrophysical gravitational waves. Several gravitational-wave observational programs, called Pulsar Timing Arrays (PTA), are being pursued around the world. Here we develop a Bayesian algorithm for measuring the stochastic gravitational-wave background (GWB) from the PTA data. Our algorithm has several strengths: (1) It analyses the data without any loss of information, (2) It trivially removes systematic errors of known functional form, including quadratic pulsar spin-down, annual modulations and jumps due to a change of equipment, (3) It measures simultaneously both the amplitude and the slope of the GWB spectrum, (4) It can deal with unevenly sampled data and coloured pulsar noise spectra. We sample the likelihood function using Markov Chain Monte Carlo (MCMC) simulations. We extensively test our approach on mock PTA datasets, and find that the algorithm has significant benefits over currently proposed counterparts. We show the importance of characterising all red noise components in pulsar timing noise by demonstrating that the presence of a red component would significantly hinder a detection of the GWB. Lastly, we explore the dependence of the signal-to-noise ratio on the duration of the experiment, number of monitored pulsars, and the magnitude of the pulsar timing noise. These parameter studies will help formulate observing strategies for the PTA experiments.

The behavior of the dark energy equation of state (EOS) is crucial in distinguishing different cosmological models. With a model independent approach, we constrain the possible evolution of the dark energy EOS. Gamma-ray bursts (GRBs) of redshifts up to $z>6$ are used, in addition to type Ia supernovae (SNe Ia). We separate the redshifts into 4 bins and assume a constant EOS parameter for dark energy in each bin. The EOS parameters are decorrelated by diagonalizing the covariance matrix. And the evolution of dark energy is estimated out of the uncorrelated EOS parameters. By including GRB luminosity data, we significantly reduce the confidence interval of the uncorrelated EOS parameter whose contribution mostly comes from the redshift bin of $0.5<z<1.8$. At high redshift where we only have GRBs, the constraints on the dark energy EOS are still very weak. However, we can see an obvious cut at about zero in the probability plot of the EOS parameter, from which we can infer that the ratio of dark energy to matter most probably continues to decrease beyond redshift 1.8. We carried out analyses with and without including the latest BAO measurements, which themselves favor a dark energy EOS of $w<-1$. If they are included, the results show some evidence of an evolving dark energy EOS. If not included, however, the results are consistent with the cosmological constant within $1 \sigma$ for redshift $0<z \lesssim 0.5$ and $2 \sigma$ for $0.5 \lesssim z<1.8$.

We explore the gamma-ray bursts' (GRBs') contributions in constraining the dark energy equation of state (EOS) at high ($1.8 < z < 7$) and at middle redshifts ($0.5 < z < 1.8$) and estimate how many GRBs are needed to get substantial constraints at high redshifts. We estimate the constraints with mock GRBs and mock type Ia supernovae (SNe Ia) for comparisons. When constraining the dark energy EOS in a certain redshift range, we allow the dark energy EOS parameter to vary only in that redshift bin and fix EOS parameters elsewhere to -1. We find that it is difficult to constrain the dark energy EOS beyond the redshifts of SNe Ia with GRBs unless some new luminosity relations for GRBs with smaller scatters are discovered. However, at middle redshifts, GRBs have comparable contributions with SNe Ia in constraining the dark energy EOS.

Aims. The high-redshift (z = 4:048) gamma-ray burst GRB 060206 showed unusual behavior, with a significant rebrightening by a factor of ~ 4 at about 3000 s after the burst. We argue that this rebrightening implies that the central engine became active again after the main burst produced by the first ejecta, then drove another more collimated jet-like ejecta with a larger viewing angle. The two ejecta both interacted with the ambient medium, giving rise to forward shocks that propagated into the ambient medium and reverse shocks that penetrated into the ejecta. The total emission was a combination of the emissions from the reverse- and forward- shocked regions. We discuss how this combined emission accounts for the observed rebrightening. Methods. We apply numerical models to calculate the light curves from the shocked regions, which include a forward shock originating in the first ejecta and a forward-reverse shock for the second ejecta. Results. We find evidence that the central engine became active again 2000 s after the main burst. The combined emission produced by interactions of these two ejecta with the ambient medium can describe the properties of the afterglow of this burst. We argue that the rapid rise in brightness at ~ 3000 s in the afterglow is due to the off-axis emission from the second ejecta. The precession of the torus or accretion disk of the central engine is a natural explanation for the departure of the second ejecta from the line of sight.

Since the inhomogeneous instrument noise can produce extra non-Gaussianity in the CMB anisotropy, its effect should be carefully subtracted in the primordial non-Gaussianity estimation. We calculate the probability distribution function of the CMB anisotropy for local type of non-Gaussianity, from which the optimal estimator in the general case (inhomogeneous noise and cut sky) is obtained. The new estimator obtained here is different from the popular one, since the inhomogeneous noise and cut sky effects are completely accounted. The CMB anisotropy in the new estimator is noise weighted. The noise weight is different from that used by WMAP Group in their 5-year data analysis. Although it is still difficult to calculate the new estimator rigorously, for the case of the slightly inhomogeneous noise, there exists a series expansion method to compute the new estimator. Each order in the series is suppressed by two factors, $(\frac{\sigma^2}{\sigma^2_{i}}-1)$ and $\frac{C_l}{C^{tot}_{l}}$, which make the method feasible. Through the Edgeworth expansion we can generalize our discussion to other types of non-Gaussianity.

The recent discovery of apparent cosmic acceleration has highlighted the depth of our ignorance of the fundamental properties of nature. It is commonly assumed that the explanation for acceleration must come from a new form of energy dominating the cosmos - dark energy - or a modification of Einstein's theory of Relativity. It is often overlooked, however, that a currently viable alternative explanation of the data is radial inhomogeneity which alters the Hubble diagram without any acceleration. This explanation is often ignored for two reasons: radial inhomogeneity significantly complicates analysis and predictions, and so the full details have not been investigated; and it is a philosophically highly controversial idea, revoking as it does the long-held Copernican Principle. To date, there has not been a general way of determining the validity if the Copernican Principle -- that we live at a typical position in the universe -- significantly weakening the foundations of cosmology as a scientific endeavour. Here we present an observational test for the Copernican assumption which can be automatically implemented while we search for dark energy in the coming decade. Our test is entirely independent of any model for dark energy or theory of gravity and thereby represents a model-independent test of the Copernican Principle.

While vector modes are usually ignored in cosmology since they are not produced during inflation they are inevitably produced from the interaction of density fluctuations of differing wavelengths. This effect may be calculated via a second-order perturbative expansion. We investigate this effect during the radiation era. We discuss the generation mechanism by investigating two scalar modes interacting, and we calculate the power of vector modes generated by a power-law spectrum of density perturbations on all scales.

We study the limits of accuracy for weak lensing maps of dark matter using diffuse 21-cm radiation from the pre-reionization epoch using simulations. We improve on previous "optimal" quadratic lensing estimators by using shear and convergence instead of deflection angles. We find that non-Gaussianity provides a limit to the accuracy of weak lensing reconstruction, even if instrumental noise is reduced to zero. The best reconstruction result is equivalent to Gaussian sources with effectively independent cell of side length 2.0/h Mpc. Using a source full map from z=10-20, this limiting sensitivity allows mapping of dark matter at a Signal-to-Noise ratio (S/N) greater than 1 out to l < 6000, which is better than any other proposed technique for large area weak lensing mapping.

Recent galaxy redshift surveys have brought in a large amount of accurate cosmological data out to redshift 0.3, and future surveys are expected to achieve a high degree of completeness out to a redshift exceeding 1. Consequently, a numerical programme for determining the metric of the universe from observational data will soon become practical; and thereby realise the ultimate application of Einstein's equations. Apart from detailing the cosmic geometry, this would allow us to verify and quantify homogeneity, rather than assuming it, as has been necessary up to now, and to do that on a metric level, and not merely at the mass distribution level. This paper is the beginning of a project aimed at such a numerical implementation. The primary observational data from our past light cone consists of galaxy redshifts, apparent luminosities, angular diameters and number densities, together with source evolution functions, absolute luminosities, true diameters and masses of sources. Here we start with the simplest case, that of spherical symmetry and a dust equation of state, and execute an algorithm that determines the unknown metric functions from this data. We discuss the challenges of turning the theoretical algorithm into a workable numerical procedure, particularly addressing the origin and the maximum in the area distance. Our numerical method is tested with several artificial data sets for homogeneous and inhomogeneous models, successfully reproducing the original models. This demonstrates the basic viability of such a scheme. Although current surveys don't have sufficient completeness or accuracy, we expect this situation to change in the near future, and in the meantime there are many refinements and generalisations to be added.

The discovery of multiband afterglows definitely shows that most $\gamma$-ray bursts are of cosmological origin. $\gamma$-ray bursts are found to be one of the most violent explosive phenomena in the Universe, in which astonishing ultra-relativistic motions are involved. In this article, the multiband observational characteristics of $\gamma$-ray bursts and their afterglows are briefly reviewed. The standard model of $\gamma$-ray bursts, i.e. the fireball model, is described. Emphasis is then put on the importance of the nonrelativistic phase of afterglows. The concept of deep Newtonian phase is elaborated. A generic dynamical model that is applicable in both the relativistic and nonrelativistic phases is introduced. Based on these elaborations, the overall afterglow behaviors, from the very early stages to the very late stages, can be conveniently calculated.

We analyze several recently detected gamma-ray bursts (GRBs) with late X-ray flares in the context of late internal shock and late external shock models. We find that the X-ray flares in GRB 050421 and GRB 050502B originate from late internal shocks, while the main X-ray flares in GRB 050406 and GRB 050607 may arise from late external shocks. Under the assumption that the central engine has two periods of activities, we get four basic types of X-ray light curves. The classification of these types depends on which period of activities produces the prompt gamma-ray emission (Type 1 and Type 2: the earlier period; Type 3 and Type 4: the late period), and on whether the late ejecta catching up with the early ejecta happens earlier than the deceleration of the early ejecta (Type 1 and Type 3) or not (Type 2 and Type 4). We find that the X-ray flare caused by a late external shock is a special case of Type 1. Our analysis reveals that the X-ray light curves of GRBs 050406, 050421, and 050607 can be classified as Type 1, while the X-ray light curve of GRB 050502B is classified as Type 2. However, the X-ray light curve of GRB 050406 is also likely to be Type 2. We also predict a long-lag short-lived X-ray flare caused by the inner external shock, which forms when a low baryon-loading long-lag late ejecta decelerates in the non-relativistic tail of an outer external shock driven by an early ejecta.

Gamma-ray burst afterglows have been observed for months or even years in a few cases. It deserves noting that at such late stages, the remnants should have entered the deep Newtonian phase, during which the majority of shock-accelerated electrons will no longer be highly relativistic. To calculate the afterglows, we must assume that the electrons obey a power-law distribution according to their kinetic energy, not simply the Lorentz factor.

Soft gamma-ray repeaters are at determined distances and their positions are known accurately. If observed, afterglows from their soft gamma-ray bursts will provide important clues to the study of the so called "classical gamma-ray bursts". On applying the popular fireball/blastwave model of classical gamma-ray bursts to soft gamma-ray repeaters, it is found that their X-ray and optical afterglows are detectable. Monitoring of the three repeaters is solicited.

Polarization behaviors of optical afterglows from two-component gamma-ray burst jets are investigated, assuming various configurations for the two components. In most cases, the observed polarization is dominated by the inner narrow component for a long period. Interestingly, it is revealed that different assumptions about the lateral expansion of the jet can lead to different evolutions of the position angle of polarization. The observed afterglow light curve and polarization behaviors of GRB 020813 can be well explained by the two-component jet model. Particularly, the model is able to explain the constancy of the observed position angle in this event, given that the line of sight is slightly outside the narrow component.

(abridged) Afterglow light curves are constructed analytically for realistic gamma-ray burst remnants decelerating in either a homogeneous interstellar medium or a stellar wind environment, taking into account the radiative loss of the blast wave, which affects the temporal behaviors significantly. Inverse Compton scattering is considered. The inverse Compton effect prolongs the fast-cooling phase markedly, during which the relativistic shock is semi-radiative and the radiation efficiency is approximately constant. It is further shown that the shock is still semi-radiative for quite a long time after it transits into the slow-cooling phase. The temporal decaying index of the X-ray afterglow light curve in this semi-radiative phase is more consistent with the observed $<\alpha_{X}> \sim 1.3$ than the commonly used adiabatic one. To manifest as a bump or even dominant in the X-ray afterglows during the relativistic stage, it is required that the density should be larger than about 1-10 cm$^{-3}$ in the interstellar medium case, or the wind parameter $A_{\ast}$ should be larger than unity in the stellar wind case.

The overall dynamical evolution and radiation mechanism of $\gamma$-ray burst jets are briefly introduced. Various interesting topics concerning beaming in $\gamma$-ray bursts are discussed, including jet structures, orphan afterglows and cylindrical jets. The possible connection between $\gamma$-ray bursts and neutron star kicks is also addressed.

Optical afterglows from two-component jets under various configurations are investigated numerically. Generally, the light curve is characterized by a rapid rebrightening when the observer is off-axis with respect to the narrow component, with the amplitude and peak time depending on detailed parameters. We further show that the optical afterglow of XRF 030723, especially its notable and rapid rebrightening, can be well explained by a typical two-component jet. This X-ray flash, together with GRB 030329, strongly hints the two-component jet model as a unified picture for X-ray flashes and gamma-ray bursts. With a narrow but ultra-relativistic inner outflow and a wide but less energetic outer ejecta, a two-component jet will be observed as a typical gamma-ray burst if our line of sight is within the angular scope of the narrow outflow. Otherwise, if the line of sight is within or slightly beyond the cone of the wide component, an X-ray flash will be detected.

It has been suggested that electromagnetic cascade of very high energy gamma-rays from gamma-ray bursts (GRBs) in the infrared/microwave background can produce delayed MeV-GeV photons. This delay could be caused by the angular spreading effect of the scattered microwave photons or deflection of the secondly pairs due to intergalactic magnetic field. Very high energy TeV photons of GRBs could be produced by a few mechanisms including the proton-synchrotron radiation and electron inverse Compton emission from GRB internal shocks as well as external shocks. We suggest that the information provided by the delayed emission could give constraints on models for TeV gamma-rays. A more accurate calculation of the delayed time caused by the angular spreading effect is presented by considering recent observations of the extragalactic infrared background and the theoretic high-redshift infrared background. We also suggest that the dependence of the maximum time delay of scattered photons on their energies, if determined by future GLAST detector, could differentiate the two mechanisms causing the time delay.

In many GRBs, afterglows have been observed for months or even years. It deserves noting that at such late stages, the remnants should have entered the deep Newtonian phase, during which the majority of shock-accelerated electrons will no longer be highly relativistic. However, a small portion of electrons are still ultra-relativistic and capable of emitting synchrotron radiation. Under the assumption that the electrons obey a power-law distribution according to their kinetic energy (not simply the Lorentz factor), we calculate optical afterglows from both isotropic fireballs and beamed ejecta, paying special attention to the late stages. In the beamed cases, it is found that the light curves are universally characterized by a flattening during the deep Newtonian phase. Implication of our results on orphan afterglows is also addressed.

The shock model has successfully explained the observed behaviors of afterglows from long gamma-ray bursts (GRBs). Here we use it to investigate the so-called early afterglows from short GRBs, which arises from blast waves that are not decelerated considerably by their surrounding medium. We consider a nearby medium loaded with $e^{\pm}$ pairs (Beloborodov 2002). The temporal behaviors show first a soft-to-hard spectral evolution, from the optical to hard X-ray, and then a usual hard-to-soft evolution after the blast waves begin to decelerate. The light curves show variability, and consist of two peaks. The first peak, due to the pair effect, can be observed in the X-ray, though too faint and too short in the optical. The second peak will be easily detected by {\it Swift}. We show that detections of the double-peak structure in the light curves of early afterglows are very helpful to determine all the shock parameters of short GRBs, including both the parameters of the relativistic source and the surroundings. Besides, from the requirement that the forward-shock emission in short GRBs should be below the BATSE detection threshold, we give a strong constraint on the shock model parameters. In particular, the initial Lorentz factor of the source is limited to be no more than $\sim 10^3$, and the ambient medium density is inferred to be low, $n\la 10^{-1}$ cm$^{-3}$.

Gamma-ray bursts (GRBs) are believed to originate from ultra-relativistic winds/fireballs to avoid the "compactness problem". However, the most energetic photons in GRBs may still suffer from $\gamma-\gamma$ absorption leading to electron/positron pair production in the winds/fireballs. We show here that in a wide range of model parameters, the resulting pairs may dominate those electrons associated with baryons. Later on, the pairs would be carried into a reverse shock so that a shocked pair-rich fireball may produce a strong flash at lower frequencies, i.e. in the IR band, in contrast with optical/UV emission from a pair-poor fireball. The IR emission would show a 5/2 spectral index due to strong self-absorption. Rapid responses to GRB triggers in the IR band would detect such strong flashes. The future detections of many IR flashes will infer that the rarity of prompt optical/UV emissions is in fact due to dust obscuration in the star formation regions.

We present the first optimal power spectrum estimation and three dimensional deprojections for the dark and luminous matter and their cross correlations. The results are obtained using a new optimal fast estimator (Pen 2003) deprojected using minimum variance and SVD techniques. We show the resulting 3-D power spectra for dark matter and galaxies, and their covariance for the VIRMOS-DESCART weak lensing shear and galaxy data. The survey is most sensitive to nonlinear scales k_NL ~ 1 h Mpc^-1. On these scales, our 3-D power spectrum of dark matter is in good agreement with the RCS 3-D power spectrum found by (Hoekstra et al 2002). Our galaxy power is similar to that found by the 2MASS survey, and larger than that of SDSS, APM and RCS, consistent with the expected difference in galaxy population. We find an average bias b=1.24+/-0.18 for the I selected galaxies, and a cross correlation coefficient r=0.75+/-0.23. Together with the power spectra, these results optimally encode the entire two point information about dark matter and galaxies, including galaxy-galaxy lensing. We address some of the implications regarding galaxy halos and mass-to-light ratios. The best fit ``halo'' parameter h=r/b=0.57+/-0.16, suggesting that dynamical masses estimated using galaxies systematically underestimate total mass. Ongoing surveys, such as the Canada-France-Hawaii-Telescope-Legacy survey will significantly improve on the dynamic range, and future photometric redshift catalogs will allow tomography along the same principles.

The optical flash accompanying GRB 990123 is believed to be powered by the reverse shock of a thin shell. With the best fitted physical parameters for GRB 990123 (Panaitescu & Kumar 2001) and the assumption that the parameters in the optical flash are the same as those in the afterglow, we show that: 1) the shell is thick but not thin, and we have provided the light curve for the thick shell case which coincides with the observation; 2) the theoretical peak flux of the optical flash accounts for only 3\times10^{-4} of the observed. In order to compensate this divergency, the physical parameters the electron energy ration and the magnetic ratio $\epsilon_e$, $\epsilon_B$ should be 0.61, 0.39 respectively, which are much different from those in the late afterglow.

The idea that gamma-ray bursts might be a kind of phenomena associated with neutron star kicks was first proposed by Dar & Plaga (1999). Here we study this mechanism in more detail and point out that the neutron star should be a high speed one (with proper motion larger than $\sim 1000$ km/s). It is shown that the model agrees well with observations in many aspects, such as the energetics, the event rate, the collimation, the bimodal distribution of durations, the narrowly clustered intrinsic energy, and the association of gamma-ray bursts with supernovae and star forming regions. We also discuss the implications of this model on the neutron star kick mechanism, and suggest that the high kick speed were probably acquired due to the electromagnetic rocket effect of a millisecond magnetar with an off-centered magnetic dipole.

The interaction of a relativistic fireball with its ambient medium is described through two shocks: a reverse shock that propagates into the fireball, and a forward shock that propagates into the medium. The observed optical flash of GRB 990123 has been considered to be the emission from such a reverse shock. The observational properties of afterglows suggest that the progenitors of some GRBs may be massive stars and their surrounding media may be stellar winds. We here study very early afterglows from the reverse and forward shocks in winds. An optical flash mainly arises from the relativistic reverse shock while a radio flare is produced by the forward shock. The peak flux densities of optical flashes are larger than 1 Jy for typical parameters, if we do not take into account some appropriate dust obscuration along the line of sight. The radio flare always has a long lasting constant flux, which will not be covered up by interstellar scintillation. The non-detections of optical flashes brighter than about 9th magnitude may constrain the GRBs isotropic energies to be no more than a few $10^{52}$ ergs and wind intensities to be relatively weak.

Large-scale, decelerating, relativistic X-ray jets due to material ejected from the black-hole candidate X-ray transient and microquasar XTE J1550-564 has been recently discovered with Chandra by Corbel et al. (2002). We find that the dynamical evolution of the eastern jet at the late time is consistent with the well-known Sedov evolutionary phase. A trans-relativistic external shock dynamic model by analogy with the evolution of gamma-ray burst remnants, is shown to be able to fit the observation data reasonably well. The inferred interstellar medium density around the source is well below the canonical value $n_ISM \sim 1 cm^{-3}$. We find that the emission from the continuously shocked interstellar medium (forward shock region) decays too slowly to be a viable mechanism for the eastern X-ray jet. However, the rapidly fading X-ray emission can be interpreted as synchrotron radiation from the non-thermal electrons in the adiabatically expanding ejecta. These electrons were accelerated by the reverse shock (moving back into the ejecta) which becomes important when the inertia of the swept external matter leads to an appreciable slowing down of the original ejecta. To ensure the dominance of the emission from the shocked ejecta over that from the forward shock region during the period of the observations, the magnetic field and electron energy fractions in the forward shock region must be far below equipartition. Future continuous, follow-up multi-wavelength observations of new ejection events from microquasars up to the significant deceleration phase should provide more valuable insight into the nature of the interaction between the jets and external medium.

The absorption feature detected in the prompt X-ray emission of GRB990705 has important consequences for its circum-burst environment and therefore on its afterglows. Here we investigate whether the circum-burst environment constrained by the absorption feature could be consistent with the observed H-band afterglow, which exhibits an earlier power law decay but a much faster decay about one day after the burst. Two possible geometries of the afterglow-emitting regions are suggested: 1) afterglow emission produced by the impact of the fireball on the surrounding torus, which serves as the absorbing material of the X-ray feature, as would be expected in the models involving that a supernova explosion precedes the gamma-ray burst by some time; 2)afterglow emission produced in the dense circum-burst medium inside the torus. In case 1), the faster decay at the later time is attributed to the disappearance of the shock due to the counter-pressure in the hot torus illuminated by the burst and afterglow photons. For case 2), the circum-burst medium density is found to be very high ($n\ga 10^4-10^5 {\rm cm^{-3}}$) if the emitting plasma is a jet or even higher if it is spherical. Future better observations of afterglows of GRBs that have absorption features might make it possible to make a more definite choice between these two scenarios.

Synchrotron self-Compton (SSC) process in the reverse shocks of gamma-ray bursts is suggested to be responsible for the observed prompt high-energy gamma-ray emissions from several gamma-ray bursts. We find that the SSC emission from the reverse shocks dominates over other emission processes in energy bands from tens of MeV to tens of GeV, for a wide range of shock parameters. This model is favorable for escape of energetic photons from the emitting regions due to a lower internal pair-production optical depth, as the characteristic size of the reverse shock region is much larger than that of internal shocks. We predict that, in this model, the prompt high-energy emissions are correlated with the prompt optical flashes, which can be test in the forthcoming GLAST era.

In previous works, a generic dynamical model has been suggested by Huang et al., which is shown to be correct for both adiabatic and radiative blastwaves, and in both ultra-relativistic and non-relativistic phases. In deriving their equations, Huang et al. have assumed that the radiative efficiency of the fireball is constant. They then applied their model directly to realistic cases where the radiative efficiency evolves with time. In this paper, we abandon the above assumption and re-derive a more accurate dynamical equation for gamma-ray burst remnants. Numerical results show that Huang et al.'s model is accurate enough in general cases.

Gamma-ray bursts may occur within pulsar wind bubbles (PWBs) under a number of scenarios, such as the supranova-like models in which the progenitor pulsar drives a powerful wind shocking against the ambient medium before it comes to death and produces a fireball. We here study the early afterglow emission from GRBs expanding into such a PWB environment. Different from the usual cold GRB external medium, the PWBs consist of a hot electron-positron medium with typical 'thermal' Lorentz factor of the order of gamma_w, the Lorentz factor of the pulsar particle wind. After GRB blast waves shock these hot electron-positron pairs, they will emit synchrotron radiation peaking at GeV bands. It is shown that GeV photons suffer negligible absorption by the soft photons radiation field in PWBs. Thus, strong GeV emissions in the early afterglow phases are expected, providing a plausible explanation for the long-duration GeV emission from GRB940217 detected by EGRET. Future GLAST may have the potential to test this GRB-PWB interaction model.

It is widely believed that multiwavelength afterglows of gamma-ray bursts (GRBs) originate from relativistic blast waves. We here show that in such blast waves, a significant fraction of the energy of shock-accelerated protons would be lost due to pion production by interactions with afterglow photons. This could lead to long-term production of $10^{16}$--$10^{18}$ eV neutrinos and sub-TeV $\gamma$-rays that accompany with usual afterglows, provided that the protons are accelerated to $10^{19}$ eV in the blast waves.

Based on a refined generic dynamical model, we investigate afterglows from jetted gamma-ray burst (GRB) remnants numerically. In the relativistic phase, the light curve break could marginally be seen. However, an obvious break does exist at the transition from the relativistic phase to the non-relativistic phase, which typically occurs at time 10 to 30 days. It is very interesting that the break is affected by many parameters, especially by the electron energy fraction (xi_e), and the magnetic energy fraction (xi_B^2). Implication of orphan afterglow surveys on GRB beaming is investigated. The possible existence of a kind of cylindrical jets is also discussed.

We generally analyze prompt high-energy emission above a few hundreds of GeV due to synchrotron self-Compton scattering in internal shocks. However, such photons cannot be detected because they may collide with cosmic infrared background photons, leading to electron/positron pair production. Inverse-Compton scattering of the resulting electron/positron pairs off cosmic microwave background photons will produce delayed MeV-GeV emission, which may be much stronger than a typical high-energy afterglow in the external shock model. We expand on the Cheng & Cheng model by deriving the emission spectrum and duration in the standard fireball shock model. A typical duration of the emission is ~ 10^3 seconds, and the time-integrated scattered photon spectrum is nu^{-(p+6)/4}, where p is the index of the electron energy distribution behind internal shocks. This is slightly harder than the synchrotron photon spectrum, nu^{-(p+2)/2}. The lower energy property of the scattered photon spectrum is dependent on the spectral energy distribution of the cosmic infrared background radiation. Therefore, future observations on such delayed MeV-GeV emission and the higher-energy spectral cutoff by the Gamma-Ray Large Area Space Telescope (GLAST) would provide a probe of the cosmic infrared background radiation.

It is believed that orphan afterglow searches can help to measure the beaming angle in gamma-ray bursts (GRBs). Great expectations have been put on this method. We point out that the method is in fact not as simple as we originally expected. Due to the baryon-rich environment that is common to almost all popular progenitor models, there should be many failed gamma-ray bursts, i.e., fireballs with Lorentz factor much less than 100 -- 1000, but still much larger than unity. In fact, the number of failed gamma-ray bursts may even be much larger than that of successful bursts. Owing to the existence of these failed gamma-ray bursts, there should be many orphan afterglows even if GRBs are due to isotropic fireballs, then the simple discovery of orphan afterglows never means that GRBs be collimated. Unfortunately, to distinguish a failed-GRB orphan and a jetted but off-axis GRB orphan is not an easy task. The major problem is that the trigger time is unknown. Some possible solutions to the problem are suggested.

Whether gamma-ray bursts are highly beamed or not is a very important question, since it has been pointed out that the beaming will lead to a sharp break in the afterglow light curves during the ultra-relativistic phase, with the breaking point determined by $\Gamma \sim 1/\theta_{0}$, where $\Gamma$ is the bulk Lorentz factor and $\theta_{0}$ is the initial half opening angle of the ejecta, and such a break is claimed to be present in the light curves of some GRBs. In this paper we will examine whether all the observed breaks in GRB afterglow light curves can be explained by jet effects. Here we present a detailed calculation of the jet evolution and emission, and have obtained a simple formula of bulk Lorentz factor evolution. We show that the light curves are very smoothly steepened by jet effect, and the shape of the light curve is determined by only one parameter -- $(E/n)^{1/8} \theta_{0}^{3/4}$, where E and n are the fireball energy and surrounding medium density. We find that for GRB990123 and GRB991216, the jet model can approximately fit their light curves, and the values of $(E_{51}/n_1)^{1/8} \theta_{0}^{3/4}$ are about 0.17 and 0.22 respectively. While for GRB990510, GRB000301c, GRB000926 and GRB010222, their light curves cannot be fitted by the jet model, which suggests that the breaks may be caused by some other reasons, jet should be not the unique reason.

We analyze in detail the hydrodynamics and afterglow emission of an ultrarelativistic blast wave when it expands in a density-jump medium. Such a medium is likely to appear in the vicinity of gamma-ray bursts (GRBs) associated with massive stars. The interaction of the blast wave with this medium is described through a reverse shock and a forward shock. We show that the reverse shock is initially relativistic if the factor of a density jump ($\alpha$) is much larger than 21, and Newtonian if $1<\alpha\ll 21$. We also calculate light curves of the afterglow emission during the interaction if the reverse shock is relativistic, and find that the optical flux density initially decays abruptly, then rises rapidly, and finally fades based on a power-law, which could be followed by an abrupt decay when the reverse shock has just crossed the originally swept-up matter. Therefore, one property of an afterglow occurring in a large-density-jump medium is an abrupt drop followed by a bump in the light curve and thus provides a probe of circumburst environments. In addition, this property could not only account for the optical afterglows of GRB 970508 and GRB 000301C but also explain the X-ray afterglow of GRB 981226.

We discuss the GRB afterglow emission from a relativistic jet that is initially in the radiative regime in which the accelerated electrons are fast cooling. We note that such a ``semiradiative'' jet decelerates faster than an adiabatic jet does. We also take into account the effect of strong inverse-Compton scattering on the cooling frequency in the synchrotron component and therefore on the light curve decay index. We find that there are two kinds of light-curve break for the jet effect. The first is an ``adiabatic break'' if the electrons become slow cooling before the jet enters a spreading phase, and the second is a ``radiative break'' which appears on the contrary case. We then show how a relativistic jet evolves dynamically and derive the overall temporal synchrotron emission in both cases, focusing on the change in light curve decay index around the break time. Finally, in view of our results, we rule out two cases for relativistic jets to account for the observed light curve breaks in a few afterglows: (i) an adiabatic jet with strong Compton cooling (Y>1) and with the cooling frequency nu_c locating in the observed energy range; (ii) a radiative jet with a significant fraction of total energy occupied by electrons (epsilon_e 1).

Sharp breaks have been observed in the afterglow light curves of several GRBs; this is generally explained by the jet model. However, there are still some uncertainties concerning this interpretation due to the unclear hydrodynamics of jet sideways expansion. Here we propose an alternative explanation to these observed breaks. If we assume that the multiwavelength spectra of GRB afterglows are not made of exact power law segments but their slope changes smoothly, i.e. $d\beta/d log\nu<0$, where $\beta$ is the spectral index, we find that this fact can very nicely explain the afterglow light curves showing breaks. Therefore we suggest that some breaks in the afterglow light curves may be caused by their curved spectra. The main feature of this interpretation is that the break time is dependent on the observed frequency, while the jet model produces achromatic breaks in the light curves. In addition, it is very important to know the position of the characteristic frequency $\nu_{c}$ in the multiwavelength spectrum at the time of the break, since it is a further discriminant between our model and the jet model. We find that although the optical light curves of seven GRB afterglows can be well fitted by the model we propose, in fact only one of them (i.e. GRB000926) can be explained in this framework, since for other ones the characteristic frequency $\nu_{c}$ is either above the optical after the break or below the optical before the break.

We calculate the spectra of inverse Compton (IC) emissions in gamma-ray burst (GRB) shocks produced when relativistic ejecta encounters the external interstellar medium, assuming a broken power-law approximation to the synchrotron seed spectrum. Four IC processes, including the synchrotron self-Compton (SSC) processes in GRB forward and reverse shocks, and two combined-IC processes (i.e. scattering of reverse shock photons on the electrons in forward shocks and forward shock photons on the electrons in reverse shocks), are considered. We find that the SSC emission from reverse shocks dominates over other emission processes in energy bands from tens of MeV to tens of GeV, for a wide range of shock parameters. This mechanism may be responsible for the prompt high energy gamma-rays detected by the Energetic Gamma Ray Experiment Telescope (EGRET). At TeV energy bands, however, the combined-IC emissions and/or the SSC emission from the forward shocks become increasingly dominant for a moderately steep distribution of shocked electrons.

The afterglow emission from a spreading jet expanding in a circumstellar cloud is discussed. Prompt X-ray radiation and a strong UV flash from the reverse shock produced by the interaction of the jet with the cloud may destroy and clear the dust out to about 30 pc within the initial solid angle of the jet. As the sideways expansion of the jet becomes significant, most of the optical radiation from the high-latitude part of the jet may be absorbed by the dust outside the initial solid angle of the jet, but only the radiation from the part within the initial solid angle can be observed. We analytically show that the flux of the observational radiation decays as $\propto t^{-(p+1)}$ (where $p$ is the power-law index of the electron distribution) in the relativistic phase. This preliminary result motivates us to perform numerical calculations. Our results show that one break in the optical afterglow ligh curve extends over a factor of $\sim 3$ in time rather than one decade in time in the previous jet model. These results may provide a way to judge whether GRBs locate in dense clouds or not. Finally, we carry out a detailed modelling for the R-band afterglow of GRB 000926.

Nearly all previous discussion on beaming effects in GRBs have assumed a conical geometry. However, more and more observations on relativistic jets in radio galaxies, active galactic nuclei, and "microquasars" in the Galaxy have shown that many of these outflows are not conical, but cylindrical, i.e., they maintain constant cross sections at large scales. Thus it is necessary to discuss the possibility that GRBs may be due to highly collimated cylindrical jets, not conical ones. Here we study the dynamical evolution of cylindrical jets and discuss their afterglows. Both analytical and numerical results are presented. It is shown that when the lateral expansion is not taken into account, a cylindrical jet typically remains to be highly relativistic for \sim 10^8 - 10^9 s. During this relativistic phase, the optical afterglow decays as \propto t^{-p/2} at first, where p is the index characterizing the power-law energy distribution of electrons. Then the light curve steepens to be \propto t^{-(p+1)/2} due to cooling of electrons. After entering the non-relativistic phase (i.e., t > 10^{11} s), the afterglow is \propto t^{-(5p-4)/6}. But if the cylindrical jet expands laterally at co-moving sound speed, then the decay becomes \propto t^{-p} and \propto t^{-(15p-21)/10} - t^{-(15p-20)/10} in the ultra-relativistic and non-relativistic phase respectively. Note that in both cases, the light curve turns flatter after the relativistic-Newtonian transition point, which differs markedly from the behaviour of a conical jet. It is suggested that some GRBs with afterglows decaying as t^{-1.1} - t^{-1.3} may be due to cylindrical jets, not necessarily isotropic fireballs.

Gamma-ray bursts (GRBs) are most intense transient gamma-ray events in the sky when they are on together with the strong evidences (i.e. the isotropic and inhomogeneous distribution of GRBs detected by BASTE) that they are located at cosmological distances, which make them the most energetic events ever known. For example, the observed radiation energies of some GRBs are equivalent to convert more than one solar mass energy into radiation completely. This is thousand times stronger than that of supernova explosion. Unconventional energy mechanisms and extremely high conversion efficiency for these mysterious events are required. The discovery of host galaxies and association with supernovae in the cosmological distances by the recently launched satellite of BeppoSAX and ground based radio and optical telescopes in GRB afterglow provides further support to the cosmological origin of GRBs and put strong constraints on central engines of GRBs. It is the aim of this article to review the possible central engines, energy mechanisms, dynamical and spectral evolution of GRBs, especially focusing on the afterglows in multi-wavebands.

There have been great and rapid progresses in the field of $\gamma$-ray bursts since BeppoSAX and other telescopes discovered their afterglows in 1997. In this talk, the main observational facts of $\gamma$-ray bursts and their afterglows, and the standard fireball shock model are reviewed briefly. And then, various post-standard effects, deviations from the standard model, are presented.

Currently popular models for progenitors of gamma-ray bursts (GRBs) are the mergers of compact objects and the explosions of massive stars. These two cases have distinctive environments for GRBs: compact object mergers occur in the interstellar medium (ISM) and the explosions of massive stars occur in the preburst stellar wind. We here discuss neutrino afterglows from reverse shocks as a result of the interaction of relativistic fireballs with their surrounding wind matter. After comparing with the analytical result of Waxman & Bahcall (2000) for the homogeneous ISM case, we find that the differential spectrum of neutrinos with energy from $\sim 3\times 10^{15}$ to $\sim 3\times 10^{17}$ eV in the wind case is softer by one power of the energy than in the ISM case. Furthermore, the expected flux of upward moving muons produced by neutrino interactions below a detector on the surface of the Earth in the wind case is $\sim 5$ events per year per km$^2$, which is about one order of magnitude larger than in the ISM case. In addition, these properties are independent of whether the fireballs are isotropic or beamed. Therefore, neutrino afterglows, if detected, may provide a way of distinguishing between GRB progenitor models based on the differential spectra of neutrinos and their event rates in a detector.

In the standard scenario of the fireball model of gamma-ray bursts(GRBs), the huge initial energy release produces a relativistic blast wave expanding into the external medium and a reverse shock moving into and heating the fireball ejecta. We calculate the high energy gamma-ray emission due to inverse Compton scattering of the synchrotron photons from relativistic electrons in the reverse shock. Under the favorable values of the physical parameters of the GRBs and the interstellar medium, our result shows that during the prompt phase, this emission dominates over the component from the forward shock at high energy gamma-ray bands. This mechanism can excellently account for the observations of the prompt high energy gamma-rays detected by EGRET, such as from GRB930131.

Gamma-ray bursts with long durations are widely thought to arise from the collapse of massive stars, where the wind environment is unavoidable. It is also believed that $\gamma$-ray bursts come from jets. Considering these two points in this paper, we calculate the evolution of a highly collimated jet that expands in a stellar wind environment and the expected afterglow from such a jet. We use a set of refined dynamical equations and a realistic lateral speed of the jet, and find: (1) There is no observable break at the time when the Lorentz factor of the jet is equal to the inverse of its initial half-opening angle. (2) No obvious break appears at the time when the blast wave transits from the relativistic to the non-relativistic phase. (3) For the wind case, there is no flattening tendency even up to $10^9$ s. (4) Compared with the homogeneous medium case, our calculated flux is weaker in the stellar wind case. Finally, we find that two kinds of GRB models (neutron star mergers and massive star collapses) may be discriminated in our numerical results.

We have constrained the intrinsic parameters, such as the magnetic energy density fraction ($\epsilon_{B}$), the electron energy density fraction ($\epsilon_e$), the initial Lorentz factor ($\Gamma_0$) and the Lorentz factor of the reverse external shock ($\Gamma_{rs}$), of GRB990123, in terms of the afterglow information (forward shock model) and the optical flash information (reverse shock model). Our result shows: 1) the inferred values of $\epsilon_e$ and $\epsilon_B$ are consistent with the suggestion that they may be universal parameters, comparing to those inferred for GRB970508; 2) the reverse external shock may have become relativistic before it passed through the ejecta shell. Other instrinsic parameters of GRB990123, such as energy contained in the forward shock $E$ and the ambient density $n$ are also determined and discussed in this paper.

There has been increasing evidence that at least some GRBs are emission beamed. The beamed GRB afterglow evolution has been discussed by several authors in the ultra-relativistic case. It has been shown that the dynamics of the blast wave will be significantly modified by the sideways expansion, and there may be a sharp break in the afterglow light curves under certain circumstances. However, this is true only when the fireball is still relativistic. Here we present an analytical approach to the evolution of the beamed GRB blast wave expanding in the surrounding medium (density $n\propto r^{-s}$) in the non-relativistic case, our purpose is to explore whether the sideways expansion will strongly affect the blast wave evolution as in the relativistic case. We find that the blast wave evolution is strongly dependent on the speed of the sideways expansion. If it expands with the sound speed, then the jet angle $\theta $ increases with time as $\theta \propto ln t$, which means that the sideways expansion has little effect on the afterglow light curves, the flux $F\propto t^{-\frac{3(5\alpha -1)}{5}}$ for $s=0$ and $F\propto t^{-\frac{7\alpha +1}{3}}$ for $s=2$.It is clear that the light curve of $s=2$ is not always steeper than that of $s=0$, as in the relativistic case.We also show that if the expansion speed is a constant, then the jet angle $\theta \propto t$, and the radius $r \propto t^{0}$, in this case the sideways expansion has the most significant effect on the blast wave evolution, the flux $F\propto t^{-(5\alpha -1)}$ independent of $s$, and we expect that there should be a smooth and gradual break in the light curve.

The fireball model of gamma-ray bursts predicted that when the energetic blast wave encountered the surrounding medium, there will be afterglow emission, and the subsequent afterglow observations appeared to confirm this prediction. In this simplest fireball model, the electrons have been accelerated to a power law energy distribution in a relativistic blast wave, then they give afterglow emission through synchrotron radiation. Up to now synchrotron radiation is believed to be the main mechanism of GRB emission, however, here we will show that under some circumstances, the inverse Compton scattering (ICS) may play an important role, and can change the light curves of GRB afterglows. Here we investigate the effects of ICS in the relativistic case (the surrounding medium density $\rho \propto r^{-2}$) and in the non-relativistic case (for both $\rho =constant$ and $\rho \propto r^{-2}$), we find that in the relativistic case the effect of ICS is usually important, while in the non-relativistic case, this effect is usually unimportant, unless the surrounding medium density is very high. We show that if ICS is important, then it can flatten and steepen the light curves of GRB afterglows, and this may provide the explanation for some afterglow observations.

Conversion from neutron stars to strange stars as a possible mechanism of cosmological gamma-ray bursts (GRBs) has been discussed in previous works, although the existence of strange stars is still an open question. On the basis of this mechanism, we here outline an explanation of the connection between supernovae (SNe) and GRBs, which has got increasing evidence recently. An asymmetric but normal SN explosion leaves a massive ($\geq1.8{\rm M_\odot}$) and rapidly rotating neutron star, which then converts to a strange star few days later, due to its rapid spindown. The accompanied fireball, which can be accelerated to ultra-relativistic velocity ($\Gamma_0\sim 100$) due to the very low baryon contamination of the strange star, flows out along the direction of the high-velocity SN jet and subsequently produces a GRB and the following low energy afterglows by interacting with the surrounding stellar wind. We will also expect a very luminous supernova like SN1998bw, if a large fraction of the conversion energy finally turns into the kinetic energy of the supernova ejecta.

GRB{980519} is characterized by its rapidly declining optical and X-ray afterglows. Explanations of this behavior include models invoking a dense medium environment which makes the shock wave evolve quickly into the sub-relativistic phase, a jet-like outflow, and a wind-shaped circumburst medium environment. Recently, Frail {et al}. (1999a) found that the latter two cases are consistent with the radio afterglow of this burst. Here, by considering the trans-relativistic shock hydrodynamics, we show that the dense medium model can also account for the radio light curve quite well. The potential virtue of the dense medium model for GRB{980519} is that it implies a smaller angular size of the afterglow, which is essential for interpreting the strong modulation of the radio light curve. Optical extinction due to the dense medium is not important if the prompt optical-UV flash accompanying the $\gamma$-ray emission can destroy dust by sublimation out to an appreciable distance. Comparisons with some other radio afterglows are also discussed.

Due to some refinements in the dynamics, we can follow the overall evolution of a realistic jet numerically till its bulk velocity being as small as $\beta c \sim 10^{-3} c$. We find no obvious break in the optical light curve during the relativistic phase itself. However, an obvious break does exist at the transition from the relativistic phase to the non-relativistic phase, which typically occurs at time $t \sim 10^6$ --- $10^{6.5}$ s (i.e., 10 --- 30 d). The break is affected by many parameters, such as the electron energy fraction $\xi_{\rm e}$, the magnetic energy fraction $\xi_{\rm B}^2$, the initial half opening angle $\theta_0$, and the medium number density $n$. Increase of any of them to a large enough value will make the break disappear. Although the break itself is parameter-dependent, afterglows from jetted GRB remnants are uniformly characterized by a quick decay during the non-relativistic phase, with power law timing index $\alpha \geq 2.1$. This is quite different from that of isotropic fireballs, and may be of fundamental importance for determining the degree of beaming in $\gamma$-ray bursts observationally.

Whether gamma-ray bursts are highly beamed or not is a very difficult but important problem that we are confronted with. Some theorists suggest that beaming effect usually leads to a sharp break in the afterglow light curve during the ultra-relativistic phase, with the breaking point determined by $\gamma = 1 / \theta_0$, where $\gamma$ is the Lorentz factor of the blastwave and $\theta_0$ is the initial half opening angle of the ejecta, but numerical studies tend to reject the suggestion. We note that previous studies are uniformly based on dynamics that is not proper for non-relativistic blastwaves. Here we investigate the problem in more detail, paying special attention to the transition from the ultra-relativistic phase to the non-relativistic phase. Due to some crucial refinements in the dynamics, we can follow the overall evolution of a realistic jet till its velocity is as small as $\beta c \sim 10^{-3} c$. We find no obvious break in the optical light curve during the relativistic phase itself. However, an obvious break does appear at the transition from the relativistic phase to the Newtonian phase if the physical parameters involved are properly assumed. Generally speaking, the Newtonian phase is characterized by a sharp decay of optical afterglows, with the power law timing index $\alpha \sim 1.8$ - 2.1. This is due to the quick lateral expansion at this stage. The quick decay of optical afterglows from GRB 970228, 980326, and 980519, and the breaks in the optical light curves of GRB 990123 and 990510 may indicate the presence of highly collimated $\gamma$-ray burst ejecta.

Recent observations on the GRB 000301c afterglow show that three breaks appear in the R-band light curve, and in particular the decay slope at late times is as steep as -3.0. This unusual afterglow is clearly inconsistent with the standard afterglow shock model. Here we propose a non-standard model for the unusual R-band afterglow of GRB 000301c. In this model, an ultra-relativistic shock in a dense medium ("dirty environment")rapidly evolved to the non-relativistic phase in initial 1 day. During such a phase, the shock happened to be caught up with by two energetic shells ejected from the central engine at two different times, and the shock was refreshed, leading to two flattenings of the light curve. After each interaction between the shock and shell, the afterglow decayed as $\propto t^{-3.0}$ if the electron distribution index of the shocked medium, $p\approx 3.4$, derived from the optical spectrum. Therefore, this model can provide an excellent explanation for the flattening and steepening features of the GRB 000301c optical afterglow light curve. We further point out that the energy injection shells ejected from the central engine at later times may be material shells (e.g., in the massive star progenitor models related to black holes) or radiation shells (e.g., in the millisecond pulsar progenitor models).

Beaming of relativistic ejecta in GRBs has been postulated by many authors in order to reduce the total GRB energy, thus it is very important to look for the observational evidence of beaming. Rhoads (1999) has pointed out that the dynamics of the blast wave, which is formed when the beamed ejecta sweeping the external medium, will be significantly modified by the sideways expansion due to the increased swept up matter. He claimed that shortly after the bulk Lorentz factor ($\Gamma $) of the blast wave drops below the inverse of the initial opening angle ($\theta_{0}$) of the beamed ejecta, there will be a sharp break in the afterglow light curves. However, some other authors have performed numerical calculations and shown that the break of the light curve is weaker and much smoother than the one analytically predicted. In this paper we reanalyse the dynamical evolution of the jet blast wave, calculate the jet emission analytically, we find that the sharp break predicted by Rhoads will actually not exist, and for most cases the afterglow light curve will almost not be affected by sideways expansion unless the beaming angle is extremely small. We demonstrate that only when $\theta_{0}<0.1$, the afterglow light curves may be steepened by sideways expansion, and in fact there cannot be two breaks as claimed before. We have also constructed a simple numerical code to verify our conclusion.

Based on the refined dynamical model proposed by us earlier for beamed $\gamma$-ray burst ejecta, we carry out detailed numerical procedure to study those $\gamma$-ray bursts with rapidly fading afterglows (i.e., $\sim t^{-2}$). It is found that optical afterglows from GRB 970228, 980326, 980519, 990123, 990510 and 991208 can be satisfactorily fitted if the $\gamma$-ray burst ejecta are highly collimated, with a universal initial half opening angle $\theta_0 \sim 0.1$. The obvious light curve break observed in GRB 990123 is due to the relativistic-Newtonian transition of the beamed ejecta, and the rapidly fading optical afterglows come from synchrotron emissions during the mildly relativistic and non-relativistic phases. We strongly suggest that the rapid fading of afterglows currently observed in some $\gamma$-ray bursts is evidence for beaming in these cases.

In a recent paper (Dai & Lu 1999), we have proposed a simple model in which the steepening in the light curve of the R-band afterglow of the gamma-ray burst (GRB) 990123 is caused by the adiabatic shock which has evolved from an ultrarelativistic phase to a nonrelativistic phase in a dense medium. We find that such a model is quite consistent with observations if the medium density is about $3\times 10^6 {\rm cm}^{-3}$. Here we discuss this model in more details. In particular, we investigate the effects of synchrotron self absorption and energy injection. A shock in a dense medium becomes nonrelativistic rapidly after a short relativistic phase. The afterglow from the shock at the nonrelativistic stage decays more rapidly than at the relativistic stage. Since some models for GRB energy sources predict that a strongly magnetic millisecond pulsar may be born during the formation of GRB, we discuss the effect of such a pulsar on the evolution of the nonrelativistic shock through magnetic dipole radiation. We find that after the energy which the shock obtains from the pulsar is much more than the initial energy of the shock, the afterglow decay will flatten significantly. When the pulsar energy input effect disappears, the decay will steepen again. These features are in excellent agreement with the afterglows of GRB 980519, GRB 990510 and GRB 980326. Furthermore, our model fits very well all the observational data of GRB 980519 including the last two detections.

There have been great and rapid progresses in the field of $\gamma$-ray bursts (denoted as GRBs) since BeppoSAX and other telescopes discovered their afterglows in 1997. Here, we will first give a brief review on the observational facts of GRBs and direct understanding from these facts, which lead to the standard fireball model. The dynamical evolution of the fireball is discussed, especially a generic model is proposed to describe the whole dynamical evolution of GRB remnant from highly radiative to adiabatic, and from ultra-relativistic to non-relativistic phase. Then, Various deviations from the standard model are discussed to give new information about GRBs and their environment. In order to relax the energy crisis, the beaming effects and their possible observational evidences are also discussed in GRB's radiations.

Conventional dynamic model of gamma-ray burst remnants is found to be incorrect for adiabatic blastwaves during the non-relativistic phase. A new model is derived, which is shown to be correct for both radiative and adiabatic blastwaves during both ultra-relativistic and non-relativistic phase. Our model also takes the evolution of the radiative efficiency into account. The importance of the transition from the ultra-relativistic phase to the non-relativistic phase is stressed.

Recent observations show that the temporal decay of the R-band afterglow from GRB 990123 steepened about 2.5 days after the burst. We here propose a possible explanation for such a steepening: a shock expanding in a dense medium has undergone the transition from a relativistic phase to a nonrelativistic phase. We find that this model is consistent with the observations if the medium density is about $3\times 10^6 {\rm cm}^{-3}$. By fitting our model to the observed optical and X-ray afterglow quantitatively, we further infer the electron and magnetic energy fractions of the shocked medium and find these two parameters are about 0.1 and $2\times 10^{-8}$ respectively. The former parameter is near the equipartition value while the latter is about six orders of magnitude smaller than inferred from the GRB 970508 afterglow. We also discuss possibilities that the dense medium can be produced.

The conventional generic model is deemed to explain the dynamics of $\gamma$-ray burst remnants very well, no matter whether they are adiabatic or highly radiative. However, we find that for adiabatic expansion, the model could not reproduce the Sedov solution in the non-relativistic phase, thus the model needs to be revised. In the present paper, a new differential equation is derived. The generic model based on this equation has been shown to be correct for both radiative and adiabatic fireballs, and in both ultra-relativistic and non-relativistic phase.

The detection of delayed X-ray, optical and radio emission, "afterglow", associated with $\gamma$-ray bursts (GRBs) is consistent with fireball models, where the emission are produced by relativistic expanding blast wave, driven by expanding fireball at cosmogical distances. The emission mechanisms of GRB afterglow have been discussed by many authors and synchrotron radiation is believed to be the main mechanism. The observations show that the optical light curves of two observed gamma-ray bursts, GRB970228 and GRB GRB970508, can be described by a simple power law, which seems to support the synchrotron radiation explanation. However, here we shall show that under some circumstances, the inverse Compton scattering (ICS) may play an important role in emission spectrum and this may influence the temporal properties of GRB afterglow. We expect that the light curves of GRB afterglow may consist of multi-components, which depends on the fireball parameters.

Millisecond pulsars with strong magnetic fields may be formed through several processes, e.g. accretion-induced collapse of magnetized white dwarfs, merger of two neutron stars. During the birth of such a pulsar, an initial fireball available for a gamma-ray burst (GRB) may occur. We here study evolution of a postburst relativistic fireball with energy injection from the pulsar through magnetic dipole radiation, and find that the magnitude of the optical afterglow from this fireball first decreases with time, subsequently flattens, and finally declines again. This may provide a natural explanation for the behavior of the lightcurve of the afterglow of GRB970228 if this burst resulted from the birth of a strongly magnetic millisecond pulsar.

We here discuss a new model of $\gamma$-ray bursts (GRBs) based on differentially rotating strange stars. Strange stars in this model and differentially rotating neutron stars in the Klu\'zniak-Ruderman model can produce extremely relativistic, variable fireballs required by GRBs and then become millisecond pulsars. The effect of such pulsars on expansion of the postburst fireballs through magnetic dipole radiation is studied. We show that these two models can explain naturally not only various features of GRBs but also light curves of afterglows.

Postburst evolution of adiabatic fireballs that produce $\gamma$-ray bursts is studied. Emphasis has been put on the transition from the highly relativistic phase to the non-relativistic phase, which, according to our calculation, should happen much earlier than previously expected. The theoretical light curves decline a little sharper at the non-relativistic stage than at early times, but still can fit the observations well. However, disagreements are obvious when $t \ge 100$ d, implying a large initial energy $E_0$, a low interstellar density $n$, or the possible existence of a persistent energy source at the center of the fireball.

In this talk, after briefly reviewing some historical remarks concerning strange stars, the achievements in physics and dynamical behavior of strange stars are discussed. Especially, various observational effects in distinguishing strange stars from neutron stars such as mechanical effects, cooling effects, phase transition and related interesting phenomena are stressed.

A GRB afterglow has been commonly thought to be due to continuous deceleration of a postburst fireball. Many analytical models have made simplifications for deceleration dynamics of the fireball and its radiation property, although they are successful at explaining the overall features of the observed afterglows. We here propose a model for a GRB afterglow in which the evolution of a postburst fireball is in an intermediate case between the adiabatic and highly radiative expansion. In our model, the afterglow is both due to the contribution of the adiabatic electrons behind the external blastwave of the fireball and due to the contribution of the radiative electrons. In addition, this model can describe evolution of the fireball from the extremely relativistic phase to the non-relativistic phase. Our calculations show that the fireball will go to the adiabatic expansion phase after about a day if the accelerated electrons are assumed to occupy the total internal energy. In all cases considered, the fireball will go to the mildly relativistic phase about $10^4$ seconds later, and to the non-relativistic phase after several days. These results imply that the relativistic adiabatic model cannot describe the deceleration dynamics of the several-days-later fireball. The comparison of the calculated light curves with the observed results at late times may imply the presence of impulsive events or energy injection with much longer durations.

The postburst evolution of fireballs that produce $\gamma$-ray bursts is studied, assuming the expansion of fireballs to be adiabatic and relativistic. Numerical results as well as an approximate analytic solution for the evolution are presented. Due to adoption of a new relation among $t$, $R$ and $\gamma$ (see the text), our results differ markedly from the previous studies. Synchrotron radiation from the shocked interstellar medium is attentively calculated, using a convenient set of equations. The observed X-ray flux of GRB afterglows can be reproduced easily. Although the optical afterglows seem much more complicated, our results can still present a rather satisfactory approach to observations. It is also found that the expansion will no longer be highly relativistic about 4 days after the main GRB. We thus suggest that the marginally relativistic phase of the expansion should be investigated so as to check the afterglows observed a week or more later.

The afterglow of a gamma-ray burst (GRB) is commonly thought to be due to continuous deceleration of a relativistically expanding fireball in the surrounding medium. Assuming that the expansion of the fireball is adiabatic and that the density of the medium is a power-law function of shock radius, viz., $n_{ext}\propto R^{-k}$, we analytically study the effects of the first-order radiative correction and the nonuniformity of the medium on a GRB afterglow. We first derive a new relation among the observed time, the shock radius and the fireball's Lorentz factor: $t_\oplus=R/4(4-k)\gamma^2c$, and also derive a new relation among the comoving time, the shock radius and the fireball's Lorentz factor: $t_{co}=2R/(5-k)\gamma c$. We next study the evolution of the fireball by using the analytic solution of Blandford and McKee (1976). The radiation losses may not significantly influence this evolution. We further derive new scaling laws both between the X-ray flux and observed time and between the optical flux and observed time. We use these scaling laws to discuss the afterglows of GRB 970228 and GRB 970616, and find that if the spectral index of the electron distribution is $p=2.5$, implied from the spectra of GRBs, the X-ray afterglow of GRB970616 is well fitted by assuming $k=2$.

The unusal hard x-ray burster GRO J1744-28 recently discovered by the Compton Gamma-ray Observatory (GRO) can be modeled as a strange star with a dipolar magnetic field $\le 10^{11} $Gauss. When the accreted mass of the star exceeds some critical mass, its crust may break, resulting in conversion of the accreted matter into strange matter and release of energy. Subsequently, a fireball may form and expand relativistically outward. The expanding fireball may interact with the surrounding interstellar medium, causing its kinetic energy to be radiated in shock waves, producing a burst of x-ray radiation. The burst energy, duration, interval and spectrum derived from such a model are consistent with the observations of GRO J1744-28.

The detection of the delayed emission in X-ray, optical and radio band, i.e. the afterglow of $\gamma$-ray bursts (GRBs), suggestes that the sources of GRBs are likely to be at cosmological distances. Here we explore the interaction of a relativistic shell with a uniform interstellar medium (ISM) and obtain the exact solution of the evolution of $\gamma$-ray burst remnants, including the radiative losses. We show that in general the evolution of bulk Lorentz factor $\gamma$ satisfies $\gamma \propto t^{-\alpha_{t}}$ when $\gamma \gg 1$, here $\alpha_{t}$ is mainly in the range $9/22 \sim 3/8$, the latter corresponds to adiabatic expansion. So it is clear that adiabatic expansion is a good approximation even when radiative loss is considered. However, in fact, $\alpha_{t}$ is slightly larger than 3/8, which may have some effects on a detailed data analysis. Synchrotron-self-absorption is also calculated and it is demonstrated that the radio emission may become optically thin during the afterglow. Our solution can also apply to the nonrelativistic case ($\gamma \sim 1$), at that time the observed flux decrease more rapidly than that in the relativistic case.