Compilation of papers presented by the JEM-EUSO Collaboration at the 36th International Cosmic Ray Conference (ICRC), held July 24 through August 1, 2019 in Madison, Wisconsin.

The JVO ALMA WebQL web service - available through the JVO ALMA FITS archive - has been upgraded to include legacy data from other telescopes, for example Nobeyama NRO45M in Japan. The updated server software has been renamed FITSWebQL. In addition, a standalone desktop version supporting Linux, macOS and Windows 10 Linux Subsystem (Bash on Windows) is also available for download from http://jvo.nao.ac.jp/~chris/ . The FITSWebQL server enables viewing of even 100GB-large FITS files in a web browser running on a PC with a limited amount of RAM. Users can interactively zoom-in to selected areas of interest with the corresponding frequency spectrum being calculated on the server in near real-time. The client (a browser) is a JavaScript application built on WebSockets, HTML5, WebGL and SVG. There are many challenges when providing a web browser-based real-time FITS data cube preview service over high-latency low-bandwidth network connections. The upgraded version tries to overcome the latency issue by predicting user mouse movements with a Kalman Filter in order to speculatively deliver the real-time spectrum data at a point where the user is likely to be looking at. The new version also allows one to view multiple FITS files simultaneously in an RGB composite mode (NRO45M FUGIN only), where each dataset is assigned one RGB channel to form a colour image. Spectra from multiple FITS cubes are shown together too. The paper briefly describes main features of FITSWebQL. We also touch on some of the recent developments, such as an experimental switch from C/C++ to Rust (see https://www.rust-lang.org/) for improved stability, better memory management and fearless concurrency, or attempts to display FITS data cubes in the form of interactive on-demand video streams in a web browser.

We develop an infinite mixture model of Ornstein-Uhlenbeck(OU) processes for describing the optical variability of QSOs based on treating the variability as a stochastic process. This enables us to get the parameters of the power spectral densities(PSDs) on their brightness variations by providing more flexible description of PSDs than the models based on single OU process(damped random walk). We apply this model to 67,507 variable objects extracted from SDSS Stripe82 photometric data and succeed in showing very high precision in identifying QSOs (~99% levels in completeness and purity) among variable objects based only on their variability, by investigating on 9,855 spectroscopically confirmed objects(7,714 QSOs and 2,141 stars) in the data of SDSS Data Release 12(DR12), with sufficient and accurate multiple measurements of their brightness. By comparing our results with the values based on other models that are used in previous research, it is revealed that our model can be used as the most effective method for selecting QSOs from variable object catalog, especially regarding completeness and purity. The main reason of improved identification rates are the ability of our model to separate clearly QSOs and stars, especially on the small fraction of QSOs with variabilities which can be described better than simple damped random walk model.

EUSO-Balloon is a pathfinder mission for the Extreme Universe Space Observatory onboard the Japanese Experiment Module (JEM-EUSO). It was launched on the moonless night of the 25$^{th}$ of August 2014 from Timmins, Canada. The flight ended successfully after maintaining the target altitude of 38 km for five hours. One part of the mission was a 2.5 hour underflight using a helicopter equipped with three UV light sources (LED, xenon flasher and laser) to perform an inflight calibration and examine the detectors capability to measure tracks moving at the speed of light. We describe the helicopter laser system and details of the underflight as well as how the laser tracks were recorded and found in the data. These are the first recorded laser tracks measured from a fluorescence detector looking down on the atmosphere. Finally, we present a first reconstruction of the direction of the laser tracks relative to the detector.

We have measured the clustering of galaxies around active galactic nuclei (AGN) for which single-epoch virial masses of the super-massive black hole (SMBH) are available to investigate the relation between the large scale environment of AGNs and the evolution of SMBHs. The AGN samples used in this work were derived from the Sloan Digital Sky Survey (SDSS) observations and the galaxy samples were from 240~deg$^{2}$ S15b data of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). The investigated redshift range is 0.6--3.0, and the masses of the SMBHs lie in the range $10^{7.5}$--$10^{10}$~M$_{\odot}$. The absolute magnitude of the galaxy samples reaches to $M_{\lambda 310}$ $\sim$ $-$18 at rest frame wavelength 310~nm for the low-redshift end of the samples. More than 70\% of the galaxies in the analysis are blue. We found a significant dependence of the cross-correlation length on redshift, which primarily reflects the brightness dependence of the galaxy clustering. At the lowest redshifts the cross-correlation length increases from 7~$h^{-1}$~Mpc around $M_{\lambda 310}$$ = $$-19$~mag to $>$10~$h^{-1}$~Mpc beyond $M_{\lambda 310}$$ = $$-20$~mag. No significant dependence of the cross-correlation length on BH mass was found for whole galaxy samples dominated by blue galaxies, while there was an indication of BH mass dependence in the cross-correlation with red galaxies. These results provides a picture of the environment of AGNs studied in this paper being enriched with blue starforming galaxies, and a fraction of the galaxies are being evolved to red galaxies along with the evolution of SMBHs in that system.

