The structure of geomagnetism is very complex, and there are still some problems left in magnetohydrodynamics (MHD) simulations. Recently, the macro spin model has been suggested. This is the idea that geomagnetism is described by interaction with many local dynamo elements (macro-spins). This model can reproduce some features of geomagnetism by solving equations of motion with only one variable. In this paper, we study this model to make more general by adding one variable. In this result, our model becomes possible for several things that are not treated in the previous study, for example, migration of the North (or South) Magnetic Pole, comparison with observed data of magnetic field distributions expressed in two directions, etc. Moreover, as a result of an application to the sun, we could reproduce periodic polarity reversals and the power index of the power spectrum, etc. In addition, we investigate the statistical properties of the pole migration.

We promote the precise gamma-ray observation project Gamma-Ray Astro-Imager with Nuclear Emulsion (GRAINE), which uses balloon-borne emulsion gamma-ray telescopes. The emulsion telescope realizes observations with high angular resolution, polarization sensitivity, and large aperture area in the 0.01--100 GeV energy region. Herein, we report the data analysis of emulsion tracks and the first demonstration of gamma-ray imaging via an emulsion telescope by using the flight data from the balloon experiment performed in 2015 (GRAINE 2015). The emulsion films were scanned by the latest read-out system for a total area of 41 m$^2$ in three months, and then the gamma-ray event selection was automatically processed. Millions of electron-pair events are accumulated in the balloon-borne emulsion telescope. The emulsion telescope detected signals from a calibration source (gamma rays from the interaction of cosmic rays with an aluminum plate) with a high significance during the balloon observation and created a gamma-ray image consistent with the source size and the expected angular resolution in the energy range of 100--300 MeV. The flight performance obtained in the GRAINE 2015 experiment proves that balloon-borne emulsion telescope experiments with larger area are feasible while maintaining expected imaging performance.

We study the accelerating present universe in terms of the time evolution of the equation of state $w(z)$ (redshift $z$) due to thawing and freezing scalar potentials in the quintessence model. The values of $dw/da$ and $d^2w/da^2$ at scale factor of $a = 1$ are associated with two parameters of each potential. For five type scalar potentials, the scalar fields $Q$ and $w$ as function of time $t$ and/or $z$ are numerically calculated under the fixed boundary condition of $w(z=0)=-1+\Delta$. The observational constraint $w_{obs}$ (Planck 2015) is imposed to test whether the numerical $w(z)$ is in $w_{obs}$. Some solutions show the thawing features in the freezing potentials. Mutually exclusive allowed regions in $dw/da$ vs. $d^2w/da^2$ diagram are obtained in order to identify the likely scalar potential and even the potential parameters by future observational tests.

We investigate how accurately phase space distribution functions (DFs) in galactic models can be reconstructed by a made-to-measure (M2M) method, which constructs $N$-particle models of stellar systems from photometric and various kinematic data. The advantage of the M2M method is that this method can be applied to various galactic models without assumption of the spatial symmetries of gravitational potentials adopted in galactic models, and furthermore, numerical calculations of the orbits of the stars cannot be severely constrained by the capacities of computer memories. The M2M method has been applied to various galactic models. However, the degree of accuracy for the recovery of DFs derived by the M2M method in galactic models has never been investigated carefully. Therefore, we show the degree of accuracy for the recovery of the DFs for the anisotropic Plummer model and the axisymmetric St\"{a}ckel model, which have analytic solutions of the DFs. Furthermore, this study provides the dependence of the degree of accuracy for the recovery of the DFs on various parameters and a procedure adopted in this paper. As a result, we find that the degree of accuracy for the recovery of the DFs derived by the M2M method for the spherical target model is a few percent, and more than ten percent for the axisymmetric target model.

