On January 14, 2019, the Major Atmospheric Gamma Imaging Cherenkov telescopes detected GRB 190114C above 0.2 TeV, recording the most energetic photons ever observed from a gamma-ray burst. We use this unique observation to probe an energy dependence of the speed of light in vacuo for photons as predicted by several quantum gravity models. Based on a set of assumptions on the possible intrinsic spectral and temporal evolution, we obtain competitive lower limits on the quadratic leading order of speed of light modification.

An optical cluster finder inevitably suffers from projection effects, where it misidentifies a superposition of galaxies in multiple halos along the line-of-sight as a single cluster. Using mock cluster catalogs built from cosmological N-body simulations, we quantify the impact of these projection effects with a particular focus on the observables of interest for cluster cosmology, namely the cluster lensing and the cluster clustering signals. We find that "observed" clusters, i.e. clusters identified by our cluster finder algorithm, exhibit lensing and clustering signals that deviate from expectations based on a statistically isotropic halo model -- while both signals agree with halo model expectations on small scales, they show unexpected boosts on large scales, by up to a factor of 1.2 or 1.4 respectively. We identify the origin of these boosts as the inherent selection bias of optical cluster finders for clusters embedded within filaments aligned with the line-of-sight, and show that a minority ($\sim 30\%$) of such clusters within the entire sample is responsible for this observed boost. We discuss the implications of our results on previous studies of optical cluster, as well as prospects for identifying and mitigating projection effects in future cluster cosmology analyses.

Extreme High-frequency BL~Lacs (EHBL) feature their synchrotron peak of the broadband spectral energy distribution (SED) at $\nu_{\rm s} \geq $10$^{17}$\,Hz. The BL~Lac object 1ES~2344+514 was included in the EHBL family because of its impressive shift of the synchrotron peak in 1996. During the following years, the source appeared to be in a low state without showing any extreme behaviours. In August 2016, 1ES~2344+514 was detected with the ground-based $\gamma$-ray telescope FACT during a high $\gamma$-ray state, triggering multi-wavelength (MWL) observations. We studied the MWL light curves of 1ES~2344+514 during the 2016 flaring state, using data from radio to VHE $\gamma$ rays taken with OVRO, KAIT, KVA, NOT, some telescopes of the GASP-WEBT collaboration at the Teide, Crimean, and St. Petersburg observatories, \textit{Swift}-UVOT, \textit{Swift}-XRT, \textit{Fermi}-LAT, FACT and MAGIC. With simultaneous observations of the flare, we built the broadband SED and studied it in the framework of a leptonic and an hadronic model. The VHE $\gamma$-ray observations show a flux level of 55\% of the Crab Nebula flux above 300\,GeV, similar to the historical maximum of 1995. The combination of MAGIC and \textit{Fermi}-LAT spectra provides an unprecedented characterization of the inverse-Compton peak for this object during a flaring episode. The $\Gamma$ index of the intrinsic spectrum in the VHE $\gamma$-ray band is $2.04\pm0.12_{\rm stat}\pm0.15_{\rm sys}$. We find the source in an extreme state with a shift of the position of the synchrotron peak to frequencies above or equal to $10^{18}$\,Hz

Context. It has become evident that one-zone synchrotron self-Compton models are not always adequate for very-high-energy (VHE) gamma-ray emitting blazars. While two-component models are performing better, they are difficult to constrain due to the large number of free parameters. Aims. In this work, we make a first attempt to take into account the observational constraints from Very Long Baseline Interferometry (VLBI) data, long-term light curves (radio, optical, and X-rays) and optical polarisation to limit the parameter space for a two-component model and test if it can still reproduce the observed spectral energy distribution (SED) of the blazars. Methods. We selected five TeV BL Lac objects based on the availability of VHE gamma-ray and optical polarisation data. We collected constraints for the jet parameters from VLBI observations. We evaluated the contributions of the two components to the optical flux by means of decomposition of long-term radio and optical light curves as well as modeling of the optical polarisation variability of the objects. We selected eight epochs for these five objects, based on the variability observed at VHE gamma rays, for which we constructed the SEDs that we then modeled with a two-component model. Results. We found parameter sets which can reproduce the broadband SED of the sources in the framework of two-component models considering all available observational constraints from VLBI observations. Moreover, the constraints obtained from the long-term behavior of the sources in the lower energy bands could be used to determine the region where the emission in each band originates. Finally, we attempted to use optical polarisation data to shed new light on the behavior of the two components in the optical band. Our observationally constrained two zone model allows explanation of the entire SED from radio to VHE with two co-located emission regions.

We provide a new fitting formula of the matter bispectrum in the nonlinear regime calibrated by high-resolution cosmological $N$-body simulations of $41$ cold dark matter ($w$CDM, $w=$ constant) models around the Planck 2015 best-fit parameters. As the parameterization in our fitting function is similar to that in Halofit, our fitting is named BiHalofit. The simulation volume is sufficiently large ($> 10 \, {\rm Gpc}^3$) to cover almost all measurable triangle bispectrum configurations in the universe. The function is also calibrated using one-loop perturbation theory at large scales ($k<0.3 \, h \, {\rm Mpc}^{-1}$). Our formula reproduced the matter bispectrum to within $10 \, (15) \, \%$ accuracy in the Planck 2015 model at wavenumber $k< 3 \, (10) \, h \, {\rm Mpc}^{-1}$ and redshifts $z=0-3$. The other $40$ $w$CDM models obtained poorer fits, with accuracy approximating $20 \, \%$ at $k<3 \, h \, {\rm Mpc}^{-1}$ and $z=0-1.5$ (the deviation includes the $10 \, \%$-level sample variance of the simulations). We also provide a fitting formula that corrects the baryonic effects such as radiative cooling and active galactic nucleus feedback, using the latest hydrodynamical simulation IllustrisTNG. We demonstrate that our new formula more accurately predicts the weak-lensing bispectrum than the existing fitting formulas. This formula will assist current and future weak-lensing surveys and cosmic microwave background lensing experiments. Numerical codes of the formula are available, written in Python, C and Fortran.

