Long-time behavior of spatial power-law correlation in one-dimensional self-gravitating system is numerically investigated. It is found that the power-law structure persists even after the system is virialized. The structure gradually disappears owing to energy exchange among particles. Energy distribution is quasi-stationary after the system is virialized.

A new phase of temporal evolution of the one-dimensional self-gravitating system is numerically discovered. Fractal structure is dynamically created from non-fractal initial conditions. Implication to astrophysics and mathematical physics is discussed.

The process of formation of fractal structure in one-dimensional self-gravitating system is examined numerically. It is clarified that structures created in small spatial scale grow up to larger scale through clustering of clusters, and form power-law correlation.

Observational and theoretical results indicate that low-redshift BL Lacertae objects are the most likely extragalactic sources to be detectable at TeV energies. In this paper we present the results of observations of 4 BL Lacertae objects (PKS0521-365, EXO0423.4-0840, PKS2005-489 and PKS2316-423) made between 1993 and 1996 with the CANGAROO 3.8m imaging Cherenkov telescope. During the period of these observations the gamma-ray energy threshold of the 3.8m telescope was ~2TeV. Searches for steady long-term emission have been made, and, inspired by the TeV flares detected from Mkn421 and Mkn501, a search on a night-by-night timescale has also been performed for each source. Comprehensive Monte Carlo simulations are used to estimate upper limits for both steady and short timescale emission.

This paper reports the first discovery of TeV gamma-ray emission from a supernova remnant made with the CANGAROO 3.8 m Telescope. TeV gamma rays were detected at the sky position and extension coincident with the north-east (NE) rim of shell-type Supernova remnant (SNR) SN1006 (Type Ia). SN1006 has been a most likely candidate for an extended TeV Gamma-ray source, since the clear synchrotron X-ray emission from the rims was recently observed by ASCA (Koyama et al. 1995), which is a strong evidence of the existence of very high energy electrons up to hundreds of TeV in the SNR. The observed TeV gamma-ray flux was $(2.4\pm 0.5(statistical) \pm 0.7(systematic)) \times 10^{-12}$ cm$^{-2}$ s$^{-1}$ ($\ge 3.0\pm 0.9$ TeV) and $ (4.6\pm 0.6 \pm 1.4) \times 10^{-12}$ cm$^{-2}$ s$^{-1}$ ($\ge 1.7\pm 0.5$ TeV) from the 1996 and 1997 observations, respectively. Also we set an upper limit on the TeV gamma-ray emission from the SW rim, estimated to be $ 1.1 \times 10^{-12}$ cm$^{-2}$ s$^{-1}$ ($\ge 1.7\pm 0.5$ TeV, 95% CL) in the 1997 data. The TeV gamma rays can be attributed to the 2.7 K cosmic background photons up-scattered by electrons of energies up to about 10$^{14}$ eV by the inverse Compton (IC) process. The observed flux of the TeV gamma rays, together with that of the non-thermal X-rays, gives firm constraints on the acceleration process in the SNR shell; a magnetic field of $6.5\pm2$ $\mu$G is inferred from both the synchrotron X-rays and inverse Compton TeV gamma-rays, which gives entirely consistent mechanisms that electrons of energies up to 10$^{14}$ eV are produced via the shock acceleration in SN1006.

Gamma rays with energies greater than 7 TeV from the Crab pulsar/nebula have been observed at large zenith angles, using the Imaging Atmospheric Technique from Woomera, South Australia. CANGAROO data taken in 1992, 1993 and 1995 indicate that the energy spectrum extends up to at least 50 TeV, without a change of the index of the power law spectrum. The observed differential spectrum is \noindent $(2.01\pm 0.36)\times 10^{-13}(E/{7 TeV})^{-2.53 \pm 0.18} TeV^{-1}cm^{-2}s^{-1}$ between 7 TeV and 50 TeV. There is no apparent cut-off. The spectrum for photon energies above $\sim$10 TeV allows the maximum particle acceleration energy to be inferred, and implies that this unpulsed emission does not originate near the light cylinder of the pulsar, but in the nebula where the magnetic field is not strong enough to allow pair creation from the TeV photons. The hard gamma-ray energy spectrum above 10 TeV also provides information about the varying role of seed photons for the inverse Compton process at these high energies, as well as a possible contribution of $\pi ^{\circ}$-gamma rays from proton collisions.

We have observed the Vela pulsar region at TeV energies using the 3.8 m imaging Cherenkov telescope near Woomera, South Australia between January 1993 and March 1995. Evidence of an unpulsed gamma-ray signal has been detected at the 5.8 sigma level. The detected gamma-ray flux is (2.9 +/- 0.5 +/- 0.4) x 10^-12 photons cm^-2 sec^-1 above 2.5 +/- 1.0 TeV and the signal is consistent with steady emission over the two years. The gamma-ray emission region is offset from the Vela pulsar position to the southeast by about 0.13 deg. No pulsed emission modulated with the pulsar period has been detected and the 95 % confidence flux upper limit to the pulsed emission from the pulsar is (3.7 +/- 0.7) x 10^-13 photons cm^-2 sec^-1 above 2.5 +/- 1.0 TeV.

