The unique combination of superb spatial resolution, wide-area coverage, and deep depth of the optical imaging from the Hyper Suprime-Cam (HSC) Subaru Strategic Program is utilized to search for dual quasar candidates. Using an automated image analysis routine on 34,476 known SDSS quasars, we identify those with two (or more) distinct optical point sources in HSC images covering 796 deg^2. We find 421 candidates out to a redshift of 4.5 of which one hundred or so are more likely after filtering out contaminating stars. Angular separations of 0.6 - 4.0" correspond to projected separations of 3 - 30 kpc, a range relatively unexplored for population studies of luminous dual quasars. Using Keck-I/LRIS and Gemini-N/NIFS, we spectroscopically confirm three dual quasar systems at z < 1, two of which are previously unknown out of eight observed, based on the presence of characteristic broad emission lines in each component, while highlighting that the continuum of one object in one of the pairs is reddened. In all cases, the [OIII]5007 emission lines have mild velocity offsets, thus the joint [OIII] line profile is not double-peaked. We find a dual quasar fraction of 0.26+/-0.18% and no evidence for evolution. A comparison with the Horizon-AGN simulation seems to support the case of no evolution in the dual quasar fraction when broadly matching the quasar selection. These results may indicate a scenario in which the frequency of the simultaneous triggering of luminous quasars is not as sensitive as expected to the cosmic evolution of the merger rate or gas content of galaxies.

We report a discovery of $z = 4.72$ radio galaxy, HSC J083913.17+011308.1, by using the Lyman break technique with the Hyper Suprime-Cam Subaru Strategic Survey (HSC-SSP) catalog for VLA FIRST radio sources. The number of known high-$z$ radio galaxies (HzRGs) at $z > 3$ is quite small to constrain the evolution of HzRGs so far. The deep and wide-area optical survey by HSC-SSP enables us to apply the Lyman break technique to a large search for HzRGs. For an HzRG candidate among pre-selected $r$-band dropouts with a radio detection, a follow-up optical spectroscopy with GMOS/Gemini has been performed. The obtained spectrum presents a clear Ly$\alpha$ emission line redshifted to $z=4.72$. The SED fitting analysis with the rest-frame UV and optical photometries suggests the massive nature of this HzRG with $\log{M_*/M_{\odot}} = 11.4$. The small equivalent width of Ly$\alpha$ and the moderately red UV colors indicate its dusty host galaxy, implying a chemically evolved and dusty system. The radio spectral index does not meet a criterion for an ultra-steep spectrum: $\alpha^{325}_{1400}$ of $-1.1$ and $\alpha^{150}_{1400}$ of $-0.9$, demonstrating that the HSC-SSP survey compensates for a sub-population of HzRGs which are missed in surveys focusing on an ultra-steep spectral index.

We present our new ALMA observations targeting [OIII]88um, [CII]158um, [NII]122um, and dust continuum emission for three Lyman break galaxies at z=6.0293-6.2037 identified in the Subaru/Hyper Suprime-Cam survey. We clearly detect [OIII] and [CII] lines from all of the galaxies at 4.3-11.8sigma levels, and identify multi-band dust continuum emission in two of the three galaxies, allowing us to estimate infrared luminosities and dust temperatures simultaneously. In conjunction with previous ALMA observations for six galaxies at z>6, we confirm that all the nine z=6-9 galaxies have high [OIII]/[CII] ratios of L[OIII]/L[CII]~3-20, ~10 times higher than z~0 galaxies. We also find a positive correlation between the [OIII]/[CII] ratio and the Lya equivalent width (EW) at the ~90% confidence level. We carefully investigate physical origins of the high [OIII]/[CII] ratios at z=6-9 using Cloudy, and find that high density of the interstellar medium, low C/O abundance ratio, and the cosmic microwave background attenuation are responsible to only a part of the z=6-9 galaxies. Instead, the observed high [OIII]/[CII] ratios are explained by 10-100 times higher ionization parameters or low photodissociation region (PDR) covering fractions of 0-10%, both of which are consistent with our [NII] observations. The latter scenario can be reproduced with a density bounded nebula with PDR deficit, which would enhance the Lya, Lyman continuum, and C+ ionizing photons escape from galaxies, consistent with the [OIII]/[CII]-Lya EW correlation we find.

We examine the absorption line spectra of a sample of 31 luminous (M_UV~-23) Lyman break galaxies at redshift z~6 using data taken with the FOCAS and OSIRIS spectrographs on the Subaru and GTC telescopes. For two of these sources we present longer exposure data taken at higher spectral resolution from ESO's X-shooter spectrograph. Using these data, we demonstrate the practicality of stacking our lower resolution data to measure the depth of various interstellar and stellar absorption lines to probe the covering fraction of low ionization gas and the gas-phase and stellar metallicities near the end of the era of cosmic reionization. From maximum absorption line depths of SiII1260 and CII1334, we infer a mean covering fraction of >0.85 +/-0.16 for our sample. This is larger than that determined using similar methods for lower luminosity galaxies at slightly lower redshifts, suggesting luminous galaxies do not play a prominent role in concluding reionization. Using various interstellar absorption lines we deduce gas-phase metallicities close to solar indicative of substantial early enrichment. Using selected stellar absorption lines, we model our spectra with a range of metallicities using techniques successfully employed at lower redshift and deduce a stellar metallicity of 0.4 +0.3/-0.1 solar, consistent with the stellar mass - stellar metallicity relation recently found at z~3-5. We discuss the implications of these metallicity estimates for the typical ages of our luminous galaxies and conclude our results imply initial star formation at redshifts z~10, consistent with independent analyses of earlier objects.

We study variability of active galactic nuclei (AGNs) by using the deep optical multiband photometry data obtained from the Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) survey in the COSMOS field. The images analyzed here were taken with 8, 10, 13, and 15 epochs over three years in the $g$, $r$, $i$, and $z$ bands, respectively. We identified 491 robust variable AGN candidates, down to $i=25$ mag and with redshift up to $4.26$. Ninety percent of the variability-selected AGNs are individually identified with the X-ray sources detected in the Chandra COSMOS Legacy survey. We investigate their properties in variability by using structure function analysis and find that the structure function for low-luminosity AGNs ($L_{\mathrm{bol}}\lesssim10^{45}$ erg s$^{-1}$) shows a positive correlation with luminosity, which is the opposite trend for the luminous quasars. This trend is likely to be caused by larger contribution of the host galaxy light for lower-luminosity AGNs. Using the model templates of galaxy spectra, we evaluate the amount of host galaxy contribution to the structure function analysis and find that dominance of the young stellar population is needed to explain the observed luminosity dependence. This suggests that low-luminosity AGNs at $0.8\lesssim z\lesssim1.8$ are predominantly hosted in star-forming galaxies. The X-ray stacking analysis reveals the significant emission from the individually X-ray undetected AGNs in our variability-selected sample. The stacked samples show very large hardness ratios in their stacked X-ray spectrum, which suggests that these optically variable sources have large soft X-ray absorption by dust-free gas.

We present ALMA [CII] line and far-infrared (FIR) continuum observations of three $z > 6$ low-luminosity quasars ($M_{\rm 1450} > -25$) discovered by our Subaru Hyper Suprime-Cam (HSC) survey. The [CII] line was detected in all three targets with luminosities of $(2.4 - 9.5) \times 10^8~L_\odot$, about one order of magnitude smaller than optically luminous ($M_{\rm 1450} \lesssim -25$) quasars. The FIR continuum luminosities range from $< 9 \times 10^{10}~L_\odot$ (3$\sigma$ limit) to $\sim 2 \times 10^{12}~L_\odot$, indicating a wide range in star formation rates in these galaxies. Most of the HSC quasars studied thus far show [CII]/FIR luminosity ratios similar to local star-forming galaxies. Using the [CII]-based dynamical mass ($M_{\rm dyn}$) as a surrogate for bulge stellar mass ($M_{\rm bulge}$), we find that a significant fraction of low-luminosity quasars are located on or even below the local $M_{\rm BH} - M_{\rm bulge}$ relation, particularly at the massive end of the galaxy mass distribution. In contrast, previous studies of optically luminous quasars have found that black holes are overmassive relative to the local relation. Given the low luminosities of our targets, we are exploring the nature of the early co-evolution of supermassive black holes and their hosts in a less biased way. Almost all of the quasars presented in this work are growing their black hole mass at much higher pace at $z \sim 6$ than the parallel growth model, in which supermassive black holes and their hosts grow simultaneously to match the local $M_{\rm BH} - M_{\rm bulge}$ relation at all redshifts. As the low-luminosity quasars appear to realize the local co-evolutionary relation even at $z \sim 6$, they should have experienced vigorous starbursts prior to the currently observed quasar phase to catch up with the relation.

This paper presents the second data release of the Hyper Suprime-Cam Subaru Strategic Program, a wide-field optical imaging survey on the 8.2 meter Subaru Telescope. The release includes data from 174 nights of observation through January 2018. The Wide layer data cover about 300 deg^2 in all five broadband filters (grizy) to the nominal survey exposure (10min in gr and 20min in izy). Partially observed areas are also included in the release; about 1100 deg^2 is observed in at least one filter and one exposure. The median seeing in the i-band is 0.6 arcsec, demonstrating the superb image quality of the survey. The Deep (26 deg^2) and UltraDeep (4 deg^2) data are jointly processed and the UltraDeep-COSMOS field reaches an unprecedented depth of i~28 at 5 sigma for point sources. In addition to the broad-bands, narrow-band data are also available in the Deep and UltraDeep fields. This release includes a major update to the processing pipeline, including improved sky subtraction, PSF modeling, object detection, and artifact rejection. The overall data quality has been improved, but this release is not without problems; there is a persistent deblender problem as well as new issues with masks around bright stars. The user is encouraged to review the issue list before utilizing the data for scientific explorations. All the image products as well as catalog products are available for download. The catalogs are also loaded to a database, which provides an easy interface for users to retrieve data for objects of interest. In addition to these main data products, detailed galaxy shape measurements withheld from the Public Data Release 1 (PDR1) are now available to the community. The shape catalog is drawn from the S16A internal release, which has a larger area than PDR1 (160 deg^2). All products are available at the data release site, https://hsc-release.mtk.nao.ac.jp/.

We report the discovery of 28 quasars and 7 luminous galaxies at 5.7 $\le$ z $\le$ 7.0. This is the tenth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The total number of spectroscopically identified objects in SHELLQs has now grown to 93 high-z quasars, 31 high-z luminous galaxies, 16 [O III] emitters at z ~ 0.8, and 65 Galactic cool dwarfs (low-mass stars and brown dwarfs). These objects were found over 900 deg2, surveyed by HSC between 2014 March and 2018 January. The full quasar sample includes 18 objects with very strong and narrow Ly alpha emission, whose stacked spectrum is clearly different from that of other quasars or galaxies. While the stacked spectrum shows N V 1240 emission and resembles that of lower-z narrow-line quasars, the small Ly alpha width may suggest a significant contribution from the host galaxies. Thus these objects may be composites of quasars and star-forming galaxies.

We present deep near-infrared spectroscopy of six quasars at 6.1<z<6.7 with VLT/X-Shooter and Gemini-N/GNIRS. Our objects, originally discovered through a wide-field optical survey with the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP), have the lowest luminosities (-25.5< M1450<-23.1 mag) of the z>5.8 quasars with measured black hole masses. From single-epoch mass measurements based on MgII2798, we find a wide range in black hole masses, from M_BH=10^7.6 to 10^9.3 Msun. The Eddington ratios L_bol/L_Edd range from 0.16 to 1.1, but the majority of the HSC quasars are powered by M_BH=10^9 Msun supermassive black holes (SMBHs) accreting at sub-Eddington rates. The Eddington ratio distribution of the HSC quasars is inclined to lower accretion rates than those of Willott et al. (2010a), who measured the black hole masses for similarly faint z=6 quasars. This suggests that the global Eddington ratio distribution is wider than has previously been thought. The presence of M_BH=10^9 Msun SMBHs at z=6 cannot be explained with constant sub-Eddington accretion from stellar remnant seed black holes. Therefore, we may be witnessing the first buildup of the most massive black holes in the first billion years of the universe, the accretion activity of which is transforming from active growth to a quiescent phase. Measurements of a larger complete sample of z>6 low-luminosity quasars, as well as deeper observations with future facilities will enable us to better understand the early SMBH growth in the reionization epoch.

