We report the confirmation and mass determination of three hot Jupiters discovered by the Transiting Exoplanet Survey Satellite (TESS) mission: HIP 65Ab (TOI-129, TIC-201248411) is an ultra-short-period Jupiter orbiting a bright (V=11.1 mag) K4-dwarf every 0.98 days. It is a massive 3.213 +/- 0.078 Mjup planet in a grazing transit configuration with impact parameter b = 1.17 +0.10/-0.08. As a result the radius is poorly constrained, 2.03 +0.61/-0.49 Rjup. We perform a full phase-curve analysis of the TESS data and detect both illumination- and ellipsoidal variations as well as Doppler boosting. HIP 65A is part of a binary stellar system, with HIP 65B separated by 269 AU (3.95 arcsec on sky). TOI-157b (TIC 140691463) is a typical hot Jupiter with a mass 1.18 +/- 0.13 Mjup and radius 1.29 +/- 0.02 Rjup. It has a period of 2.08 days, which corresponds to a separation of just 0.03 AU. This makes TOI-157 an interesting system, as the host star is an evolved G9 sub-giant star (V=12.7). TOI-169b (TIC 183120439) is a bloated Jupiter orbiting a V=12.4 G-type star. It has a mass of 0.79 +/- 0.06 Mjup and radius 1.09 +0.08/-0.05 Rjup. Despite having the longest orbital period (P=2.26 days) of the three planets, TOI-169b receives the most irradiation and is situated on the edge of the Neptune desert. All three host stars are metal rich with Fe/H ranging from 0.18 - 0.24.

WIMPs are promising dark matter candidates. A WIMP occasionally collides with a mirror equipped with interferometric gravitational wave detectors such as LIGO, Virgo, KAGRA and the Einstein Telescope (ET). When WIMPs collide with a mirror of an interferometer, we expect that characteristic motions of the pendulum and mirror are excited, and those signals could be extracted by highly sophisticated sensors developed for gravitational wave detection. We analyze the motions of the pendulum and mirror, and estimate the detectability of these motions. For the "Thin-ET" detector, the signal-to-noise ratio may be $ 1.7 \left( \frac{ m_{\rm{DM}} }{ 100 \rm{GeV} } \right) $, where $ m_{\rm{DM}} $ is the mass of a WIMP. We may set a more strict upper limit on the cross section between a WIMP and a nucleon than the limits obtained by other experiments so far when $ m_{\rm{DM}} $ is approximately lower than 0.2 GeV. We find an order-of-magnitude improvement in the upper limit around $ m_{\rm{DM}} = 0.2 {\rm{GeV}} $.

We present the discovery and characterisation of two transiting planets observed by the Transiting Exoplanet Survey Satellite (TESS) orbiting the nearby (d ~ 22 pc), bright (J ~ 9 mag) M3.5 dwarf LP 729-54 (TOI-732). We confirm both planets and their association with LP 729-54 via ground-based photometry and determine their masses using precise radial velocities measured with the CARMENES spectrograph. Precise stellar parameters determined from CARMENES high resolution spectra confirm that LP 729-54 is a mid-M dwarf with an effective temperature of T_eff = 3360 +\- 51 K, a surface gravity of log(g) = 4.81 +/- 0.04 (cgs), and an iron abundance of [Fe/H] = 0.09 +/- 0.16 dex, with an inferred mass of M_star = 0.379 +/- 0.016 M_sun and a radius of R_star = 0.382 +/- 0.012 R_sun. The ultra-short period planet LP 729-54 b (P_b = 0.77 d) with a radius of 1.35^{+0.06}_{-0.06} R_earth, a mass of 2.34^{+0.24}_{-0.23} M_earth, and a bulk density of 5.24^{+0.94}_{-0.81} g cm^{-3} joins the population of Earth-size planets with rocky, terrestrial composition. The outer planet, LP 729-54 c, with an orbital period of 12.25 d, radius of 2.42^{+0.10}_{-0.10} R_earth, mass of 6.29^{+0.63}_{-0.61} M_earth, and mean density of 2.45^{+0.44}_{-0.37} g cm^{-3} belongs to the population of dense sub-Neptunes. With the two planets located on opposite sides of the radius gap, this planetary system is an excellent target for testing planetary formation, evolution and atmospheric models. In particular, LP 729-54 c is an ideal object for atmospheric studies with the James Webb Space Telescope.

We investigate the distribution of metals in the cosmological volume at $z\sim3$, in particular, provided by massive population III (Pop III) stars using a cosmological $N$-body simulation in which a model of Pop III star formation is implemented. Owing to the simulation, we can choose minihaloes where Pop III star formation occurs at $z>10$ and obtain the spatial distribution of the metals at lower-redshifts. To evaluate the amount of heavy elements provided by Pop III stars, we consider metal yield of pair-instability or core-collapse supernovae (SNe) explosions of massive stars. By comparing our results to the Illustris-1 simulation, we find that heavy elements provided by Pop III stars often dominate those from galaxies in low density regions. The median value of the volume averaged metallicity is $Z\sim 10^{-4.5 - -2} Z_{\odot}$ at the regions. Spectroscopic observations with the next generation telescopes are expected to detect the metals imprinted on quasar spectra.

We report a search for astronomical neutrinos in the energy region from several GeV to TeV in the direction of the blazar TXS0506+056 using the Super-Kamiokande detector following the detection of a 100 TeV neutrino from the same location by the IceCube collaboration. Using Super-Kamiokande neutrino data across several data samples observed from April 1996 to February 2018 we have searched for both a total excess above known backgrounds across the entire period as well as localized excesses on smaller time scales in that interval. No significant excess nor significant variation in the observed event rate are found in the blazar direction. Upper limits are placed on the electron and muon neutrino fluxes at 90\% confidence level as $6.03 \times 10^{-7}$ and $4.52 \times 10^{-7}$ to $9.26 \times 10^{-10}$ [${\rm erg}/{\rm cm}^2/{\rm s}$], respectively.

We report the discovery of TOI 263.01 (TIC 120916706), a transiting substellar object (R = 0.87 RJup) orbiting a faint M3.5~V dwarf (V=18.97) on a 0.56~d orbit. We set out to determine the nature of the TESS planet candidate TOI 263.01 using ground-based multicolour transit photometry. The host star is faint, which makes RV confirmation challenging, but the large transit depth makes the candidate suitable for validation through multicolour photometry. Our analysis combines three transits observed simultaneously in r', i', and z_s bands using the MuSCAT2 multicolour imager, three LCOGT-observed transit light curves in g, r', and i' bands, a TESS light curve from Sector 3, and a low-resolution spectrum for stellar characterisation observed with the ALFOSC spectrograph. We model the light curves with PyTransit using a transit model that includes a physics-based light contamination component that allows us to estimate the contamination from unresolved sources from the multicolour photometry. This allows us to derive the true planet-star radius ratio marginalised over the contamination allowed by the photometry, and, combined with the stellar radius, gives us a reliable estimate of the object's absolute radius. The ground-based photometry excludes contamination from unresolved sources with a significant colour difference to TOI 263. Further, contamination from sources of same stellar type as the host is constrained to levels where the true radius ratio posterior has a median of 0.217. The median radius ratio corresponds to an absolute planet radius of 0.87 RJup, which confirms the substellar nature of the planet candidate. The object is either a giant planet or a brown dwarf (BD) located deep inside the so-called "brown dwarf desert". Both possibilities offer a challenge to current planet/BD formation models and makes 263.01 an object deserving of in-depth follow-up studies.

We report the analysis of additional multiband photometry and spectroscopy and new adaptive optics (AO) imaging of the nearby planetary microlensing event TCP J05074264+2447555 (Kojima-1), which was discovered toward the Galactic anticenter in 2017 (Nucita et al.). We confirm the planetary nature of the light-curve anomaly around the peak while finding no additional planetary feature in this event. We also confirm the presence of apparent blending flux and the absence of significant parallax signal reported in the literature. The AO image reveals no contaminating sources, making it most likely that the blending flux comes from the lens star. The measured multiband lens flux, combined with a constraint from the microlensing model, allows us to narrow down the previously unresolved mass and distance of the lens system. We find that the primary lens is a dwarf on the K/M boundary (0.581 \pm 0.033 M_sun) located at 505 \pm 47 pc and the companion (Kojima-1Lb) is a Neptune-mass planet (20.0 \pm 2.0 M_earth) with a semi-major axis of 1.08 ^{+0.62}_{-0.18} au. This orbit is a few times smaller than those of typical microlensing planets and is comparable to the snow-line location at young ages. We calculate that the a priori detection probability of Kojima-1Lb is only \sim 35%, which may imply that Neptunes are common around the snow line, as recently suggested by the transit and radial velocity techniques. The host star is the brightest among the microlensing planetary systems (Ks = 13.7), offering a great opportunity to spectroscopically characterize this system, even with current facilities.

Supernova detection is a major objective of the Super-Kamiokande (SK) experiment. In the next stage of SK (SK-Gd), gadolinium (Gd) sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons. A core-collapse supernova will be preceded by an increasing flux of neutrinos and anti-neutrinos, from thermal and weak nuclear processes in the star, over a timescale of hours; some of which may be detected at SK-Gd. This could provide an early warning of an imminent core-collapse supernova, hours earlier than the detection of the neutrinos from core collapse. Electron anti-neutrino detection will rely on inverse beta decay events below the usual analysis energy threshold of SK, so Gd loading is vital to reduce backgrounds while maximising detection efficiency. Assuming normal neutrino mass ordering, more than 200 events could be detected in the final 12 hours before core collapse for a 15-25 solar mass star at around 200 pc, which is representative of the nearest red supergiant to Earth, $\mathrm{\alpha}$Ori (Betelgeuse). At a statistical false alarm rate of 1 per century, detection could be up to 10 hours before core collapse, and a pre-supernova star could be detected by SK-Gd up to 600 pc away. A pre-supernova alert could be provided to the astrophysics community following gadolinium loading.

Supernova neutrinos are crucially important to probe the final phases of massive star evolution. As is well known from observations of SN1987A, neutrinos provide information on the physical conditions responsible for neutron star formation and on the supernova explosion mechanism. However, there is still no complete understanding of the long-term evolution of neutrino emission in supernova explosions, although there are a number of modern simulations of neutrino radiation hydrodynamics, which study neutrino emission at times less than one second after the bounce. In the present work we systematically calculate the number of neutrinos that can be observed in Super-Kamiokande over periods longer than ten seconds using the database of Nakazato et al. (2013) anticipating that neutrinos from a Galactic supernova can be detected for several tens of seconds. We find that for a supernova at a distance of 10 kpc, neutrinos remain observable for longer than 30 s for a low-mass neutron star ($1.20M_\odot$ gravitational mass) and even longer than 100 s for a high-mass neutron star ($2.05M_\odot$). These scenarios are much longer than the observations of SN1987A and longer than the duration of existing numerical simulations. We propose a new analysis method based on the cumulative neutrino event distribution as a function of reverse time from the last observed event, as a useful probe of the neutron star mass. Our result demonstrates the importance of complete modeling of neutrino light curves in order to extract physical quantities essential for understanding supernova explosion mechanisms, such as the mass and radius of the resulting neutron star.

Context. The TESS and PLATO missions are expected to find vast numbers of new transiting planet candidates. However, only a fraction of these candidates will be legitimate planets, and the candidate validation will require a significant amount of follow-up resources. Radial velocity follow-up can be carried out only for the most promising candidates around bright, slowly rotating, stars. Thus, before devoting RV resources to candidates, they need to be vetted using cheaper methods, and, in the cases for which an RV confirmation is not feasible, the candidate's true nature needs to be determined based on these alternative methods alone. Aims. We study the applicability of multicolour transit photometry in the validation of transiting planet candidates when the candidate signal arises from a real astrophysical source. We seek to answer how securely can we estimate the true uncontaminated star-planet radius ratio when the light curve may contain contamination from unresolved light sources inside the photometry aperture when combining multicolour transit observations with a physics-based contamination model. Methods. The study is based on simulations and ground-based transit observations. The analyses are carried out with a contamination model integrated into the PyTransit v2 transit modelling package, and the observations are carried out with the MuSCAT2 multicolour imager installed in the 1.5 m TCS in the Teide Observatory. Results. We show that multicolour transit photometry can be used to estimate the amount of flux contamination and the true radius ratio. Combining the true radius ratio with an estimate for the stellar radius yields the true absolute radius of the transiting object, which is a valuable quantity in statistical candidate validation, and enough in itself to validate a candidate whose radius falls below the theoretical lower limit for a brown dwarf.

A sensitive survey of the MeV gamma-ray sky is needed to understand important astrophysical problems such as gamma-ray bursts in the early universe, progenitors of Type Ia supernovae, and the nature of dark matter. However, the study has not progressed remarkably since the limited survey by COMPTEL onboard CGRO in the 1990s. Tanimori et al. have developed a Compton camera that tracks the trajectory of each recoil electron in addition to the information obtained by the conventional Compton cameras, leading to superior imaging. This Electron Tracking Compton Camera (ETCC) facilitates accurate reconstruction of the incoming direction of each MeV photon from a wide sky at ~degree angular resolution and with minimized particle background using trajectory information. The latest ETCC model, SMILE-2+, made successful astronomical observations during a day balloon flight in 2018 April and detected diffuse continuum and 511 keV annihilation line emission from the Galactic Center region at a high significance in ~2.5 hours. We believe that MeV observations from space with upgraded ETCCs will dramatically improve our knowledge of the MeV universe. We advocate for a space-based all-sky survey mission with multiple ETCCs onboard and detail its expected benefits.