The time variation of the equation of state $w$ for quintessence scenario with a scalar field as dark energy is studied up to the third derivative ($d^3w/da^3$) with respect to the scale factor $a$, in order to predict the future observations and specify the scalar potential parameters with the observables. The third derivative of $w$ for general potential $V$ is derived and applied to several types of potentials. They are the inverse power-law ($V=M^{4+\alpha}/Q^{\alpha}$), the exponential ($V=M^4\exp{(\beta M/Q)}$), the cosine ($V=M^4(\cos (Q/f)+1)$) and the Gaussian types ($V=M^4\exp(-Q^2/\sigma^2)$), which are prototypical potentials for the freezing and thawing models. If the parameter number for a potential form is $ n$, it is necessary to find at least for $n+2$ independent observations to identify the potential form and the evolution of the scalar field ($Q$ and $ \dot{Q} $). Such observations would be the values of $ \Omega_Q, w, dw/da. \cdots $, and $ dw^n/da^n$. Since four of the above mentioned potentials have two parameters, it is necessary to calculate the third derivative of $w$ for them to estimate the predict values. If they are tested observationally, it will be understood whether the dark energy could be described by the scalar field with this potential. Numerical analysis for $d^3w/da^3$ are made under some specified parameters in the investigated potentials. It becomes possible to distinguish the freezing and thawing modes by the accurate observing $dw/da$ and $d^2w/da^2$ in some parameters.

The time evolution of the equation of state $w$ for quintessence models with a scalar field as dark energy is studied up to the third derivative ($d^3w/da^3$) with respect to the scale factor $a$, in order to predict the future observations and specify the scalar potential parameters with the observables. The third derivative of $w$ for general potential $V$ is derived and applied to several types of potentials. They are the inverse power-law ($V=M^{4+\alpha}/Q^{\alpha}$), the exponential ($V=M^4\exp{(\beta M/Q)}$), the mixed ( $V=M^{4+\gamma}\exp{(\beta M/Q)}/Q^{\gamma}$), the cosine ($V=M^4(\cos (Q/f)+1)$) and the Gaussian types ($V=M^4\exp(-Q^2/\sigma^2)$), which are prototypical potentials for the freezing and thawing models. If the parameter number for a potential form is $ n$, it is necessary to find at least for $n+2$ independent observations to identify the potential form and the evolution of the scalar field ($Q$ and $ \dot{Q} $). Such observations would be the values of $ \Omega_Q, w, dw/da. \cdots $, and $ dw^n/da^n$. From these specific potentials, we can predict the $ n+1 $ and higher derivative of $w$ ; $ dw^{n+1}/da^{n+1}, \cdots$. Since four of the above mentioned potentials have two parameters, it is necessary to calculate the third derivative of $w$ for them to estimate the predict values. If they are tested observationally, it will be understood whether the dark energy could be described by the scalar field with this potential. At least it will satisfy the necessary conditions. Numerical analysis for $d^3w/da^3$ are made under some specified parameters in the investigated potentials, except the mixed one. It becomes possible to distinguish the potentials by the accurate observing $dw/da$ and $d^2w/da^2$ in some parameters.

The time variation of the equation of state ($w_Q$) for the dark energy is analyzed by the current values of parameters $\Omega_Q $, $w_Q $ and their time derivatives. In the future, detailed feature of the dark energy could be observed, so we have considered the second derivatives of $w_Q$ for two types of potential: One is an inverse power-law type ($V=M^{4+\alpha}/Q^{\alpha}$) and the other is an exponential one ($V=M^4\exp{(\beta M/Q)}$). The first derivative $dw_Q/da$ and the second derivative $d^2 w_Q/da^2$ for both potentials are derived. The first derivative is estimated by the observed two parameters $\Delta=w_Q+1$ and $\Omega_Q$, with the assuming for $Q_0$. In the limit $\Delta \rightarrow 0$, the first derivative is null and, under the tracker approximation, the second derivative also becomes null. For the inverse power potential $V=M^{4+\alpha}/Q^{\alpha}$, the observed first and second derivatives are used to determine the potential parameter $M$ and $\alpha$. For the potential of $V=M^4\exp{(\beta M/Q)}$, the second derivative is estimated by the observed parameters $\Delta$, $\Omega_Q$ and $dw_Q/da$.