We have obtained three-dimensional maps of the universe in $\sim200\times200\times80$ comoving Mpc$^3$ (cMpc$^3$) volumes each at $z=5.7$ and $6.6$ based on a spectroscopic sample of 179 galaxies that achieves $\gtrsim80$\% completeness down to the Ly$\alpha$ luminosity of $\log(L_{\rm Ly\alpha}/[\mathrm{erg\ s^{-1}}])=43.0$, based on our Keck and Gemini observations and the literature. The maps reveal filamentary large-scale structures and two remarkable overdensities made out of at least 44 and 12 galaxies at $z=5.692$ (z57OD) and $z=6.585$ (z66OD), respectively, making z66OD the most distant overdensity spectroscopically confirmed to date with $>10$ spectroscopically confirmed galaxies. We compare spatial distributions of submillimeter galaxies at $z\simeq 4-6$ with our $z=5.7$ galaxies forming the large-scale structures, and detect a $99.97\%$ signal of cross correlation, indicative of a clear coincidence of dusty star-forming galaxy and dust unobscured galaxy formation at this early epoch. The galaxies in z57OD and z66OD are actively forming stars with star formation rates (SFRs) $\gtrsim5$ times higher than the main sequence, and particularly the SFR density in z57OD is 10 times higher than the cosmic average at the redshift (a.k.a. the Madau-Lilly plot). Comparisons with numerical simulations suggest that z57OD and z66OD are protoclusters that are progenitors of the present-day clusters with halo masses of $\sim10^{14}\ \mathrm{M_\odot}$.

We present results of the cross-correlation Lya intensity mapping with Subaru/Hyper Suprime-Cam (HSC) ultra-deep narrowband images and Lya emitters (LAEs) at $z=5.7$ and $6.6$ in a total area of $4$ deg$^2$. Although overwhelming amount of data quality controls have been performed for the narrowband images and the LAE samples, we further conduct extensive analysis evaluating systematics of large-scale point spread function wings, sky subtractions, and unknown errors on the basis of physically uncorrelated signals and sources found in real HSC images and object catalogs, respectively. Removing the systematics, we carefully calculate cross-correlations between Lya intensity of the narrowband images and the LAEs. We identify very diffuse Lya emission with the $3\sigma$ ($2\sigma$) significance at $>150$ comoving kpc (ckpc) far from the LAEs at $z = 5.7$ ($6.6$), beyond a virial radius of star-forming galaxies with $M_\mathrm{h} \sim 10^{11} M_\odot$. The diffuse Lya emission possibly extends up to $1,000$ ckpc with the surface brightness of $10^{-20}-10^{-19}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. We confirm that the small-scale ($<150$ ckpc) Lya radial profiles of LAEs in our Lya intensity maps are consistent with those obtained by recent MUSE observations (Leclercq et al. 2017). Comparisons with numerical simulations suggest that the large-scale ($\sim 150-1,000$ ckpc) Lya emission are not explained by unresolved faint sources of neighboring galaxies including satellites, but by a combination of Lya photons emitted from the central LAE and other unknown sources, such as a cold-gas stream and galactic outflow. We find no evolution in the Lya radial profiles of our LAEs from $z=5.7$ to $6.6$, where theoretical models predict a flattening of the profile slope made by cosmic reionization, albeit with our moderately large observational errors.

We present physical properties of radio galaxies (RGs) with $f_{\rm 1.4 GHz} >$ 1 mJy discovered by Subaru Hyper Supreme-Cam (HSC) and VLA Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) survey. For 1056 FIRST RGs at $0 < z \leq 1.7$ with HSC counterparts in about 100 deg$^2$, we compiled multi-wavelength data of optical, near-infrared (IR), mid-IR, far-IR, and radio (150 MHz). We derived their color excess ($E (B-V)_{*}$), stellar mass, star formation rate (SFR), IR luminosity, the ratio of IR and radio luminosity ($q_{\rm IR}$), and radio spectral index ($\alpha_{\rm radio}$) that are derived from the SED fitting with CIGALE. We also estimated Eddington ratio based on stellar mass and integration of the best-fit SEDs of AGN component. We found that $E (B-V)_{*}$, SFR, and IR luminosity clearly depend on redshift while stellar mass, $q_{\rm IR}$, and $\alpha_{\rm radio}$ do not significantly depend on redshift. Since optically-faint ($i_{\rm AB} \geq 21.3$) RGs that are newly discovered by our RG survey tend to be high redshift, they tend to not only have a large dust extinction and low stellar mass but also have high SFR and AGN luminosity, high IR luminosity, and high Eddington ratio compared to optically-bright ones. The physical properties of a fraction of RGs in our sample seem to differ from a classical view of RGs with massive stellar mass, low SFR, and low Eddington ratio, demonstrating that our RG survey with HSC and FIRST provides us curious RGs among entire RG population.

We report on the optical properties of infrared (IR)-bright dust-obscured galaxies (DOGs) that are defined as (i-[22])_AB > 7.0. Because supermassive black holes (SMBHs) in IR-bright DOGs are expected to be rapidly growing in the major merger scenario, they provide useful clues for understanding the co-evolution of SMBHs and their host galaxies. However, the optical properties of IR-bright DOGs remain unclear because the optical emission of a DOG is very faint. By combining ~ 105 deg^2 images of the optical, near-IR, and mid-IR data obtained from the Subaru Hyper Suprime-Cam (HSC) survey, the VISTA VIKING survey, and the WISE all-sky survey, respectively, 571 IR-bright DOGs were selected. We found that IR-bright DOGs show a redder (g-z)_AB color than other populations of dusty galaxies, such as ultra-luminous IR galaxies (ULIRGs) at a similar redshift, with a significantly large dispersion. Among the selected DOGs, star-formation (SF) dominated DOGs show a relatively red color, while active galactic nucleus (AGN) dominated DOGs show a rather blue color in optical. This result is consistent with the idea that the relative AGN contribution in the optical emission becomes more significant at a later stage in the major merger scenario. We discovered eight IR-bright DOGs showing a significant blue excess in blue HSC bands (BluDOGs). This blue excess can be interpreted as a leaked AGN emission that is either a directly leaking or a scattered AGN emission, as proposed for some blue-excess Hot DOGs in earlier studies.

We report the discovery of a quasar at z = 7.07, which was selected from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. This quasar, HSC J124353.93+010038.5, has an order of magnitude lower luminosity than do the other known quasars at z > 7. The rest-frame ultraviolet absolute magnitude is M1450 = -24.13 +/- 0.08 mag and the bolometric luminosity is Lbol = (1.4 +/- 0.1) x 10^{46} erg/s. Its spectrum in the optical to near-infrared shows strong emission lines, and shows evidence for a fast gas outflow, as the C IV line is blueshifted and there is indication of broad absorption lines. The Mg II-based black hole mass is Mbh = (3.3 +/- 2.0) x 10^8 Msun, thus indicating a moderate mass accretion rate with an Eddington ratio 0.34 +/- 0.20. It is the first z > 7 quasar with sub-Eddington accretion, besides being the third most distant quasar, known to date. The luminosity and black hole mass are comparable to, or even lower than, those measured for the majority of low-z quasars discovered by the Sloan Digital Sky Survey, and thus this quasar likely represents a z > 7 counterpart to quasars commonly observed in the low-z universe.

Using the data release 1 of the Hyper Suprime-Cam Subaru Strategic Program covering about 130 square degrees at high galactic latitudes, we have obtained L dwarf counts based on the selection criteria on colors, limiting magnitude and PSF morphology using $i$, $z$, and $y$ bands. 3665 L dwarfs brighter than $z=24$ have been detected by these criteria. The surface number counts obtained differentially in $z$ magnitude are compared with predictions of an exponential disk model to estimate the thin-disk scale height in the vicinity of the Sun. In the exponential disk model, we first fix the local luminosity function (LLF) to the mean LLF of Cruz et al. (2007) and derive the best fit scale height of 260 pc. However this fit appears to be poor. We then allow the LLF to vary along with the scale height. We use the LLF of Cruz et al. as a starting point. The best-fit model is found for the vertical scale height of 380 pc. However the $\chi^2$ minimum is rather broad and the 90% confidence interval is between 320 and 520 pc. We investigate another model by varying the scale height and the density of the brightest magnitude bin, while other magnitude bins are fixed to the mean LLF of Cruz et al.. We find an equally good fit with the two free parameters and the best-fit scale height is again 380 pc, but the 90% confidence interval is between 340 and 420 pc.