In this paper, we present the analysis and results of a direct measurement of the cosmic-ray proton spectrum with the CALET instrument onboard the International Space Station, including the detailed assessment of systematic uncertainties. The observation period used in this analysis is from October 13, 2015 to August 31, 2018 (1054 days). We have achieved the very wide energy range necessary to carry out measurements of the spectrum from 50 GeV to 10 TeV covering, for the first time in space, with a single instrument the whole energy interval previously investigated in most cases in separate subranges by magnetic spectrometers (BESS-TeV, PAMELA, and AMS-02) and calorimetric instruments (ATIC, CREAM, and NUCLEON). The observed spectrum is consistent with AMS-02 but extends to nearly an order of magnitude higher energy, showing a very smooth transition of the power-law spectral index from -2.81 +- 0.03 (50--500 GeV) neglecting solar modulation effects (or -2.87 +- 0.06 including solar modulation effects in the lower energy region) to -2.56 +- 0.04 (1--10 TeV), thereby confirming the existence of spectral hardening and providing evidence of a deviation from a single power law by more than 3 sigma.

The CALorimetric Electron Telescope (CALET) is a high-energy astroparticle physics space experiment installed on the International Space Station (ISS), developed and operated by Japan in collaboration with Italy and the United States. The CALET mission goals include the investigation of possible nearby sources of high-energy electrons, of the details of galactic particle acceleration and propagation, and of potential signatures of dark matter. CALET measures the cosmic-ray electron + positron flux up to 20 TeV, gamma-rays up to 10 TeV, and nuclei with Z=1 to 40 up to 1,000 TeV for the more abundant elements during a long-term observation aboard the ISS. Starting science operation in mid-October 2015, CALET performed continuous observation without major interruption with close to 20 million triggered events over 10 GeV per month. Based on the data taken during the first two-years, we present an overview of CALET observations: uses w/o major interruption 1) Electron + positron energy spectrum, 2) Nuclei analysis, 3) Gamma-ray observation including a characterization of on-orbit performance. Results of the electromagnetic counterpart search for LIGO/Virgo gravitational wave events are discussed as well.

We report the Transiting Exoplanet Survey Satellite ($TESS$) detection of a multi-planet system orbiting the $V=10.9$ K0 dwarf TOI 125. We find evidence for up to five planets, with varying confidence. Three high signal-to-noise transit signals correspond to sub-Neptune-sized planets ($2.76$, $2.79$, and $2.94\ R_{\oplus}$), and we statistically validate the planetary nature of the two inner planets ($P_b = 4.65$ days, $P_c = 9.15$ days). With only two transits observed, we report the outer object ($P_{.03} = 19.98$ days) as a high signal-to-noise ratio planet candidate. We also detect a candidate transiting super-Earth ($1.4\ R_{\oplus}$) with an orbital period of only $12.7$ hours and a candidate Neptune-sized planet ($4.2\ R_{\oplus}$) with a period of $13.28$ days, both at low signal-to-noise. This system is amenable to mass determination via radial velocities and transit timing variations, and provides an opportunity to study planets of similar size while controlling for age and environment. The ratio of orbital periods between TOI 125 b and c ($P_c/P_b = 1.97$) is slightly smaller than an exact 2:1 commensurability and is atypical of multiple planet systems from $Kepler$, which show a preference for period ratios just $wide$ of first-order period ratios. A dynamical analysis refines the allowed parameter space through stability arguments and suggests that, despite the nearly commensurate periods, the system is unlikely to be in resonance.

We report the development of a 4-color simultaneous camera for the 1.52~m Telescopio Carlos S\'anchez (TCS) in the Teide Observatory, Canaries, Spain. The new instrument, named MuSCAT2, has a capability of 4-color simultaneous imaging in $g$ (400--550 nm), $r$ (550--700 nm), $i$ (700--820 nm), and $z_s$ (820--920 nm) bands. MuSCAT2 equips four 1024$\times$1024 pixel CCDs, having a field of view of 7.4$\times$7.4 arcmin$^2$ with a pixel scale of 0.44 arcsec per pixel. The principal purpose of MuSCAT2 is to perform high-precision multi-color exoplanet transit photometry. We have demonstrated photometric precisions of 0.057%, 0.050%, 0.060%, and 0.076% as root-mean-square residuals of 60~s binning in $g$, $r$, $i$ and $z_s$ bands, respectively, for a G0 V star WASP-12 ($V=11.57\pm0.16$). MuSCAT2 has started science operations since January 2018, with over 250 telescope nights per year. MuSCAT2 is expected to become a reference tool for exoplanet transit observations, and will substantially contribute to the follow-up of the TESS and PLATO space missions.

We present results on searches for gamma-ray counterparts of the LIGO/Virgo gravitational-wave events using CALorimetric Electron Telescope ({\sl CALET}) observations. The main instrument of {\sl CALET}, CALorimeter (CAL), observes gamma-rays from $\sim1$ GeV up to 10 TeV with a field of view of nearly 2 sr. In addition, the {\sl CALET} gamma-ray burst monitor (CGBM) views $\sim$3 sr and $\sim2\pi$ sr of the sky in the 7 keV -- 1 MeV and the 40 keV -- 20 MeV bands, respectively, by using two different crystal scintillators. The {\sl CALET} observations on the International Space Station started in October 2015, and here we report analyses of events associated with the following gravitational wave events: GW151226, GW170104, GW170608, GW170814 and GW170817. Although only upper limits on gamma-ray emission are obtained, they correspond to a luminosity of $10^{49}\sim10^{53}$ erg s$^{-1}$ in the GeV energy band depending on the distance and the assumed time duration of each event, which is approximately the order of luminosity of typical short gamma-ray bursts. This implies there will be a favorable opportunity to detect high-energy gamma-ray emission in further observations if additional gravitational wave events with favorable geometry will occur within our field-of-view. We also show the sensitivity of {\sl CALET} for gamma-ray transient events which is the order of $10^{-7}$~erg\,cm$^{-2}$\,s$^{-1}$ for an observation of 100~s duration.

Extended results on the cosmic-ray electron + positron spectrum from 11 GeV to 4.8 TeV are presented based on observations with the Calorimetric Electron Telescope (CALET) on the International Space Station utilizing the data up to November 2017. The analysis uses the full detector acceptance at high energies, approximately doubling the statistics compared to the previous result. CALET is an all-calorimetric instrument with a total thickness of 30 $X_0$ at normal incidence and fine imaging capability, designed to achieve large proton rejection and excellent energy resolution well into the TeV energy region. The observed energy spectrum in the region below 1 TeV shows good agreement with Alpha Magnetic Spectrometer (AMS-02) data. In the energy region below $\sim$300 GeV, CALET's spectral index is found to be consistent with the AMS-02, Fermi Large Area Telescope (Fermi-LAT) and Dark Matter Particle Explorer (DAMPE), while from 300 to 600 GeV the spectrum is significantly softer than the spectra from the latter two experiments. The absolute flux of CALET is consistent with other experiments at around a few tens of GeV. However, it is lower than those of DAMPE and Fermi-LAT with the difference increasing up to several hundred GeV. The observed energy spectrum above $\sim$1 TeV suggests a flux suppression consistent within the errors with the results of DAMPE, while CALET does not observe any significant evidence for a narrow spectral feature in the energy region around 1.4 TeV. Our measured all-electron flux, including statistical errors and a detailed breakdown of the systematic errors, is tabulated in the Supplemental Material in order to allow more refined spectral analyses based on our data.

We report the results of a neutrino search in Super-Kamiokande for coincident signals with the first detected gravitational wave produced by a binary neutron star merger, GW170817, which was followed by a short gamma-ray burst, GRB170817A, and a kilonova/macronova. We searched for coincident neutrino events in the range from 3.5 MeV to $\sim$100 PeV, in a time window $\pm$500 seconds around the gravitational wave detection time, as well as during a 14-day period after the detection. No significant neutrino signal was observed for either time window. We calculated 90% confidence level upper limits on the neutrino fluence for GW170817. From the upward-going-muon events in the energy region above 1.6 GeV, the neutrino fluence limit is $16.0^{+0.7}_{-0.6}$ ($21.3^{+1.1}_{-0.8}$) cm$^{-2}$ for muon neutrinos (muon antineutrinos), with an error range of $\pm5^{\circ}$ around the zenith angle of NGC4993, and the energy spectrum is under the assumption of an index of $-2$. The fluence limit for neutrino energies less than 100 MeV, for which the emission mechanism would be different than for higher-energy neutrinos, is also calculated. It is $6.6 \times 10^7$ cm$^{-2}$ for anti-electron neutrinos under the assumption of a Fermi-Dirac spectrum with average energy of 20 MeV.

The CALorimetric Electron Telescope (CALET), launched for installation on the International Space Station (ISS) in August, 2015, has been accumulating scientific data since October, 2015. CALET is intended to perform long-duration observations of high-energy cosmic rays onboard the ISS. CALET directly measures the cosmic-ray electron spectrum in the energy range of 1 GeV to 20 TeV with a 2% energy resolution above 30 GeV. In addition, the instrument can measure the spectrum of gamma rays well into the TeV range, and the spectra of protons and nuclei up to a PeV. In order to operate the CALET onboard ISS, JAXA Ground Support Equipment (JAXA-GSE) and the Waseda CALET Operations Center (WCOC) have been established. Scientific operations using CALET are planned at WCOC, taking into account orbital variations of geomagnetic rigidity cutoff. Scheduled command sequences are used to control the CALET observation modes on orbit. Calibration data acquisition by, for example, recording pedestal and penetrating particle events, a low-energy electron trigger mode operating at high geomagnetic latitude, a low-energy gamma-ray trigger mode operating at low geomagnetic latitude, and an ultra heavy trigger mode, are scheduled around the ISS orbit while maintaining maximum exposure to high-energy electrons and other high-energy shower events by always having the high-energy trigger mode active. The WCOC also prepares and distributes CALET flight data to collaborators in Italy and the United States. As of August 31, 2017, the total observation time is 689 days with a live time fraction of the total time of approximately 84%. Nearly 450 million events are collected with a high-energy (E>10 GeV) trigger. By combining all operation modes with the excellent-quality on-orbit data collected thus far, it is expected that a five-year observation period will provide a wealth of new and interesting results.

We present wide and deep photometry of the northwest part of the halo of the Andromeda galaxy (M31) using Hyper Suprime-Cam on the Subaru Telescope. The survey covers 9.2 deg$^{2}$ field in the $g$, $i$, and $NB515$ bands and shows a clear red giant branch (RGB) of M31's halo stars and a pronounced red clump (RC) feature. The spatial distribution of RC stars shows a prominent stream feature, the North Western (NW) Stream, and a diffuse substructure in the south part of our survey field. We estimate the distances based on the RC method and obtain $(m-M)$ = 24.63$\pm 0.191$(random)$\pm0.057$(systematic) and 24.29$\pm 0.211$(random)$\pm0.057$(systematic) mag for the NW stream and diffuse substructure, respectively, implying that the NW Stream is located behind M31, whereas the diffuse substructure is located in front. We also estimate line-of-sight distances along the NW Stream and find that the south part of the stream is $\sim$20 kpc closer to us relative to the north part. The distance to the NW Stream inferred from the isochrone fitting to the color-magnitude diagram favors the RC-based distance, but the TRGB-based distance estimated for $NB515$-selected RGB stars does not agree with it. The surface number density distribution of RC stars across the NW Stream is found to be approximately Gaussian with a FWHM of $\sim$25 arcmin (5.7 kpc), with a slight skew to the south-west side. That along the NW Stream shows a complicated structure including variations in number density and a significant gap in the stream.

We present the detection of a blackbody component in GRB 160107A emission by using the combined spectral data of the CALET Gamma-ray Burst Monitor (CGBM) and the MAXI Gas Slit Camera (GSC). The MAXI/GSC detected the emission $\sim$45 s prior to the main burst episode observed by the CGBM. The MAXI/GSC and the CGBM spectrum of this prior emission period is well fit by a blackbody with the temperature of $1.0^{+0.3}_{-0.2}$ keV plus a power-law with the photon index of $-1.6 \pm 0.3$. We discuss the radius to the photospheric emission and the main burst emission based on the observational properties. We stress the importance of the coordinated observations via various instruments collecting the high quality data over a broad energy coverage in order to understand the GRB prompt emission mechanism.

First results of a cosmic-ray electron + positron spectrum, from 10 GeV to 3 TeV, is presented based upon observations with the CALET instrument on the ISS starting in October, 2015. Nearly a half million electron + positron events are included in the analysis. CALET is an all-calorimetric instrument with total vertical thickness of 30 $X_0$ and a fine imaging capability designed to achieve a large proton rejection and excellent energy resolution well into the TeV energy region. The observed energy spectrum over 30 GeV can be fit with a single power law with a spectral index of -3.152 $\pm$ 0.016 (stat.+ syst.). Possible structure observed above 100 GeV requires further investigation with increased statistics and refined data analysis.