A phase-space distribution function of the steady state in galaxy models that admits regular orbits overall in the phase-space can be represented by a function of three action variables. This type of distribution function in Galactic models is often constructed theoretically for comparison of the Galactic models with observational data as a test of the models. On the other hand, observations give Cartesian phase-space coordinates of stars. Therefore it is necessary to relate action variables and Cartesian coordinates in investigating whether the distribution function constructed in galaxy models can explain observational data. Generating functions are very useful in practice for this purpose, because calculations of relations between action variables and Cartesian coordinates by generating functions do not require a lot of computational time or computer memory in comparison with direct numerical integration calculations of stellar orbits. Here, we propose a new method called a torus-fitting method, by which a generating function is derived numerically for models of the Galactic potential in which almost all orbits are regular. We confirmed the torus-fitting method can be applied to major orbit families (box and loop orbits) in some two-dimensional potentials. Furthermore, the torus-fitting method is still applicable to resonant orbit families, besides major orbit families. Hence the torus-fitting method is useful for analyzing real Galactic systems in which a lot of resonant orbit families might exist.

The CALET Gamma-ray Burst Monitor (CGBM) is the secondary scientific instrument of the CALET mission on the International Space Station (ISS), which is scheduled for launch by H-IIB/HTV in 2014. The CGBM provides a broadband energy coverage from 7 keV to 20 MeV, and simultaneous observations with the primary instrument Calorimeter (CAL) in the GeV - TeV gamma-ray range and Advanced Star Camera (ASC) in the optical for gamma-ray bursts (GRBs) and other X-gamma-ray transients. The CGBM consists of two kinds of scintillators: two LaBr$_3$(Ce) (7 keV - 1 MeV) and one BGO (100 keV - 20 MeV) each read by a single photomultiplier. The LaBr$_3$(Ce) crystal, used in space for the first time here for celestial gamma-ray observations, enables GRB observations over a broad energy range from low energy X-ray emissions to gamma rays. The detector performance and structures have been verified using the bread-board model (BBM) via vibration and thermal vacuum tests. The CALET is currently in the development phase of the proto-flight model (PFM) and the pre-flight calibration of the CGBM is planned for August 2013. In this paper, we report on the current status and expected performance of CALET for GRB observations.

The probability is investigated that the meteorites originating on Earth are transferred to other planets in our Solar System and to extra solar planets. We take the collisional Chicxulub crater event, and material that was ejected as an example of Earth-origin meteors. If we assume the appropriate size of the meteorites as 1cm in diameter, the number of meteorites to reach the exoplanet system (further than 20 ly) would be much greater than one. We have followed the ejection and capture rates estimated by Melosh (2003) and the discussion by Wallis and Wickramasinghe (2004). If we consider the possibility that the fragmented ejecta (smaller than 1cm) are accreted to comets and other icy bodies, then buried fertile material could make the interstellar journey throughout Galaxy. If life forms inside remain viable, this would be evidence of life from Earth seeding other planets. We also estimate the transfer velocity of the micro-organisms in the interstellar space. In some assumptions, it could be estimated that, if life has originated $10^{10}$\ years ago anywhere in our Galaxy as theorized by Joseph and Schild (2010a, b), it will have since propagated throughout our Galaxy and could have arrived on Earth by 4.6 billion years ago. Organisms disperse.

Recently, it has been observed the extreme metal-poor stars in the Galactic halo, which must be formed just after Pop III objects. On the other hand, the first gas clouds of mass $\sim 10^6 M_{\odot}$ are supposed to be formed at $ z \sim $ 10, 20, and 30 for the $1\sigma$, $2\sigma $ and $3\sigma$, where the density perturbations are assumed of the standard $\Lambda$CDM cosmology. If we could apply this gaussian distribution to the extreme small probability, the gas clouds would be formed at $ z \sim $40, 60, and 80 for the $4\sigma$, $6\sigma$, and $8\sigma$. The first gas clouds within our galaxy must be formed around $z\sim 40$. Even if the gas cloud is metal poor, there is a lot of possibility to form the planets around such stars. The first planetary systems could be formed within $\sim 6\times 10^7$ years after the Big Bang in the universe. Even in our galaxies, it could be formed within $\sim 1.7\times 10^8$ years. It is interesting to wait the observations of planets around metal-poor stars. For the panspermia theory, the origin of life could be expected in such systems.