We present new measurements of the quasar luminosity function (LF) at $z \sim 6$, over an unprecedentedly wide range of the rest-frame ultraviolet luminosity $M_{1450}$ from $-30$ to $-22$ mag. This is the fifth in a series of publications from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The LF was calculated with a complete sample of 110 quasars at $5.7 \le z \le 6.5$, which includes 48 SHELLQs quasars discovered over 650 deg$^2$, and 63 brighter quasars discovered by the Sloan Digital Sky Survey and the Canada-France-Hawaii Quasar Survey (including one overlapping object). This is the largest sample of $z \sim 6$ quasars with a well-defined selection function constructed to date, and has allowed us to detect significant flattening of the LF at its faint end. A double power-law function fit to the sample yields a faint-end slope $\alpha = -1.23^{+0.44}_{-0.34}$, a bright-end slope $\beta = -2.73^{+0.23}_{-0.31}$, a break magnitude $M_{1450}^* = -24.90^{+0.75}_{-0.90}$, and a characteristic space density $\Phi^* = 10.9^{+10.0}_{-6.8}$ Gpc$^{-3}$ mag$^{-1}$. Integrating this best-fit model over the range $-18 < M_{1450} < -30$ mag, quasars emit ionizing photons at the rate of $\dot{n}_{\rm ion} = 10^{48.8 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$ at $z = 6.0$. This is less than 10 % of the critical rate necessary to keep the intergalactic medium ionized, which indicates that quasars are not a major contributor to cosmic reionization.

We report the result of optical identifications of FIRST radio sources with the Hyper Suprime-Cam Subaru Strategic Program survey (HSC-SSP). The positional cross-match within 1" between the FIRST and HSC-SSP catalogs (i ~< 26) produced more than 3600 optical counterparts in the 156 deg^2 of the HSC-SSP field. The matched counterparts account for more than 50% of the FIRST sources in the search field, which substantially exceed previously reported fractions of SDSS counterparts (i ~< 22) of ~30%. Among the matched sample, 9% are optically unresolved sources such as radio-loud quasars. The optically faint (i > 21) radio galaxies (RGs) show that the fitting linear function of the 1.4 GHz source counts has a slope that is flatter than that of the bright RGs, while optically faint radio quasars show a slope steeper than that of bright radio quasars. The optically faint RGs show a flat slope in the i-band number counts down to 24 mag, implying either less massive or distant radio-active galactic nuclei (AGNs) beyond 24 mag. The photometric redshift and the comparison of colors with the galaxy models show that most of the matched RGs are distributed at redshifts from 0 to 1.5. The optically faint sample includes the high radio-loudness sources that are not seen in the optically bright sample. Such sources are located at redshift z > 1. This study gives ~1500 radio AGNs lying at the optically faint end and high-redshift regime not probed by previous searches.

We reanalyze the Imaging Photopolarimeter data from Pioneer 10 to study the zodiacal light in the B and R bands beyond Earth orbit, applying an improved method to subtract integrated star light (ISL) and diffuse Galactic light (DGL). We found that there exists a significant instrumental offset, making it difficult to examine the absolute sky brightness. Instead, we analyzed the differential brightness, i.e., the difference in sky brightness from the average at high ecliptic latitude, and compared with that expected from the model zodiacal light. At a heliocentric distance of r<2 au, we found a fairly good correlation between the J-band model zodiacal light and the residual sky brightness after subtracting the ISL and DGL. The reflectances of the interplanetary dust derived from the correlation study are marginally consistent with previous works. The zodiacal light is not significantly detectable at r>3 au, as previously reported. However, a clear discrepancy from the model is found at r=2.94 au which indicates the existence of a local dust cloud produced by the collision of asteroids or dust trail from active asteroids (or main-belt comets). Our result confirms that the main component of the zodiacal light (smooth cloud) is consistent with the model even beyond the earth orbit, which justifies the detection of the extragalactic background light after subtracting the zodiacal light based on the model.

We present measurements of the mean and scatter of the IGM Lyman-{\alpha} opacity at 4.9 < z < 6.1 along the lines of sight of 62 quasars at z > 5.7, the largest sample assembled at these redshifts to date by a factor of two. The sample size enables us to sample cosmic variance at these redshifts more robustly than ever before. The spectra used here were obtained by the SDSS, DES-VHS and SHELLQs collaborations, drawn from the ESI and X-Shooter archives, reused from previous studies or observed specifically for this work. We measure the effective optical depth of Lyman-{\alpha} in bins of 10, 30, 50 and 70 cMpc h-1, construct cumulative distribution functions under two treatments of upper limits on flux and explore an empirical analytic fit to residual Lyman-{\alpha} transmission. We verify the consistency of our results with those of previous studies via bootstrap re-sampling and confirm the existence of tails towards high values in the opacity distributions, which may persist down to z = 5.2. Comparing our results with predictions from cosmological simulations, we find further strong evidence against models that include a spatially uniform ionizing background and temperature-density relation. We also compare to IGM models that include either a fluctuating UVB dominated by rare quasars or temperature fluctuations due to patchy reionization. Although both models produce better agreement with the observations, neither fully captures the observed scatter in IGM opacity. Our sample of 62 z > 5.7 quasar spectra opens many avenues for future study of the reionisation epoch.

We report discovery of 41 new high-z quasars and luminous galaxies, which were spectroscopically identified at 5.7 < z < 6.9. This is the fourth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, based on the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We selected the photometric candidates by a Bayesian probabilistic algorithm, and then carried out follow-up spectroscopy with the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous papers, we have now spectroscopically identified 137 extremely-red HSC sources over about 650 deg2, which include 64 high-z quasars, 24 high-z luminous galaxies, 6 [O III] emitters at z ~ 0.8, and 43 Galactic cool dwarfs (low-mass stars and brown dwarfs). The new quasars span the luminosity range from M1450 ~ -26 to -22 mag, and continue to populate a few magnitude lower luminosities than have been probed by previous wide-field surveys. In a companion paper, we derive the quasar luminosity function at z ~ 6 over an unprecedentedly wide range of M1450 ~ -28 to -21 mag, exploiting the SHELLQs and other survey outcomes.

We investigate Lya, [OIII]5007, Ha, and [CII]158um emission from 1124 galaxies at z=4.9-7.0. Our sample is composed of 1092 Lya emitters (LAEs) at z=4.9, 5.7, 6.6, and 7.0 identified by Subaru/Hyper Suprime-Cam (HSC) narrowband surveys covered by Spitzer large area survey with Subaru/HSC (SPLASH) and 34 galaxies at z=5.148-7.508 with deep ALMA [CII]158um data in the literature. Fluxes of strong rest-frame optical lines of [OIII] and Ha (Hb) are constrained by significant excesses found in the SPLASH 3.6 and 4.5um photometry. At z=4.9, we find that the rest-frame Ha equivalent width and the Lya escape fraction f_Lya positively correlate with the rest-frame Lya equivalent width EW^0_Lya. The f_Lya-EW^0_Lya correlation is similarly found at z~0-2, suggesting no evolution of the correlation over z~0-5. The typical ionizing photon production efficiency of LAEs is logxi_ion/[Hz erg^-1]~25.5 significantly (60-100%) higher than those of LBGs at a given UV magnitude. At z=5.7-7.0, there exists an interesting turn-over trend that the [OIII]/Ha flux ratio increases in EW^0_Lya~0-30 A, and then decreases out to EW^0_Lya~130 A. We also identify an anti-correlation between a [CII] luminosity to star-formation rate ratio (L_[CII]/SFR) and EW^0_Lya at the >99% confidence level. We carefully investigate physical origins of the correlations with stellar-synthesis and photoionization models, and find that a simple anti-correlation between EW_Lya^0 and metallicity explains self-consistently all of the correlations of Lya, Ha, [OIII]/Ha, and [CII] identified in our study, indicating detections of metal-poor (~0.03 Zo) galaxies with EW^0_Lya~200 A.

Narrow-line regions excited by active galactic nuclei (AGN) are important for studying AGN photoionization and feedback. Their strong [O III] lines can be detected with broadband images, allowing morphological studies of these systems with large-area imaging surveys. We develop a new technique to reconstruct the [O III] images using the Subaru Hyper Suprime-Cam (HSC) Survey aided with spectra from the Sloan Digital Sky Survey (SDSS). The technique involves a careful subtraction of the galactic continuum to isolate emission from the [O III]$\lambda$5007 and [O III]$\lambda$4959 lines. Compared to traditional targeted observations, this technique is more efficient at covering larger samples with less dedicated observational resources. We apply this technique to an SDSS spectroscopically selected sample of 300 obscured AGN at redshifts 0.1 - 0.7, uncovering extended emission-line region candidates with sizes up to tens of kpc. With the largest sample of uniformly derived narrow-line region sizes, we revisit the narrow-line region size-luminosity relation. The area and radii of the [O III] emission-line regions are strongly correlated with the AGN luminosity inferred from the mid-infrared (15 $\mu$m rest-frame) with a power-law slope of $0.62^{+0.05}_{-0.06} \pm 0.10$ (statistical and systemic errors), consistent with previous spectroscopic findings. We discuss the implications for the physics of AGN emission-line region and future applications of this technique, which should be useful for current and next-generation imaging surveys to study AGN photoionization and feedback with large statistical samples.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

We present our ALMA Cycle 4 measurements of the [CII] emission line and the underlying far-infrared (FIR) continuum emission from four optically low-luminosity ($M_{\rm 1450} > -25$) quasars at $z \gtrsim 6$ discovered by the Subaru Hyper Suprime Cam (HSC) survey. The [CII] line and FIR continuum luminosities lie in the ranges $L_{\rm [CII]} = (3.8-10.2) \times 10^8~L_\odot$ and $L_{\rm FIR} = (1.2-2.0) \times 10^{11}~L_\odot$, which are at least one order of magnitude smaller than those of optically-luminous quasars at $z \gtrsim 6$. We estimate the star formation rates (SFR) of our targets as $\simeq 23-40~M_\odot ~{\rm yr}^{-1}$. Their line and continuum-emitting regions are marginally resolved, and found to be comparable in size to those of optically luminous quasars, indicating that their SFR or likely gas mass surface densities (key controlling parameter of mass accretion) are accordingly different. The $L_{\rm [CII]}/L_{\rm FIR}$ ratios of the hosts, $\simeq (2.2-8.7) \times 10^{-3}$, are fully consistent with local star-forming galaxies. Using the [CII] dynamics, we derived their dynamical masses within a radius of 1.5-2.5 kpc as $\simeq (1.4-8.2) \times 10^{10}~M_\odot$. By interpreting these masses as stellar ones, we suggest that these faint quasar hosts are on or even below the star-forming main sequence at $z \sim 6$, i.e., they appear to be transforming into quiescent galaxies. This is in contrast to the optically luminous quasars at those redshifts, which show starburst-like properties. Finally, we find that the ratios of black hole mass to host galaxy dynamical mass of the most of low-luminosity quasars including the HSC ones are consistent with the local value. The mass ratios of the HSC quasars can be reproduced by a semi-analytical model that assumes merger-induced black hole-host galaxy evolution.

Collisions and interactions between gas-rich galaxies are thought to be pivotal stages in their formation and evolution, causing the rapid production of new stars, and possibly serving as a mechanism for fueling supermassive black holes (BH). Harnessing the exquisite spatial resolution (~0.5 arcsec) afforded by the first ~170 deg^2 of the Hyper Suprime-Cam (HSC) Survey, we present our new constraints on the importance of galaxy-galaxy major mergers (1:4) in growing BHs throughout the last ~8 Gyrs. Utilizing mid-infrared observations in the WISE All-Sky survey, we robustly select active galactic nuclei (AGN) and mass-matched control galaxy samples, totaling ~140,000 spectroscopically confirmed systems at i<22 mag. We identify galaxy interaction signatures using a novel machine-learning random forest decision tree technique allowing us to select statistically significant samples of major-mergers, minor-mergers/irregular-systems, and non-interacting galaxies. We use these samples to show that galaxies undergoing mergers are a factor ~2-7 more likely to contain luminous obscured AGN than non-interacting galaxies, and this is independent of both stellar mass and redshift to z < 0.9. Furthermore, based on our comparison of AGN fractions in mass-matched samples, we determine that the most luminous AGN population (L_AGN > 10^45 erg/s) systematically reside in merging systems over non-interacting galaxies. Our findings show that galaxy-galaxy interactions do, on average, trigger luminous AGN activity substantially more often than in secularly evolving non-interacting galaxies, and we further suggest that the BH growth rate may be closely tied to the dynamical time of the merger system.

We constrain the quasar contribution to cosmic reionization based on our deep optical survey of z~6 quasars down to z_R=24.15 using Subaru/Suprime-Cam in three UKIDSS-DXS fields covering 6.5 deg^2. In Kashikawa et al. (2015), we select 17 quasar candidates and report our initial discovery of two low-luminosity quasars (M_1450~ -23) from seven targets, one of which might be a Lyman alpha emitting galaxy. From an additional optical spectroscopy, none of the four candidates out of the remaining ten turn out to be genuine quasars. Moreover, the deeper optical photometry provided by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) shows that, unlike the two already-known quasars, the i-z and z-y colors of the last six candidates are consistent with M- or L-type brown dwarfs. Therefore, the quasar luminosity function (QLF) in the previous paper is confirmed. Compiling QLF measurements from the literature over a wide magnitude range, including an extremely faint AGN candidate from Parsa et al. (2017}, to fit them with a double power-law, we find that the best-fit faint-end slope is alpha=-2.04^+0.33_-0.18 (-1.98^+0.48_-0.21) and characteristic magnitude is M_1450^*=-25.8^+1.1_-1.9 (-25.7^+1.0_-1.8) in the case of two (one) quasar detection. Our result suggests that, if the QLF is integrated down to M_1450=-18, quasars produce ~1-12% of the ionizing photons required to ionize the whole universe at z~6 with 2sigma confidence level, assuming that the escape fraction is f_esc=1 and the IGM clumpy factor is C=3. Even when the systematic uncertainties are taken into account, our result supports the scenario that quasars are the minor contributors of reionization.

We study the UV luminosity functions (LFs) at $z\sim 4$, $5$, $6,$ and $7$ based on the deep large-area optical images taken by the Hyper Suprime-Cam (HSC) Subaru strategic program (SSP). On the 100 deg$^2$ sky of the HSC SSP data available to date, we make enormous samples consisting of a total of 579,565 dropout candidates at $z\sim 4-7$ by the standard color selection technique, 358 out of which are spectroscopically confirmed by our follow-up spectroscopy and other studies. We obtain UV LFs at $z \sim 4-7$ that span a very wide UV luminosity range of $\sim 0.002 - 100 \, L_{\rm UV}^\ast$ ($-26 < M_{\rm UV} < -14$ mag) by combining LFs from our program and the ultra-deep Hubble Space Telescope legacy surveys. We derive three parameters of the best-fit Schechter function, $\phi^\ast$, $M_{\rm UV}^\ast$, and $\alpha$, of the UV LFs in the magnitude range where the AGN contribution is negligible, and find that $\alpha$ and $\phi^\ast$ decrease from $z\sim 4$ to $7$ with no significant evolution of $M_{\rm UV}^\ast$. Because our HSC SSP data bridge the LFs of galaxies and AGNs with great statistical accuracy, we carefully investigate the bright end of the galaxy UV LFs that are estimated by the subtraction of the AGN contribution either aided with spectroscopy or the best-fit AGN UV LFs. We find that the bright end of the galaxy UV LFs cannot be explained by the Schechter function fits at $> 2 \sigma$ significance, and require either double power-law functions or modified Schechter functions that consider a magnification bias due to gravitational lensing.

We investigate the galaxy overdensity around proto-cluster scale quasar pairs at high (z>3) and low (z~1) redshift based on the unprecedentedly wide and deep optical survey of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Using the first-year survey data covering effectively ~121 deg^2 with the 5sigma depth of i~26.4 and the SDSS DR12Q catalog, we find two luminous pairs at z~3.3 and 3.6 which reside in >5sigma overdense regions of g-dropout galaxies at i<25. The projected separations of the two pairs are R_perp=1.75 and 1.04 proper Mpc, and their velocity offsets are Delta V=692 and 1448 km s^{-1}, respectively. This result is in clear contrast to the average z~4 quasar environments as discussed in Uchiyama et al. (2017) and implies that the quasar activities of the pair members are triggered via major mergers in proto-clusters, unlike the vast majority of isolated quasars in general fields that may turn on via non-merger events such as bar and disk instabilities. At z~1, we find 37 pairs with R_perp<2 pMpc and Delta V<2300 km s^{-1} in the current HSC-Wide coverage, including four from Hennawi et al. (2006). The distribution of the peak overdensity significance within two arcminutes around the pairs has a long tail toward high density (>4sigma) regions. Thanks to the large sample size, we find a statistical evidence that this excess is unique to the pair environments when compared to single quasar and randomly selected galaxy environments at the same redshift range. Moreover, there are nine small-scale (R_perp<1 pMpc) pairs, two of which are found to reside in cluster fields. Our results demonstrate that <2 pMpc-scale quasar pairs at both redshift range tend to occur in massive haloes, although perhaps not the most massive ones, and that they are useful to search for rare density peaks.

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/.