In August 2015, the CALorimetric Electron Telescope (CALET), designed for long exposure observations of high energy cosmic rays, docked with the International Space Station (ISS) and shortly thereafter began tocollect data. CALET will measure the cosmic ray electron spectrum over the energy range of 1 GeV to 20 TeV with a very high resolution of 2% above 100 GeV, based on a dedicated instrument incorporating an exceptionally thick 30 radiation-length calorimeter with both total absorption and imaging (TASC and IMC) units. Each TASC readout channel must be carefully calibrated over the extremely wide dynamic range of CALET that spans six orders of magnitude in order to obtain a degree of calibration accuracy matching the resolution of energy measurements. These calibrations consist of calculating the conversion factors between ADC units and energy deposits, ensuring linearity over each gain range, and providing a seamless transition between neighboring gain ranges. This paper describes these calibration methods in detail, along with the resulting data and associated accuracies. The results presented in this paper show that a sufficient accuracy was achieved for the calibrations of each channel in order to obtain a suitable resolution over the entire dynamic range of the electron spectrum measurement.

The theory of universal extra dimensions involves Kaluza-Klein (KK) particles. The lightest KK particle (LKP) is one of the good candidates for cold dark matter. Annihilation of LKP dark matter in the Galactic halo produces high-energy gamma-rays. The gamma-ray spectrum shows a characteristic peak structure around the LKP mass. This paper investigates the observability of this peak structure by near-future detectors taking account of their energy resolution, and calculates the expected energy spectrum of the gamma-ray signal. Then, by using the HESS data, we set some constraints on the boost factor, which is a product of the annihilation cross section relative to the thermal one and an uncertain factor dependent on the substructure of the LKP distribution in the Galactic halo, for each LKP mass. The resulting upper limit on the boost factor is in the range from 1 to 30. The constraints can be regarded as comparable with the results of previous work for gamma-ray and electron-positron observation. However, the observational data for the TeV or higher energy region are still limited, and the possible LKP signal is not conclusive. Thus, we expect near-future missions with better sensitivity will clarify whether the LKP dark matter should exist or not.

We present the SILVERRUSH program strategy and clustering properties investigated with $\sim 2,000$ Ly$\alpha$ emitters at $z=5.7$ and $6.6$ found in the early data of the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey exploiting the carefully designed narrowband filters. We derive angular correlation functions with the unprecedentedly large samples of LAEs at $z=6-7$ over the large total area of $14-21$ deg$^2$ corresponding to $0.3-0.5$ comoving Gpc$^2$. We obtain the average large-scale bias values of $b_{\rm avg}=4.1\pm 0.2$ ($4.5\pm 0.6$) at $z=5.7$ ($z=6.6$) for $\gtrsim L^*$ LAEs, indicating the weak evolution of LAE clustering from $z=5.7$ to $6.6$. We compare the LAE clustering results with two independent theoretical models that suggest an increase of an LAE clustering signal by the patchy ionized bubbles at the epoch of reionization (EoR), and estimate the neutral hydrogen fraction to be $x_{\rm HI}=0.15^{+0.15}_{-0.15}$ at $z=6.6$. Based on the halo occupation distribution models, we find that the $\gtrsim L^*$ LAEs are hosted by the dark-matter halos with the average mass of $\log (\left < M_{\rm h} \right >/M_\odot) =11.1^{+0.2}_{-0.4}$ ($10.8^{+0.3}_{-0.5}$) at $z=5.7$ ($6.6$) with a Ly$\alpha$ duty cycle of 1 % or less, where the results of $z=6.6$ LAEs may be slightly biased, due to the increase of the clustering signal at the EoR. Our clustering analysis reveals the low-mass nature of $\gtrsim L^*$ LAEs at $z=6-7$, and that these LAEs probably evolve into massive super-$L^*$ galaxies in the present-day universe.

Polytropic transonic solutions of spherically symmetric and steady galactic winds in the gravitational potential of a dark matter halo (DMH) with a supermassive black hole (SMBH) are studied. The solutions are classified in terms of their topological features, and the gravitational potential of the SMBH adds a new branch to the transonic solutions generated by the gravity of the DMH. The topological types of the transonic solutions depend on the mass distribution, the amount of supplied energy, the polytropic index $\gamma$, and the slope $\alpha$ of the DMH mass distribution. When $\alpha$ becomes larger than a critical value $\alpha_\mathrm{c}$, the transonic solution types change dramatically. Further, our model predicts that it is possible for a slowly accelerating outflow to exist, even in quiescent galaxies with small $\gamma$. This slowly accelerating outflow differs from those considered in many of the previous studies focusing on supersonic outflows in active star-forming galaxies. In addition, our model indicates that outflows in active star-forming galaxies have only one transonic point in the inner region ($\sim$ 0.01 kpc). The locus of this transonic point does not strongly depend on $\gamma$. We apply the polytropic model incorporating mass flux supplied by stellar components to the Sombrero Galaxy, and conclude that it can reproduce the observed gas density and the temperature distribution well. This result differs significantly from the isothermal model, which requires an unrealistically large mass flux (Igarashi et al. 2014). Thus, we conclude that the polytropic model is more realistic than the isothermal model, and that the Sombrero Galaxy can have a slowly accelerating outflow.

We examine the nature, possible orbits and physical properties of the progenitor of the North-western stellar stream (NWS) in the halo of the Andromeda galaxy (M31). The progenitor is assumed to be an accreting dwarf galaxy with globular clusters (GCs). It is, in general, difficult to determine the progenitor's orbit precisely because of many necessary parameters. Recently, Veljanoski et al. 2014 reported five GCs whose positions and radial velocities suggest an association with the stream. We use this data to constrain the orbital motions of the progenitor using test-particle simulations. Our simulations split the orbit solutions into two branches according to whether the stream ends up in the foreground or in the background of M31. Upcoming observations that will determine the distance to the NWS will be able to reject one of the two branches. In either case, the solutions require that the pericentric radius of any possible orbit be over 2 kpc. We estimate the efficiency of the tidal disruption and confirm the consistency with the assumption for the progenitor being a dwarf galaxy. The progenitor requires the mass $\ga 2\times10^6 M_{\sun}$ and half-light radius $\ga 30$ pc. In addition, $N$-body simulations successfully reproduce the basic observed features of the NWS and the GCs' line-of-sight velocities.

Motivated by a recently found interesting property of the dark halo surface density within a radius, $r_{\rm max}$, giving the maximum circular velocity, $V_{\rm max}$, we investigate it for dark halos of the Milky Way's and Andromeda's dwarf satellites based on cosmological simulations. We select and analyze the simulated subhalos associated with Milky Way-sized dark halos and find that the values of their surface densities, $\Sigma_{V_{\rm max}}$, are in good agreement with those for the observed dwarf spheroidal satellites even without employing any fitting procedures. This implies that this surface density would not be largely affected by any baryonic feedbacks and thus universal. Moreover, all subhalos on the small scales of dwarf satellites are expected to obey the relation $\Sigma_{V_{\rm max}}\propto V_{\rm max}$, irrespective of differences in their orbital evolutions, host halo properties, and observed redshifts. Therefore, we find that the universal scaling relation for dark halos on dwarf galaxy mass scales surely exists and provides us important clues to understanding fundamental properties of dark halos. We also investigate orbital and dynamical evolutions of subhalos to understand the origin of this universal dark halo relation and find that most of subhalos evolve generally along the $r_{\rm max}\propto V_{\rm max}$ sequence, even though these subhalos have undergone different histories of mass assembly and tidal stripping. This sequence, therefore, should be the key feature to understand the nature of the universality of $\Sigma_{V_{\rm max}}$.

The lightest Kaluza-Klein particle (LKP), which appears in the theory of universal extra dimensions, is one of good candidates for cold dark matter (CDM). When LKP pairs annihilate around the center of the Galaxy where CDM is concentrated, there are some modes which produce electrons and positrons as final products, and we categorize them into two components. One of them is the "Line" component, which directly annihilates into electron--positron pair. Another one is the "Continuum" component, which consists of secondarily produced electrons and positrons via some decay modes. Before reaching Earth, directions of electrons and positrons are randomized by the Galactic magnetic field, and their energies are reduced by energy loss mechanisms. We assume the LKP is in the mass range from 300 GeV to 1500 GeV. We calculate the electron plus positron spectrum after propagation in the Galactic halo to Earth, and we analyze the resulting spectrum and positron fraction. We also point out that the energy dependence of observed positron fraction is well reproduced by the mixture of "line" and "continuum" components. We can fit the electron plus positron spectrum and the positron fraction by assuming appropriate boost factors describing dark matter concentration in the Galactic halo. However, it is difficult to explain both the electron plus positron spectrum and the positron fraction by a single boost factor, if we take account of observational data obtained by AMS-02 only.

We report the result from observations conducted with the Atacama Large Millimeter/submillimeter Array (ALMA) to detect [CII] 158 um fine structure line emission from galaxies embedded in one of the most spectacular Lyman-alpha blobs (LABs) at z=3.1, SSA22-LAB1. Of three dusty star-forming galaxies previously discovered by ALMA 860 um dust continuum survey toward SSA22-LAB1, we detected the [CII] line from one, LAB1-ALMA3 at z=3.0993+/-0.0004. No line emission was detected, associated with the other ALMA continuum sources or from three rest-frame UV/optical selected z_spec~3.1 galaxies within the field of view. For LAB1-ALMA3, we find relatively bright [CII] emission compared to the infrared luminosity (L_[CII]/L_[CII]) and an extremely high [CII] 158 um and [NII] 205 um emission line ratio (L_[CII]/L_[NII]>55). The relatively strong [CII] emission may be caused by abundant photodissociation regions and sub-solar metallicity, or by shock heating. The origin of the unusually strong [CII] emission could be causally related to the location within the giant LAB, although the relationship between extended Lyman-alpha emission and ISM conditions of associated galaxies is yet to be understand.

We focus on the evidence of a past minor merger discovered in the halo of the Andromeda galaxy (M31). Previous N-body studies have enjoyed moderate success in producing the observed giant stellar stream (GSS) and stellar shells in M31's halo. The observed distribution of stars in the halo of M31 shows an asymmetric surface brightness profile across the GSS; however, the effect of the morphology of the progenitor galaxy on the internal structure of the GSS requires further investigation in theoretical studies. To investigate the physical connection between the characteristic surface brightness in the GSS and the morphology of the progenitor dwarf galaxy, we systematically vary the thickness, rotation velocity and initial inclination of the disc dwarf galaxy in N-body simulations. The formation of the observed structures appears to be dominated by the progenitor's rotation. Besides reproducing the observed GSS and two shells in detail, we predict additional structures for further observations. We predict the detectability of the progenitor's stellar core in the phase-space density distribution, azimuthal metallicity gradient of the western shell-like structure and an additional extended shell in the north-western direction that may constrain the properties of the progenitor galaxy.

We present upper limits in the hard X-ray and gamma-ray bands at the time of the LIGO gravitational-wave event GW 151226 derived from the CALorimetric Electron Telescope (CALET) observation. The main instrument of CALET, CALorimeter (CAL), observes gamma-rays from ~1 GeV up to 10 TeV with a field of view of ~2 sr. The CALET gamma-ray burst monitor (CGBM) views ~3 sr and ~2pi sr of the sky in the 7 keV - 1 MeV and the 40 keV - 20 MeV bands, respectively, by using two different scintillator-based instruments. The CGBM covered 32.5% and 49.1% of the GW 151226 sky localization probability in the 7 keV - 1 MeV and 40 keV - 20 MeV bands respectively. We place a 90% upper limit of 2 x 10^{-7} erg cm-2 s-1 in the 1 - 100 GeV band where CAL reaches 15% of the integrated LIGO probability (~1.1 sr). The CGBM 7 sigma upper limits are 1.0 x 10^{-6} erg cm-2 s-1 (7-500 keV) and 1.8 x 10^{-6} erg cm-2 s-1 (50-1000 keV) for one second exposure. Those upper limits correspond to the luminosity of 3-5 x 10^{49} erg s-1 which is significantly lower than typical short GRBs.

To unveil a progenitor of the Andromeda Giant Stellar Stream, we investigate the interaction between an accreting satellite galaxy and the Andromeda Galaxy using an $N$-body simulation. A comprehensive parameter study with 247 models is performed by varying size and mass distribution of the progenitor dwarf galaxy. We show that the binding energy of the progenitor is the crucial parameter in reproducing the Andromeda Giant Stellar Stream and the shell-like structures surrounding the Andromeda Galaxy. As a result of the simulations, the progenitor must satisfy a simple scaling relation between the core radius, the total mass and the tidal radius. Using this relation, we successfully constrain the physical properties of the progenitors to have mass ranging from $5\times10^8 M_\odot$ to $5\times10^9 M_\odot$ and central surface density around $10^3\, M_\odot\, \mathrm{pc}^{-2}$. A detailed comparison between our result and the nearby observed galaxies indicates that possible progenitors of the Andromeda Giant Stellar Stream include a dwarf elliptical galaxy, a dwarf irregular galaxy, and a small spiral galaxy.