Detecting the first electron pairs with nuclear emulsion allows a precise measurement of the direction of incident gamma-rays as well as their polarization. With recent innovations in emulsion scanning, emulsion analyzing capability is becoming increasingly powerful. Presently, we are developing a balloon-borne gamma-ray telescope using nuclear emulsion. An overview and a status of our telescope is given.

We report the detection, with the CANGAROO-III imaging atmospheric Cherenkov telescope array, of a very high energy gamma-ray signal from the unidentified gamma-ray source HESS J1614-518, which was discovered in the H.E.S.S. Galactic plane survey. Diffuse gamma-ray emission was detected above 760 GeV at the 8.9 sigma level during an effective exposure of 54 hr from 2008 May to August. The spectrum can be represented by a power-law: 8.2+-2.2_{stat}+-2.5_{sys}x10^{-12}x (E/1TeV)^{-Gamma} cm^{-2} s^{-1} TeV^{-1} with a photon index Gamma of 2.4+-0.3_{stat}+-0.2_{sys}, which is compatible with that of the H.E.S.S. observations. By combining our result with multi-wavelength data, we discuss the possible counterparts for HESS J1614-518 and consider radiation mechanisms based on hadronic and leptonic processes for a supernova remnant, stellar winds from massive stars, and a pulsar wind nebula. Although a leptonic origin from a pulsar wind nebula driven by an unknown pulsar remains possible, hadronic-origin emission from an unknown supernova remnant is preferred.

SS433, located at the center of the supernova remnant W50, is a close proximity binary system consisting of a compact star and a normal star. Jets of material are directed outwards from the vicinity of the compact star symmetrically to the east and west. Non-thermal hard X-ray emission is detected from lobes lying on both sides. Shock accelerated electrons are expected to generate sub-TeV gamma rays through the inverse-Compton process in the lobes. Observations of the western X-ray lobe region of SS433/W50 system have been performed to detect sub-TeV gamma-rays using the 10m CANGAROO-II telescope in August and September, 2001, and July and September, 2002. The total observation times are 85.2 hours for ON source, and 80.8 hours for OFF source data. No significant excess of sub-TeV gamma rays has been found at 3 regions of the western X-ray lobe of SS433/W50 system. We have derived 99% confidence level upper limits to the fluxes of gamma rays and have set constraints on the strengths of the magnetic fields assuming the synchrotron/inverse-Compton model for the wide energy range of photon spectrum from radio to TeV. The derived lower limits are 4.3 microgauss for the center of the brightest X-ray emission region and 6.3 microgauss for the far end from SS433 in the western X-ray lobe. In addition, we suggest that the spot-like X-ray emission may provide a major contribution to the hardest X-ray spectrum in the lobe.

Because accretion and merger shocks in clusters of galaxies may accelerate particles to high energies, clusters are candidate sites for the origin of ultra-high-energy (UHE) cosmic-rays. A prediction was presented for gamma-ray emission from a cluster of galaxies at a detectable level with the current generation of imaging atmospheric Cherenkov telescopes. The gamma-ray emission was produced via inverse Compton upscattering of cosmic microwave background (CMB) photons by electron-positron pairs generated by collisions of UHE cosmic rays in the cluster. We observed two clusters of galaxies, Abell 3667 and Abell 4038, searching for very-high-energy gamma-ray emission with the CANGAROO-III atmospheric Cherenkov telescope system in 2006. The analysis showed no significant excess around these clusters, yielding upper limits on the gamma-ray emission. From a comparison of the upper limit for the north-west radio relic region of Abell 3667 with a model prediction, we derive a lower limit for the magnetic field of the region of ~0.1 micro G. This shows the potential of gamma-ray observations in studies of the cluster environment. We also discuss the flux upper limit from cluster center regions using a model of gamma-ray emission from neutral pions produced in hadronic collisions of cosmic-ray protons with the intra-cluster medium (ICM). The derived upper limit of the cosmic-ray energy density within this framework is an order of magnitude higher than that of our Galaxy.