We present results of the clustering analysis between active galactic nuclei (AGNs) and galaxies at redshift 0.1-1.0 for investigating properties of galaxies associated with the AGNs, revealing the nature of fueling mechanism of supermassive black holes (SMBHs). We used 8059 SDSS AGNs/QSOs for which virial masses of individual SMBHs were measured, and divided them into four mass groups. Cross-correlation analysis was performed and bias for each mass group was derived. The averaged color and luminosity distributions of galaxies around the AGNs/QSOs were also derived for each mass group. The galaxy color was estimated for SED constructed from a merged SDSS and UKIDSS catalog. The distributions of color and luminosity were derived by the subtraction method, which does not require redshift information of galaxies. The main results of this work are: (1) a bias increases by a factor two from the lower mass group to the highest mass group; (2) the environment around AGNs with the most massive SMBH (Mbh > 10^9 Msun) is dominated by red sequence galaxies; (3) marginal indication of decline in luminosity function at dimmer side of M > -19.5 mag is found for galaxies around AGNs with Mbh = 10^8.2 - 10^9 Msun and nearest redshift group (z=0.1-0.3). These results indicate that AGNs with the most massive SMBHs reside in haloes where large fraction of galaxies have been transited to the red sequence. The accretion of hot halo gas as well as recycled gas from evolving stars can be the most plausible mechanism to fuel the SMBHs above ~10^9 Msun.

The JVO ALMA Archive provides users one of the easiest ways to access the ALMA archival data. The users can have a quick look at a 3 or 4-dimensional data cube without downloading multiple huge tarballs from a science portal of ALMA Regional Centers (ARCs). Since we just synchronize all datasets with those of ARCs, the metadata are identical to the upstream, including ``target name'' for each dataset. The name is not necessarily a common one like NGC numbers, but sometimes one of sequential numbers assigned in an observation proposal. Compilation of these artificial names into astronomical ones could provide users more flexible and powerful search interfaces; for instance, with the knowledge of the redshift for each source, the users can easily find the datasets which observed their interested emission/absorption lines at not the observer frame but the rest frame, fitting well with theoretical studies. To implement this functionality, cross-identification of all the sources in our archive with those in some other astronomical databases such as NED and SIMBAD is required. We developed a tiny Java application named ``Blade Runner'' for this purpose. The program works as a crawler for both the JVO ALMA Archive and SIMBAD, storing all information onto a SQLite-based database file; this portable design enables us to communicate results to each other even under different computing environments. In this paper, we introduce its software design and our recent work on the application, and report a preliminary result on the source identification in our archive.

We present results of the cross-correlation analysis between active galactic nuclei (AGNs) and galaxies at redshift 0.1-1. We obtain data of ~10,000 SDSS AGNs in which their virial masses with a super-massive black hole (SMBH) were estimated. The UKIDSS galaxy samples around the AGNs were obtained using the virtual observatory. The scale length of AGN-galaxy cross-correlation for all of the samples is measured to be r_0= 5.8+0.8-0.6 h^-1Mpc (for the fixed slope parameter \gamma=1.8). We also derived a dependence of r_0 on the BH mass, M_BH, and found an indication of an increasing trend of r_0 at M_BH > 10^8 Msun. It is suggested that the growth of SMBH is mainly driven by interactions with the surrounding environment for M_BH > 10^8 Msun. On the other hand, at M_BH <~10^8 Msun, we did not find the BH mass dependence. This would imply that, for less massive BHs, the mass growth process can be different from that for massive BHs.

We evaluate the exposure during nadir observations with JEM-EUSO, the Extreme Universe Space Observatory, on-board the Japanese Experiment Module of the International Space Station. Designed as a mission to explore the extreme energy Universe from space, JEM-EUSO will monitor the Earth's nighttime atmosphere to record the ultraviolet light from tracks generated by extensive air showers initiated by ultra-high energy cosmic rays. In the present work, we discuss the particularities of space-based observation and we compute the annual exposure in nadir observation. The results are based on studies of the expected trigger aperture and observational duty cycle, as well as, on the investigations of the effects of clouds and different types of background light. We show that the annual exposure is about one order of magnitude higher than those of the presently operating ground-based observatories.