We present statistical properties of diffuse Lyman-alpha halos (LAHs) around high-$z$ star-forming galaxies with large Subaru samples of Lyman-alpha emitters (LAEs) at $z=2.2$. We make subsamples defined by the physical quantities of LAEs' central Lyman-alpha luminosities, UV magnitudes, Lyman-alpha equivalent widths, and UV slopes, and investigate LAHs' radial surface brightness (SB) profiles and scale lengths $r_n$ as a function of these physical quantities. We find that there exist prominent LAHs around LAEs with faint Lyman-alpha luminosities, bright UV luminosities, and small Lyman-alpha equivalent widths in cumulative radial Lyman-alpha SB profiles. We confirm this trend with the anti-correlation between $r_n$ and Lyman-alpha luminosities (equivalent widths) based on the Spearman's rank correlation coefficient that is $\rho=-0.9$ ($-0.7$) corresponding to the $96\%$ ($93\%$) confidence level, although the correlation between $r_n$ and UV magnitudes is not clearly found in the rank correlation coefficient. Our results suggest that LAEs with properties similar to typical Lyman-break galaxies (with faint Lyman-alpha luminosities and small equivalent widths) possess more prominent LAHs. We investigate scenarios for the major physical origins of LAHs with our results, and find that the cold stream scenario is not preferred, due to the relatively small equivalent widths up to $77$\AA\ in LAHs that include LAEs' central components. There remain two possible scenarios of Lyman-alpha scattering in circum-galactic medium and satellite galaxies that cannot be tested with our observational data.

The lightest Kaluza-Klein particle (LKP), which appears in the theory of universal extra dimensions, is one of the good candidates for cold dark matter. The gamma-ray spectrum from annihilation of LKP dark matter shows a characteristic peak structure around the LKP mass. We investigate the detectability of this peak structure by considering energy resolution of near-future detectors, and calculate the expected count spectrum of the gamma-ray signal. In order to judge whether the count spectrum contains the LKP signal, the {\chi} squared test is employed. If the signal is not detected, we set some constraints on the boost factor that is an uncertain factor dependent on the substructure of the LKP distribution in the galactic halo. Detecting such peak structure would be conclusive evidence that dark matter is made of LKP.

During a search for gamma-ray emission from NGC 3628 (Arp 317), two new unidentified gamma-ray sources, Fermi J1049.7+0435 and J1103.2+1145 have been discovered \cite{ATel}. The detections are made in data from the Large Area Telescope (LAT), on board the Fermi Gamma-ray Space Telescope, in the 100\,MeV to 300\,GeV band during the period between 2008 August 5 and 2012 October 27. Neither is coincident with any source listed in the 2FGL catalogue \cite{Nolan2012}. Fermi J1049.7+0435 is at Galactic coordinates $(l,b) = (245.34^\circ, 53.27^\circ)$, $(\alpha_{J2000}, \delta_{J2000}) = (162.43^\circ, 4.60^\circ)$. Fermi J1103.2+1145 is at Galactic coordinates $(l,b) = (238.85^\circ, 60.33^\circ)$, $(\alpha_{J2000},\delta_{J2000})= (165.81^\circ, 11.75^\circ)$. Possible radio counterparts are found for both sources, which show flat radio spectra similar to other Fermi LAT detected AGN, and their identifications are discussed. These identification have been supoorted by snap-shot observations with the Australia Telescope Compact Array at several epochs in 2013 and 2014,

This study investigates the dynamical response of dark matter (DM) halos to recurrent starbursts in forming less-massive galaxies to solve the core-cusp problem. The gas, which is heated by supernova feedback after a starburst, expands and the star formation then terminates. This expanding gas loses energy by radiative cooling and then falls back toward the galactic center. Subsequently, the starburst is enhanced again. This cycle of expansion and contraction of the interstellar gas leads to a repetitive change in the gravitational potential of the gas. The resonance between DM particles and the density wave excited by the oscillating potential plays a key role in understanding the physical mechanism of the cusp-core transition of DM halos. DM halos effectively gain kinetic energy from the baryon potential through the energy transfer driven by the resonance between the particles and density waves. We determine that the critical condition for the cusp-core transition is such that the oscillation period of the gas potential is approximately the same as the local dynamical time of DM halos. We present the resultant core radius of a DM halo after the cusp-core transition induced by the resonance by using the conventional mass density profile predicted by the cold dark matter models. Moreover, we verify the analytical model by using $N$-body simulations, and the results validate the resonance model.

The cold dark matter (CDM) cosmology, which is the standard theory of the structure formation in the universe, predicts that the outer density profile of dark matter halos decreases with the cube of distance from the center. However, so far not much effort has examined this hypothesis. In the halo of the Andromeda galaxy (M31), large-scale stellar structures detected by the recent observations provide a potentially suitable window to investigate the mass--density distribution of the dark matter halo. We explore the density structure of the dark matter halo in M31 using an N-body simulation of the interaction between an accreting satellite galaxy and M31. To reproduce the Andromeda Giant Southern Stream and the stellar shells at the east and west sides of M31, we find the sufficient condition for the power-law index $\alpha$ of the outer density distribution of the dark matter halo. The best-fit parameter is $\alpha=-3.7$, which is steeper than the CDM prediction.

We have classified possible transonic solutions of galactic outflows in the gravitational potential of the dark matter halo (DMH) and super massive black hole (SMBH) under the assumptions of isothermal, spherically symmetric and steady state. It is clarified that the gravity of SMBH adds a new branch of transonic solutions with the transonic point in very close proximity to the centre in addition to the outer transonic point generated by the gravity of DMH. Because these two transonic solutions have substantially different mass fluxes and starting points, these solutions may have different influences on the evolution of galaxies and the release of metals into intergalactic space. We have applied our model to the Sombrero galaxy and obtained a new type of galactic outflow: a slowly accelerated transonic outflow through the transonic point at very distant region ($\simeq 126$\ kpc). In this galaxy, previous works reported that although the trace of the galactic outflow is observed by X-ray, the gas density distribution is consistent with the hydrostatic state. We have clarified that the slowly accelerating outflow has a gas density profile quite similar to that of the hydrostatic solution in the widely spread subsonic region. Thus, the slowly accelerating transonic solution cannot be distinguished from the hydrostatic solution in the observed region ($\leq 25$\ kpc) even if slow transonic flow exists. Our model provides a new perspective of galactic outflows and is applicable even to quiescent galaxies with inactive star formation.

Galaxies and massive black holes (BHs) presumably grow via galactic merging events and subsequent BH coalescence. As a case study, we investigate the merging event between the Andromeda galaxy (M31) and a satellite galaxy. We compute the expected observational appearance of the massive BH that was at the center of the satellite galaxy prior to the merger, and is currently wandering in the M31 halo. We demonstrate that a radiatively inefficient accretion flow with a bolometric luminosity of a few tens of solar luminosities develops when Hoyle-Lyttleton accretion onto the BH is assumed. We compute the associated broadband spectrum and show that the radio band (observable with EVLA, ALMA and SKA) is the best frequency range to detect the emission. We also evaluate the mass and the luminosity of the stars bound by the wandering BH and find that such a star cluster is sufficiently luminous that it could correspond to one of the star clusters found by the PAndAS survey. The discovery of a relic massive BH wandering in a galactic halo will provide a direct means to investigate in detail the coevolution of galaxies and BHs. It also means a new population of BHs (off-center massive BHs), and offers targets for clean BH imaging that avoids strong interstellar scattering in the center of galaxies.

We have categorized possible transonic solutions of galactic outflows in the gravitational potential of DMH and SMBH using the isothermal, spherically symmetric and steady model. We conclude that the gravitational potential of SMBH generates a new transonic branch while Tsuchiya et al. (2013) concluded that the gravitational potential of DMH forms one transonic solution. Because these two transonic solutions have different mass fluxes and starting points, these solutions will make different influences to the star formation rate, the evolution of galaxies, and the chemical evolution of the intergalactic medium. Therefore, we conclude that the influence of galactic outflows to the intergalactic medium depends not only on the mass distribution but also on the selected transonic solution. In addition, we have estimated range of parameters (KDMH; KBH) for actual galaxies. Moreover, it may be possible to estimate the galactic mass distributions of DMH and SMBH applying the model to the observed profile of the outflow velocity. Although it is difficult to determine the velocity of hot gas in the galactic halos from the current X-ray observations, but the next-generation X-ray observatory will be able to detect the detailed profiles of outflow velocities.

We study fundamental properties of isothermal, steady and spherically symmetric galactic outflow in the gravitational potential of a cold dark matter halo and a central super-massive black hole. We find that there are two transonic solutions having different properties: each solution is mainly produced by the dark matter halo and the super-massive black hole, respectively. Furthermore, we apply our model to the Sombrero galaxy. In this galaxy, Chandra X-ray observatory detected the diffuse hot gas as the trace of galactic outflows while the star-formation rate is low and the observed gas density distribution presumably indicates the hydrostatic equilibrium. To solve this discrepancy, we propose a solution that this galaxy has a transonic outflow, however, the transonic point forms in a very distant region from the galactic center (?$\sim$ 127 kpc). In this slowly accelerated transonic outflow, the outflow velocity is less than the sound velocity for most of the galactic halo. Since the gas density distribution in this subsonic region is similar to the hydrostatic one, it is difficult to distinguish the wide subsonic region from hydrostatic state. Such galactic outflows are dfferent from the conventional supersonic outflows observed in star-forming galaxies.

We present a statistical study of velocities of Lya, interstellar (IS) absorption, and nebular lines and gas covering fraction for Lya emitters (LAEs) at z~2. We make a sample of 22 LAEs with a large Lya equivalent width (EW) of > 50A based on our deep Keck/LRIS observations, in conjunction with spectroscopic data from the Subaru/FMOS program and the literature. We estimate the average velocity offset of Lya from a systemic redshift determined with nebular lines to be dv_Lya=234+-9 km s-1. Using a Kolmogorv-Smirnov test, we confirm the previous claim of Hashimoto et al. (2013) that the average dv_Lya of LAEs is smaller than that of LBGs. Our LRIS data successfully identify blue-shifted multiple IS absorption lines in the UV continua of four LAEs on an individual basis. The average velocity offset of IS absorption lines from a systemic redshift is dv_IS=204+-27 km s-1, indicating LAE's gas outflow with a velocity comparable to typical LBGs. Thus, the ratio, R^Lya_ IS = dv_Lya/dv_IS of LAEs, is around unity, suggestive of low impacts on Lya transmission by resonant scattering of neutral hydrogen in the IS medium. We find an anti-correlation between Lya EW and the covering fraction, f_c, estimated from the depth of absorption lines, where f_c is an indicator of average neutral hydrogen column density, N_HI. The results of our study support the idea that N_HI is a key quantity determining Lya emissivity.

We present diffuse Lyman-alpha halos (LAHs) identified in the composite Subaru narrowband images of 100-3600 Lyman-alpha emitters (LAEs) at z=2.2, 3.1, 3.7, 5.7, and 6.6. First, we carefully examine potential artifacts mimicking LAHs that include a large-scale point-spread function (PSF) made by instrumental and atmospheric effects. Based on our critical test with composite images of non-LAE samples whose narrowband-magnitude and source-size distributions are the same as our LAE samples, we confirm that no artifacts can produce a diffuse extended feature similar to our LAHs. After this test, we measure the scale lengths of exponential profile for the LAHs estimated from our z=2.2-6.6 LAE samples of L(Lyman-alpha) > 2 x 10^42 erg s^-1. We obtain the scale lengths of ~ 5-10 kpc at z=2.2-5.7, and find no evolution of scale lengths in this redshift range beyond our measurement uncertainties. Combining this result and the previously-known UV-continuum size evolution, we infer that the ratio of LAH to UV-continuum sizes is nearly constant at z=2.2-5.7. The scale length of our z=6.6 LAH is larger than 5-10 kpc just beyond the error bar, which is a hint that the scale lengths of LAHs would increase from z=5.7 to 6.6. If this increase is confirmed by future large surveys with significant improvements of statistical and systematical errors, this scale length change at z > 6 would be a signature of increasing fraction of neutral hydrogen scattering Lyman-alpha photons, due to cosmic reionization.

We present the results of structure analyses for a large sample of 426 Lya emitters (LAEs) at z~2.2 that are observed with HST/ACS and WFC3-IR by deep extra-galactic legacy surveys. We confirm that the merger fraction and the average ellipticity of LAE's stellar component are 10-30 % and 0.4-0.6, respectively, that are comparable with previous study results. We successfully identify that some LAEs have a spatial offset between Lya and stellar-continuum emission peaks, d_Lya, by ~2.5-4 kpc beyond our statistical errors. To uncover the physical origin of strong Lya emission found in LAEs, we investigate Lya equivalent width (EW) dependences of these three structural parameters, merger fraction, d_Lya, and ellipticity of stellar distribution in the range of EW(Lya)=20-250A. Contrary to expectations, we find that merger fraction does not significantly increase with Lya EW. We reveal an anti-correlation between d_Lya and EW(Lya) by Kolmogorov-Smirnov (KS) test. There is a trend that the LAEs with a large Lya EW have a small ellipticity. This is consistent with the recent theoretical claims that Lya photons can more easily escape from face-on disks having a small ellipticity, due to less inter-stellar gas along the line of sight, although our KS test indicates that this trend is not statistically significant. Our results of Lya-EW dependence generally support the idea that an HI column density is a key quantity determining Lya emissivity.