We have reanalyzed data from observations of PSR B1706-44, SN 1006, and the Vela pulsar region made with the CANGAROO 3.8 m imaging atmospheric Cherenkov telescope between 1993 and 1998 in response to the results reported for these sources by the H.E.S.S. collaboration. In our reanalysis, in which gamma-ray selection criteria have been determined exclusively using gamma-ray simulations and OFF-source data as background samples, no significant TeV gamma-ray signals have been detected from compact regions around PSR B1706-44 or within the northeast rim of SN 1006. We discuss reasons why the original analyses gave the source detections. The reanalysis did result in a TeV gamma-ray signal from the Vela pulsar region at the 4.5 sigma level using 1993, 1994, and 1995 data. The excess was located at the same position, 0.13 deg. to the southeast of the Vela pulsar, as that reported in the original analysis. We have investigated the effect of the acceptance distribution in the field of view of the 3.8 m telescope, which rapidly decreases toward the edge of the field of the camera, on the detected gamma-ray morphology. The expected excess distribution for the 3.8 m telescope has been obtained by reweighting the distribution of HESS J0835-455 measured by H.E.S.S. with the acceptance of the 3.8 m telescope. The result is morphologically comparable to the CANGAROO excess distribution, although the profile of the acceptance-reweighted H.E.S.S. distribution is more diffuse than that of CANGAROO. The integral gamma-ray flux from HESS J0835-455 has been estimated for the same region as defined by H.E.S.S. from the 1993-1995 data of CANGAROO to be F(> 4.0 +/- 1.6 TeV) = (3.28 +/- 0.92) x 10^{-12} photons cm^{-2} s^{-1}, which is statistically consistent with the integral flux obtained by H.E.S.S.

Observation by the CANGAROO-III stereoscopic system of the Imaging Cherenkov Telescope has detected extended emission of TeV gamma rays in the vicinity of the pulsar PSR B1706$-$44. The strength of the signal observed as gamma-ray-like events varies when we apply different ways of emulating background events. The reason for such uncertainties is argued in relevance to gamma-rays embedded in the "off-source data", that is, unknown sources and diffuse emission in the Galactic plane, namely, the existence of a complex structure of TeV gamma-ray emission around PSR B1706$-$44.

We have observed Markarian 421 in January and March 2001 with the CANGAROO-II imaging Cherenkov telescope during an extraordinarily high state at TeV energies. From 14 hours observations at very large zenith angles, $\sim70^\circ$, a signal of 298 $\pm$ 52 gamma-ray--like events (5.7 $\sigma$) was detected at $E>10$ TeV, where a higher sensitivity is achieved than those of usual observations near the zenith, owing to a greatly increased collecting area. Under the assumption of an intrinsic power-law spectrum, we derived a differential energy spectrum $dN/dE = (3.3 \pm 0.9_{stat.} \pm 0.3_{syst.})\times10^{-13} (E/10 {Te V})^{-(4.0 ^{+0.9}_{-0.6}_{stat.} \pm 0.3_{syst.})}$ ph./cm$^2$/sec/TeV, which is steeper than those previously measured around 1 TeV, and supports the evidence for a cutoff in the spectrum of Markarian 421. However, the 4 $\sigma$ excess at energies greater than 20 TeV in our data favors a cutoff energy of $\sim$8 TeV, at the upper end of the range previously reported from measurements at TeV energies.

A preliminary result from Markarian 421 observations in the energy region above 10 TeV with the CANGAROO-II 10 m telescope is presented. In January 2001, the HEGRA group reported that Markarian 421 had become very active, with flux levels up to 4 times that of the Crab Nebula. As a result, we observed Mkn 421 during six nights from January 24th to February 1st, and four nights from March 1st to 4th. Observations were carried out using the very large zenith angle technique ($\sim$70 degree) and the energy threshold is estimated from Monte Carlo simulations to be around 10 TeV. We have detected gamma-ray emission in this energy range.