ALMA is estimated to generate TB scale data during only one observation; astronomers manage to identify which part of the data they are really interested in. Now we have been developing new GUI software for this purpose utilizing the VO interface: ALMA Web Quick Look System (ALMAWebQL) and ALMA Desktop Application (Vissage). The former is written in JavaScript and HTML5 generated from Java codes by Google Web Toolkit, and the latter is in pure Java. An essential point of our approach is how to reduce network traffic: we prepare, in advance, "compressed" FITS files of 2x2x1 (horizontal, vertical, and spectral directions, respectively) binning, 2x2x2 binning, 4x4x2 binning data, and so on. These files are hidden from users, and Web QL automatically choose proper one by each user operation. Through this work, we find that network traffic in our system is still a bottleneck towards TB scale data distribution. Hence we have to develop alternative data containers for much faster data processing. In this paper, I introduce our data analysis systems, and describe what we learned through the development.

Contributions of the JEM-EUSO Collaboration to the 32nd International Cosmic Ray Conference, Beijing, August, 2011.

We have searched for very high energy (VHE) gamma rays from four blazars using the CANGAROO-III imaging atmospheric Cherenkov telescope. We report the results of the observations of H 2356-309, PKS 2155-304, PKS 0537-441, and 3C 279, performed from 2005 to 2009, applying a new analysis to suppress the effects of the position dependence of Cherenkov images in the field of view. No significant VHE gamma ray emission was detected from any of the four blazars. The GeV gamma-ray spectra of these objects were obtained by analyzing Fermi/LAT archival data. Non-simultaneous wide range (radio to VHE gamma-ray bands) spectral energy distributions (SEDs) including CANGAROO-III upper limits, GeV gamma-ray spectra, and archival data are discussed using a one-zone synchrotron self-Compton (SSC) model in combination with a external Compton (EC) radiation. The HBLs (H 2356-309 and PKS 2155-304) can be explained by a simple SSC model, and PKS 0537-441 and 3C 279 are well modeled by a combination of SSC and EC model. We find a consistency with the blazar sequence in terms of strength of magnetic field and component size.

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.

We present the result of projected cross correlation analysis of AGNs and galaxies at redshifts from 0.3 to 3.0. The Japanese Virtual Observatory (JVO) was used to obtain the Subaru Suprime-Cam images and UKIDSS catalog data around AGNs. We investigated 1,809 AGNs, which is about ten times larger a sample than that used in previous studies on AGN-galaxy clustering at redshifts larger than 0.6. 90% of the AGN samples are optically-selected AGN from the SDSS and 2dF catalogs. The galaxy samples at low redshift includes many redder objects from UKIDSS survey, while at higher redshift they are mainly blue galaxies from Suprime-Cam. We found significant excess of galaxies around the AGNs at redshifts from 0.3 to 1.8. For the low redshift samples ($z<0.9$), we obtained correlation length of $r_{0} = $5--6 $h^{-1}$Mpc ($\gamma = 1.8$), which indicates that the AGNs at this redshift range reside in a similar environment around typical local galaxies. We also found that AGNs at higher redshift ranges reside in a denser environment than lower redshift AGNs; For $z=1.3 \sim 1.8$ AGNs, the cross correlation length was measured as 11$^{+6}_{-3}$ $h^{-1}$Mpc ($\gamma=1.8$). Considering that our galaxies sample is based on optical observations with Suprime-Cam at the redshift range, it is expected that blue star-forming galaxies comprise the majority of objects that are observed to be clustered around the AGNs. It is successfully demonstrated that the use of the archive through the Virtual Observatory system can provide a powerful tool for investigating the small scale environment of the intermediate redshift AGNs.