An accurate theoretical template of the redshift-space galaxy power spectrum, if applicable out to nonlinear scales, enables us to extract more stringent and robust constraints on cosmological parameters from the measured galaxy clustering. In this work we develop a simulation-based template, so-called emulator, for the redshift-space power spectrum of dark matter halos. Using the redshift-space halo power spectra measured from the Dark Quest $N$-body simulation suite that covers 101 flat-geometry $w$CDM cosmologies around the Planck $\Lambda$CDM model, we feed these data into a feed-forward neural network to build the fast and accurate emulation of the power spectrum from the linear to nonlinear scales up to $k \simeq 0.6 \, h \, {\rm Mpc}^{-1}$. Our emulator achieves about 1 and 5% fractional accuracies in predicting the monopole and quadrupole moments of the power spectrum, respectively, for halos of $\sim 10^{13}h^{-1}M_\odot$ that correspond to host halos of the SDSS LOWZ- and CMASS-like galaxies, where the achieved accuracies are sufficient compared to the statistical errors of SDSS volume. The validation and performance of the emulator are given by the comparison of the emulator predictions with the power spectra directly measured from the simulations for validation sets that are not used in the training. We demonstrate that the emulator outputs can be used to make model predictions for the redshift-space power spectrum of galaxies by employing user-fed models for the halo-galaxy connection, such as the halo occupation distribution. The emulator allows us to easily incorporate the Finger-of-God effect due to the virial motions of galaxies and the Alcock--Paczy\'{n}ski distortions. Our code can compute the redshift-space galaxy power spectrum in a CPU subseconds, and is ready to perform the emulator-based cosmological analysis for the exiting and upcoming galaxy redshift surveys.

We use a suite of $N$-body simulations to study intrinsic alignments (IA) of halo shapes with the surrounding large-scale structure in the $\Lambda$CDM model. For this purpose, we develop a novel method to measure multipole moments of the three-dimensional power spectrum of the $E$-mode field of halo shapes with the matter/halo distribution, $P_{\delta E}^{(\ell)}(k)$ (or $P^{(\ell)}_{{\rm h}E}$), and those of the auto-power spectrum of the $E$ mode, $P^{(\ell)}_{EE}(k)$, based on the $E$/$B$-mode decomposition. The IA power spectra have non-vanishing amplitudes over the linear to nonlinear scales, and the large-scale amplitudes at $k\lesssim 0.1~h~{\rm Mpc}^{-1}$ are related to the matter power spectrum via a constant coefficient ($A_{\rm IA}$), similar to the linear bias parameter. We find that the cross- and auto-power spectra $P_{\delta E}$ and $P_{EE}$ at nonlinear scales, $k\gtrsim 0.1~h~{\rm Mpc}^{-1}$, show different $k$-dependences relative to the matter power spectrum, suggesting a violation of the nonlinear alignment model commonly used to model contaminations of cosmic shear signals. The IA power spectra exhibit baryon acoustic oscillations, and vary with halo samples of different masses, redshifts and cosmological parameters ($\Omega_{\rm m}, S_8$). The cumulative signal-to-noise ratio for the IA power spectra is about 60\% of that for the halo density power spectrum, where the super-sample covariance is found to give a significant contribution to the total covariance. Our results demonstrate that the IA power spectra of galaxy shapes, measured from imaging and spectroscopic surveys for an overlapping area of the sky, can be powerful tools to probe the underlying matter power spectrum, the primordial curvature perturbations, and cosmological parameters, in addition to the standard galaxy density power spectrum.

We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer.

Recent developments of transition-edge sensors (TESs), based on extensive experience in ground-based experiments, have been making the sensor techniques mature enough for their application on future satellite CMB polarization experiments. LiteBIRD is in the most advanced phase among such future satellites, targeting its launch in Japanese Fiscal Year 2027 (2027FY) with JAXA's H3 rocket. It will accommodate more than 4000 TESs in focal planes of reflective low-frequency and refractive medium-and-high-frequency telescopes in order to detect a signature imprinted on the cosmic microwave background (CMB) by the primordial gravitational waves predicted in cosmic inflation. The total wide frequency coverage between 34GHz and 448GHz enables us to extract such weak spiral polarization patterns through the precise subtraction of our Galaxy's foreground emission by using spectral differences among CMB and foreground signals. Telescopes are cooled down to 5Kelvin for suppressing thermal noise and contain polarization modulators with transmissive half-wave plates at individual apertures for separating sky polarization signals from artificial polarization and for mitigating from instrumental 1/f noise. Passive cooling by using V-grooves supports active cooling with mechanical coolers as well as adiabatic demagnetization refrigerators. Sky observations from the second Sun-Earth Lagrangian point, L2, are planned for three years. An international collaboration between Japan, USA, Canada, and Europe is sharing various roles. In May 2019, the Institute of Space and Astronautical Science (ISAS), JAXA selected LiteBIRD as the strategic large mission No. 2.

Clustering properties and peculiar velocities of halos in large-scale structure carry a wealth of cosmological information over a wide range of scales from linear to nonlinear scales. We use halo catalogs in a suite of high-resolution $N$-body simulations to construct mock catalogs of galaxies that resemble the SDSS-like luminous early-type galaxies at three redshift bins in the range $0.15\le z\le 0.7$. To do this we include 10 nuisance parameters to model variations in halo-galaxy connections for each redshift bin; the halo occupation distribution, and the spatial and velocity distributions of galaxies in the host halos. We evaluate the Fisher information matrix for the redshift-space power spectrum of SDSS-like galaxies using different sets of the mock catalogs that are generated from changes in each of model parameters; cosmological parameters ($\sigma_8$ and $\Omega_{\rm m}$), the halo-galaxy connection parameters, and the cosmological distances ($D_{\rm A}$ and $H$ parameters at each redshift bin) for modeling an apparent geometrical distortion of the redshift-space power spectrum (the Alcock-Paczynski effect). We show that combining the monopole and quadrupole power spectra of galaxies allows for accurate estimations of the cosmological parameters and the cosmological distances, even after marginalization over the halo-galaxy parameters, by lifting the parameter degeneracies that are otherwise inevitable if either of the two spectra alone is used. When including the galaxy power spectrum information up to $k=0.3\,h\,$Mpc$^{-1}$, we find about factor of 6 gain in the cosmological information content of ($\sigma_8,\Omega_{\rm m},D_{\rm A}$'s and $H$'s) compared to $k=0.2\,h\,$Mpc$^{-1}$.