We present a new analysis of the stellar mass function and morphology of recently-quenched galaxies (RQGs), whose star formation has been recently quenched for some reason. The COSMOS2015 catalog was exploited to select those galaxies at 0.2 < z < 4.8, over 1.5 deg2 of the Cosmic Evolution Survey (COSMOS) UltraVISTA field. This is the first time that RQGs are consistently selected and studied in such a wide range of redshift. We find increasing number density of RQGs with time in a broad mass range at z>1, while low-mass RQGs start to grow very rapidly at z < 1. We also demonstrate that the migration of RQGs may largely drive the evolution of the stellar mass function of passive galaxies. Moreover, we find that the morphological type distribution of RQGs are intermediate between those of star-forming and passive galaxies. These results indicate that RQGs represent a major transitional phase of galaxy evolution, in which star-forming galaxies turn into passive galaxies, accompanied by the build up of spheroidal component.

We present spectroscopic identification of 32 new quasars and luminous galaxies discovered at 5.7 < z < 6.8. This is the second in a series of papers presenting the results of the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The photometric candidates were selected by a Bayesian probabilistic algorithm, and then observed with spectrographs on the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous paper, we have now identified 64 HSC sources over about 430 deg2, which include 33 high-z quasars, 14 high-z luminous galaxies, 2 [O III] emitters at z ~ 0.8, and 15 Galactic brown dwarfs. The new quasars have considerably lower luminosity (M1450 ~ -25 to -22 mag) than most of the previously known high-z quasars. Several of these quasars have luminous (> 10^(43) erg/s) and narrow (< 500 km/s) Ly alpha lines, and also a possible mini broad absorption line system of N V 1240 in the composite spectrum, which clearly separate them from typical quasars. On the other hand, the high-z galaxies have extremely high luminosity (M1450 ~ -24 to -22 mag) compared to other galaxies found at similar redshift. With the discovery of these new classes of objects, we are opening up new parameter spaces in the high-z Universe. Further survey observations and follow-up studies of the identified objects, including the construction of the quasar luminosity function at z ~ 6, are ongoing.

Active galactic nuclei (AGNs) in the high-redshift Universe are thought to reside in overdense environments. However, recent works provide controversial results partly due to the use of different techniques and possible suppression of nearby galaxy formation by AGN feedback. We conducted deep and wide-field imaging observations with the Suprime-Cam on the Subaru telescope and searched for Lyman-alpha emitters (LAEs) around two QSOs (quasi-stellar objects) at $z\sim4.9$ and a radio galaxy at $z\sim4.5$ by using narrow-band filters to address these issues more robustly. In the QSO fields, we obtained additional broad-band images to select Lyman-break galaxies (LBGs) at $z\sim5$ for comparison. We constructed a photometric sample of 301 LAEs and 170 LBGs in total. A wide field of view (34arcmin$\times$27arcmin, corresponding to 80$\times$60 comoving Mpc$^2$) of the Suprime-Cam enabled us to probe galaxies in the immediate vicinities of the AGNs and in the blank fields simultaneously and compare various properties of them in a consistent manner. The two QSOs are located near local density peaks ($<2\sigma$) and one of the QSOs has a close companion LAE with projected separation of 80 physical kpc. The radio galaxy is found to be near a void of LAEs. The number densities of LAEs/LGBs in a larger spatial scale around the AGNs are not significantly different from those in blank fields. No sign of feedback is found down to $L_{Ly\alpha}\sim10^{41.8}\mathrm{~erg~s^{-1}}$. Our results suggest that high-redshift AGNs are not associated with extreme galaxy overdensity and that this cannot be attributed to the effect of AGN feedback.