We present optical and near infrared observations of GRB 080325 classified as a "Dark GRB". Near-infrared observations with Subaru/MOIRCS provided a clear detection of afterglow in Ks band, although no optical counterpart was reported. The flux ratio of rest-wavelength optical to X-ray bands of the afterglow indicates that the dust extinction along the line of sight to the afterglow is Av = 2.7 - 10 mag. This large extinction is probably the major reason for optical faintness of GRB 080325. The J - Ks color of the host galaxy, (J - Ks = 1.3 in AB magnitude), is significantly redder than those for typical GRB hosts previously identified. In addition to J and Ks bands, optical images in B, Rc, i', and z' bands with Subaru/Suprime-Cam were obtained at about one year after the burst, and a photometric redshift of the host is estimated to be z_{photo} = 1.9. The host luminosity is comparable to L^{*} at z \sim 2 in contrast to the sub-L^{*} property of typical GRB hosts at lower redshifts. The best-fit stellar population synthesis model for the host shows that a large dust extinction (Av = 0.8 mag) attributes to the red nature of the host and that the host galaxy is massive (M_{*} = 7.0 \times 10^{10} Msun) which is one of the most massive GRB hosts previously identified. By assuming that the mass-metallicity relation for star-forming galaxies at z \sim 2 is applicable for the GRB host, this large stellar mass suggests the high metallicity environment around GRB 080325, consistent with inferred large extinction.

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 present multi-waveband optical imaging data obtained from observations of the Subaru/XMM-Newton Deep Survey (SXDS). The survey field, centered at R.A.=02:18:00, decl.=-05:00:00, has been the focus of a wide range of multi-wavelength observing programs spanning from X-ray to radio wavelengths. A large part of the optical imaging observations are carried out with Suprime-Cam on Subaru Telescope at Mauna Kea in the course of Subaru Telescope Observatory Projects. This paper describes our optical observations, data reduction and analysis procedures employed, and the characteristics of the data products. A total area of 1.22 sqdeg is covered in five contiguous sub-fields, each of which corresponds to a single Suprime-Cam field of view (34'x27'), in five broad-band filters B, V, Rc, i', z' to the depths of B=28.4, V=27.8, Rc=27.7, i'=27.7 and z'=26.6 (AB, 3-sigma, 2-arcsec aperture). The data are reduced and compiled into five multi-waveband photometric catalogs, separately for each Suprime-Cam pointing. The i'-band catalogs contain about 900,000 objects, making the SXDS catalogs one of the largest multi-waveband catalogs in corresponding depth and area coverage. The SXDS catalogs can be used for an extensive range of astronomical applications such as the number density of the Galactic halo stars to the large scale structures at the distant universe. The number counts of galaxies are derived and compared with those of existing deep extragalactic surveys. The optical data, the source catalogs, and configuration files used to create the catalogs are publicly available via the SXDS web page (http://www.naoj.org/Science/SubaruProject/SXDS/index.html)

We present our activity on making the Subaru Data Archive accessible through the Japanese Virtual Observatory (JVO) system. There are a lot of demand to use the archived Subaru data from various fields of astronomers. To be used by those who are not familiar with the way to reduce the Subaru data, the data reduction should be made before providing for them, or at least it should be easily done without precise knowledge about instrument's characteristic and environment where data are taken. For those purposes, data quality assessment system NAQATA is developed, which is presented in this meeting by Nakata et al. (2006), and the science-ready image data are provided for some of the SuprimeCam data at SMOKA data service which is presented by Enoki et al (2006). JVO portal will provide the way to access the reduce data, and also provides the way to reduce from raw data with very few efforts through the user-friendly web browser I/F. To provide such a CPU-intensive service, we have developed a GRID computing system. The architecture of this Subaru Data and Analysis system are discussed.

We present the design and implementation of the Japanese Virtual Observatory (JVO) system. JVO is a portal site to various kinds of astronomical resources distributed all over the world. We have developed five components for constructing the portal: (1) registry, (2) data service, (3) workflow system, (4) data analysis service (5) portal GUI. Registry services are used for publishing and searching data services in the VO, and they are constructed using an OAI-PMH metadata harvesting protocol and a SOAP web service protocol so that VO standard architecture is applied. Data services are developed based on the Astronomical Data Query Language (ADQL) which is an international VO standard and an extension of the standard SQL. The toolkit for building the ADQL-based service is released to the public on the JVO web site. The toolkit also provides the protocol translation from a Simple Image Access Protocol (SIAP) to ADQL protocol, so that both the VO standard service can be constructed using our toolkit. In order to federate the distributed databases and analysis services, we have designed a workflow language which is described in XML and developed execution system of the workflow. We have succeeded to connect to a hundred of data resources of the world as of April 2006. We have applied this system to the study of QSO environment by federating a QSO database, a Subaru Suprim-Cam database, and some analysis services such a SExtractor and HyperZ web services. These experiences are described is this paper.