We present the results of a dust-reverberation survey of quasars at redshifts z<0.6. We found a delayed response of the K-band flux variation after the optical flux variation in 25 out of 31 targets, and obtained the lag time between them for 22 targets. Combined with the results for nearby Seyfert galaxies, we provide the largest homogeneous collection of K-band dust-reverberation data for 36 type 1 active galactic nuclei (AGNs). This doubles the sample and includes the most distant AGN and the largest lag so far measured. We estimated the optical luminosity of the AGN component of each target using three different methods: spectral decomposition, the flux-variation-gradient method, and image decomposition. We found a strong correlation between the reverberation radius for the innermost dust torus and the optical luminosity over a range of approximately four orders of magnitude in luminosity, as is already known for Seyfert galaxies. We estimated the luminosity distances of the AGNs based on their dust-reverberation lags, and found that the data in the redshift-distance diagram are consistent with the current standard estimates of the cosmological parameters. We also present the radius-luminosity relations for isotropic luminosity indicators such as the hard X-ray (14--195 keV), [OIV] 25.89 um, and mid-infrared (12 um) continuum luminosities, which are applicable to obscured AGNs.

We perform an ensemble of $N$-body simulations with $2048^3$ particles for 101 flat $w$CDM cosmological models sampled based on a maximin-distance Sliced Latin Hypercube Design. By using the halo catalogs extracted at multiple redshifts in the range of $z=[0,1.48]$, we develop Dark Emulator, which enables fast and accurate computations of the halo mass function, halo-matter cross-correlation, and halo auto-correlation as a function of halo masses, redshift, separations and cosmological models, based on the Principal Component Analysis and the Gaussian Process Regression for the large-dimensional input and output data vector. We assess the performance of the emulator using a validation set of $N$-body simulations that are not used in training the emulator. We show that, for typical halos hosting CMASS galaxies in the Sloan Digital Sky Survey, the emulator predicts the halo-matter cross correlation, relevant for galaxy-galaxy weak lensing, with an accuracy better than $2\%$ and the halo auto-correlation, relevant for galaxy clustering correlation, with an accuracy better than $4\%$. We give several demonstrations of the emulator. It can be used to study properties of halo mass density profiles such as the mass-concentration relation and splashback radius for different cosmologies. The emulator outputs can be combined with an analytical prescription of halo-galaxy connection such as the halo occupation distribution at the equation level, instead of using the mock catalogs, to make accurate predictions of galaxy clustering statistics such as the galaxy-galaxy weak lensing and the projected correlation function for any model within the $w$CDM cosmologies, in a few CPU seconds.

The Cherenkov Telescope Array (CTA) is the major next-generation observatory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possible by using tens of imaging Cherenkov telescopes of three successive sizes. They will be arranged into two arrays, one per hemisphere, located on the La Palma island (Spain) and in Paranal (Chile). We present here the optimised and final telescope arrays for both CTA sites, as well as their foreseen performance, resulting from the analysis of three different large-scale Monte Carlo productions.