We report on the status of the construction of an array of four 10 m atmospheric Cherenkov telescopes for gamma-ray astronomy, near Woomera, in South Australia -- the CANGAROO-III project. The first telescope of this array is the upgraded version of the CANGAROO-II 7 m telescope and has been in operation since March 2000. The second telescope, an improved version of the first, is being constructed for installation in late 2001. Stereoscopic observation of sub TeV gamma-rays with the two 10 m telescopes will begin in 2002 and the full array will be operational in 2004.

The atmospheric Cerenkov imaging technique has been used to search for point-like and diffuse TeV gamma-ray emission from the southern supernova remnant, W28, and surrounding region. The search, made with the CANGAROO 3.8m telescope, encompasses a number of interesting features, the supernova remnant itself, the EGRET source 3EG J1800-2338, the pulsar PSR J1801-23, strong 1720 MHz OH masers and molecular clouds on the north and east boundaries of the remnant. An analysis tailored to extended and off-axis point sources was used, and no evidence for TeV gamma-ray emission from any of the features described above was found in data taken over the 1994 and 1995 seasons. Our upper limit (E>1.5 TeV) for a diffuse source of radius 0.25deg encompassing both molecular clouds was calculated at 6.64e-12 photons cm^-2 s^-1 (from 1994 data), and interpreted within the framework of a model predicting TeV gamma-rays from shocked-accelerated hadrons. Our upper limit suggests the need for some cutoff in the parent spectrum of accelerated hadrons and/or slightly steeper parent spectra than that used here (-2.1). As to the nature of 3EG J1800-2338, it possibly does not result entirely from pi-zero decay, a conclusion also consistent with its location in relation to W28.

The gamma-ray pulsar PSR B1509-58 and its surrounding nebulae have been observed with the CANGAROO 3.8m imaging atmospheric Cherenkov telescope. The observations were performed from 1996 to 1998 in Woomera, South Australia, under different instrumental conditions with estimated threshold energies of 4.5 TeV (1996), 1.9 TeV (1997) and 2.5 TeV (1998) at zenith angles of ~30 deg. Although no strong evidence of the gamma-ray emission was found, the lowest energy threshold data of 1997 showed a marginal excess of gamma-ray--like events at the 4.1 sigma significance level. The corresponding gamma-ray flux is calculated to be (2.9 +/- 0.7) * 10^{-12}cm^{-2}s^{-1} above 1.9 TeV. The observations of 1996 and 1998 yielded only upper limits (99.5% confidence level) of 1.9 * 10^{-12}cm^{-2}s^{-1} above 4.5 TeV and 2.0 * 10^{-12}cm^{-2}s^{-1} above 2.5 TeV, respectively. Assuming that the 1997 excess is due to Very High-Energy (VHE) gamma-ray emission from the pulsar nebula, our result, when combined with the X-ray observations, leads to a value of the magnetic field strength ~5 micro G. This is consistent with the equipartition value previously estimated in the X-ray nebula surrounding the pulsar. No significant periodicity at the 150ms pulsar period has been found in any of the three years' data. The flux upper limits set from our observations are one order of magnitude below previously reported detections of pulsed TeV emission.

We report the results of TeV gamma-ray observations of the shell type SNR RXJ1713.7-3946 (G347.3-0.5). The discovery of strong non-thermal X-ray emission from the northwest part of the remnant strongly suggests the existence of electrons with energies up to 100 TeV in the remnant, making the SNR a good candidate TeV gamma-ray source. We observed RXJ1713.7-3946 from May to August 1998 with the CANGAROO 3.8m atmospheric imaging Cerenkov telescope and obtained evidence for TeV gamma-ray emission from the NW rim of the remnant with the significance of 5.6 sigma. The observed TeV gamma-ray flux from the NW rim region was estimated to be (5.3 +/- 0.9[statistical] +/- 1.6[systematic]) * 10^{-12} photons cm^{-2} s^{-1} at energies >= 1.8 +/- 0.9 TeV. The data indicate that the emitting region is much broader than the point spread function of our telescope. The extent of the emission is consistent with that of hard X-rays observed by ASCA. This TeV gamma-ray emission can be attributed to the Inverse Compton scattering of the Cosmic Microwave Background Radiation by shock accelerated ultra-relativistic electrons. Under this assumption, a rather low magnetic field of 11 micro gauss is deduced for the remnant from our observation.