Gamma-ray bursts (GRBs) and their afterglows have been proposed as an excellent probe to study the evolution of cosmic star formation, the reionization of the intergalactic medium, and the metal enrichment history of the universe, since the prompt gamma-ray emission of GRBs should be detectable out to distances z>10. Hitherto, the highest measured redshift for a GRB has been z=4.50. Here we report the optical spectrum of the afterglow of GRB 050904 obtained 3.4 days after the burst. The spectrum shows a clear continuum at the long wavelength end of the spectrum with a sharp cutoff at around 9000 A due to Ly alpha absorption at a redshift of 6.3 with a damping wing. Little flux is present in the waveband shortward of the Ly alpha break. A system of absorption lines of heavy elements at redshift z=6.295 +- 0.002 were also detected, yielding a precise measurement of the largest known redshift of a GRB. Analysis of the Si II fine structure lines suggest a dense metal-enriched environment around the GRB progenitor, providing unique information on the properties of the gas in a galaxy when the universe was younger than one billion years.

During the period when the Sun was intensely active on October-November 2003, two remarkable solar neutron events were observed by the ground-based neutron monitors. On October 28, 2003, in association with an X17.2 large flare, solar neutrons were detected with high statistical significance (6.4 sigma) by the neutron monitor at Tsumeb, Namibia. On November 4, 2003, in association with an X28 class flare, relativistic solar neutrons were observed by the neutron monitors at Haleakala in Hawaii and Mexico City, and by the solar neutron telescope at Mauna Kea in Hawaii simultaneously. Clear excesses were observed at the same time by these detectors, with the significance calculated as 7.5 sigma for Haleakala, and 5.2 sigma for Mexico City. The detector onboard the INTEGRAL satellite observed a high flux of hard X-rays and gamma-rays at the same time in these events. By using the time profiles of the gamma-ray lines, we can explain the time profile of the neutron monitor. It appears that neutrons were produced at the same time as the gamma-ray emission.

Currently two query languages are defined as standards for the Virtual Observatory (VO). Astronomical Data Query Language (ADQL) is used for catalog data query and Simple Image Access Protocol (SIAP) is for image data query. As a result, when we query each data service, we need to know in advance which language is supported and then construct a query language accordingly. The construct of SIAP is simple, but they have a limited capability. For example, there is no way to specify multiple regions in one query, and it is difficult to specify complex query conditions. In this paper, we propose a unified query language for any kind of astronomical database on the basis of SQL99. SQL is a query language optimized for a table data, so to apply the SQL to the image and spectrum data set, the data structure need to be mapped to a table like structure. We present specification of this query language and an example of the architecture for the database system.

The Subaru Deep Field (SDF) project is a program of Subaru Observatory to carry out a deep galaxy survey over a blank field as large as 34'x27'. The program consists of very deep multi-band optical imaging, near infrared imaging for smaller portions of the field and follow-up optical spectroscopy. Major scientific goals of the project are to construct large samples of Lyman-break galaxies at z~4-5 and Lyman alpha emitters at z~5.7 and 6.6, and to make detailed studies these very high-redshift galaxy populations. In this paper, we describe the optical imaging observations and data reduction, presenting mosaicked images and object catalogs in seven bandpasses.The optical imaging was made through five broad-band filters, B, V, R, i', z', and two narrow-band filters, NB816 (lambda_c=8150A) and NB921 (lambda_c=9196A) with almost 10 hours long integrations for each band. The limiting magnitudes measured at 3-sigma on a 2" aperture are B=28.45, V=27.74, R=27.80, i'=27.43, z'=26.62, NB816=26.63, and NB921=26.54 in the AB system. The object catalog constructed for each of the seven bands contains more than 10^5 objects. The galaxy number counts corrected for detection incompleteness and star count contribution are found to be consistent with previous results in the literature. The mosaicked images and catalogs of all the bands have been made open to the public on Oct. 1, 2004 on the SDF project website at http://soaps.naoj.org/sdf/.