The CANGAROO experiment has observed gamma-ray above 7TeV from the Crab pulsar/nebula at large zenith angle in Woomera, South Australia. We report the CANGAROO data taken in 1992, 1993 and 1995, from which it appears that the energy spectrum extends at least up to 50 TeV. The observed integral spectrum is (8.4+-1.0) x 10^{-13}(E/7 TeV)^(-1.53+-0.15)cm^{-2}s^{-1} between 7 TeV and 50 TeV. In November 1996, the 3.8m mirror was recoated in Australia, and its reflectivity was improved to be about 90% as twice as before. Due to this recoating, the threshold energy of ~4 TeV for gamma rays has been attained in the observation of the Crab at large zenith angle. Here we also report the preliminary result taken in 1996.

We report recent results of the CANGAROO Collaboration on very high energy gamma ray emission from pulsars, their nebulae, SNR and AGN in the southern sky. Observations are made in South Australia using the imaging technique of detecting atmospheric Cerenkov light from gamma rays higher than about 1 TeV. The detected gamma rays are most likely produced by the inverse Compton process by electrons which also radiate synchrotron X-rays. Together with information from longer wavelengths, our results can be used to infer the strength of magnetic field in the emission region of gamma rays as well as the energy of the progenitor electrons. A description of the CANGAROO project is also given, as well as details of the new telescope of 7 m diameter which is scheduled to be in operation within two years.

The data for VHE (~TeV) gamma rays from young gamma-ray pulsar PSR1509-58 observed in 1996 with the CANGAROO 3.8m Cerenkov imaging telescope are presented, as well as the additional data from March to June of 1997. The high spin-down luminosity of the pulsar and the plerionic feature around the pulsar observed with radio and X-rays suggest that VHE gamma-ray emission is quite likely above the sensitivity of the CANGAROO telescope. The CANGAROO results on other pulsars, such as PSR1259-63, are also presented. PSR1259-63 is a highly eccentric X-ray binary system, which includes a high mass Be companion star, and a preliminary analysis on the data taken 4 months after the periastron in 1994 suggests emission of VHE gamma rays.

We have observed the Vela pulsar region at TeV energies using the 3.8 m imaging Cherenkov telescope near Woomera, South Australia every year since 1992. This is the first concerted search for pulsed and unpulsed emission from the Vela region, and the imaging technique also allows the location of the emission within the field of view to be examined. A significant excess of gamma-ray-like events is found offset from the Vela pulsar to the southeast by about 0.13deg. The excess shows the behavior expected of gamma-ray images when the asymmetry cut is applied to the data. There is no evidence for the emission being modulated with the pulsar period -- in contrast to earlier claims of signals from the Vela pulsar direction.

Since 1992 the CANGAROO 3.8m imaging telescope has been used to search for sources of TeV gamma-rays. Results are presented here for observations of four Southern Hemisphere BL-Lacs - PKS0521-365, PKS2316-423, PKS2005-489 and EXO0423-084. In addition to testing for steady DC emission, a night by night burst excess search has been performed for each source.

This is the third paper of the series of our studies of the one-dimensional self-gravitating many-body systems. In this paper, we thus study the transition phenomena after the first transition from a quasiequilibrium. We found that irrespective of the initial conditions, the system wanders between many states which yield the same properties of the quasiequilibrium as the water-bag. This itinerancy is most prominent in the time scale $t\sim 10^5 \sim 10^6 t_c$. In the midway between two succeeding quasiequilibria, the system experiences a transient state, where one particle keeps exclusive high energy and its motion decouples with the others. Though the transient state is not the thermally relaxed equilibrium, its distribution quite resembles the isothermal distribution. Thus the macroscopic relaxation we discussed in the previous papers corresponds to the transition from one of quasiequilibria to a transient state. Averaging the behavior over a time scale much longer than the macroscopic relaxation time gives the isothermal distribution. Distribution of lifetime $\tau$ of the transient states yields a power-law distribution of $\tau^{-2}$. The transient state gives a clear example of chaotic itinerancy in conserved dynamical systems. The mechanism of the onset of itinerancy is examined.

Though one dimensional self-gravitating $N$-body systems have been studied for three decade, the nature of relaxation was still unclear. There were inconsistent results about relaxation time; some initial state relaxed in the time scale $T\sim N\,t_c$, but another state did not relax even after $T\sim N^2\,t_c$, where $t_c$ is the crossing time. The water-bag distribution was believed not to relax after $T\sim N^2\,t_c$. In our previous paper, however, we found there are two different relaxation times in the water-bag distribution;in the faster relaxation ( microscopic relaxation ) the equipartition of energy distribution is attains but the macroscopic distribution turns into the isothermal distribution in the later relaxation (macroscopic relaxation). In this paper, we investigated the properties of the two relaxation. We found that the microscopic relaxation time is $T\sim N\,t_c$, and the macroscopic relaxation time is proportional to $N\,t_c$, thus the water-bag does relax. We can see the inconsistency about the relaxation times is resolved as that we see the two different aspect of relaxations. Further, the physical mechanisms of the relaxations are presented.

The microscopic dynamics of one-dimensional self-gravitating many-body systems is studied. We examine two courses of the evolution which has the isothermal and stationary water-bag distribution as initial conditions. We investigate the evolution of the systems toward thermal equilibrium. It is found that when the number of degrees of freedom of the system is increased, the water-bag distribution becomes a quasi-equilibrium, and also the stochasticity of the system reduces. This results suggest that the phase space of the system is effectively not ergodic and the system with large degreees of freedom approaches to the near-integrable one.