Observations have revealed interesting universal properties of dark matter (DM) halos especially around low-mass galaxies. Strigari et al. (2008) showed that DM halos have common enclosed masses within 300pc (Strigari relation). Kormendy & Freeman (2004) reported DM halos having almost identical central surface densities (the $\mu_{\rm 0D}$ relation). In addition, there exists a core-cusp problem, a discrepancy of the central density distribution between simulated halos and observations. We investigate whether a scenario where cuspy halos transform into cores by some dynamical processes can also explain their universal structural properties. It is shown that a cusp-to-core transformation model naturally reproduces the $\mu_{\rm 0D}$ relation and that Strigari relation follows from the $\mu_{\rm 0D}$ relation for dwarf galaxies. We also show that the central densities of cored dark halos provide valuable information about their formation redshifts.

In the hierarchical structure formation scenario, galaxies enlarge through multiple merging events with less massive galaxies. In addition, the Magorrian relation indicates that almost all galaxies are occupied by a central supermassive black hole (SMBH) of mass $10^{-3}$ of its spheroidal component. Consequently, SMBHs are expected to wander in the halos of their host galaxies following a galaxy collision, although evidence of this activity is currently lacking. We investigate a current plausible location of an SMBH wandering in the halo of the Andromeda galaxy (M31). According to theoretical studies of $N$-body simulations, some of the many substructures in the M31 halo are remnants of a minor merger occurring about 1 Gyr ago. First, to evaluate the possible parameter space of the infalling orbit of the progenitor, we perform numerous parameter studies using a Graphics Processing Unit (GPU) cluster. To reduce uncertainties in the predicted position of the expected SMBH, we then calculate the time evolution of the SMBH in the progenitor dwarf galaxy from $N$-body simulations using the plausible parameter sets. Our results show that the SMBH lies within the halo ($\sim$20--50 kpc from the M31 center), closer to the Milky Way than the M31 disk. Furthermore, the predicted current positions of the SMBH were restricted to an observational field of $0\degr.6 \times 0\degr.7$ in the northeast region of the M31 halo. We also discuss the origin of the infalling orbit of the satellite galaxy and its relationships with the recently discovered vast thin disk plane of satellite galaxies around M31.

The Second Fermi LAT source catalog (Nolan et al. 2012) includes as many as 1,873 sources, but initial attempts to identify counterparts at other wavelengths resulted in 575 sources remaining unidentified. The 2FGL catalog is based on the first 24 months of LAT observation since its launch in 2008, but the LAT has now accumulated more than 5 years of high-energy gamma-ray data almost flawlessly, presenting the possibility of finding new sources which were too faint to be detected in the first two years of data or showed flaring activity after the catalog was created. In this paper we report on two new gamma-ray sources serendipitously discovered in the constellation Leo and discuss possible counterparts based on radio observations.

We present the first systematic search for extended metal-line [OII]{\lambda}{\lambda}3726,3729 nebulae, or [OII] blobs (OIIBs), at z=1.2 using deep narrowband imaging with a survey volume of 1.9x10^5 Mpc^3 on the 0.62 deg^2 sky of Subaru-XMM Deep Survey (SXDS) field. We discover a giant OIIB, dubbed 'OIIB 1', with a spatial extent over ~75 kpc at a spectroscopic redshift of z=1.18, and also identify a total of twelve OIIBs with a size of >30 kpc. Our optical spectrum of OIIB 1 presents [NeV]{\lambda}3426 line at the 6{\sigma} level, indicating that this object harbors an obscured type-2 AGN. The presence of gas outflows in this object is suggested by two marginal detections of FeII{\lambda}2587 absorption and FeII*{\lambda}2613 emission lines both of which are blueshifted at as large as 500-600 km/s, indicating that the heating source of OIIB 1 is AGN or associated shock excitation rather than supernovae produced by starbursts. The number density of OIIB 1-type giant blobs is estimated to be ~5x10^{-6} Mpc^{-3} at z~1.2, which is comparable with that of AGNs driving outflow at a similar redshift, suggesting that giant OIIBs are produced only by AGN activity. On the other hand, the number density of small OIIBs, 6x10^{-5} Mpc^{-3}, compared to that of z~1 galaxies in the blue cloud in the same M_B range, may imply that 3% of star-forming galaxies at z~1 are quenching star formation through outflows involving extended [OII] emission.

We have observed the lightcurves of 13 V-type asteroids ((1933) Tinchen, (2011) Veteraniya, (2508) Alupka, (3657) Ermolova, (3900) Knezevic, (4005) Dyagilev, (4383) Suruga, (4434) Nikulin, (4796) Lewis, (6331) 1992 $\mathrm{FZ_{1}}$, (8645) 1998 TN, (10285) Renemichelsen, and (10320) Reiland). Using these observations we determined the rotational rates of the asteroids, with the exception of Nikulin and Renemichelsen. The distribution of rotational rates of 59 V-type asteroids in the inner main belt, including 29 members of the Vesta family that are regarded as ejecta from the asteroid (4) Vesta, is inconsistent with the best-fit Maxwellian distribution. This inconsistency may be due to the effect of thermal radiation Yarkovsky--O'Keefe--Radzievskii--Paddack (YORP) torques, and implies that the collision event that formed V-type asteroids is sub-billion to several billion years in age.

Radiation from the blazar class of of active galactic nuclei (AGN) exhibits fast time variability which is usually ascribed to instabilities in the emission region near the central supermassive black hole. The variability time scale is generally faster in higher energy region, and data recently provided by the {\it Fermi} Gamma-ray Space Telescope in the GeV energy band enable a detailed study of the temporal behavior of AGN. Due to its wide field-of-view in the scanning mode, most sky regions are observed for several hours per day and daily light curves of many AGN have been accumulated for more than 4 years. In this paper we investigate the time variability of 15 well-detected AGNs by studying the normalized power spectrum density (NPSD) of their light curves in the GeV energy band. One source, 3C 454.3, shows a specific time scale of $6.8\times10^5$ s, and this value suggests, assuming the internal shock model, a mass for the central black hole of $(10^8$--$10^{10})M_\odot$ which is consistent with other estimates. It also indicates the typical time interval of ejected blobs is $(7$--$70)$ times the light crossing time of the Schwarzschild radius.

Lower bounds are derived on the amplitude B of intergalactic magnetic fields (IGMFs) in the region between Galaxy and the blazar Mrk 421, from constraints on the delayed GeV pair-echo flux that are emitted by secondary e^-e^+ produced in \gamma\gamma interactions between primary TeV gamma-rays and the cosmic infrared background. The distribution of galaxies mapped by the Sloan Digital Sky Survey shows that this region is dominated by a large intergalactic void. We utilize data from long-term, simultaneous GeV-TeV observations by the Fermi Large Area Telescope and the ARGO-YBJ experiment extending over 850 days. For an assumed value of B, we evaluate the daily GeV pair-echo flux expected from the TeV data, select the dates where this exceeds the Fermi 2-\sigma sensitivity, compute the probability that this flux is excluded by the Fermi data for each date, and then combine the probabilities using the inverse normal method. Consequently, we exclude B < 10^{-20.5} G for a field coherence length of 1 kpc at 4- \sigma level, as long as plasma instabilities are unimportant for cooling of the pair beam. This is much more significant than the 2-\sigma bounds we obtained previously from observations of Mrk 501, by virtue of more extensive data from the ARGO-YBJ, as well as improved statistical analysis. Compared with most other studies of IGMF bounds, the evidence we present here for a non-zero IGMF is more robust as it does not rely on unproven assumptions on the primary TeV emission during unobserved periods.

Large-scale faint structure detected by the recent observations in the halo of the Andromeda galaxy (M31) provides an attractive window to explore the structure of outer cold dark matter (CDM) halo in M31. Using an N-body simulation of the interaction between an accreting satellite galaxy and M31, we investigate the mass density profile of the CDM halo. We find the sufficient condition of the outer density profile of CDM halo in M31 to reproduce the Andromeda giant stream and the shells at the east and west sides of M31. The result indicates that the density profile of the outer dark matter halo of M31 is a steeper than the prediction of the theory of the structure formation based on the CDM model.

We study fundamental properties of steady, spherically symmetric, isothermal galactic outflow in appropriate gravitational potential models. We aim at constructing a universal scale free theory not only for galactic winds, but also for winds from clusters/groups of galaxies. In particular, we consider effects of mass-density distribution on the formation of transonic galactic outflows under several models of the density distribution profile predicted by cosmological simulations of structure formation based on the cold dark matter (CDM) scenario. In this study, we have clarified that there exists two types of transonic solutions: outflows from the central region and from distant region with a finite radius, depending upon the density distribution of the system. The system with sufficiently steep density gradient at the center is allowed to have the transonic outflows from the center. The resultant criterion intriguingly indicates that the density gradient at the center must be steeper than that of the prediction of conventional CDM model including Navarro, Frenk & White (1997) and Moore et al. (1999). This result suggests that an additional steeper density distribution originated by baryonic systems such as the stellar component and/or the central massive black hole is required to realize transonic outflow from the central region. On the other hand, we predict the outflow, which is started at the outskirts of the galactic center and is slowly-accelerated without any drastic energy injection like starburst events. These transonic outflows may contribute secularly to the metal enrichment of the intergalactic medium.

GeV gamma-ray emission from two gamma-ray binary candidates, HESS J0632+057 and AGL J2241+4454, which were recently reported by H.E.S.S. and AGILE, respectively, have been searched for using the Fermi-LAT archival dataset. Spatial and temporal distribution of gamma-ray events are studied, but there was no evidence for GeV gamma-ray signal from either sources.

A heterogeneous CPU-GPU node is getting popular in HPC clusters. We need to rethink algorithms and optimization techniques for such system depending on the relative performance of CPU vs. GPU. In this paper, we report a performance optimized particle simulation code "OTOO", that is based on the octree method, for heterogenous systems. Main applications of OTOO are astrophysical simulations such as N-body models and the evolution of a violent merger of stars. We propose optimal task split between CPU and GPU where GPU is only used to compute the calculation of the particle force. Also, we describe optimization techniques such as control of the force accuracy, vectorized tree walk, and work partitioning among multiple GPUs. We used OTOO for modeling a merger of two white dwarf stars and found that OTOO is powerful and practical to simulate the fate of the process.

Using stacks of Ly-a images of 2128 Ly-a emitters (LAEs) and 24 protocluster UV-selected galaxies (LBGs) at z=3.1, we examine the surface brightness profiles of Ly-a haloes around high-z galaxies as a function of environment and UV luminosity. We find that the slopes of the Ly-a radial profiles become flatter as the Mpc-scale LAE surface densities increase, but they are almost independent of the central UV luminosities. The characteristic exponential scale lengths of the Ly-a haloes appear to be proportional to the square of the LAE surface densities (r(Lya) \propto Sigma(LAE)^2). Including the diffuse, extended Ly-a haloes, the rest-frame Ly-a equivalent width of the LAEs in the densest regions approaches EW_0(Lya) ~ 200 A, the maximum value expected for young (< 10^7 yr) galaxies. This suggests that Ly-a photons formed via shock compression by gas outflows or cooling radiation by gravitational gas inflows may partly contribute to illuminate the Ly-a haloes; however, most of their Ly-a luminosity can be explained by photo-ionisation by ionising photons or scattering of Ly-a photons produced in HII regions in and around the central galaxies. Regardless of the source of Ly-a photons, if the Ly-a haloes trace the overall gaseous structure following the dark matter distributions, it is not surprising that the Ly-a spatial extents depend more strongly on the surrounding Mpc-scale environment than on the activities of the central galaxies.

Temporal changes of X-ray to very-high-energy gamma-ray emissions from the pulsar-Be star binary PSR B1259-63/LS 2883 are studied based on 3-D SPH simulations of pulsar wind interaction with Be-disk and wind. We focus on the periastron passage of the binary and calculate the variation of the synchrotron and inverse-Compton emissions using the simulated shock geometry and pressure distribution of the pulsar wind. The characteristic double-peaked X-ray light curve from observations is reproduced by our simulation under a dense Be disk condition (base density ~10^{-9} g cm^{-3}). We interpret the pre- and post-periastron peaks as being due to a significant increase in the conversion efficiency from pulsar spin down power to the shock-accelerated particle energy at orbital phases when the pulsar crosses the disk before periastron passage, and when the pulsar wind creates a cavity in the disk gas after periastron passage, respectively. On the contrary, in the model TeV light curve, which also shows a double peak feature, the first peak appears around the periastron phase. The possible effects of cooling processes on the TeV light curve are briefly discussed.

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.