We report on the discovery of a star-forming galaxy at z=6.33 in the Subaru Deep Field. This object is selected as a candidate of an i'-dropout, high-redshift galaxy around z=6 because of its red i'-z' color in our deep optical imaging survey in the Subaru Deep Field. Our follow up optical spectroscopy reveals that this object is a strong Ly-alpha emitter with only very faint ultraviolet continuum. The rest-frame equivalent width of the detected Ly-alpha emission is as much as 130 A. Thus the light detected in our z' image is largely attributed to the Ly-alpha emission, i.e., ~40% of the z'-band flux is the strong Ly-alpha emission, giving a very red i'-z' color. This is consistent with the photometric property of this object because the narrow-band data obtained with the NB921 filter shows a significant depression, z'-NB921 = -0.54 mag. By using the photometric data, we show that some other objects among the 48 i'-dropout high-redshift galaxy candidates found in the Subaru Deep Field also show a significant NB921 depression. We briefly discuss the nature of these NB921-depressed objects.

Spectroscopic and spectropolarimetric observations of SN 2003dh/GRB 030329 obtained in 2003 May using the Subaru 8.2 m telescope are presented. The properties of the SN are investigated through a comparison with spectra of the Type Ic hypernovae SNe 1997ef and 1998bw. (Hypernovae being a tentatively defined class of SNe with very broad absorption features: these features suggest a large velocity of the ejected material and possibly a large explosion kinetic energy.) Comparison with spectra of other hypernovae shows that the spectrum of SN 2003dh obtained on 2003 May 8 and 9, i.e., 34-35 rest-frame days after the GRB (for z=0.1685), are similar to those of SN 1997ef obtained ~34-42 days after the fiducial time of explosion of that SN. The match with SN 1998bw spectra is not as good (at rest 7300-8000 A, but again spectra obtained ~33-43 days after GRB 980425 are preferred. This indicates that the SN may have intermediate properties between SNe 1997ef and 1998bw. Based on the analogy with the other hypernovae, the time of explosion of SN 2003dh is then constrained to be between -8 and +2 days of the GRB. The Si and O P-Cygni lines of SN 2003dh seem comparable to those of SN 1997ef, which suggests that the ejected mass in SN 2003dh may match that in SN 1997ef. Polarization was marginally detected at optical wavelengths. This is consistent with measurements of the late afterglow, implying that it mostly originated in the interstellar medium of the host galaxy.

It has been suggested that cosmic strings produced at a phase transition in the early universe can be the origin of the extremely high energy cosmic rays (EHCR) observed by AGASA above 10^20 eV. Superheavy cosmic strings with linear mass density of 10^22 g/cm can be indirectly observed through the gravitational lensing effect the distant galaxies. The lensing effect by a long straight object can be characterized by a line of double galaxies or quasars with angular separation of about 5 arcsec. We have searched for aligned double objects from the archived data taken by the Subaru Prime Focus Camera (Suprime-Cam). The Suprime-Cam has a great advantage in observing the wide field of view (30x30 arcmin^2) with high sensitivity (R<26 400s exposure), so it is suitable for this research. In this paper, we describe the result of simulation study for developing the method of searching the objects lensed by cosmic strings, and present the observational result obtained by this method.

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 have performed a long-slit K band spectroscopic observation of the luminous infrared galaxy NGC 6240. The peak position of the H_2 v=1-0 S(1) emission in the slit is located ~0.3" - 0.4" north of the southern nucleus. It is almost the midpoint between the southern nucleus and the peak position of the ^12CO J=1-0 emission. Based on the line-ratio analyses, we suggest the excitation mechanism of H_2 is pure thermal at most positions. In the southern region we find the following three velocity components in the H_2 emission: the blueshifted shell component (~-250 km s^-1 with respect to V_sys) which is recognized as a distinct C-shape distortion in the velocity field around the southern nucleus, the high-velocity blueshifted ``wing'' component (~-1000 km s^-1 with respect to V_sys), and the component indicating possible line splitting of ~500 km s^-1. The latter two components are extended to the south from the southern nucleus. We show that the kinematic properties of these three components can be reproduced by expanding motion of a shell-like structure around the southern nucleus. The offset peak position of the H_2 emission can be understood if we assume that the shell expanding to the north interacts with the extragalactic molecular gas. At the interface between the shell and the molecular gas concentration the cloud-crushing mechanism proposed by Cowie et al. (1981) may work efficiently, and the intense H_2 emission is thus expected there. All these findings lead us to propose a model that the most H_2 emission is attributed to the shock excitation driven by the superwind activity of the southern nucleus.

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.

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

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.

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.

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.

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.

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.

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.

The Telescope Array prototype detectors were installed at Akeno Observatory and at the Utah Fly's Eye site. Using these detectors, we have observed the Crab Nebula and AGN's since the end of 1995. The successful detections of TeV gamma rays from Crab Nebula and Mkn501 are reported.