We examine the clustering of quasars over a wide luminosity range, by utilizing 901 quasars at $\overline{z}_{\rm phot}\sim3.8$ with $-24.73<M_{\rm 1450}<-22.23$ photometrically selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) S16A Wide2 date release and 342 more luminous quasars at $3.4<z_{\rm spec}<4.6$ having $-28.0<M_{\rm 1450}<-23.95$ from the Sloan Digital Sky Survey (SDSS) that fall in the HSC survey fields. We measure the bias factors of two quasar samples by evaluating the cross-correlation functions (CCFs) between the quasar samples and 25790 bright $z\sim4$ Lyman Break Galaxies (LBGs) in $M_{\rm 1450}<-21.25$ photometrically selected from the HSC dataset. Over an angular scale of \timeform{10.0"} to \timeform{1000.0"}, the bias factors are $5.93^{+1.34}_{-1.43}$ and $2.73^{+2.44}_{-2.55}$ for the low and high luminosity quasars, respectively, indicating no luminosity dependence of quasar clustering at $z\sim4$. It is noted that the bias factor of the luminous quasars estimated by the CCF is smaller than that estimated by the auto-correlation function (ACF) over a similar redshift range, especially on scales below \timeform{40.0"}. Moreover, the bias factor of the less-luminous quasars implies the minimal mass of their host dark matter halos (DMHs) is $0.3$-$2\times10^{12}h^{-1}M_{\odot}$, corresponding to a quasar duty cycle of $0.001$-$0.06$.

We present the luminosity function of z=4 quasars based on the Hyper Suprime-Cam Subaru Strategic Program Wide layer imaging data in the g, r, i, z, and y bands covering 339.8 deg^2. From stellar objects, 1666 z~4 quasar candidates are selected by the g-dropout selection down to i=24.0 mag. Their photometric redshifts cover the redshift range between 3.6 and 4.3 with an average of 3.9. In combination with the quasar sample from the Sloan Digital Sky Survey in the same redshift range, the quasar luminosity function covering the wide luminosity range of M1450=-22 to -29 mag is constructed. It is well described by a double power-law model with a knee at M1450=-25.36+-0.13 mag and a flat faint-end slope with a power-law index of -1.30+-0.05. The knee and faint-end slope show no clear evidence of redshift evolution from those at z~2. The flat slope implies that the UV luminosity density of the quasar population is dominated by the quasars around the knee, and does not support the steeper faint-end slope at higher redshifts reported at z>5. If we convert the M1450 luminosity function to the hard X-ray 2-10keV luminosity function using the relation between UV and X-ray luminosity of quasars and its scatter, the number density of UV-selected quasars matches well with that of the X-ray-selected AGNs above the knee of the luminosity function. Below the knee, the UV-selected quasars show a deficiency compared to the hard X-ray luminosity function. The deficiency can be explained by the lack of obscured AGNs among the UV-selected quasars.

X-ray and gamma-ray polarimetry is a promising tool to study the geometry and the magnetic configuration of various celestial objects, such as binary black holes or gamma-ray bursts (GRBs). However, statistically significant polarizations have been detected in few of the brightest objects. Even though future polarimeters using X-ray telescopes are expected to observe weak persistent sources, there are no effective approaches to survey transient and serendipitous sources with a wide field of view (FoV). Here we present an electron-tracking Compton camera (ETCC) as a highly-sensitive gamma-ray imaging polarimeter. The ETCC provides powerful background rejection and a high modulation factor over a FoV of up to 2$\pi$ sr thanks to its excellent imaging based on a well-defined point spread function. Importantly, we demonstrated for the first time the stability of the modulation factor under realistic conditions of off-axis incidence and huge backgrounds using the SPring-8 polarized X-ray beam. The measured modulation factor of the ETCC was 0.65 $\pm$ 0.01 at 150 keV for an off-axis incidence with an oblique angle of 30$^\circ$ and was not degraded compared to the 0.58 $\pm$ 0.02 at 130 keV for on-axis incidence. These measured results are consistent with the simulation results. Consequently, we found that the satellite-ETCC proposed in Tanimori et al. (2015) would provide all-sky surveys of weak persistent sources of 13 mCrab with 10% polarization for a 10$^{7}$ s exposure and over 20 GRBs down to a $6\times10^{-6}$ erg cm$^{-2}$ fluence and 10% polarization during a one-year observation.

We have developed an Electron Tracking Compton Camera (ETCC), which provides a well-defined Point Spread Function (PSF) by reconstructing a direction of each gamma as a point and realizes simultaneous measurement of brightness and spectrum of MeV gamma-rays for the first time. Here, we present the results of our on-site pilot gamma-imaging-spectroscopy with ETCC at three contaminated locations in the vicinity of the Fukushima Daiichi Nuclear Power Plants in Japan in 2014. The obtained distribution of brightness (or emissivity) with remote-sensing observations is unambiguously converted into the dose distribution. We confirm that the dose distribution is consistent with the one taken by conventional mapping measurements with a dosimeter physically placed at each grid point. Furthermore, its imaging spectroscopy, boosted by Compton-edge-free spectra, reveals complex radioactive features in a quantitative manner around each individual target point in the background-dominated environment. Notably, we successfully identify a "micro hot spot" of residual caesium contamination even in an already decontaminated area. These results show that the ETCC performs exactly as the geometrical optics predicts, demonstrates its versatility in the field radiation measurement, and reveals potentials for application in many fields, including the nuclear industry, medical field, and astronomy.

Since the discovery of nuclear gamma-rays, its imaging has been limited to pseudo imaging, such as Compton Camera (CC) and coded mask. Pseudo imaging does not keep physical information (intensity, or brightness in Optics) along a ray, and thus is capable of no more than qualitative imaging of bright objects. To attain quantitative imaging, cameras that realize geometrical optics is essential, which would be, for nuclear MeV gammas, only possible via complete reconstruction of the Compton process. Recently we have revealed that "Electron Tracking Compton Camera" (ETCC) provides a well-defined Point Spread Function (PSF). The information of an incoming gamma is kept along a ray with the PSF and that is equivalent to geometrical optics. Here we present an imaging-spectroscopic measurement with the ETCC. Our results highlight the intrinsic difficulty with CCs in performing accurate imaging, and show that the ETCC surmounts this problem. The imaging capability also helps the ETCC suppress the noise level dramatically by ~3 orders of magnitude without a shielding structure. Furthermore, full reconstruction of Compton process with the ETCC provides spectra free of Compton edges. These results mark the first proper imaging of nuclear gammas based on the genuine geometrical optics.

We report the serendipitous discovery of a quadruply lensed source at $z_{\rm s}=3.76$, HSC~J115252+004733, from the Hyper Suprime-Cam (HSC) Survey. The source is lensed by an early-type galaxy at $z_{\rm l}=0.466$ and a satellite galaxy. Here, we investigate the properties of the source by studying its size and luminosity from the imaging and the luminosity and velocity width of the Ly-$\alpha$ line from the spectrum. Our analyses suggest that the source is most probably a low-luminosity active galactic nucleus (LLAGN) but the possibility of it being a compact bright galaxy (e.g., a Lyman-$\alpha$ emitter or Lyman Break Galaxy) cannot be excluded. The brighter pair of lensed images appears point-like except in the HSC $i$-band (with a seeing $\sim0.5"$). The extended emission in the $i$-band image could be due to the host galaxy underneath the AGN, or alternatively, due to a highly compact lensed galaxy (without AGN) which appears point-like in all bands except in $i$-band. We also find that the flux ratio of the brighter pair of images is different in the Ks-band compared to optical wavelengths. Phenomena such as differential extinction and intrinsic variability cannot explain this chromatic variation. While microlensing from stars in the foreground galaxy is less likely to be the cause, it cannot be ruled out completely. If the galaxy hosts an AGN, then this represents the highest redshift quadruply imaged AGN known to date, enabling study of a distant LLAGN. Discovery of this unusually compact and faint source demonstrates the potential of the HSC survey.

We present an analysis of the blank sky spectra observed with the Faint Object Spectrograph on board the Hubble Space Telescope. We study the diffuse sky emission from ultraviolet to optical wavelengths, which is composed of the zodiacal light (ZL), diffuse Galactic light (DGL), and residual emission. The observations were performed toward 54 fields distributed widely over the sky, with the spectral coverage from 0.2 to 0.7 um. In order to avoid contaminating light from the earthshine, we use the data collected only in orbital nighttime. The observed intensity is decomposed into the ZL, DGL, and residual emission, in eight photometric bands spanning our spectral coverage. We found that the derived ZL reflectance spectrum is flat in the optical, which indicates major contribution of C-type asteroids to the interplanetary dust (IPD). In addition, the ZL reflectance spectrum has an absorption feature at ~0.3 um. The shape of the DGL spectrum is consistent with those found in earlier measurements and model predictions. While the residual emission contains a contribution from the extragalactic background light, we found that the spectral shape of the residual looks similar to the ZL spectrum. Moreover, its optical intensity is much higher than that measured from beyond the IPD cloud by Pioneer10/11, and also than that of the integrated galaxy light. These findings may indicate the presence of an isotropic ZL component, which is missed in the conventional ZL models.

We present measurements of the clustering properties of a sample of infrared (IR) bright dust-obscured galaxies (DOGs). Combining 125 deg$^2$ of wide and deep optical images obtained with the Hyper Suprime-Cam on the Subaru Telescope and all-sky mid-IR (MIR) images taken with Wide-Field Infrared Survey Explorer, we have discovered 4,367 IR-bright DOGs with $(i - [22])_{\rm AB}$ $>$ 7.0 and flux density at 22 $\mu$m $>$ 1.0 mJy. We calculate the angular autocorrelation function (ACF) for a uniform subsample of 1411 DOGs with 3.0 mJy $<$ flux (22 $mu$m) $<$ 5.0 mJy and $i_{\rm AB}$ $<$ 24.0. The ACF of our DOG subsample is well-fit with a single power-law, $\omega (\theta)$ = (0.010 $\pm$ 0.003) $\theta^{-0.9}$, where $\theta$ in degrees. The correlation amplitude of IR-bright DOGs is larger than that of IR-faint DOGs, which reflects a flux-dependence of the DOG clustering, as suggested by Brodwin et al. (2008). We assume that the redshift distribution for our DOG sample is Gaussian, and consider 2 cases: (1) the redshift distribution is the same as IR-faint DOGs with flux at 22 $\mu$m $<$ 1.0 mJy, mean and sigma $z$ = 1.99 $\pm$ 0.45, and (2) $z$ = 1.19 $\pm$ 0.30, as inferred from their photometric redshifts. The inferred correlation length of IR-bright DOGs is $r_0$ = 12.0 $\pm$ 2.0 and 10.3 $\pm$ 1.7 $h^{-1}$ Mpc, respectively. IR-bright DOGs reside in massive dark matter halos with a mass of $\log [\langle M_{\mathrm{h}} \rangle / (h^{-1} M_{\odot})]$ = 13.57$_{-0.55}^{+0.50}$ and 13.65$_{-0.52}^{+0.45}$ in the two cases, respectively.

We report the discovery of 15 quasars and bright galaxies at 5.7 < z < 6.9. This is the initial result from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the exquisite multiband imaging data produced by the Subaru Hyper Suprime-Cam (HSC) Strategic Program survey. The candidate selection is performed by combining several photometric approaches including a Bayesian probabilistic algorithm to reject stars and dwarfs. The spectroscopic identification was carried out with the Gran Telescopio Canarias and the Subaru Telescope for the first 80 deg2 of the survey footprint. The success rate of our photometric selection is quite high, approaching 100 % at the brighter magnitudes (zAB < 23.5 mag). Our selection also recovered all the known high-z quasars on the HSC images. Among the 15 discovered objects, six are likely quasars, while the other six with interstellar absorption lines and in some cases narrow emission lines are likely bright Lyman-break galaxies. The remaining three objects have weak continua and very strong and narrow Ly alpha lines, which may be excited by ultraviolet light from both young stars and quasars. These results indicate that we are starting to see the steep rise of the luminosity function of z > 6 galaxies, compared with that of quasars, at magnitudes fainter than M1450 ~ -22 mag or zAB ~24 mag. Follow-up studies of the discovered objects as well as further survey observations are ongoing.

Recent large-scale infrared surveys have been revealing stellar populations in the inner Galaxy seen through strong interstellar extinction in the disk. In particular, classical Cepheids with their period-luminosity and period-age relations are useful tracers of Galactic structure and evolution. Interesting groups of Cepheids reported recently include four Cepheids in the Nuclear Stellar Disk (NSD), about 200 pc around the Galactic Centre, found by Matsunaga et al. and those spread across the inner part of the disk reported by Dekany and collaborators. We here report our discovery of nearly thirty classical Cepheids towards the bulge region, some of which are common with Dekany et al., and discuss the large impact of the reddening correction on distance estimates for these objects. Assuming that the four Cepheids in the NSD are located at the distance of the Galactic Centre and that the near-infrared extinction law, i.e. wavelength dependency of the interstellar extinction, is not systematically different between the NSD and other bulge lines-of-sight, most of the other Cepheids presented here are located significantly further than the Galactic Centre. This suggests a lack of Cepheids in the inner 2.5 kpc region of the Galactic disk except the NSD. Recent radio observations show a similar distribution of star-forming regions.

We explore the effect of varying the mass of the seed black hole on the resulting black hole mass - bulge mass relation at z ~ 0, using a semi-analytic model of galaxy formation combined with large cosmological N-body simulations. We constrain our model by requiring the observed properties of galaxies at z ~ 0 are reproduced. In keeping with previous semi-analytic models, we place a seed black hole immediately after a galaxy forms. When the mass of the seed is set at 10^5 M_sun, we find that the model results become inconsistent with recent observational results of the black hole mass - bulge mass relation for dwarf galaxies. In particular, the model predicts that bulges with ~ 10^9 M_sun harbour larger black holes than observed. On the other hand, when we employ seed black holes with 10^3 M_sun, or randomly select their mass within a 10^(3-5) M_sun range, the resulting relation is consistent with observation estimates, including the observed dispersion. We find that to obtain stronger constraints on the mass of seed black holes, observations of less massive bulges at z ~ 0 are a more powerful comparison than the relations at higher redshifts.

Using all-sky maps obtained from COBE/DIRBE at 3.5 and 4.9 um, we present a reanalysis of diffuse sky emissions such as zodiacal light (ZL), diffuse Galactic light (DGL), integrated starlight (ISL), and isotropic residual emission including the extragalactic background light (EBL). Our new analysis, which includes an improved estimate of ISL using the Wide-field Infrared Survey Explorer (WISE) data, enabled us to find the DGL signal in a direct linear correlation between diffuse near-infrared and 100 um emission at high Galactic latitudes (|b| > 35 degree). At 3.5um, the high-latitude DGL result is comparable to the low-latitude value derived from the previous DIRBE analysis. In comparison with models of the DGL spectrum assuming a size distribution of dust grains composed of amorphous silicate, graphite, and polycyclic aromatic hydrocarbon (PAH), the measured DGL values at 3.5 and 4.9 um constrain the mass fraction of PAH particles in the total dust species to be more than ~ 2%. This was consistent with the results of Spitzer/IRAC toward the lower Galactic latitude regions. The derived residual emission of 8.9 +/- 3.4 nW m^{-2} sr^{-1} at 3.5 um is marginally consistent with the level of integrated galaxy light and the EBL constraints from the gamma-ray observations. The residual emission at 4.9 um is not significantly detected due to the large uncertainty in the ZL subtraction, same as previous studies. Combined with our reanalysis of the DIRBE data at 1.25 and 2.2 um, the residual emission in the near-infrared exhibits the Rayleigh-Jeans spectrum.

For MeV gamma-ray astronomy, we have developed an electron-tracking Compton camera (ETCC) as a MeV gamma-ray telescope capable of rejecting the radiation background and attaining the high sensitivity of near 1 mCrab in space. Our ETCC comprises a gaseous time-projection chamber (TPC) with a micro pattern gas detector for tracking recoil electrons and a position-sensitive scintillation camera for detecting scattered gamma rays. After the success of a first balloon experiment in 2006 with a small ETCC (using a 10$\times$10$\times$15 cm$^3$ TPC) for measuring diffuse cosmic and atmospheric sub-MeV gamma rays (Sub-MeV gamma-ray Imaging Loaded-on-balloon Experiment I; SMILE-I), a (30 cm)$^{3}$ medium-sized ETCC was developed to measure MeV gamma-ray spectra from celestial sources, such as the Crab Nebula, with single-day balloon flights (SMILE-II). To achieve this goal, a 100-times-larger detection area compared with that of SMILE-I is required without changing the weight or power consumption of the detector system. In addition, the event rate is also expected to dramatically increase during observation. Here, we describe both the concept and the performance of the new data-acquisition system with this (30 cm)$^{3}$ ETCC to manage 100 times more data while satisfying the severe restrictions regarding the weight and power consumption imposed by a balloon-borne observation. In particular, to improve the detection efficiency of the fine tracks in the TPC from $\sim$10\% to $\sim$100\%, we introduce a new data-handling algorithm in the TPC. Therefore, for efficient management of such large amounts of data, we developed a data-acquisition system with parallel data flow.

Quasar host galaxies are key for understanding the relation between galaxies and the supermassive black holes (SMBHs) at their centers. We present a study of 191 broad-line quasars and their host galaxies at z < 1, using high signal-to-noise ratio (SNR) spectra produced by the Sloan Digital Sky Survey Reverberation Mapping project. Clear detection of stellar absorption lines allows a reliable decomposition of the observed spectra into nuclear and host components, using spectral models of quasar and stellar radiations as well as emission lines from the interstellar medium. We estimate age, mass (M*), and velocity dispersion (sigma*) of the host stars, the star formation rate (SFR), quasar luminosity, and SMBH mass (Mbh), for each object. The quasars are preferentially hosted by massive galaxies with M* ~ 10^{11} Msun characterized by stellar ages around a billion years, which coincides with the transition phase of normal galaxies from the blue cloud to the red sequence. The host galaxies have relatively low SFRs and fall below the main sequence of star-forming galaxies at similar redshifts. These facts suggest that the hosts have experienced an episode of major star formation sometime in the past billion years, which was subsequently quenched or suppressed. The derived Mbh - sigma* and Mbh - M* relations agree with our past measurements and are consistent with no evolution from the local Universe. The present analysis demonstrates that reliable measurements of stellar properties of quasar host galaxies are possible with high-SNR fiber spectra, which will be acquired in large numbers with future powerful instruments such as the Subaru Prime Focus Spectrograph.

Using all-sky maps obtained with COBE/DIRBE, we reanalyzed the diffuse sky brightness at 1.25 and 2.2 um, which consists of zodiacal light, diffuse Galactic light (DGL), integrated starlight (ISL), and isotropic emission including the extragalactic background light. Our new analysis including an improved estimate of the DGL and the ISL with the 2MASS data showed that deviations of the isotropic emission from isotropy were less than 10% in the entire sky at high Galactic latitude (|b|>35). The result of our analysis revealed a significantly large isotropic component at 1.25 and 2.2 um with intensities of 60.15 +/- 16.14 and 27.68 +/- 6.21 nWm-2sr-1, respectively. This intensity is larger than the integrated galaxy light, upper limits from gamma-ray observation, and potential contribution from exotic sources (i.e., Population III stars, intrahalo light, direct collapse black holes, and dark stars). We therefore conclude that the excess light may originate from the local universe; the Milky Way and/or the solar system.

An electron-tracking Compton camera (ETCC) is a detector that can determine the arrival direction and energy of incident sub-MeV/MeV gamma-ray events on an event-by-event basis. It is a hybrid detector consisting of a gaseous time projection chamber (TPC), that is the Compton-scattering target and the tracker of recoil electrons, and a position-sensitive scintillation camera that absorbs of the scattered gamma rays, to measure gamma rays in the environment from contaminated soil. To measure of environmental gamma rays from soil contaminated with radioactive cesium (Cs), we developed a portable battery-powered ETCC system with a compact readout circuit and data-acquisition system for the SMILE-II experiment. We checked the gamma-ray imaging ability and ETCC performance in the laboratory by using several gamma-ray point sources. The performance test indicates that the field of view (FoV) of the detector is about 1$\;$sr and that the detection efficiency and angular resolution for 662$\;$keV gamma rays from the center of the FoV is $(9.31 \pm 0.95) \times 10^{^-5}$ and $5.9^{\circ} \pm 0.6^{\circ}$, respectively. Furthermore, the ETCC can detect 0.15$\;\mu\rm{Sv/h}$ from a $^{137}$Cs gamma-ray source with a significance of 5$\sigma$ in 13 min in the laboratory. In this paper, we report the specifications of the ETCC and the results of the performance tests. Furthermore, we discuss its potential use for environmental gamma-ray measurements.

Photon imaging for MeV gammas has serious difficulties due to huge backgrounds and unclearness in images, which are originated from incompleteness in determining the physical parameters of Compton scattering in detection, e.g., lack of the directional information of the recoil electrons. The recent major mission/instrument in the MeV band, Compton Gamma Ray Observatory/COMPTEL, which was Compton Camera (CC), detected mere $\sim30$ persistent sources. It is in stark contrast with $\sim$2000 sources in the GeV band. Here we report the performance of an Electron-Tracking Compton Camera (ETCC), and prove that it has a good potential to break through this stagnation in MeV gamma-ray astronomy. The ETCC provides all the parameters of Compton-scattering by measuring 3-D recoil electron tracks; then the Scatter Plane Deviation (SPD) lost in CCs is recovered. The energy loss rate (dE/dx), which CCs cannot measure, is also obtained, and is found to be indeed helpful to reduce the background under conditions similar to space. Accordingly the significance in gamma detection is improved severalfold. On the other hand, SPD is essential to determine the point-spread function (PSF) quantitatively. The SPD resolution is improved close to the theoretical limit for multiple scattering of recoil electrons. With such a well-determined PSF, we demonstrate for the first time that it is possible to provide reliable sensitivity in Compton imaging without utilizing an optimization algorithm. As such, this study highlights the fundamental weak-points of CCs. In contrast we demonstrate the possibility of ETCC reaching the sensitivity below $1\times10^{-12}$ erg cm$^{-2}$ s$^{-1}$ at 1 MeV.

We present the photometric properties of a sample of infrared (IR) bright dust obscured galaxies (DOGs). Combining wide and deep optical images obtained with the Hyper Suprime-Cam (HSC) on the Subaru Telescope and all-sky mid-IR (MIR) images taken with Wide-Field Infrared Survey Explorer (WISE), we discovered 48 DOGs with $i - K_\mathrm{s} > 1.2$ and $i - [22] > 7.0$, where $i$, $K_\mathrm{s}$, and [22] represent AB magnitude in the $i$-band, $K_\mathrm{s}$-band, and 22 $\mu$m, respectively, in the GAMA 14hr field ($\sim$ 9 deg$^2$). Among these objects, 31 ($\sim$ 65 %) show power-law spectral energy distributions (SEDs) in the near-IR (NIR) and MIR regime, while the remainder show a NIR bump in their SEDs. Assuming that the redshift distribution for our DOGs sample is Gaussian, with mean and sigma $z$ = 1.99 $\pm$ 0.45, we calculated their total IR luminosity using an empirical relation between 22 $\mu$m luminosity and total IR luminosity. The average value of the total IR luminosity is (3.5 $\pm$ 1.1) $\times$ $10^{13}$ L$_{\odot}$, which classifies them as hyper-luminous infrared galaxies (HyLIRGs). We also derived the total IR luminosity function (LF) and IR luminosity density (LD) for a flux-limited subsample of 18 DOGs with 22 $\mu$m flux greater than 3.0 mJy and with $i$-band magnitude brighter than 24 AB magnitude. The derived space density for this subsample is log $\phi$ = -6.59 $\pm$ 0.11 [Mpc$^{-3}$]. The IR LF for DOGs including data obtained from the literature is well fitted by a double-power law. The derived lower limit for the IR LD for our sample is $\rho_{\mathrm{IR}}$ $\sim$ 3.8 $\times$ 10$^7$ [L$_{\odot}$ Mpc$^{-3}$] and its contributions to the total IR LD, IR LD of all ultra-luminous infrared galaxies (ULIRGs), and that of all DOGs are $>$ 3 %, $>$ 9 %, and $>$ 15 %, respectively.

We have developed an electron-tracking Compton camera (ETCC) for use in next-generation MeV gamma ray telescopes. An ETCC consists of a gaseous time projection chamber (TPC) and pixel scintillator arrays (PSAs). Since the TPC measures the three dimensional tracks of Compton-recoil electrons, the ETCC can completely reconstruct the incident gamma rays. Moreover, the ETCC demonstrates efficient background rejection power in Compton-kinematics tests, identifies particle from the energy deposit rate (dE/dX) registered in the TPC, and provides high quality imaging by completely reconstructing the Compton scattering process. We are planning the "Sub-MeV gamma ray Imaging Loaded-on-balloon Experiment" (SMILE) for our proposed all-sky survey satellite. Performance tests of a mid-sized 30 cm-cubic ETCC, constructed for observing the Crab nebula, are ongoing. However, observations at balloon altitudes or satellite orbits are obstructed by radiation background from the atmosphere and the detector itself. The background rejection power was checked using proton accelerator experiments conducted at the Research Center for Nuclear Physics, Osaka University. To create the intense radiation fields encountered in space, which comprise gamma rays, neutrons, protons, and other energetic entities, we irradiated a water target with a 140 MeV proton beam and placed a SMILE-II ETCC near the target. In this situation, the counting rate was five times than that expected at the balloon altitude. Nonetheless, the ETCC stably operated and identified particles sufficiently to obtain a clear gamma ray image of the checking source. Here, we report the performance of our detector and demonstrate its effective background rejection based in electron tracking experiments.

To explore the sub-MeV/MeV gamma-ray window for astronomy, we have developed the Electron-Tracking Compton Camera (ETCC), and carried out the first performance test at room condition using several gamma-ray sources in the sub-MeV energy band. Using a simple track analysis for a quick first test of the performance, the gamma-ray imaging capability was demonstrated by clear images and 5.3 degrees of angular resolution measure (ARM) measured at 662 keV. As the greatest impact of this work, a gamma-ray detection efficiency on the order of $10^{-4}$ was achieved at the sub-MeV gamma-ray band, which is one order of magnitude higher than our previous experiment. This angular resolution and detection efficiency enables us to detect the Crab Nebula at the 5 sigma level with several hours observation at balloon altitude in middle latitude. Furthermore, good consistency of efficiencies between this performance test and simulation including only physical processes has a large importance; it means we achieve nearly 100% detection of Compton recoil electrons. Thus, our estimation of enhancements by upgrades of the detector is more dependable. We are planning to confirm the imaging capability of the ETCC by observation of celestial objects in the SMILE-II (Sub-MeV gamma ray Imaging Loaded-on-balloon Experiment II). The SMILE-II and following SMILE-III project will be an important key of sub-MeV/MeV gamma-ray astronomy.

The stellar properties of about 800 galaxies hosting optically luminous, unobscured quasars at z < 0.6 are analyzed. Deep co-added Sloan Digital Sky Survey (SDSS) images of the quasars on Stripe 82 are decomposed into nucleus and host galaxy using point spread function and Sersic models. The systematic errors in the measured galaxy absolute magnitudes and colors are estimated to be less than 0.5 mag and 0.1 mag, respectively, with simulated quasar images. The effect of quasar light scattered by the interstellar medium is also carefully addressed. The measured quasar-to-galaxy ratio in total flux decreases toward longer wavelengths, from ~8 in the u band to ~1 in the i and z bands. We find that the SDSS quasars are hosted exclusively by massive galaxies (stellar mass Mstar > 10^{10} Msun), which is consistent with previous results for less luminous narrow-line (obscured) active galactic nuclei (AGNs). The quasar hosts are very blue and almost absent on the red sequence, showing stark contrast to the color-magnitude distribution of normal galaxies. The fact that more powerful AGNs reside in galaxies with higher star-formation efficiency may indicate that negative AGN feedback, if it exists, is not concurrent with the most luminous phase of AGNs. We also find positive correlation between the mass of supermassive black holes (SMBHs; Mbh) and host stellar mass, but the Mbh - Mstar relation is offset toward large Mbh or small Mstar compared to the local relation. While this could indicate that SMBHs grow earlier than do their host galaxies, such an argument is not conclusive, as the effect may be dominated by observational biases.

We have conducted B, g, V, and R-band imaging in a 45x40 arcmin^2 field containing part of the high Galactic latitude translucent cloud MBM32, and correlated the intensity of diffuse optical light S_\nu(\lambda) with that of 100 micron emission S_\nu(100um). A chi^2 minimum analysis is applied to fit a linear function to the measured correlation and derive the slope parameter b(\lambda)= \Delta S_\nu(\lambda) / \Delta S_\nu(100um) of the best-fit linear function. Compiling a sample by combining our b(\lambda) and published ones, we show that the b(\lambda) strength varies from cloud to cloud by a factor of 4. Finding that b(\lambda) decreases as S_\nu(100um) increases in the sample, we suggest that a non-linear correlation including a quadratic term of S_\nu(100um)^2 should be fitted to the measured correlation. The variation of optical depth, which is A_V = 0.16 - 2.0 in the sample, can change b(\lambda) by a factor of 2 - 3. There would be some contribution to the large b(\lambda) variation from the forward-scattering characteristic of dust grains which is coupled to the non-isotropic interstellar radiation field (ISRF). Models of the scattering of diffuse Galactic light (DGL) underestimate the b(\lambda) values by a factor of 2. This could be reconciled by deficiency in UV photons in the ISRF or by a moderate increase in dust albedo. Our b(\lambda) spectrum favors a contribution from extended red emission (ERE) to the diffuse optical light; b(\lambda) rises from B to V faster than the models, seems to peak around 6000 \AA, and decreases towards long wavelengths. Such a characteristic is expected from the models in which the DGL is combined with ERE.

Merger processes play an important role in galaxy formation and evolution. To study the influence of merger processes on the evolution of dust properties and cosmic star formation rate, we investigate a local sample of major merger galaxies and a control sample of isolated galaxies using GALEX ultraviolet (UV) and Spitzer infrared (IR) images. Through a statistical study, we find that dust attenuation in merger galaxies is enhanced with respect to isolated galaxies. We find this enhancement is contributed mainly by spiral galaxies in spiral-spiral (S-S) pairs, and increases with the increasing stellar mass of a galaxy. Combining the IR and UV parts of star formation rates (SFRs), we then calculated the total SFRs and specific star formation rates (SSFRs). We find the SSFRs to be enhanced in merger galaxies. This enhancement depends on galaxy stellar mass and the companion's morphology, but depends little on whether the galaxy is a primary or secondary component or on the separation between two components. These results are consistent with a previous study based only on IR images. In addition, we investigate the nuclear contributions to SFRs. SFRs in paired galaxies are more concentrated in the central part of the galaxies than in isolate galaxies. Our studies of dust attenuation show that the nuclear parts of pairs most resemble ULIRGs. Including UV data in the present work not only provides reliable information on dust attenuation, but also refines analyses of SFRs.

The unification of active galactic nuclei (AGN) is a model that has been difficult to test due to the lack of knowledge on the intrinsic luminosities of the objects. We present a test were we probe the model by statistical investigation of the neighbours to AGN at redshifts 0.03 < z < 0.2 within a projected distance of 350 kpc and |\Delta z|<0.001, 0.006, 0.012 and 0.03 between AGN and neighbour. 1658 Type-1 (broad-line) AGN-galaxy pairs and 5698 Type-2 AGN-galaxy pairs with spectroscopic redshifts from the Data Release 7 of Sloan Digital Sky Survey were used together with a complementary set of pairs with photometric redshifts on the neighbour galaxies (13519 Type-1 AGN-galaxy and 58743 Type-2 AGN-galaxy pairs). Morphologies for the AGN host galaxies were derived from the Galaxy Zoo project. Our results suggest that broad-line AGN and narrow-line AGN reside in widely different environments where the neighbours to Type-2 AGN are more star-forming and bluer than those of Type-1 AGN. There is a colour-dependency only detectable in the neighbours with photometric redshifts for the Type-2 AGN. We see that the ratio between Type-1/Type-2 neighbours to Type-2 AGN decreases steadily at short separations with a statistical significance of 4.5 sigma. The lack of change in the morphology of the Type-2 AGN hosts having a close companion (contrary to the case of Type-1 AGN hosts) suggests that the innate state of Type-2 AGN is extremely short-lived and is not preserved in subsequent mergers. Finally, we perform a hypothetical luminosity test to investigate whether a mass bias in our selection could explain the observed differences in our samples. Our conclusion is that AGN unification is consistently not supported by the environment of the two types of AGN, but that an evolutionary connection between them might exist.

HI 1225+01 is an intergalactic gas cloud located on the outskirts of Virgo cluster. Its main components are two large clumps of comparable HI masses (M_HI ~ 10^9 Msun) separated by about 100 kpc. One of the clumps hosts a blue low-surface-brightness galaxy J1227+0136, while the other has no identified stellar emission and is sometimes referred to as a promising candidate of a "dark galaxy", an optically invisible massive intergalactic system. We present a deep optical image covering the whole HI 1225+01 structure for the first time, as well as a collection of archival data from ultraviolet to far-infrared (IR) spectral region of the brightest knot "R1" in J1227+0136. We find that R1 has a young stellar population of age 10-100 Myr and mass ~ 10^6 Msun, near-IR excess brightness which may point to the presence of hot dust with color temperature ~ 600 K, and relatively faint mid- to far-IR fluxes corresponding to the dust mass of up to ~ 100 Msun. Overall, it seems to share the general properties with low-metallicity blue compact dwarf galaxies. On the other hand, no optical counterpart to the other clump is found in our deepest-ever image. Now the limiting surface brightness reaches down to R_AB > 28 mag/arcsec2 for any emission extended over 10" (comparable to R1), which is more than one hundred times fainter than the brightest part of the companion galaxy J1227+0136.

We present the results of our follow-up observation program of gamma-ray sources detected by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. 26 blazars and 39 sources unidentified at other wavelengths were targeted at IRSF 1.4 m telescope equipped with the SIRIUS/SIRPOL imager and polarimeter. H-band magnitudes of the blazars at the epoch of 2010 Dec -- 2011 Feb are presented, which reveal clear flux variation since the Two Micron All Sky Survey observations and can be useful data for variation analyses of these objects in longer periods. We also find that nearly half of the gamma-ray blazars are highly (>10 %) polarized in near-infrared wavelengths. Combining the polarization and variation properties, most (~90 %) of the blazars are clearly distinguished from all other types of objects at high Galactic latitudes. On the other hand, we find only one highly polarized and/or variable object in the fields of unidentified sources. This object is a counterpart of the optical variable source PQV1 J131553.00-073302.0 and the radio source NVSS J131552-073301, and is a promising candidate of new gamma-ray blazars. From the measured polarization and variation statistics, we conclude that most of the Fermi/LAT unidentified sources are not likely similar types of objects to the known gamma-ray blazars.

A comprehensive analysis of extended emission line region (EELR) around quasars is presented. New Subaru/Suprime-Cam observation is combined with literature search, resulting in a compilation of 81 EELR measurements for type-1 and type-2 quasars with associated active galactic nucleus (AGN) and host galaxy properties. It is found that EELR phenomenon shows clear correlation with Eddington ratio, which links EELR to the constituents of the principal component 1 (PC 1), or eigenvector 1, of the AGN emission correlations. We also find that EELR is preferentially associated with gas-rich, massive blue galaxies. It supports the idea that the primary determinant of EELR creation is the gas availability and that the gas may be brought in by galaxy merger triggering the current star formation as well as AGN activity, and also gives an explanation for the fact that most luminous EELR is found around radio-loud sources with low Eddington ratio. By combining all the observations, it is suggested that EELR quasars occupy the massive blue corner of the green valley, the AGN realm, on the galaxy color - stellar mass diagram. Once a galaxy is pushed to this corner, activated AGN would create EELR by the energy injection into the interstellar gas and eventually blow it away, leading to star-formation quenching. The results presented here provide a piece of evidence for the presence of such AGN feedback process, which may be playing a leading role in the co-evolution of galaxies and central super-massive black holes.

We present the results of the near-infrared (IR) spectroscopy of nine IRAS galaxies (NGC 1266, NGC 1320, NGC 2633, NGC 2903, NGC 3034, Mrk 33, NGC 7331, NGC 7625, NGC 7714) with the ISLE imager and spectrograph mounted on the Okayama Astrophysical Observatory 1.88 m telescope. [Fe II] 1.257 um and Pa beta emission lines were observed for the whole sample while H2 2.121 um and Br gamma lines were additionally obtained for two sources, whose flux ratios are used as a diagnostic tool of dominant energy sources of the galaxies. We find that the nucleus of NGC 1266 is most likely a low ionization nuclear emission-line region (LINER), while NGC 2633 and NGC 2903 possibly harbor active galactic nuclei (AGNs). No AGN or LINER signal is found for other objects. In addition, we find the spectral features which is indicative of some unusual phenomena occurring in the galaxies, such as the large [Fe II] line widths compared to the local escape velocity in NGC 1266. The present work shows the potential ability of the ISLE to shed new light on the nature of infrared galaxies, either through a statistical survey of galaxies or an exploration of spectral features found in individual objects.

We present the new constraints on the cosmic optical background (COB) obtained from an analysis of the Pioneer 10/11 Imaging Photopolarimeter (IPP) data. After careful examination of data quality, the usable measurements free from the zodiacal light are integrated into sky maps at the blue (~0.44 um) and red (~0.64 um) bands. Accurate starlight subtraction is achieved by referring to all-sky star catalogs and a Galactic stellar population synthesis model down to 32.0 mag. We find that the residual light is separated into two components: one component shows a clear correlation with thermal 100 um brightness, while another betrays a constant level in the lowest 100 um brightness region. Presence of the second component is significant after all the uncertainties and possible residual light in the Galaxy are taken into account, thus it most likely has the extragalactic origin (i.e., the COB). The derived COB brightness is (1.8 +/- 0.9) x 10^(-9) and (1.2 +/- 0.9) x 10^(-9) erg/s/cm2/sr/A at the blue and red band, respectively, or 7.9 +/- 4.0 and 7.7 +/- 5.8 nW/m2/sr. Based on a comparison with the integrated brightness of galaxies, we conclude that the bulk of the COB is comprised of normal galaxies which have already been resolved by the current deepest observations. There seems to be little room for contributions of other populations including "first stars" at these wavelengths. On the other hand, the first component of the IPP residual light represents the diffuse Galactic light (DGL) - scattered starlight by the interstellar dust. We derive the mean DGL-to-100 um brightness ratios of 2.1 x 10^(-3) and 4.6 x 10^(-3) at the two bands, which are roughly consistent with the previous observations toward denser dust regions. Extended red emission in the diffuse interstellar medium is also confirmed.

We report the discovery of two field brown dwarfs, ULAS J0128-0041 and ULAS J0321+0051, and the rediscovery of ULAS J0226+0051 (IfA 0230-Z1), in the Sloan Digital Sky Survey (SDSS) southern equatorial stripe. They are found in the course of our follow-up observation program of 1 um excess sources in the United Kingdom Infrared Telescope Infrared Deep Sky Survey. The Gemini Multi-Object Spectrographs spectra at red optical wavelengths (6500-10500 A) are presented, which reveal that they are early-T dwarfs. The classification is also supported by their optical to near-infrared colors. It is noted that ULAS J0321+0051 is one of the faintest currently known T dwarfs. The estimated distances to the three objects are 50-110 pc, thus they are among the most distant field T dwarfs known. Dense temporal coverage of the target fields achieved by the SDSS-II Supernova Survey allows us to perform a simple time-series analysis, which leads to the finding of significant proper motions of 150-290 mas/yr or the transverse velocities of 40-100 km/s for ULAS J0128-0041 and ULAS J0226+0051. We also find that there are no detectable, long-term (a-few-year) brightness variations above a few times 0.1 mag for the two brown dwarfs.

We present the discovery and mass measurement of the cold, low-mass planet MOA-2009-BLG-266Lb, made with the gravitational microlensing method. This planet has a mass of m_p = 10.4 +- 1.7 Earth masses and orbits a star of mass M_* = 0.56 +- 0.09 Solar masses at a semi-major axis of a = 3.2 (+1.9 -0.5) AU and an orbital period of P = 7.6 (+7.7 -1.5} yrs. The planet and host star mass measurements are enabled by the measurement of the microlensing parallax effect, which is seen primarily in the light curve distortion due to the orbital motion of the Earth. But, the analysis also demonstrates the capability to measure microlensing parallax with the Deep Impact (or EPOXI) spacecraft in a Heliocentric orbit. The planet mass and orbital distance are similar to predictions for the critical core mass needed to accrete a substantial gaseous envelope, and thus may indicate that this planet is a "failed" gas giant. This and future microlensing detections will test planet formation theory predictions regarding the prevalence and masses of such planets.

We apply the supernova(SN) extinction curves to reproduce the observed properties of SST J1604+4304 which is a young infrared (IR) galaxy at z = 1. The SN extinction curves used in this work were obtained from models of unmixed ejecta of type II supernovae(SNe II) for the Salpeter initial mass function (IMF) with a mass range from 8 to 30 M_sun or 8 to 40 M_sun. The effect of dust distributions on the attenuation of starlight is investigated by performing the chi-square fitting method against various dust distributions. These are the commonly used uniform dust screen, the clumpy dust screen, and the internal dust geometry. We add to these geometries three scattering properties, namely, no-scattering, isotropic scattering, and forward-only scattering. Judging from the chi-square values, we find that the uniform screen models with any scattering property provide good approximations to the real dust geometry. Internal dust is inefficient to attenuate starlight and thus cannot be the dominant source of the extinction. We show that the SN extinction curves reproduce the data of SST J1604+4304 comparable to or better than the Calzetti extinction curve. The Milky Way extinction curve is not in satisfactory agreement with the data unless several dusty clumps are in the line of sight. This trend may be explained by the abundance of SN-origin dust in these galaxies; SN dust is the most abundant in the young IR galaxy at z = 1, abundant in local starbursts, and less abundant in the Galaxy. If dust in SST J1604+4304 is dominated by SN dust, the dust production rate is about 0.1 M_sun per SN.

We investigate FeII emission in Broad Line Region (BLR) of AGNs by analyzing the FeII(UV), FeII(4570) and MgII emission lines in 884 quasars in the Sloan Digital Sky Survey (SDSS) Quasar catalog in a redshift range of 0.727 < z < 0.804. FeII(4570)/FeII(UV) is used to infer the column density of FeII-emitting clouds and explore the excitation mechanism of FeII emission lines. As suggested before in various works, the classical photoionization models fail to account for FeII(4570)/FeII(UV) by a factor of 10, which may suggest anisotropy of UV FeII emission; otherwise, an alternative heating mechanism like shock is working. The column density distribution derived from FeII(4570)/FeII(UV) indicates that radiation pressure plays an important role in BLR gas dynamics. We find a positive correlation between FeII(4570)/FeII(UV) and the Eddington ratio. We also find that almost all FeII-emitting clouds are to be under super-Eddington conditions unless ionizing photon fraction is much smaller than that previously suggested. Finally we propose a physical interpretation of a striking set of correlations between various emission-line properties, known as ``Eigenvector 1''.

We present a clustering analysis of ~60,000 massive (stellar mass Mstar > 10^{11} Msun) galaxies out to z = 1 drawn from 55.2 deg2 of the UKIRT Infrared Deep Sky Survey (UKIDSS) and the Sloan Digital Sky Survey (SDSS) II Supernova Survey. Strong clustering is detected for all the subsamples of massive galaxies characterized by different stellar masses (Mstar = 10^{11.0-11.5} Msun, 10^{11.5-12.0} Msun) or rest-frame colors (blue: U - V < 1.0, red: U - V > 1.0). We find that more mature (more massive or redder) galaxies are more clustered, which implies that more mature galaxies have started stellar-mass assembly earlier within the highly-biased region where the structure formation has also started earlier. By means of halo occupation distribution (HOD) models fitted to the observed angular correlation function, we infer the properties of the underlying host dark halos. We find that the estimated bias factors and host halo masses are systematically larger for galaxies with larger stellar masses, which is consistent with the general agreement that the capability of hosting massive galaxies depends strongly on halo mass. The estimated effective halo masses are ~10^{14} Msun, which gives the stellar-mass to halo-mass ratios of ~0.003. The observed evolution of bias factors indicates rapid evolution of spatial distributions of cold dark matter relative to those traced by the massive galaxies, while the transition of host halo masses might imply that the fractional mass growth rate of halos is less than those of stellar systems. The inferred halo masses and high fractions of central galaxies indicate that the massive galaxies in the current sample are possibly equivalent to central galaxies of galaxy clusters.

We present an analysis of ~60 000 massive (stellar mass M_star > 10^{11} M_sun) galaxies out to z = 1 drawn from 55.2 deg2 of the United Kingdom Infrared Telescope (UKIRT) Infrared Deep Sky Survey (UKIDSS) and the Sloan Digital Sky Survey (SDSS) II Supernova Survey. This is by far the largest survey of massive galaxies with robust mass estimates, based on infrared (K-band) photometry, reaching to the Universe at about half its present age. We find that the most massive (M_star > 10^{11.5} M_sun) galaxies have experienced rapid growth in number since z = 1, while the number densities of the less massive systems show rather mild evolution. Such a hierarchical trend of evolution is consistent with the predictions of the current semi-analytic galaxy formation model based on Lambda CDM theory. While the majority of massive galaxies are red-sequence populations, we find that a considerable fraction of galaxies are blue star-forming galaxies. The blue fraction is smaller in more massive systems and decreases toward the local Universe, leaving the red, most massive galaxies at low redshifts, which would support the idea of active 'bottom-up' formation of these populations during 0 < z < 1.

We demonstrate the capability of AKARI for mapping diffuse far-infrared emission and achieved reliability of all-sky diffuse map. We have conducted an all-sky survey for more than 94 % of the whole sky during cold phase of AKARI observation in 2006 Feb. -- 2007 Aug. The survey in far-infrared waveband covers 50 um -- 180 um with four bands centered at 65 um, 90 um, 140 um, and 160 um and spatial resolution of 3000 -- 4000 (FWHM).This survey has allowed us to make a revolutionary improvement compared to the IRAS survey that was conducted in 1983 in both spatial resolution and sensitivity after more than a quarter of a century. Additionally, it will provide us the first all-sky survey data with high-spatial resolution beyond 100 um. Considering its extreme importance of the AKARI far-infrared diffuse emission map, we are now investigating carefully the quality of the data for possible release of the archival data. Critical subjects in making image of diffuse emission from detected signal are the transient response and long-term stability of the far-infrared detectors. Quantitative evaluation of these characteristics is the key to achieve sensitivity comparable to or better than that for point sources (< 20 -- 95 [MJy/sr]). We describe current activities and progress that are focused on making high quality all-sky survey images of the diffuse far-infrared emission.

We present the detailed optical to far-infrared observations of SST J1604+4304, an ULIRG at z = 1.135. Analyzing the stellar absorption lines, namely, the CaII H & K and Balmer H lines in the optical spectrum, we derive the upper limits of an age for the stellar population. Given this constraint, the minimum {chi}^2 method is used to fit the stellar population models to the observed SED from 0.44 to 5.8um. We find the following properties. The stellar population has an age 40 - 200 Myr with a metallicity 2.5 Z_{sun}. The starlight is reddened by E(B-V) = 0.8. The reddening is caused by the foreground dust screen, indicating that dust is depleted in the starburst site and the starburst site is surrounded by a dust shell. The infrared (8-1000um) luminosity is L_{ir} = 1.78 +/- 0.63 * 10^{12} L_{sun}. This is two times greater than that expected from the observed starlight, suggesting either that 1/2 of the starburst site is completely obscured at UV-optical wavelengths, or that 1/2 of L_{ir} comes from AGN emission. The inferred dust mass is 2.0 +/- 1.0 * 10^8 M_{sun}. This is sufficient to form a shell surrounding the galaxy with an optical depth E(B-V) = 0.8. From our best stellar population model - an instantaneous starburst with an age 40 Myr, we infer the rate of 19 supernovae(SNe) per year. Simply analytical models imply that 2.5 Z_{sun} in stars was reached when the gas mass reduced to 30% of the galaxy mass. The gas metallcity is 4.8 Z_{sun} at this point. The gas-to-dust mass ratio is then 120 +/- 73. The inferred dust production rate is 0.24 +/- 0.12 M_{sun} per SN. If 1/2 of L_{ir} comes from AGN emission, the rate is 0.48 +/- 0.24 M_{sun} per SN. We discuss the evolutionary link of SST J1604+4304 to other galaxy populations in terms of the stellar masses and the galactic winds.

We present rest-frame optical/near-infrared spectra of the gravitationally lense d quasar APM 08279+5255 at $z=3.91$ that has been taken using the Infrared Camera (IRC) onboard the AKARI infrared satellite. The observed continuum consists of two components; a power-law component dominating optical wavelengths which is the direct light from the central source and thermal emission dominating near-infrared wavelengths which is attributed to the emission from hot dust in the circumnuclear region. The thermal emission well represents optically thick emission by hot dust at T ~ 1300K with tau(2micron)> 2 and apparent mass, M(hot) >10 M_sun. Thus, our observations directly detected the optically thick region of hot dust in APM 08279+5255. HI recombination lines of H_alpha(0.656micron), Pa_alpha(1.875micron), and Pa_beta(1.282micron) are clearly detected at 3.2, 6.3, and 9.3 micron. Simulations with the photoionization models suggest that APM 08279+5255 has BLR(Broad Line Region) clouds characterized by log n_{H} ~ 12- 14 for the gas density, log U ~ -2 - -6$ for the ionization parameter, and E(B-V) ~ 0.3-0.6$for the broad line region. Thus, optically thick emission of hot dust support an idea on non-spherical distribution of dust near the central source, consistent with the Active Galactic Nuclei model with the dust torus. The temperature of hot dust and flux ratios of these HI lines are similar to those observed in low-redshift quasars. There are significant time-variations in the HI lines, which are probably caused by variations in the brightness of the central source.

We present spectra of six luminous quasars at z ~ 2, covering rest wavelengths 1600-3200 A. The fluxes of the UV Fe II emission lines and Mg II 2798 doublet, the line widths of Mg II, and the 3000 A luminosity were obtained from the spectra. These quantities were compared with those of low-redshift quasars at z = 0.06 - 0.55 studied by Tsuzuki et al. In a plot of the Fe II(UV)/Mg II flux ratio as a function of the cental black hole mass, Fe II(UV)/Mg II in our z ~ 2 quasars is systematically greater than in the low-redshift quasars. We confermed that luminosity is not responsible for this excess. It is unclear whether this excess is caused by rich Fe abundance at z ~ 2 over low-redshift or by non-abundance effects such as high gas density, strong radiation field, and high microturbulent velocity.

We present the optical to near-infrared (IR) spectrum of the galaxy TSPS J1329-0957, a red and bright member of the class of extremely red objects (EROs) at z = 1.26. This galaxy was found in the course of the Tokyo-Stromlo Photometry Survey (TSPS) which we are conducting in the southern sky. The spectroscopic observations were carried out with the Gemini Multi-Object Spectrograph (GMOS) and the Gemini Near Infra-Red Spectrograph (GNIRS) mounted on the Gemini-South telescope. The wide wavelength coverage of 0.6 - 2.3 um provides useful clues as to the nature of EROs while most published spectra are limited to a narrower spectral range which is dictated by the need for efficient redshift determination in a large survey. We compare our spectrum with several optical composite spectra obtained in recent large surveys, and with stellar population synthesis models. The effectiveness of using near-IR broad-band data, instead of the spectral data, in deriving the galaxy properties are also investigated. We find that TSPS J1329-0957 formed when the universe was 2 - 3 Gyr old, and subsequently evolved passively to become one of the most massive galaxies found in the z = 1 - 2 universe. Its early type and estimated stellar mass of M* = 10^{11.5} Msun clearly point to this galaxy being a direct ancestor of the brightest elliptical and spheroidal galaxies in the local universe.

We have compiled the emission-line fluxes of O I 8446, O I 11287, and the near-IR Ca II triplet (8579) observed in 11 quasars. These lines are considered to emerge from the same gas as do the Fe II lines in the low-ionized portion of the broad emission line region (BELR). The compiled quasars are distributed over wide ranges of redshift (0.06 < z < 1.08) and of luminosity (-29.8 < M_B < -22.1), thus representing a useful sample to investigate the line-emitting gas properties in various quasar environments. The measured line strengths and velocities, as functions of the quasar properties, are analyzed using photoionization model calculations. We found that the flux ratio between Ca II and O I 8446 is hardly dependent on the redshift or luminosity, indicating similar gas density in the emission region from quasar to quasar. On the other hand, a scatter of the O I 11287/8446 ratios appears to imply the diversity of the ionization parameter. These facts invoke a picture of the line-emitting gases in quasars that have similar densities and are located at regions exposed to various ionizing radiation fluxes. The observed O I line widths are found to be remarkably similar over more than 3 orders of magnitude in luminosity, which indicates a kinematically determined location of the emission region and is in clear contrast to the well-studied case of H I lines. We also argue about the dust presence in the emission region since the region is suggested to be located near the dust sublimation point at the outer edge of the BELR.

We present results of the near-infrared (IR) spectroscopy of six quasars whose redshifts range from 0.158 to 1.084. Combined with the satellite ultraviolet data, the relative line strengths of OI 1304, OI 8446, OI 11287, and the near-IR CaII triplet are given. In addition, the corresponding OI line strengths measured in normal Seyfert 1s and narrow-line Seyfert 1s are collected from the literature. These lines are thought to emerge from the same gas as do the FeII lines, so they are good tracers of the FeII emission region within a broad emission line region (BELR) in active galactic nuclei (AGNs). In order to reveal the physical condition within the relevant emission region, we performed photoionized model calculations and compared them to the observations. It suggests that a rather dense gas with density nH ~ 10^(11.5) cm-3 is present at an outer portion of the BELR, illuminated by the ionizing radiation corresponding to an ionization parameter U ~ 10^(-2.5) and is primarily responsible for the observed OI, CaII, and FeII lines, based on the resemblance of their profiles. The three OI lines are proven to be formed through Ly beta fluorescence and collisional excitation. We also show that the lambda1304 bump typically observed in AGN spectra consists of the comparable contributions of OI and SiII multiplets, and we discuss the origin of such a strong SiII emission. The results are interpreted in the context of the locally optimally emitting cloud (LOC) scenario to find the plausible gas distribution within the BELR as a function of distance from the central source and density.

We present infrared spectroscopy of the recurrent nova RS Ophiuchi, obtained 11.81, 20.75 and 55.71 days following its 2006 eruption. The spectra are dominated by hydrogen recombination lines, together with HeI, OI and OII lines; the electron temperature of ~10^4 K implied by the recombination spectrum suggests that we are seeing primarily the wind of the red giant, ionized by the ultraviolet flash when RS Oph erupted. However, strong coronal emission lines (i.e. emission from fine structure transitions in ions having high ionization potential) are present in the last spectrum. These imply a temperature of 930000K for the coronal gas; this is in line with x-ray observations of the 2006 eruption. The emission line widths decrease with time in a way that is consistent with the shock model for the x-ray emission.

We present the spectra of 14 quasars with a wide coverage of rest wavelengths from 1000 to 7300 A. The redshift ranges from z = 0.061 to 0.555 and the luminosity from M_{B} = -22.69 to -26.32. We describe the procedure of generating the template spectrum of Fe II line emission from the spectrum of a narrow-line Seyfert 1 galaxy I Zw 1 that covers two wavelength regions of 2200-3500 A and 4200-5600 A. Our template Fe II spectrum is semi-empirical in the sense that the synthetic spectrum calculated with the CLOUDY photoionization code is used to separate the Fe II emission from the Mg II line. The procedure of measuring the strengths of Fe II emission lines is twofold; (1) subtracting the continuum components by fitting models of the power-law and Balmer continua in the continuum windows which are relatively free from line emissions, and (2) fitting models of the Fe II emission based on the Fe II template to the continuum-subtracted spectra. From 14 quasars, we obtained the Fe II fluxes in five wavelength bands, the total flux of Balmer continuum, and the fluxes of Mg II, Halpha, and other emission lines, together with the full width at half maxima (FWHMs) of these lines. Regression analysis was performed by assuming a linear relation between any two of these quantities. Eight correlations were found with a confidence level higher than 99%. The fact that six of these eight are related to FWHM or M_{BH} may imply that M_{BH} is a fundamental quantity that controls Gamma or the spectral energy distribution (SED) of the incident continuum, which in turn controls the Fe II emission. Furthermore, it is worthy of noting that Fe II(O1)/Fe II(U1) is found to tightly correlate with Fe II(O1)/Mg II, but not with Fe II(U1)/Mg II.

By observing the near-infrared spectrum of the quasar PG 1116+215 at z = 0.176 and combining with the HST/FOS spectrum, we obtained the relative strengths of three permitted OI lines ({lambda}1304, {lambda}8446, and {lambda}11287) in a quasar for the first time. The photon flux ratios of the OI lines of the quasar were compared with those previously measured in a Seyfert 1 and six narrow-line Seyfert 1s. No significant differences were found in the OI line flux ratios between the quasar and the other Seyferts, suggesting that the gas density in the OI and FeII line-emitting regions in the quasar is of the same order as those in low-luminosity AGNs. It was also found that the line width of OI {lambda}11287 is significantly narrower than that of Ly{alpha}, which is consistent with OI and FeII emission occurring in the partly ionized regions at the outermost portion of the broad-line region where velocities are small.

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

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

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