Based on a three-dimensional model of an early star-forming galaxy, we explore the evolution of the sub-millimeter brightness. The model galaxy is employed from an ultra-high-resolution chemodynamic simulation of a primordial galaxy by Mori & Umemura, where the SFR is ~ 10 Msun/yr at t<0.3 Gyr and several Msun/yr at t>0.3 Gyr. The former phase well reproduces the observed properties of LAEs and the latter does LBGs. We solve the three-dimensional radiative transfer in the clumpy interstellar media in this model galaxy, taking the size distributions of dust grains into account, and calculate the dust temperature as a function of galactic evolutionary time. We find that the clumpiness of interstellar media plays an important role for the sub-millimeter brightness. In the LAE phase, dust grains are concentrated on clumpy star-forming regions that are distributed all over the galaxy, and the grains can effectively absorb UV radiation from stars. As a result, the dust is heated up to T>35 K. In the LBG phase, the continuous supernovae drive dust grains far away from star-forming regions. Then, the grains cannot absorb much radiation from stars, and becomes into a cold state close to the CMB temperature. Consequently, the dust temperature decreases with the evolutionary time, where the mass-weighted mean temperature is T=26 K at t=0.1 Gyr and T=21 K at t=1.0 Gyr. By this analysis, it turns out that the sub-millimeter brightness is higher in the LAE phase than that in the LBG phase, although the dust-to-gas ratio increases monotonically as a function of time. We derive the spectral energy distributions by placing the model galaxy at a given redshift. The peak flux at 850 micron is found to be S_850 ~ 0.2 - 0.9 mJy if the model galaxy is placed at 6>z>2. This means that ALMA can detect an early star-forming galaxy with SFR of ~ 10 Msun/yr by less than one hour integration with 16 antennas.

PSR B1259-63/SS2883 is a binary system which consists of a 48-ms radio pulsar and a massive star in a highly eccentric orbit with a period of about 3.4 years. Non-pulsed and non-thermal emissions from this binary have been reported in the radio, X-ray and TeV gamma-ray energy ranges. Light curves in the radio and X-ray bands showed characteristic double-peaked features which can be attibuted to the interactions of the pulsar wind and the Be disk during the crossings by the pulsar. The TeV light curves around periastron differ between 2004 and 2007 observations, and the feature is not conclusive. We report a detection of GeV gamma-ray emission around the periastron passage in December 2010 with Fermi-LAT.

We derive lower bounds on intergalactic magnetic fields (IGMFs) from upper limits on the pair echo emission from the blazar Mrk 501, that is, delayed GeV emission from secondary $e^{-}e^{+}$ pairs produced via interactions of primary TeV gamma rays with the cosmic infrared background. Utilizing only simultaneous GeV-TeV light curves observed by VERITAS, MAGIC and {\it Fermi}-LAT during a multiwavelength campaign in 2009 that included a TeV flare, bounds are deduced on the IGMF strength of $B \gtrsim 10^{-20} {\rm G}$ at 90% confidence level for a field coherence length of 1 kpc. Since our analysis is based firmly on the observational data alone and nearly free of assumptions concerning the primary TeV flux in unobserved periods or spectral bands, our evaluation of the pair echo flux is conservative and the evidence for a non-zero IGMF is more robust compared to previous studies.

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.

According to recent models, the accretion disk and black hole in active galactic nuclei are surrounded by a clumpy torus. We investigate the NIR flux variation of the torus in response to a UV flash for various geometries. Anisotropic illumination by the disk and the torus self-occultation contrast our study with earlier works. Both the waning effect of each clump and the torus self-occultation selectively reduce the emission from the region with a short delay. Therefore, the NIR delay depends on the viewing angle (where a more inclined angle leads to a longer delay) and the time response shows an asymmetric profile with a negative skewness, opposing to the results for optically thin tori. The range of the computed delay coincides with the observed one, suggesting that the viewing angle is primarily responsible for the scatter of the observed delay. We also propose that the red NIR-to-optical color of type-1.8/1.9 objects is caused by not only the dust extinction but also the intrinsically red color. Compared with the modest torus thickness, both a thick and a thin tori display the weaker NIR emission. A selection bias is thus expected such that NIR-selected AGNs tend to possess moderately thick tori. A thicker torus shows a narrower and more heavily skewed time profile, while a thin torus produces a rapid response. A super-Eddington accretion rate leads to a much weaker NIR emission due to the disk self-occultation and the disk truncation by the self-gravity. A long delay is expected from an optically thin and/or a largely misaligned torus. A very weak NIR emission, such as in hot-dust-poor active nuclei, can arise from a geometrically thin torus, a super-Eddington accretion rate or a slightly misaligned torus.

The core-cusp problem remains as one of the unsolved discrepancies between observations and theories predicted by the standard paradigm of cold dark matter (CDM) cosmology. To solve this problem, we perform N-body simulations to study the nonlinear response of CDM halos to the variance of the gravitational potential induced by gas removal from galaxy centers. In this study, we focus on the timescale of the gas ejection, which is strongly correlated with stellar activities, and demonstrate that it is one of the key factors in determining the dynamical response of CDM halos. The results of simulations show that the power-low index of the mass-density profile of the dark matter halo correlated with the timescale of the mass loss, and it is flatter when the mass loss occurs over a short time than when it occurs over a long time. However, it is still larger than typical observational values; in other words, the central cusp remains for any mass loss model in the simulations. Moreover, for the slow mass-loss case, the final density profile of the dark matter halo recovers the universal density profiles predicted by the CDM cosmology. Therefore, mass loss driven by stellar feedback may not be an effective mechanism to flatten the central cusp.

A sublimation process governs the innermost region of the dusty torus of active galactic nuclei. However, the observed inner radius of the torus is systematically smaller than the expected radius by a factor of ~ 1/3. We show that the anisotropy of the emission from accretion disks resolves this conflict naturally and quantitatively. An accretion disk emits lesser radiation in the direction closer to its equatorial plane (i.e., to the torus). We find that the anisotropy makes the torus inner region closer to the central black hole and concave. Moreover, the innermost edge of the torus may connect with the outermost edge of the disk continuously. Considering the anisotropic emission of each clump in the torus, we calculate the near-infrared flux variation in response to a UV flash. For an observer at the polar angle theta_obs = 25 deg, the centroid of the time delay is found to be 37% of the delay expected in the case of isotropic illumination, which explains the observed systematic deviation.

We present results of a survey for giant Ly-alpha nebulae (LABs) at z=3 with Subaru/Suprime-Cam. We obtained Ly-alpha imaging at z=3.09+-0.03 around the SSA22 protocluster and in several blank fields. The total survey area is 2.1 square degrees, corresponding to a comoving volume of 1.6 x 10^6 Mpc^3. Using a uniform detection threshold of 1.4 x 10^{-18} erg s^{-1} cm^{-2} arcsec^{-2} for the Ly-alpha images, we construct a sample of 14 LAB candidates with major-axis diameters larger than 100 kpc, including five previously known blobs and two known quasars. This survey triples the number of known LABs over 100 kpc. The giant LAB sample shows a possible "morphology-density relation": filamentary LABs reside in average density environments as derived from compact Ly-alpha emitters, while circular LABs reside in both average density and overdense environments. Although it is hard to examine the formation mechanisms of LABs only from the Ly-alpha morphologies, more filamentary LABs may relate to cold gas accretion from the surrounding inter-galactic medium (IGM) and more circular LABs may relate to large-scale gas outflows, which are driven by intense starbursts and/or by AGN activities. Our survey highlights the potential usefulness of giant LABs to investigate the interactions between galaxies and the surrounding IGM from the field to overdense environments at high-redshift.

The prospects for future blazar surveys by next-generation very high energy (VHE) gamma-ray telescopes, such as Advanced Gamma-ray Imaging System (AGIS) and Cherenkov Telescope Array (CTA), are investigated using the latest model of blazar luminosity function and its evolution which is in good agreement with the flux and redshift distribution of observed blazars as well as the extragalactic gamma-ray background. We extend and improve the template of spectral energy distributions (SEDs) based on the blazar SED sequence paradigm, to make it reliable also in the VHE bands (above 100 GeV) by comparing with the existing VHE blazar data. Assuming the planned CTA sensitivities, a blind survey using a total survey time of ~100 hrs could detect ~3 VHE blazars, with larger expected numbers for wider/shallower surveys. We also discuss following-up of Fermi blazars. Detectability of VHE blazars in the plane of Fermi flux and redshift is presented, which would be useful for future survey planning. Prospects and strategies are discussed to constrain the extragalactic background light (EBL) by using the absorption feature of brightest blazar spectra, as well as cut-offs in the redshift distribution. We will be able to get useful constraints on EBL by VHE blazars at different redshifts ranging 0.3-1 TeV corresponding to z=0.10-0.36.

A cluster of galaxies is a huge system bounded by gravitation, and cosmic rays are thought to be confined in the system, thus it should contain much non-thermal components. Many theories predict significant gamma-ray emission that could be detectable by state-of-the-art gamma-ray telescopes. Some clusters have already been studied by using Fermi gamma-ray space telescope in the GeV band and Cherenkov telescopes in the TeV band, but most clusters are not studied in both energy bands. Here I present results on GeV gamma-ray emission from clusters of galaxies which have been given upper limits by Cherenkov telescopes using Fermi archival data.

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.

In order to assess the contribution of Lyman break galaxies (LBGs) and Lyman alpha emitters (LAEs) at redshifts 3<z<7 to the ionization of intergalactic medium (IGM), we investigate the escape fractions of ionizing photons from supernova-dominated primordial galaxies by solving the three-dimensional radiative transfer. The model galaxy is employed from an ultra-high-resolution chemodynamic simulation of a primordial galaxy by Mori & Umemura (2006),which well reproduces the observed properties of LAEs and LBGs. The total mass of model galaxy is 10^{11}M_sun. We solve not only photo-ionization but also collisional ionization by shocks. In addition, according to the chemical enrichment, we incorporate the effect of dust extinction, taking the size distributions of dust into account. As a result, we find that dust extinction reduces the escape fractions by a factor 1.5-8.5 in the LAE phase and by a factor 2.5-11 in the LBG phase, while the collisional ionization by shocks increases the escape fractions by a factor ~ 2. The resultant escape fractions are 0.07-0.47 in the LAE phase and 0.06-0.17 in the LBG phase. These results are well concordant with the recent estimations derived from the flux ratio at 1500 angstrom to 900 angstrom of LAEs and LBGs. Combining the resultant escape fractions with the luminosity functions of LAEs and LBGs, we find that high-z LAEs and LBGs can ionize the IGM at z=3-5. However, ionizing radiation from LAEs as well as LBGs falls short to ionize the IGM at z>6. That implies that additional ionization sources may required at z>6.

A new calculation of nuclear enhancement factor, used in estimation of Galactic diffuse gamma-ray flux from proton-proton interaction in order to take account of heavy nuclei included in cosmic-rays and interstellar matter, is presented by use of a Monte Carlo simulator, DPMJET-3. A new value of 1.8-2.0 in the energy range of 6-1000 GeV/nucleon, slightly increasing with kinetic energy of projectile cosmic rays, is about 20% larger than previous estimates.

Black-hole accretion systems are known to possess several distinct modes (or spectral states), such as low/hard state, high/soft state, and so on. Since the dynamics of the corresponding flows is distinct, theoretical models were separately discussed for each state. We here propose a unified model based on our new, global, two-dimensional radiation-magnetohydrodynamic simulations. By controlling a density normalization we could for the first time reproduce three distinct modes of accretion flow and outflow with one numerical code. When the density is large (model A), a geometrically thick, very luminous disk forms, in which photon trapping takes place. When the density is moderate (model B), the accreting gas can effectively cool by emitting radiation, thus generating a thin disk, i.e., the soft-state disk. When the density is too low for radiative cooling to be important (model C), a disk becomes hot, thick, and faint; i.e., the hard-state disk. The magnetic energy is amplified within the disk up to about twice, 30%, and 20% of the gas energy in models A, B, and C, respectively. Notably, the disk outflows with helical magnetic fields, which are driven either by radiation pressure force or magnetic pressure force, are ubiquitous in any accretion modes. Finally, our simulations are consistent with the phenomenological alpha-viscosity prescription, that is, the disk viscosity is proportional to the pressure.

With the advent of imaging atmospheric Cherenkov telescopes in late 1980's, ground-based observation of TeV gamma-rays came into reality after struggling trials by pioneers for twenty years, and the number of gamma-ray sources detected at TeV energies has increased to be over seventy now. In this review, recent findings from ground-based very-high-energy gamma-ray observations are summarized (as of 2008 March), and up-to-date problems in this research field are presented.

The halo of M31 shows a wealth of substructures that are consistent with satellite accretion. Here we report on kinematic and abundance results from Keck/DEIMOS spectroscopy in the calcium triplet region of over 3500 red giant star candidates along the minor axis and in off-axis spheroid fields of M31. Our data reach out to large radial distances of 160 kpc. The derived velocity distributions show a kinematically cold substructure at 17 kpc that has been reported before. We devise an improved method to measure accurate metallicities from the calcium triplet in low signal-to-noise spectra using a coaddition of the individual lines. The resulting distribution leads us to note an even stronger gradient in the abundance distribution along M31's minor axis than previously detected. The mean metallicity in the outer halo reaches below -2 dex, with individual values as low as -2.6 dex. In the inner spheroid, at 17-19 kpc, we find a sharp decline of ~0.5 dex in metallicity, which roughly coincides with the edge of an extended disk, previously detected from star count maps. A large fraction of red giants in the most distant fields are likely members of M33's overlapping halo. A comparison of our velocities with those predicted by new N-body simulations argues that the event responsible for the giant Stream is most likely not responsible for the full population of the inner halo; we show further that the abundance distribution of the Stream is different from that of the inner halo, from which it becomes evident that the merger event that formed the outer halo cannot have contributed any significant material to the inner spheroid. All this evidence of severe structure changes in the halo suggests a high degree of infall and stochastic abundance accretion governing the build-up of M31's halo. (Abridged)