CANGAROO group constructed an imaging air Cherenkov telescope (CANGAROO-II) in March 1999 atWoomera, South Australia to observe celestial gamma-rays in hundreds GeV region. It has a 7m parabolic mirror consisting of 60 small plastic spherical mirrors, and the prime focus is equipped with a multi-pixel camera of 512 PMTs covering the field of view of 3 degrees. We report initial performance of the telescope.

The southern supernova remnant (SNR)W28 was observed in 1994 and 1995 by the CANGAROO 3.8m telescope in a search formulti-TeV gamma ray emission, using the Cerenkov imaging technique. We obtained upper limits for a variety of point-like and extended features within a +-1 degree-region and briefly discuss these results, together with that of EGRET within the framework of a shock acceleration model of the W28 SNR.

The CANGAROO-III project,which consists of an array of four 10 m imaging Cherenkov telescopes,has just started being constructed in Woomera, South Australia,in a collaboration between Australia and Japan. The first stereoscopic observation of celestial high-energy gamma-rays in the 100 GeV region with two telescopes will start in 2002,and the four telescope array will be completed in 2004. The concept of the project and the expected performance are discussed.

We have been successful in developing light and durable mirrors made of CFRP (Carbon Fiber Reinforced Plastic) laminates for the reflector of the new CANGAROO-II 7 m telescope. The reflector has a parabolic shape (F/1.1) with a 30 m^2 effective area which consists of 60 small spherical mirrors of CFRP laminates. The orientation of each mirror can be remotely adjusted by stepping motors. After the first adjustment work, the reflector offers a point image of about $0.^\circ 14$ (FWHM) on the optic axis.

The shell type SNR RXJ1713.7-3946 is a new SNR discovered by the ROSAT all sky survey. Recently, strong non-thermal X-ray emission from the northwest part of the remnant was detected by the ASCA satellite. This synchrotron X-ray emission strongly suggests the existence of electrons with energies up to hundreds of TeV in the remnant. This SNR is, therefore, a good candidate TeV gamma ray source, due to the Inverse Compton scattering of the Cosmic Microwave Background Radiation by the shock accelerated ultra-relativistic electrons, as seen in SN1006. In this paper, we report a preliminary result of TeV gamma-ray observations of the SNR RXJ1713.7-3946 by the CANGAROO 3.8m telescope at Woomera, South Australia.

CANGAROO group has constructed the new large imaging Air \v Cerenkov telescope to exploit hundred GeV region gamma-ray astronomy in March 1999 at Woomera, South Australia. It has a 7m parabolic mirror consisting of 60 small plastic spherical mirrors, and a fine imaging camera with 512 PMTs covering the field of view of 3 degree. Observation will start from July 1999.

The data acquisition system for the new CANGAROO-II 7m telescope is described.

We have searched for TeV gamma-rays from Centaurus A and surrounding region out to +/- 1.0 deg using the CANGAROO 3.8m telescope. No evidence for TeV gamma-ray emission was observed from the search region, which includes a number of interesting features located away from the tracking centre of our data. The 3 sigma upper limit to the flux of gamma-rays above 1.5 TeV from an extended source of radius 14' centred on Centaurus A is 1.28e-11 photons cm^-2 s^-1.

Despite extensive surveys of extragalactic TeV gamma-ray candidates only 3 sources have so far been detected. All three are northern hemisphere objects and all three are low-redshift X-ray selected BL Lacs (XBLs). In this paper we present the results of observations of the three nearest southern hemisphere XBLs (PKS0548-322, PKS2005-489 and PKS2155-304) with the CANGAROO 3.8m imaging telescope. During the period of observation we estimate that the threshold of the 3.8m telescope was around 1.5TeV. Searches for both steady and short timescale emission have been performed for each source. Additionally, we are able to monitor the X-ray state of each source on a daily basis and we have made contemporaneous measurements of optical activity for PKS0548-322 and PKS2155-304.

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 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.

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.

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.

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.

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.