The FOREST Unbiased Galactic plane Imaging survey with the Nobeyama 45-m telescope (FUGIN) project is one of the legacy projects using the new multi-beam FOREST receiver installed on the Nobeyama 45-m telescope. This project aims to investigate the distribution, kinematics, and physical properties of both diffuse and dense molecular gas in the Galaxy at once by observing 12CO, 13CO, and C18O J=1-0 lines simultaneously. The mapping regions are a part of the 1st quadrant (10d < l < 50d, |b| < 1d) and the 3rd quadrant (198d < l <236d, |b| < 1d) of the Galaxy, where spiral arms, bar structure, and the molecular gas ring are included. This survey achieves the highest angular resolution to date (~20") for the Galactic plane survey in the CO J=1-0 lines, which makes it possible to find dense clumps located farther away than the previous surveys. FUGIN will provide us with an invaluable dataset for investigating the physics of the galactic interstellar medium (ISM), particularly the evolution of interstellar gas covering galactic scale structures to the internal structures of giant molecular clouds, such as small filament/clump/core. We present an overview of the FUGIN project, observation plan, and initial results, which reveal wide-field and detailed structures of molecular clouds, such as entangled filaments that have not been obvious in previous surveys, and large-scale kinematics of molecular gas such as spiral arms.

In Yoshii et al. (2014), we described a new method for measuring extragalactic distances based on dust reverberation in active galactic nuclei (AGNs), and we validated our new method with Cepheid variable stars. In this paper, we validate our new method with Type Ia supernovae (SNe Ia) which occurred in two of the AGN host galaxies during our AGN monitoring program: SN 2004bd in NGC 3786 and SN 2008ec in NGC 7469. Their multicolor light curves were observed and analyzed using two widely accepted methods for measuring SN distances, and the distance moduli derived are $\mu=33.47\pm 0.15$ for SN 2004bd and $33.83\pm 0.07$ for SN 2008ec. These results are used to obtain independently the distance measurement calibration factor, $g$. The $g$ value obtained from the SN Ia discussed in this paper is $g_{\rm SN} = 10.61\pm 0.50$ which matches, within the range of 1$\sigma$ uncertainty, $g_{\rm DUST} = 10.60$, previously calculated ab initio in Yoshii et al. (2014). Having validated our new method for measuring extragalactic distances, we use our new method to calibrate reverberation distances derived from variations of H$\beta$ emission in the AGN broad line region (BLR), extending the Hubble diagram to $z\approx 0.3$ where distinguishing between cosmologies is becoming possible.

We present the results of a dust reverberation survey for 17 nearby Seyfert 1 galaxies, which provides the largest homogeneous data collection for the radius of the innermost dust torus. A delayed response of the K-band light curve after the V-band light curve was found for all targets, and 49 measurements of lag times between the flux variation of the dust emission in the K band and that of the optical continuum emission in the V band were obtained. The lag times strongly correlated with the optical luminosity in the luminosity range of M_V=-16 to -22 mag, and the regression analysis was performed to obtain the correlation log $\Delta t$ (days) = -2.11 -0.2 M_V assuming $\Delta t \propto L^{0.5}$, which was theoretically expected. We discuss the possible origins of the intrinsic scatter of the dust lag-luminosity correlation, which was estimated to be about 0.13 dex, and we find that the difference of internal extinction and delayed response of changes in lag times to the flux variations could have partly contributed to intrinsic scatter. However, we could not detect any systematic change of the correlation with the subclass of the Seyfert type or the Eddington ratio. Finally, we compare the dust reverberation radius with the near-infrared interferometric radius of the dust torus and the reverberation radius of broad Balmer emission lines. The interferometric radius in the K band was found to be systematically larger than the dust reverberation radius in the same band by about a factor of two, which could be interpreted by the difference between the flux-weighted radius and the response-weighted radius of the innermost dust torus. The reverberation radius of the broad Balmer emission lines was found to be systematically smaller than the dust reverberation radius by about a factor of 4-5, which strongly supports the unified scheme of the Seyfert type of active galactic nuclei. (Abridged)

We have pioneered a new method for the measurement of extragalactic distances. This method uses the time-lag between variations in the short wavelength and long wavelength light from an active galactic nucleus (AGN), based on a quantitative physical model of dust reverberation that relates the time-lag to the absolute luminosity of the AGN. We use the large homogeneous data set from intensive monitoring observations in optical and near-infrared wavelength bands with the dedicated 2-m MAGNUM telescope to obtain the distances to 17 AGNs in the redshift range z=0.0024 to z=0.0353. These distance measurements are compared with distances measured using Cepheid variable stars, and are used to infer that H_0= 73 +- 3 (random) km/s/Mpc. The systematic error in H_0 is examined, and the uncertainty in the size distribution of dust grains is the largest source of the systematic error, which is much reduced for a sample of AGNs for which their parameter values in the model of dust reverberation are individually measured. This AGN time-lag method can be used beyond 30 Mpc, the farthest distance reached by extragalactic Cepheids, and can be extended to high-redshift quasi-stellar objects.

We present new data for five under-luminous type II-plateau supernovae (SNe IIP), namely SN 1999gn, SN 2002gd, SN 2003Z, SN 2004eg and SN 2006ov. This new sample of low-luminosity SNe IIP (LL SNe IIP) is analyzed together with similar objects studied in the past. All of them show a flat light curve plateau lasting about 100 days, an under luminous late-time exponential tail, intrinsic colours that are unusually red, and spectra showing prominent and narrow P-Cygni lines. A velocity of the ejected material below 10^3 km/s is inferred from measurements at the end of the plateau. The 56Ni masses ejected in the explosion are very small (less than 10^-2 solar masses). We investigate the correlations among 56Ni mass, expansion velocity of the ejecta and absolute magnitude in the middle of the plateau, confirming the main findings of Hamuy (2003), according to which events showing brighter plateau and larger expansion velocities are expected to produce more 56Ni. We propose that these faint objects represent the low luminosity tail of a continuous distribution in parameters space of SNe IIP. The physical properties of the progenitors at the explosion are estimated through the hydrodynamical modeling of the observables for two representative events of this class, namely SN 2005cs and SN 2008in. We find that the majority of LL SNe IIP, and quite possibly all, originate in the core-collapse of intermediate mass stars, in the mass range 10-15 solar masses.

We search for ultra-high energy photons by analyzing geometrical properties of shower fronts of events registered by the Telescope Array surface detector. By making use of an event-by-event statistical method, we derive upper limits on the absolute flux of primary photons with energies above 10^19, 10^19.5 and 10^20 eV based on the first three years of data taken.

We search for correlations between positions of extragalactic objects and arrival directions of Ultra-High Energy Cosmic Rays (UHECRs) with primary energy $E \ge 40$ EeV as observed by the surface detector array of the Telescope Array (TA) experiment during the first 40 months of operation. We examined several public astronomical object catalogs, including the Veron-Cetty and Veron catalog of active galactic nuclei. We counted the number of TA events correlated with objects in each catalog as a function of three parameters: the maximum angular separation between a TA event and an object, the minimum energy of the events, and the maximum redshift of the objects. We determine combinations of these parameters which maximize the correlations, and calculate the chance probabilities of having the same levels of correlations from an isotropic distribution of UHECR arrival directions. No statistically significant correlations are found when penalties for scanning over the above parameters and for searching in several catalogs are taken into account.

We measure the spectrum of cosmic rays with energies greater than $10^{18.2}$ eV with the Fluorescence Detectors (FDs) and the Surface Detectors (SDs) of the Telescope Array Experiment using the data taken in our first 2.3-year observation from May 27 2008 to September 7 2010. A hybrid air shower reconstruction technique is employed to improve accuracies in determination of arrival directions and primary energies of cosmic rays using both FD and SD data. The energy spectrum presented here is in agreement with our previously published spectra and the HiRes results.

We present a measurement of the energy spectrum of ultra-high-energy cosmic rays performed by the Telescope Array experiment using monocular observations from its two new FADC-based fluorescence detectors. After a short description of the experiment, we describe the data analysis and event reconstruction procedures. Since the aperture of the experiment must be calculated by Monte Carlo simulation, we describe this calculation and the comparisons of simulated and real data used to verify the validity of the aperture calculation. Finally, we present the energy spectrum calculated from the merged monocular data sets of the two FADC-based detectors, and also the combination of this merged spectrum with an independent, previously published monocular spectrum measurement performed by Telescope Array's third fluorescence detector (Abu-Zayyad {\it et al.}, {Astropart. Phys.} 39 (2012), 109). This combined spectrum corroborates the recently published Telescope Array surface detector spectrum (Abu-Zayyad {\it et al.}, {Astrophys. Journ.} 768 (2013), L1) with independent systematic uncertainties.

The Telescope Array (TA) collaboration has measured the energy spectrum of ultra-high energy cosmic rays with primary energies above 1.6 x 10^(18) eV. This measurement is based upon four years of observation by the surface detector component of TA. The spectrum shows a dip at an energy of 4.6 x 10^(18) eV and a steepening at 5.4 x 10^(19) eV which is consistent with the expectation from the GZK cutoff. We present the results of a technique, new to the analysis of ultra-high energy cosmic ray surface detector data, that involves generating a complete simulation of ultra-high energy cosmic rays striking the TA surface detector. The procedure starts with shower simulations using the CORSIKA Monte Carlo program where we have solved the problems caused by use of the "thinning" approximation. This simulation method allows us to make an accurate calculation of the acceptance of the detector for the energies concerned.

We study the anisotropy of Ultra-High Energy Cosmic Ray (UHECR) events collected by the Telescope Array (TA) detector in the first 40 months of operation. Following earlier studies, we examine event sets with energy thresholds of 10 EeV, 40 EeV, and 57 EeV. We find that the distributions of the events in right ascension and declination are compatible with an isotropic distribution in all three sets. We then compare with previously reported clustering of the UHECR events at small angular scales. No significant clustering is found in the TA data. We then check the events with E>57 EeV for correlations with nearby active galactic nuclei. No significant correlation is found. Finally, we examine all three sets for correlations with the large-scale structure of the Universe. We find that the two higher-energy sets are compatible with both an isotropic distribution and the hypothesis that UHECR sources follow the matter distribution of the Universe (the LSS hypothesis), while the event set with E>10 EeV is compatible with isotropy and is not compatible with the LSS hypothesis at 95% CL unless large deflection angles are also assumed. We show that accounting for UHECR deflections in a realistic model of the Galactic magnetic field can make this set compatible with the LSS hypothesis.

The Telescope Array (TA) experiment, located in the western desert of Utah,USA, is designed for observation of extensive air showers from extremely high energy cosmic rays. The experiment has a surface detector array surrounded by three fluorescence detectors to enable simultaneous detection of shower particles at ground level and fluorescence photons along the shower track. The TA surface detectors and fluorescence detectors started full hybrid observation in March, 2008. In this article we describe the design and technical features of the TA surface detector.

The Telescope Array's Middle Drum fluorescence detector was instrumented with telescopes refurbished from the High Resolution Fly's Eye's HiRes-1 site. The data observed by Middle Drum in monocular mode was analyzed via the HiRes-1 profile-constrained geometry reconstruction technique and utilized the same calibration techniques enabling a direct comparison of the energy spectra and energy scales between the two experiments. The spectrum measured using the Middle Drum telescopes is based on a three-year exposure collected between December 16, 2007 and December 16, 2010. The calculated difference between the spectrum of the Middle Drum observations and the published spectrum obtained by the data collected by the HiRes-1 site allows the HiRes-1 energy scale to be transferred to Middle Drum. The HiRes energy scale is applied to the entire Telescope Array by making a comparison between Middle Drum monocular events and hybrid events that triggered both Middle Drum and the Telescope Array's scintillator Ground Array.

Since 2007, the Telescope Array (TA) experiment, based in Utah, USA, has been observing ultra high energy cosmic rays to understand their origins. The experiment involves a surface detector (SD) array and three fluorescence detector (FD) stations. FD stations, installed surrounding the SD array, measure the air fluorescence light emitted from extensive air showers (EASs) for precise determination of their energies and species. The detectors employed at one of the three FD stations were relocated from the High Resolution Fly's Eye experiment. At the other two stations, newly designed detectors were constructed for the TA experiment. An FD consists of a primary mirror and a camera equipped with photomultiplier tubes. To obtain the EAS parameters with high accuracies, understanding the FD optical characteristics is important. In this paper, we report the characteristics and installation of new FDs and the performances of the FD components. The results of the monitored mirror reflectance during the observation time are also described in this report.

We examine whether the spectral energy distribution of UV continuum emission of active galactic nuclei changes during flux variation. We used multi-epoch photometric data of QSOs in the Stripe 82 observed by the SDSS Legacy Survey and selected 10 bright QSOs observed with high photometric accuracies, in the redshift range of z = 1.0-2.4 where strong broad emission lines such as Ly\alpha and CIV do not contaminate SDSS filters, to examine spectral variation of the UV continuum emission with broadband photometries. All target QSOs showed clear flux variations during the monitoring period 1998-2007, and the multi-epoch flux data in two different bands obtained on the same night showed a linear flux-to-flux relationship for all target QSOs. Assigning the flux in the longer wavelength to the x-axis in the flux-to-flux diagram, the x-intercept of the best-fit linear regression line was positive for most targets, which means that their colors in the observing bands become bluer as they become brighter. Then, the host-galaxy flux was estimated on the basis of the correlation between the stellar mass of the bulge of the host galaxy and the central black hole mass. We found that the longer-wavelength flux of the host galaxy was systematically smaller than that of the fainter extension of the best-fit regression line at the same shorter-wavelength flux for most targets. This result strongly indicates that the spectral shape of the continuum emission of QSOs in the UV region usually becomes bluer as it becomes brighter. We found that the multi-epoch flux-to-flux plots could be fitted well with the standard accretion disk model changing the mass accretion rate with a constant black hole mass for most targets. This finding strongly supports the standard accretion disk model for UV continuum emission of QSOs.

We examine whether the spectral energy distribution of optical continuum emission of active galactic nuclei (AGNs) changes during flux variation, based on accurate and frequent monitoring observations of 11 nearby Seyfert galaxies and QSOs carried out in the B, V, and I bands for seven years by the MAGNUM telescope. The multi-epoch flux data in any two different bands obtained on the same night show a very tight linear flux to flux relationship for all target AGNs. The flux of the host galaxy within the photometric aperture is carefully estimated by surface brightness fitting to available high-resolution HST images and MAGNUM images. The flux of narrow emission lines in the photometric bands is also estimated from available spectroscopic data. We find that the non-variable component of the host galaxy plus narrow emission lines for all target AGNs is located on the fainter extension of the linear regression line of multi-epoch flux data in the flux to flux diagram. This result strongly indicates that the spectral shape of AGN continuum emission in the optical region does not systematically change during flux variation. The trend of spectral hardening that optical continuum emission becomes bluer as it becomes brighter, which has been reported by many studies, is therefore interpreted as the domination of the variable component of the nearly constant spectral shape of an AGN as it brightens over the non-variable component of the host galaxy plus narrow lines, which is usually redder than AGN continuum emission.

The long-term optical and near infrared monitoring observations for a type 1 act ive galactic nucleus NGC 4151 were carried out for six years from 2001 to 2006 b y using the MAGNUM telescope, and delayed response of flux variations in the $K(2.2\mu m)$ band to those in the $V(0.55\mu m)$ band was clearly detected. Based on cross correlation analysis, we precisely measured a lag time $\Delta t$ for eight separate periods, and we found that $\Delta t$ is not constant changing be tween 30 and 70 days during the monitoring period. Since $\Delta t$ is the ligh t travel time from the central energy source out to the surrounding dust torus, this is the first convincing evidence that the inner radius of dust torus did ch ange in an individual AGN. In order to relate such a change of $\Delta t$ with a change of AGN luminosity $L$, we presented a method of taking an average of th e observed $V$-band fluxes that corresponds to the measured value of $\Delta t$, and we found that the time-changing track of NGC 4151 in the $\Delta t$ versus $L$ diagram during the monitoring period deviates from the relation of $\Delta t \propto L^{0.5}$ expected from dust reverberation. This result, combined with t he elapsed time from period to period for which $\Delta t$ was measured, indicat es that the timescale of dust formation is about one year, which should be taken into account as a new constraint in future studies of dust evolution in AGNs.

We searched for evidence of line emission around 4keV from the northwestern rim of the supernova remnant RX J0852.0-4622 using Suzaku XIS data. Several papers have reported the detection of an emission line around 4.1keV from this region of the sky. This line would arise from K-band fluorescence by Sc, the immediate decay product of 44Ti. We performed spectral analysis for the entire portion of the NW rim of the remnant within the XIS field of view, as well as various regions corresponding to regions of published claims of line emission. We found no line emission around 4.1keV anywhere, and are able to set a restrictive upper limit to the X-ray flux: 1.1x10^-6 s^-1 cm^-2 for the entire field. For every region, our flux upper limit falls below that of the previously claimed detection. Therefore, we conclude that, to date, no definite X-ray line feature from Sc-K emission has been detected in the NW rim of RX J0852.0-4622. Our negative-detection supports the recent claim that RX J0852-4622 is neither young (1700--4000 yr) nor nearby(~750 pc).

We present the most successful infrared cloud monitor for a robotic telescope. This system was originally developed for the MAGNUM 2-m telescope, which has been achieving unmanned and automated monitoring observation of active galactic nuclei at Haleakala on the Hawaiian island of Maui since 2001. Using a thermal imager and two aspherical mirrors, it at once sees almost the whole sky at a wavelength of $\lambda\sim 10\mu{\rm m}$. Its outdoor part is weather-proof and is totally maintenance-free. The images obtained every one or two minutes are analysed immediately into several ranks of weather condition, from which our automated observing system not only decides to open or close the dome, but also selects what types of observations should be done. The whole-sky data accumulated over four years show that 50$-$60 % of all nights are photometric, and about 75 % are observable with respect to cloud condition at Haleakala. Many copies of this system are now used all over the world such as Mauna Kea in Hawaii, Atacama in Chile, and Okayama and Kiso in Japan.

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

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.

We propose a new method of analysing a variable component for type 1 active galactic nuclei (AGNs) in the near-infrared wavelength region. This analysis uses a multiple regression technique and divides the variable component into two components originating in the accretion disk at the center of AGNs and from the dust torus that far surrounds the disk. Applying this analysis to the long-term $VHK$ monitoring data of MCG+08-11-011 that were obtained by the MAGNUM project, we found that the $(H-K)$-color temperature of the dust component is $T = 1635$K $\pm20$K, which agrees with the sublimation temperature of dust grains, and that the time delay of $K$ to $H$ variations is $\Delta t\approx 6$ days, which indicates the existence of a radial temperature gradient in the dust torus. As for the disk component, we found that the power-law spectrum of $f_\nu \propto \nu^\alpha$ in the $V$ to near-infrared $HK$ bands varies with a fixed index of $\alpha\approx -0.1$ -- +0.4, which is broadly consistent with the irradiated standard disk model. The outer part of the disk therefore extends out to a radial distance where the temperature decreases to radiate the light in the near-infrared.

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

We carried out a one-night optical V and near-infrared JHK monitoring observation of the least luminous Seyfert 1 galaxy, NGC4395, on 2004 May 1, and detected for the first time the intraday flux variations in the J and H bands, while such variation was not clearly seen for the K band. The detected J and H variations are synchronized with the flux variation in the V band, which indicates that the intraday-variable component of near-infrared continuum emission of the NGC4395 nucleus is an extension of power-law continuum emission to the near-infrared and originates in an outer region of the central accretion disk. On the other hand, from our regular program of long-term optical BVI and near-infrared JHK monitoring observation of NGC4395 from 2004 February 12 until 2005 January 22, we found large flux variations in all the bands on time scales of days to months. The optical BVI variations are almost synchronized with each other, but not completely with the near-infrared JHK variations. The color temperature of the near-infrared variable component is estimated to be T=1320-1710 K, in agreement with thermal emission from hot dust tori in active galactic nuclei (AGNs). We therefore conclude that the near-infrared variation consists of two components having different time scales, so that a small K-flux variation on a time scale of a few hours would possibly be veiled by large variation of thermal dust emission on a time scale of days.

Late-time BVRIJHK photometry of the peculiar Type Ic SN 2002ap, taken between 2002 June 12 and 2003 August 29 with the MAGNUM telescope, is presented. The light curve decline rate is derived in each band and the color evolution is studied through comparison with nebular spectra and with SN 1998bw. Using the photometry, the OIR bolometric light curve is built, extending from before light maximum to day 580 after explosion. The light curve has a late-time shape strikingly similar to that of the hypernova SN 1998bw. The decline rate changes from 0.018 mag/day between day 130 and 230 to 0.014 mag/day between day 270 and 580. To reproduce the late-time light curve, a dense core must be added to the 1-D hypernova model that best fits the early-time observations, bringing the ejecta mass from 2.5 Msun to 3 Msun without much change in the kinetic energy, which is 4 times 10^51 ergs. This is similar to the case of other hypernovae and suggests asymmetry. A large H-band bump developed in the spectral energy distribution after about day 300, probably caused by strong [Si I] 1.646 micron and 1.608 micron emissions. The near-infrared flux contribution increased simultaneously from <30% to >50% at day 580. The near-infrared light curves were compared with those of other Type Ib/c supernovae, among which SN 1983I seems similar to SN 2002ap both in the near-infrared and in the optical.

The most intense monitoring observations yet made in the optical (UBV) and near-infrared (JHK) wave bands were carried out for nearby Seyfert1 galaxies of NGC 5548, NGC 4051, NGC 3227, and NGC 7469. Over three years of observations with MAGNUM telescope since early 2001, clear time-delayed response of the K-band flux variations to the V-band flux variations was detected for all of these galaxies. Their H-K color temperature was estimated to be 1500-1800 K from the observed flux variation gradients, which supports a view that the bulk of the K flux should originate in the thermal radiation of hot dust that surrounds the central engine. Cross-correlation analysis was performed to quantify the lag time corresponding to the light-travel distance of the hot dust region from the central engine. The measured lag time is 47-53 days for NGC 5548, 11-18 days for NGC 4051, about 20 days for NGC 3227, and 65-87 days for NGC 7469. We found that the lag time is tightly correlated with the optical luminosity as expected from dust reverberation ($\Delta t \propto L^{0.5}$), while only weakly with the central virial mass, which suggests that an inner radius of the dust torus around the active nucleus has a one-to-one correspondence to central luminosity. In the lag time versus central luminosity diagram, the K-band lag times place an upper boundary on the similar lag times of broad-emission lines in the literature. This not only supports the unified scheme of AGNs, but also implies a physical transition from the BLR out to the dust torus that encircles the BLR. Furthermore, our V-band flux variations of NGC 5548 on timescales of up to 10 days are found to correlate with X-ray variations and delay behind them by one or two days, indicating the thermal reprocessing of X-ray emission by the central accretion flow.

The most intense monitoring observations yet made were carried out on the Seyfert 1 galaxy NGC 4151 in the optical and near-infrared wave-bands. A lag from the optical light curve to the near-infrared light curve was measured. The lag-time between the V and K light curves at the flux minimum in 2001 was precisely 48+2-3 days, as determined by a cross-correlation analysis. The correlation between the optical luminosity of an active galactic nucleus (AGN) and the lag-time between the UV/optical and the near-infrared light curves is presented for NGC 4151 in combination with previous lag-time measurements of NGC 4151 and other AGNs in the literature. This correlation is interpreted as thermal dust reverberation in an AGN, where the near-infrared emission from an AGN is expected to be the thermal re-radiation from hot dust surrounding the central engine at a radius where the temperature equals to that of the dust sublimation temperature. We find that the inner radius of the dust torus in NGC 4151 is $\sim $ 0.04 pc corresponding to the measured lag-time, well outside the broad line region (BLR) determined by other reverberation studies of the emission lines.

The light curves of 'hypernovae', i.e. very energetic supernovae with $E_{51} \equiv E/10^{51}$ergs $\gsim 5-10$ are characterized at epochs of a few months by a phase of linear decline. Classical, one-dimensional explosion models fail to simultaneously reproduce the light curve near peak and at the linear decline phase. The evolution of these light curves may however be explained by a simple model consisting of two concentric components. The outer component is responsible for the early part of the light curve and for the broad absorption features observed in the early spectra of hypernovae, similar to the one-dimensional models. In addition, a very dense inner component is added, which reproduces the linear decline phase in the observed magnitude-versus-time relation for SNe 1998bw, 1997ef, and 2002ap. This simple approach does contain one of the main features of jet-driven, asymmetric explosion models, namely the presence of a dense core. Although the total masses and energies derived with the two-component model are similar to those obtained in previous studies which also adopted spherical symmetry, this study suggests that the ejecta are aspherical, and thus the real energies and masses may deviate from those derived assuming spherical symmetry. The supernovae which were modeled are divided into two groups, according to the prominence of the inner component: the inner component of SN 1997ef is denser and more $^{56}$Ni-rich, relative to the outer component, than the corresponding inner components of SNe 1998bw and 2002ap. These latter objects have a similar inner-to-outer component ratio, although they have very different global values of mass and energy.

Supernova (SN) 2002ap in M74 was observed in the $UBVRIJHK$ bands for the first 40 days following its discovery (2002 January 29) until it disappeared because of solar conjunction, and then in June after it reappeared. The magnitudes and dates of peak brightness in each band were determined. While the rate of increase of the brightness before the peak is almost independent of wavelength, the subsequent rate of decrease becomes smaller with wavelength from the $U$ to the $R$ band, and is constant at wavelengths beyond $I$. The photometric evolution is faster than in the well-known ``hypernovae'' SNe~1998bw and 1997ef, indicating that SN 2002ap ejected less mass. The bolometric light curve of SN 2002ap for the full period of observations was constructed. The absolute magnitude is found to be much fainter than that of SN 1998bw, but is similar to that of SN 1997ef, which lies at the faint end of the hypernova population. The bolometric light curve at the early epochs was best reproduced with the explosion of a C+O star that ejects $2.5~M_\sun$ with kinetic energy $E_{\rm K}=4\times 10^{51}~{\rm ergs}$. A comparison of the predicted brightness of SN 2002ap with that observed after solar conjunction may imply that $\gamma$-ray deposition at the later epochs was more efficient than in the model. This may be due to an asymmetric explosion.

Photometric and spectroscopic data of the energetic Type Ic supernova (SN) 2002ap are presented, and the properties of the SN are investigated through models of its spectral evolution and its light curve. The SN is spectroscopically similar to the "hypernova" SN 1997ef. However, its kinetic energy [$\sim (4-10) \times 10^{51}$ erg] and the mass ejected (2.5-5 $M_{\odot}$) are smaller, resulting in a faster-evolving light curve. The SN synthesized $\sim 0.07 M_{\odot}$ of $^{56}$Ni, and its peak luminosity was similar to that of normal SNe. Brightness alone should not be used to define a hypernova, whose defining character, namely very broad spectral features, is the result of a high kinetic energy. The likely main-sequence mass of the progenitor star was 20-25 $M_{\odot}$, which is also lower than that of both hypernovae SNe 1997ef and 1998bw. SN 2002ap appears to lie at the low-energy and low-mass end of the hypernova sequence as it is known so far. Observations of the nebular spectrum, which is expected to dominate by summer 2002, are necessary to confirm these values.

Observations of 226 AGNs in the near-infrared $J$, $H$, and $K'$ bands are presented along with the analysis of the observations for variability. Our sample consists mainly of Seyfert 1 AGNs and QSOs. About a quarter of the objects in each category are radio loud. The AGNs in the entire sample have the redshifts spanning the range from $z=0$ to 1, and the absolute magnitudes from $M_B=-29$ to -18. All the objects were observed twice and their variability was measured by differential photometry. A reduction method of differential photometry, optimized to the analysis of extended images, has been developed. The systematic error in variability arising from AGNs of highly extended images is estimated to be less than 0.01 mag in each of the $J$, $H$, and $K'$ bands. The systematic error arising from the flat fielding is negligible for most AGNs, although it is more than 0.1 mag for some particular cases. The overall average flat fielding error is 0.03 mag for the image pairs. We find that these systematic errors are superseded by statistical errors, and the overall average total systematic and statistical errors amounts to 0.05 mag in the measured variability in each band. We find that 58% of all the AGNs in the entire sample show variability of more than $2\sigma$, and 44% of more than $3\sigma$. This result holds independent of the $J$, $H$, and $K'$ bands. The detection rate of variability is higher for a subsample of higher photometric accuracy, and there appears no limit to this tendency. In particular, when we consider a subsample with small photometric errors of $\sigma<0.03$ mag, the rate of $2\sigma$ detection is 80%, and 64% for $3\sigma$ detection. This suggests that most AGNs are variable in the near-infrared.

Variability of 226 AGNs in the near-infrared $J$, $H$, and $K'$ bands is analyzed and discussed. An ensemble average for measured variabilities was obtained for various samples of the AGNs divided by absolute $B$-magnitude $M_B$, redshift $z$, and radio strength. All the samples in the $J$, $H$, and $K'$ bands are found to give significant ensemble variability, but no significant wavelength dependence is found. The ensemble variability in the entire sample combining the $J$, $H$, and $K'$ samples is $\Delta m\approx 0.22$ mag, while $\Delta m\approx 0.18$ mag for the radio-quiet AGNs and $\Delta m\approx 0.26$ mag for radio-loud AGNs. The ensemble variability for the radio-quiet AGNs shows no significant $M_B$-dependence, while showing positive $M_B$-dependence for the radio-loud AGNs. In any samples the measured variability shows positive correlation among different passbands, with the correlation coefficients of $r_{JH}$, $r_{HK'}$, and $r_{JK'}$ ranging from 0.6 to 0.9. For radio-quiet AGNs, the coefficient $r_{HK'}$ in a redshift range of $0.1<z<0.3$ is significantly higher than $r_{JH}$ or $r_{JK'}$. The coefficient for the radio-loud AGNs with $0.6<z<1.0$ is as high as 0.95, irrespective of the passband. However, for the radio-quiet AGNs with $z>0.3$ and radio-loud AGNs with $z<0.3$, we cannot confirm such strong correlation among different passbands. All the features of near-infrared variability for the radio-quiet AGNs are consistent with a simple dust reverberation model of the central regions of AGNs. However, the features for the radio-loud AGNs are not fully explained by such a model, and a non-thermal variable component is suggested as a viable candidate for causing their large and fast variability in the near-infrared region.

Near-infrared $JHK'$ imaging photometry was obtained of 331 AGNs consisting mainly of Seyfert 1 AGNs and quasars (QSOs). This sample was selected to cover a range of radio emission strength, redshift from $z=0$ to 1, and absolute $B$-magnitude from $M_B=-29$ mag to -18 mag. Among low-$z$ AGNs with $z<0.3$, Seyfert $1-1.5$ AGNs are distributed over a region from a location typical of ``galaxies'' to a location typical of ``QSOs'' in the two-color $J-H$ to $H-K'$ diagram, but Seyfert $1.8-2$ AGNs are distributed around the location of ``galaxies''. Moreover, bright AGNs with respect to absolute $B$-magnitude are distributed near the location of ``QSOs'', while faint AGNs are near the location of ``galaxies''. The distribution of such low-$z$ AGNs in this diagram was found to have little dependence on their 6 cm radio flux. The near-infrared colors of the AGNs observed with an aperture of 7 pixels ($7.49''$) are more QSO-like than those observed with larger apertures up to 15 pixels ($16.1''$). This aperture effect may be explained by contamination from the light of host galaxies within larger apertures. This effect is more prominent for less luminous AGNs.

We present new deep $K'$ counts of stars at the South Galactic Pole (SGP) taken with the NAOJ PICNIC camera to $K'=17.25$. Star-galaxy separation to $K'=17.5$ was accomplished effectively using image profiles because the pixel size we used is 0.509 arcsec. We interpret these counts using the SKY (Cohen 1994) model of the Galactic point source sky and determine the relative normalization of halo-to-disk populations, and the location of the Sun relative to the Galactic plane, within the context of this model. The observed star counts constrain these parameters to be: halo/disk $\sim$ 1/900 and z$_\odot$=16.5$\pm$2.5 pc. These values have been used to correct our SGP galaxy counts for contamination by the point source Galactic foreground.

We present new K-band galaxy number counts from K=13 to 20.5 obtained from $K'$-band surveys in the south galactic pole region, which cover 180.8 arcmin$^2$ to a limiting magnitude of K=19, and 2.21 arcmin$^2$ to K=21. These are currently the most precise K-band galaxy counts at $17.5<K<19.0$ because the area of coverage is largest among the existing surveys for this magnitude range. The completeness and photometry corrections are estimated from the recovery of simulated galaxy and stellar profiles added to the obtained field image. Many simulations were carried out to construct a probability matrix which corrects the galaxy counts at the faint-end magnitudes of the surveys so the corrected counts can be compared with other observations. The K-band star counts in the south galactic pole region to $K=17.25$ are also presented for use to constrain the vertical structure of the Galaxy.

We report the discovery of a group of galaxies at redshift 2.38. We imaged about 10% of a claimed supercluster of QSO absorption-lines at z=2.38 (Francis & Hewett 1993). In this small field (2 arcmin radius) we detect two Ly-alpha emitting galaxies. The discovery of two such galaxies in our tiny field supports Francis & Hewett's interpretation of the absorption-line supercluster as a high redshift "Great Wall". One of the Ly-alpha galaxies lies 22 arcsec from a background QSO, and may be associated with a multi-component Ly-alpha absorption complex seen in the QSO spectrum. This galaxy has an extended (50kpc) lumpy Ly-alpha morphology, surrounding a compact IR-bright nucleus. The nucleus shows a pronounced break in its optical-UV colors at about 4000 A (rest-frame), consistent with a stellar population of mass about 7E11 solar masses, an age of more than 500 Myr, and little on-going star-formation. C IV emission is detected, suggesting that a concealed AGN is present. Extended H-alpha emission is also detected; the ratio of Ly-alpha flux to H-alpha is abnormally low (about 0.7), probable evidence for extended dust. This galaxy is surrounded by a number of very red (B-K>5) objects, some of which have colors suggesting that they too are at z=2.38. We hypothesize that this galaxy, its neighbors and a surrounding lumpy gas cloud may be a giant elliptical galaxy in the act of bottom-up formation.