The interaction between an accreting satellite and the Andromeda galaxy (M31) has been studied analytically and numerically, using a high-resolution N-body simulation with $4\times10^7$ particles. For the first time, we show the self-gravitating response of the disk, the bulge, and the dark matter halo of M31 to an accreting satellite. We reproduce the Stream and the shells at the East and West side of M31, by following the evolution of the collision 4 Gyr into the future, and show that recently discovered diffuse arcs on the south side of the minor axis are the remnants of a similar collision that occurred 3-4 Gyr earlier than the Stream event. The present day integrity of the M31 disk constrains the satellite mass to be $M \leq 5\times10^9 M_\odot$. The stars that were originally in the center of the satellite are now in the east shell. Therefore, observations in this region might reveal additional clues about the nature of satellite, such as the central core and any metallicity gradient.

The current understanding of galaxy formation is that it proceeds in a 'bottom up' way, with the formation of small clumps of gas and stars that merge hierarchically until giant galaxies are built up. The baryonic gas loses the thermal energy by radiative cooling and falls towards the centres of the new galaxies, while supernovae (SNe) blow gas out. Any realistic model therefore requires a proper treatment of these processes, but hitherto this has been far from satisfactory. Here we report an ultra-high-resolution simulation that follows evolution from the earliest stages of galaxy formation through the period of dynamical relaxation. The bubble structures of gas revealed in our simulation ($< 3\times10^8$ years) resemble closely the high-redshift Lyman $\alpha$ emitters (LAEs). After $10^9$ years these bodies are dominated by stellar continuum radiation and look like the Lyman break galaxies (LBGs) known as the high-redshift star-forming galaxies at which point the abundance of elements heavier than helium ("metallicity") appears to be solar. After $1.3\times10^{10}$ years, these galaxies resemble present-day ellipticals.

Development of MITSuME is reported. Two 50-cm optical telescopes have been built at Akeno in Yamanashi prefecture and at Okayama Astrophysical Observatory (OAO) in Okayama prefecture. Three CCD cameras for simultaneous g'RcIc photometry are to be mounted on each focal plane, covering a wide FOV of about 30" x 30". The limiting magnitude at V is fainter than 18. In addition to these two optical telescopes, a 91-cm IR telescope with a 1 deg x 1 deg field of view is being built at OAO, which performs photometry in YJHK bands. These robotic telescopes can start the observation of counterparts of a GRB within a minute from an alert. We aim to obtain photometric redshifts exceeding 10 with these telescopes. The performance and the current construction status of the telescopes are presented.

We report an ultra-high-resolution simulation that follows evolution from the earliest stages of galaxy formation through the period of dynamical relaxation. The bubble structures of gas revealed in our simulation ($< 3\times 10^8$ years) resemble closely the high-redshift Lyman $\alpha$ emitters (LAEs). After $10^9$ years these bodies are dominated by stellar continuum radiation and look like the Lyman break galaxies (LBGs) known as the high-redshift star-forming galaxies at which point the abundance of elements appears to be solar. After $1.3\times10^{10}$ years, these galaxies resemble present-day ellipticals. The comparisons of simulation results with the observations of elliptical galaxies allow us to conclude that LAEs and LBGs are infants of elliptical galaxies or bulge systems in the nearby universe.

The Hard X-ray Detector (HXD) on board Suzaku covers a wide energy range from 10 keV to 600 keV by combination of silicon PIN diodes and GSO scintillators. The HXD is designed to achieve an extremely low in-orbit back ground based on a combination of new techniques, including the concept of well-type active shield counter. With an effective area of 142 cm^2 at 20 keV and 273 cm2 at 150 keV, the background level at the sea level reached ~1x10^{-5} cts s^{-1} cm^{-2} keV^{-1} at 30 keV for the PI N diodes, and ~2x10^{-5} cts s^{-1} cm^{-2} keV^{-1} at 100 keV, and ~7x10^{-6} cts s^{-1} cm^{-2} keV^{-1} at 200 keV for the phoswich counter. Tight active shielding of the HXD results in a large array of guard counters surrounding the main detector parts. These anti-coincidence counters, made of ~4 cm thick BGO crystals, have a large effective area for sub-MeV to MeV gamma-rays. They work as an excellent gamma-ray burst monitor with limited angular resolution (~5 degree). The on-board signal-processing system and the data transmitted to the ground are also described.

WIDGET is a robotic telescope for monitoring the HETE-2 field-of-view to detect Gamma-ray Burst optical flashes or possible optical precursors. The system has 62degx62deg wide field-of-view which covers about 80% of HETE-2 one with a 2kx2k Apogee U10 CCD camera and a Canon EF 24mm f/1.4 wide-angle lens without a bandpass filter. WIDGET has been in operation since June 2004 at Akeno observing site where is about 200 km apart from Tokyo. Typical limiting magnitude with S/N=3 at the site is V=10mag for 5 seconds exposure and V=11mag for 30 seconds exposure. We had already six coincident observations with HETE-2 position alerts. It was, however, cloudy for all cases due to rainy season in Japan. Expected number of coincident observations under clear sky is about 5 events per year. We will extend the system in early 2005 for Swift era to monitor optical transients in wider field-of-view, multi-color or polarization modes.

he formation and growth of supermassive black holes (SMBHs) physically linked with bulges are considered. We focus on the radiation hydrodynamic process for the growth of SMBH in the optically thick starburst phase, where radiation from bulge stars drives the mass accretion on to a galactic center through radiation drag effect. In the present scenario, the AGN luminosity-dominant phase (QSO phase) is preceded by the host luminosity-dominat phase, which is called ``proto-QSO phase''. In this phase, there exists the massive dusty disks within younger bulges. Also, the proto-QSO phase is anticipated by an optically-thick ultraluminous infrared galaxy (ULIRG) phase. Furthermore, such radiation hydrodynamic model has been also applied to disk galaxies. It turns out that the mass of a SMBH primarily correlates with a bulge component even in a disk galaxy. Thus, by analogy to proto-QSOs, the BH growing phase in disk galaxies may have massive dusty disks within younger bulges.

The quasi-steady structure of super-critical accretion flows around a black hole is studied based on the two-dimensional radiation-hydrodynamical (2D-RHD) simulations. The super-critical flow is composed of two parts: the disk region and the outflow regions above and below the disk. Within the disk region the circular motion as well as the patchy density structure are observed, which is caused by Kelvin-Helmholtz instability and probably by convection. The mass-accretion rate decreases inward, roughly in proportion to the radius, and the remaining part of the disk material leaves the disk to form outflow because of strong radiation pressure force. We confirm that photon trapping plays an important role within the disk. Thus, matter can fall onto the black hole at a rate exceeding the Eddington rate. The emission is highly anisotropic and moderately collimated so that the apparent luminosity can exceed the Eddington luminosity by a factor of a few in the face-on view. The mass-accretion rate onto the black hole increases with increase of the absorption opacity (metalicity) of the accreting matter. This implies that the black hole tends to grow up faster in the metal rich regions as in starburst galaxies or star-forming regions.

We have detected gamma-ray emission at the 6sigma level at energies greater than 500GeV from the supernova remnant RX J0852.0-4622 (G266.2-1.2) using the CANGAROO-II Imaging Atmospheric Cherenkov Telescope (IACT). The flux was 0.12 times of that of Crab at 1TeV. The signal centroid is consistent with the peak of the X-ray emission in the north-west rim of the remnant.

We consider a forming galaxy undergoing multitudinous supernova (SN) explosions, as a possible model of \lya blobs (LABs). For this purpose, an ultra--high resolution hydrodynamic simulation is performed using $1024^3$ grid points, where SN remnants are resolved with sufficient accuracy. It is found that multiple SN explosions produce kpc--size expanding hot bubbles, which drive cool, dense shells by strong shock. The colliding high--density cooling shells radiate intensive \lya emission, resulting in a high \lya luminosity of $\sim10^{43}$ erg s$^{-1}$, comparable to the observed level in LABs. Also, recently discovered bubbly features in some LABs are quite similar to the structure predicted in the present simulation. Furthermore, the result demonstrates that LABs are representative of evolving primordial galaxies; they could hold direct information on the early chemical enrichment of galaxies, contrary to present--day galaxies which have undergone intense recycling of interstellar matter, thus erasing most of the early chemical history. It turns out that the metal mixing proceeds in a very inhomogeneous fashion, so that there appears a large spread of metallicity, that is,$[{\rm Fe/H}] \approx 0 {\rm ~to} -5$ or $[{\rm O/H}] \approx 1 {\rm ~to} -4$. Hence, the early galactic chemical evolution may have proceeded in a different manner from that hitherto considered in one--zone models.

Gamma-ray spectra from cosmic-ray proton and electron interactions with dense gas clouds have been calculated using a Monte Carlo event simulation code, GEANT4. Such clouds are postulated as a possible form of baryonic dark matter in the Universe. The simulation fully tracks the cascade and transport processes that are important in a dense medium, and the resulting gamma-ray spectra are computed as a function of cloud column-density. These calculations are used for predicting the Galactic diffuse gamma-ray spectrum that may be contributed by baryonic dark matter; the results are compared with data from the EGRET instrument, and used to constrain the fraction of Galactic dark matter that may be in the form of dense gas clouds. In agreement with previous authors, we find useful constraints on the fraction of Galactic dark matter that may be in the form of low column-density clouds ($\Sigma \la 10 g cm^{-2}$). However, this fraction rises steeply in the region $\Sigma \sim 10^2 g cm^{-2}$, and for $\Sigma \ga 200 g cm^{-2}$ we find that baryonic dark matter models are virtually unconstrained by the existing gamma-ray data.

We have detected sub-TeV gamma-ray emission from the direction of the Galactic Center (GC) using the CANGAROO-II Imaging Atmospheric Cherenkov Telescope (IACT). We detected a statistically signicant excess at energies greater than 250GeV. The flux was one order of magnitude lower than that of Crab at 1 TeV with a soft spectrum E^{-4.6+-0.5}. The signal centroid is consistent with the GC direction and the observed prole is consistent with a point-like source. Our data suggests that the GeV source 3EG J1746-2851 is identical with this TeV source and we study the combined spectra to determine the possible origin of the gamma-ray emission. We also obtain an upper limit on the cold dark-matter density in the Galactic halo.

Models in which a large fraction of the Galactic dark matter takes the form of cold gas clouds imply that there is thermal microwave emission from the Galactic dark halo. Such models can therefore be directly constrained by data on the microwave sky, and in particular the very sensitive observations of microwave anisotropies which are now being made. To this end we have computed the anisotropy power-spectrum expected for a Galactic dark halo made of cold, dense gas clouds, including the effects of clustering with a CDM-like mass spectrum of mini-halo substructure. The power-spectrum displays two peaks: one, at l~50, is the Poisson noise for the mini-halos, and the second, much larger and at much higher l, is the Poisson noise of the individual clouds. Because it appears on small (milli-arcsecond) angular scales, where the instrumental sensitivity is inevitably very poor, the latter signal is not directly detectable. By contrast, clusters of cold gas clouds may contribute significantly to the observed anisotropies if their emission has a grey-body spectrum. In this case the peak fluctuation, at l~50, amounts to 4/|sin b| micro-K in the Rayleigh-Jeans limit, and is the dominant Galactic foreground between 40 and 80 GHz. It will be possible to constrain this foreground component using low-latitude data from the MAP satellite, providing that its spectrum conforms to a grey-body. If the spectrum is ``dusty'' there will be relatively little power at frequencies below the thermal peak, and in this case the predicted anisotropies are shown to be negligible.

Cold, dense gas clouds have been proposed as a major component of the Galactic dark matter; such clouds can be revealed by the gamma-ray emission which arises from cosmic-rays interacting with the gas. If this dark matter component is clustered then highly luminous GeV sources result, preferentially at low to mid Galactic latitudes where they lie within the cosmic-ray disk. The predicted emission for such clusters is steady, continuum emission with peak power emerging at several hundred MeV. These sources would not have obvious counterparts at other wavelengths and are thus of interest in connection with the UnIDentified (UID) EGRET sources. Here we present a Monte Carlo simulation of the gamma-ray source population due to cold gas clouds, assuming a Cold-Dark-Matter-like mass spectrum for the clustering. We find that ~280 EGRET sources are predicted by this model, with a median Galactic latitude of 12 degrees for the population, and a median angular size of 2.2 degrees. The latitude and size distributions are consistent with the UID EGRET source data, but the source counts are clearly overpredicted. On the basis of these results we propose that clusters of cold gas clouds comprise the majority population of the observed UID EGRET sources. Our interpretation implies that there should be microwave counterparts to most of the UID sources, and will thus be strongly constrained by data from the present generation of microwave anisotropy experiments.

An 80 ksec Chandra ACIS observation of the radio galaxy 3C 452 is reported. A diffuse X-ray emission associated with the lobes has been detected with high statistical significance, together with the X-ray nucleus of the host galaxy. The 0.5--5 keV ACIS spectrum of the diffuse emission is described by a two-component model, consisting of a soft thermal plasma emission from the host galaxy halo and a hard non-thermal power-law component. The hard component is ascribed to the inverse Comptonization of cosmic microwave background photons by the synchrotron emitting electrons in the lobes, because its spectral energy index, 0.68+-0.28, is consistent with the radio synchrotron index, 0.78. These results reveal a significant electron dominance in the lobes. The electrons are inferred to have a relatively uniform distribution, while the magnetic field is compressed toward the lobe periphery.

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.

Based on the radiation hydrodynamical model for the black hole(BH) growth,incorporated with the chemical evolution of the early-type host galaxy, we construct the coevolution model of a QSO BH and the host galaxy. As a result, it is found that after a galactic wind epoch,the luminosity is shifted from the host-dominant phase to the AGN-dominant phase (QSO phase) in the timescale of a few $10^{8}$ years.The former phase corresponds to the early stage of growing BH, and can be regarded as a ``proto-QSO'' phase. It has observable characteristic properties (detail inthis paper).By comparing these predictions with recent observations, radio galaxies are a possible candidate for proto-QSOs.Also, it is anticipated that the proto-QSO phase is preceded by an optically thick phase, which may correspond to ULIRGs.In this phase, $M_{\rm BH}/M_{\rm bulge}$ is predicted to be much less than $10^{-3}$ and grow with metallicity.Moreover, as precursors of ULIRGs, optically-thin star-forming galaxies are predicted. These may be in the assembly phase of Lyman break galaxies (LBGs) or Ly$\alpha$ emitters.

We investigate the photon-trapping effects in the super-critical black hole accretion flows by solving radiation transfer as well as the energy equations of radiation and gas. It is found that the slim-disk model generally overestimates the luminosity of the disk at around the Eddington luminosity (L_E) and is not accurate in describing the effective temperature profile, since it neglects time delay between energy generation at deeper inside the disk and energy release at the surface. Especially, the photon-trapping effects are appreciable even below L ~ L_E, while they appear above ~ 3L_E according to the slim disk. Through the photon-trapping effects, the luminosity is reduced and the effective temperature profile becomes flatter than r^{-3/4} as in the standard disk. In the case that the viscous heating is effective only around the equatorial plane, the luminosity is kept around the Eddington luminosity even at very large mass accretion rate, Mdot>>L_E/c^2. The effective temperature profile is almost flat, and the maximum temperature decreases in accordance with rise in the mass accretion rate. Thus, the most luminous radius shifts to the outer region when Mdot/(L_E/c^2) >> 10^2. In the case that the energy is dissipated equally at any heights, the resultant luminosity is somewhat larger than in the former case, but the energy-conversion efficiency still decreases with increase of the mass accretion rate, as well. The most luminous radius stays around the inner edge of the disk in the latter case. Hence, the effective temperature profile is sensitive to the vertical distribution of energy production rates, so is the spectral shape. Future observations of high L/L_E objects will be able to test our model.

We present results from three-dimensional numerical simulations of the dynamics of SN-driven bubbles as they propagate through and escape the grasp of subgalactic halos with masses M ~ 10^8 Msun at redshift z=9. The hydrodynamic simulations use a nested grid method to follow the evolution of explosive multi-SN events operating on the characteristic timescale of a few x 10^7 yr, the lifetime of massive stars. The results confirm that, if the star formation efficiency of subgalactic halos is close to 10%, a significant fraction of the halo gas will be lifted out of the potential well (`blow-away'), shock the intergalactic medium, and pollute it with metal-enriched material, a scenario recently advocated by Madau, Ferrara & Rees (2001). The volume filling factor of the ejecta is close to unity. Depending on the stellar distribution, we find that less than 30% of the available SN energy is converted into kinetic energy of the blown away material, the remainder being radiated away. It appears that mechanical feedback is less efficient than expected from simple energetic arguments, as off-nuclear SN explosions drive inward-propagating shocks that tend to collect and pile up cold gas in the central regions of the host halo. Low-mass galaxies at early epochs may survive multiple SN events and continue forming stars.

Recent results of GeV and TeV observations of gamma-rays from the Universe are briefly reviewed. Topics include observational technique, diffuse gamma-rays, pulsars, unidentified sources, plerions, supernova remnants and AGNs.

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.

Preliminary results from CANGAROO-II, a 10 m imaging Cherenkov telescope, in Woomera, South Australia are presented. They include the confirmation of detections of TeV gamma-ray sources we have reported using a 3.8 m telescope, CANGAROO-I. Also the status of the construction of an array of four 10 m telescopes, called CANGAROO-III, is reported. The first telescope of the array was upgraded from the 7 m telescope and the second one is being constructed in this year. The full array will be operating in 2004.

The CANGAROO-III project, a system of four 10 m telecopes dedicated for gamma-ray astrophysics, has started in 1999 and is expected to complete in 2004 in Woomera, South Australia. We report the construction work during the first year, which includes the completion of the first 10 m telescope built as a 7 m telescope in 1999, and the work in progess to increase the performance by constructing three more similar telescopes based on the experience gained during the construction of the first telescope.

A new imaging atmospheric Cherenkov telescope (CANGAROO-II) with a light-weight reflector has been constructed. Light, robust, and durable mirror facets of containing CFRP (Carbon Fiber Reinforced Plastic) laminates were developed for the telescope. The attitude of each facet can be adjusted by stepping motors. In this paper, we describe the design, manufacturing, alignment procedure, and the performance of the CANGAROO-II optical reflector system.

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 examine various physical processes associated with the formation of globular clusters by using the three-dimensional Smoothed Particle Hydrodynamics (SPH) code. Our code includes radiative cooling of gases, star formation, energy feedback from stars including stellar winds and supernovae, and chemical enrichment by stars. We assume that, in the collapsing galaxy, isothermal cold clouds form through thermal condensations and become proto-globular clouds. We calculate the size of proto-globular clouds by solving the linearized equations for perturbation. We compute the evolution of the inner region of the proto-cloud with our SPH code for various initial radius and initial composition of gases. When the initial gases contain no heavy elements, the evolution of proto-clouds sensitively depends on the initial radius. For a smaller initial radius, the initial star burst is so intense that the subsequent star formation occurs in the central regions to form a dense star cluster as massive as the globular cluster. When the initial gases contain some heavy elements, the metallicity of gases affects the evolution and the final stellar mass. If the initial radius of the proto-globular clouds was relatively large, the formation of a star cluster as massive as the globular clusters requires the initial metallicity as high as [Fe/H] $\geq -2$. The self-enrichment of heavy elements in the star cluster does not occur in all cases.

Many clusters of galaxies contain an appreciable amount of hot gas, the intracluster medium. As a consequence, gas will be stripped from galaxies that move through the inter cluster medium, if the ram pressure exceeds the internal gravitational force. We study the interaction between the intracluster medium and an extended gas component of dwarf galaxies confined by a surrounding cold dark matter halo analytically and numerically, using axisymmetric two-dimensional hydrodynamical simulations at high resolution. The results show that the gas within the dark matter halo is totally stripped in a typical galactic cluster. The process of ram pressure stripping therefore must have played an important role during the chemo-dynamical evolution of dwarf galaxies in galactic clusters. Our results predict that most of the chemical evolution of dwarf galaxies must have occurred at high redshift, before the virialized cluster had formed.

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.

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

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

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.

We have developed a new parallel tree method which will be called the forest method hereafter. This new method uses the sectional Voronoi tessellation (SVT) for the domain decomposition. The SVT decomposes a whole space into polyhedra and allows their flat borders to move by assigning different weights. The forest method determines these weights based on the load balancing among processors by means of the over-load diffusion (OLD). Moreover, since all the borders are flat, before receiving the data from other processors, each processor can collect enough data to calculate the gravity force with precision. Both the SVT and the OLD are coded in a highly vectorizable manner to accommodate on vector parallel processors. The parallel code based on the forest method with the Message Passing Interface is run on various platforms so that a wide portability is guaranteed. Extensive calculations with 15 processors of Fujitsu VPP300/16R indicate that the code can calculate the gravity force exerted on 10^5 particles in each second for some ideal dark halo. This code is found to enable an N-body simulation with 10^7 or more particles for a wide dynamic range and is therefore a very powerful tool for the study of galaxy formation and large-scale structure in the universe.

We simulate the dynamical, chemical, and spectro-photometric evolution of dwarf and normal elliptical galaxies embedded in a dark matter halo, using a three-dimensional N-body/SPH simulation code. For the forming dwarf elliptical galaxies, supernova-driven winds propagating outwards from inside the system collide with the infalling gas and produce the super-shell in which stars are formed. The resulting stellar system forms a loosely bound virialized system due to the significant mass loss and has a large velocity dispersion and a large core. Consequently the surface brightness distribution shows an exponential profile and the color distribution shows a positive gradient such that the colors become redder away from the galaxy center in agreement with observations. On the other hand, for the normal elliptical galaxies with deep gravitational potential, the mass loss out of the system does not have a significant dynamical effect on the stellar system. The resulting surface brightness distribution has a large central concentration like de Vaucouleurs' $r^{1/4}$-profile and the color distribution shows a negative gradient as observed. Our simulation shows that different features between dwarf and normal elliptical galaxies stem from different cooling efficiencies for their respective protogalaxies in a standard CDM universe.

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.

Several weeks after the explosion of supernova (SN) 1987A, the UV flash of the SN illuminated a ring-like structure in the circumstellar material at about 0.65 ly from the SN. The interaction between the stellar winds from the SN progenitor is considered to be the candidate for the formation of the circumstellar structure. In the case that the stellar winds are spherically symmetric, the interaction should result in a shell-like structure. However, in this paper we show that the magnetic field in the winds causes an anisotropy which leads to the formation of a ring-like structure. When the fast wind of the blue supergiant phase of the progenitor sweeps up the surrounding slow wind of the red-supergiant phase, the magnetic field as well as the wind material are piled up in the interaction region. Since the magnetic energy increases in proportion to the square of the amplitude, the magnetic field exhibits its effect prominently at the interaction region; due to the magnetic pressure force the material at lower latitudes is compressed into a ring-like structure. It is suggested that this magnetic process can also explain the newly observed pair of rings of the SN 1987A nebula.

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.

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.

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.

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.

We have developed Remote-GRAPE, a subroutine library to use the special purpose computer GRAPE-3A. The GRAPE-3A can efficiently calculate gravitational force between particles, and construct neighbor lists. All other calculations are performed on the host workstation (WS) which is directly connected to GRAPE. The use of GRAPE for Smoothed Particle Hydrodynamics (GRAPE-SPH) can in principle greatly speed up the calculations on WS. However the current bottleneck of GRAPE-SPH is that its performance is limited by the speed of the host WS. To solve this problem, we implement Remote-GRAPE; it allows us to run applications which use GRAPE-3A hardware on the significantly faster computers than the physical host WS. Thus, we can take advantage of the fast computers even though they can not physically be connected to GRAPE. The Remote-GRAPE system is implemented on the Parallel Virtual Machine (PVM). In this paper, the detail of implementation and performance analysis are described.

We simulate the dynamical and chemical evolution of a dwarf galaxy embedded in a dark matter halo, using a three-dimensional N-body/SPH simulation code combined with stellar population synthesis. The initial condition is adopted in accord with a $10^{10}M_\odot$ virialized sphere in a 1-sigma CDM perturbation which contains 10% baryonic mass. A supersonic spherical outflow is driven by the first star burst near the center of the galaxy and produces an expanding super shell in which stars are subsequently formed. Consecutive formation of stars in the expanding shell makes the stellar system settled with the exponential brightness profile, the positive metallicity gradient, and the inverse color gradient in agreement with observed features of dwarf galaxies. We therefore propose that the energy feedback via stellar winds and supernovae is a decisive mechanism for formation of less compact, small systems like dwarf galaxies.

A new calculation of the Galactic diffuse gamma-ray spectrum from the decay of secondary particles produced by interactions of cosmic-ray protons with interstellar matter is presented. The calculation utilizes the modern Monte Carlo event generators, Hadrin, Fritiof and Pythia, which simulate high-energy proton-proton collisions and are widely used in studies of nuclear and particle physics, in addition to scaling calculation. This study is motivated by the result on the Galactic diffuse gamma-ray flux observed by the EGRET detector on the Compton Gamma-ray Observatory, which indicates an excess above about 1 GeV of the observed intensity compared with a model prediction. The prediction is based on cosmic-ray interactions with interstellar matter, in which secondary pion productions are treated by a simple model. With the improved interaction model used here, however, the diffuse gamma-ray flux agrees rather well with previous calculations within uncertainties, which mainly come from the unobservable demodulated cosmic-ray spectrum in interstellar space. As a possible solution to the excess flux, flatter spectra of cosmic-ray protons have been tested and we found that the power-law spectrum with an index of about $-(2.4\sim2.5)$ gives a better fit to the EGRET data, though the spectrum is not explained completely.