Search

We present the optical and infrared properties of 39 extremely radio-loud galaxies discovered by cross-matching the Subaru/Hyper Suprime-Cam (HSC) deep optical imaging survey and VLA/FIRST 1.4 GHz radio survey. The recent Subaru/HSC strategic survey revealed optically-faint radio galaxies (RG) down to $g_\mathrm{AB} \sim 26$, opening a new parameter space of extremely radio-loud galaxies (ERGs) with radio-loudness parameter of $\log \mathcal{R}_\mathrm{rest} = \log (f_{1.4 \mathrm{GHz,rest}}/f_{g,\mathrm{rest}}) >4$. Because of their optical faintness and small number density of $\sim1~$deg$^{-2}$, such ERGs were difficult to find in the previous wide but shallow, or deep but small area optical surveys. ERGs show intriguing properties that are different from the conventional RGs: (1) most ERGs reside above or on the star-forming main-sequence, and some of them might be low-mass galaxies with $\log (M_\star/M_\odot) < 10$. (2) ERGs exhibit a high specific black hole accretion rate, reaching the order of the Eddington limit. The intrinsic radio-loudness ($\mathcal{R}_\mathrm{int}$), defined by the ratio of jet power over bolometric radiation luminosity, is one order of magnitude higher than that of radio quasars. This suggests that ERGs harbor a unique type of active galactic nuclei (AGN) that show both powerful radiations and jets. Therefore, ERGs are prominent candidates of very rapidly growing black holes reaching Eddington-limited accretion just before the onset of intensive AGN feedback., 20 pages, 10 figures, accepted for publication in ApJ

Studies of white dwarfs have greatly benefited from time-domain surveys and subsequent follow-up observations. However, with the avalanche of alerts delivered by ZTF and LSST and the limited resources for follow-up, we will need brokers to select intriguing alerts that warrant follow-up in a timely manner. At the University of Arizona and NSF’s OIR Lab, we are developing the Arizona-NOAO Temporal Analysis and Response to Events System, to hunt for the rarest of the rare events in the time-domain. In this work, we provide an overview of the ANTARES system, how we use ZTF as a training set, and the way forwards to LSST.

We present ALMA [CII] line and far-infrared (FIR) continuum observations of seven z > 6 low-luminosity quasars (M1450 > −25 mag) discovered by our on-going Subaru Hyper Suprime-Cam survey. The [CII] line was detected in all targets with luminosities of ∼(2−10) × 108 L⊙, about one order of magnitude smaller than optically luminous quasars. Also found was a wide scatter of FIR continuum luminosity, ranging from LFIR < 1011L⊙ to ∼2 × 1012L⊙. With the [CII]-based dynamical mass, we suggest that a significant fraction of low-luminosity quasars are located on or even below the local Magorrian relation, particularly at the massive end of the galaxy mass distribution. This is a clear contrast to the previous finding that luminous quasars tend to have overmassive black holes relative to the relation. Our result is expected to show a less-biased nature of the early co-evolution of black holes and their host galaxies.

Izumi, T., et al., We present Atacama Large Millimeter/submillimeter Array [C ii] 158 μm line and underlying far-infrared (FIR) continuum emission observations (0.″70 × 0.″56 resolution) toward HSC J124353.93+010038.5 (J1243+0100) at z = 7.07, the only low-luminosity (M 1450 > -25 mag) quasar currently known at z > 7. The FIR continuum is bright (1.52 mJy) and resolved with a total luminosity of L FIR = 3.5 × 1012 L o˙. The spatially extended component is responsible for ∼40% of the emission. The area-integrated [C ii] spectrum shows a broad wing (FWHM = 997 km s-1, L [C ii] = 1.2 × 109 L o˙), as well as a bright core (FWHM = 235 km s-1, L [C ii] = 1.9 × 109 L o˙). This wing is the first detection of a galactic-scale quasar-driven outflow (atomic outflow rate >447 M o˙ yr-1) at z > 7. The estimated large mass-loading factor of the total outflow (e.g., ⪆9 relative to the [C ii]-based star formation rate) suggests that this outflow will soon quench the star formation of the host. The core gas dynamics are governed by rotation, with a rotation curve suggestive of a compact bulge (∼3.3 × 1010 M o˙), although it is not yet spatially resolved. Finally, we found that J1243+0100 has a black hole mass-to-dynamical mass (and -to-bulge mass) ratio of ∼0.4% (∼1%), consistent with the local value within the uncertainties. Our results therefore suggest that the black hole-host coevolution relation is already in place at z ∼ 7 for this object., Funding was contributed by the FIRST program from the Japanese Cabinet Office, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), the Japan Society for the Promotion of Science (JSPS), the Japan Science and Technology Agency (JST), the Toray Science Foundation, NAOJ, Kavli IPMU, KEK, ASIAA, and Princeton University. T.H. was supported by th

Some supernovae such as pair-instability supernovae are predicted to have the duration of more than a year in the observer frame. To constrain the rates of supernovae lasting for more than a year, we conducted a long-term deep transient survey using Hyper Suprime-Cam (HSC) on the 8.2m Subaru telescope. HSC is a wide-field (a 1.75 deg2 field-of-view) camera and it can efficiently conduct transient surveys. We observed the same 1.75 deg2 field repeatedly using the g, r, i, and z band filters with the typical depth of 26 mag for 4 seasons (from late 2016 to early 2020). Using these data, we searched for transients lasting for more than a year. Two supernovae were detected in 2 continuous seasons, one supernova was detected in 3 continuous seasons, but no transients lasted for all 4 seasons searched. The discovery rate of supernovae lasting for more than a year with the typical limiting magnitudes of 26 mag is constrained to be 1.4^{+1.3}_{-0.7}(stat.)^{+0.2}_{-0.3}(sys.) events deg-2 yr-1. All the long-lasting supernovae we found are likely Type IIn supernovae and our results indicate that about 40% of Type IIn supernovae have long-lasting light curves. No plausible pair-instability supernova candidates lasting for more than a year are discovered. By comparing the survey results and survey simulations, we constrain the luminous pair-instability supernova rate up to z ~ 3 should be of the order of 100 Gpc-3 yr-1 at most, which is 0.01 - 0.1 per cent of the core-collapse supernova rate., 18 pages, 9 figures, 6 tables, accepted by The Astrophysical Journal

We describe the Arizona-NOIRLab Temporal Analysis and Response to Events System (ANTARES), a software instrument designed to process large-scale streams of astronomical time-domain alerts. With the advent of large-format CCDs on wide-field imaging telescopes, time-domain surveys now routinely discover tens of thousands of new events each night, more than can be evaluated by astronomers alone. The ANTARES event broker will process alerts, annotating them with catalog associations and filtering them to distinguish customizable subsets of events. We describe the data model of the system, the overall architecture, annotation, implementation of filters, system outputs, provenance tracking, system performance, and the user interface., 24 Pages, 8 figures, Accepted by AJ

Galaxies in dense environments, such as groups and clusters, experience various processes by which galaxies gain and lose gas. Using data from the SDSS-IV MaNGA survey, we previously reported the discovery of a giant (6 -- 8 kpc in diameter) H$\alpha$ blob, Totoro, about 8 kpc away from a pair of galaxies (Satsuki and Mei) residing in a galaxy group which is experiencing a group-group merger. Here, we combine interferometric $^{12}$CO(1--0) molecular gas data, new wide-field H$\alpha$, $u$-band data, and published X-ray data to determine the origin of the blob. Several scenarios are discussed to account for its multi-wavelength properties, including (1) H$\alpha$ gas being stripped from galaxy Satsuki by ram-pressure; (2) a separated low-surface-brightness galaxy; (3) gas being ejected or ionized by an active galactic nucleus (AGN); and (4) a cooling intra-group medium (IGM). Scenarios (1) and (2) are less favored by the present data. Scenario (3) is also less likely as there is no evidence for an active ongoing AGN in the host galaxy. We find that the CO (cold) and H$\alpha$ (warm) gas coexist with X-ray (hot) structures; moreover, the derived cooling time is within the regime where molecular and H$\alpha$ gas are expected. The coexistence of gas with different temperatures also agrees with that of cooling gas in other systems. Our multi-wavelength results strongly suggest that the CO and H$\alpha$ gas are the product of cooling from the IGM at its current location, i.e., cooling has occurred, and may be ongoing, well outside the host-galaxy core., Comment: 23 pages, 10 figures, Accepted for publication in ApJ

The ongoing Zwicky Transient Facility (ZTF) survey is generating a massive alert rate from a variety of optical transients and variable stars, which are being filtered down to subsets meeting user-specified criteria by broker systems such as ANTARES. In a beta implementation of the algorithm of Soraisam et al. (2020) on ANTARES, we flagged AT 2020iko from the ZTF real-time alert stream as an anomalous source. This source is located close to a red extended SDSS source. In the first few epochs of detection, it exhibited a V-shaped brightness profile, preceded by non-detections both in ZTF and in ASASSN extending to 2014. Its full light curve shows a precursor event, followed by a main superoutburst and at least two rebrightenings. A low-resolution spectrum of this source points to a dwarf nova (DN) nature. Although some of the features of AT 2020iko indicate an SU UMa-type DN, its large amplitude, presence of rebrightenings, and inferred supercycle period of > 6 yr are in favor of AT 2020iko being a new WZ Sge-type dwarf nova candidate, a subset of rare DNe consisting of extreme mass-ratio (< 0.1) binaries with orbital period around the period minimum. AT 2020iko's precursor event brightened by 6.5 mag, while its decay spanned 3-5 mag. We speculate this superoutburst is associated with a less expanded accretion disk than in typical superoutbursts in WZ Sge systems, with the large depth of the precursor decay implying an extremely small mass-ratio. To the best of our knowledge, such a precursor event has not been recorded for any DN. This result serves to demonstrate the efficacy of our real-time anomaly search algorithm., 11 pages, 4 figures; accepted to AJ

The Faint Object Camera and Spectrograph (FOCAS) is an optical imaging and spectroscopy instrument for the Subaru Telescope. It has been a workhorse instrument since the first-light phase of the telescope. We describe an integral field unit (IFU) that has recently been installed in FOCAS. The IFU utilizes an image slicer that divides a ${13{^{\prime \prime }_{.}}5}$ × ${10{^{\prime \prime }_{.}}0}$ field of view into 23 stripes, with a width of ${0{^{\prime \prime }_{.}}435}$. A sky spectrum separated from an object by approximately ${5{^{\prime }_{.}}2}$ can be obtained at the same time as an object spectrum. Test observations confirmed that the image quality of the IFU does not degrade the ${0{^{\prime \prime }_{.}}435}$ sampling, and that slice width and length are consistent with the design. Highly reflective multilayer dielectric coatings were coated on all the mirrors in the IFU, thereby offering a high mean IFU throughput of ∼85% over the field. However, the outer part of the field showed throughput degradation, which was mainly caused by vignetting as a result of misalignment. The flat-fielding accuracy was degraded by the vignetting, with the variation depending on the direction of the telescope.

We have conducted a comprehensive survey of emission-line galaxies at $z\lesssim1.6$ based on narrowband (NB) imaging data taken with Hyper Suprime-Cam (HSC) on the Subaru telescope. In this paper, we update the catalogs of H$\alpha$, [OIII], and [OII] emission-line galaxies using the data from the second Public Data Release (PDR2) of Subaru Strategic Program (SSP) of HSC and Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS) survey along with the spectroscopic redshifts for 2,019 emission-line galaxies selected with the PDR1 data. The wider effective coverage of NB816 and NB921, 16.3 deg$^2$ and 16.9 deg$^2$ respectively, are available in the Deep and UltraDeep layers of HSC-SSP from the PDR2. The CHORUS survey provides us with data with additional three NBs (NB527, NB718, and NB973) in the COSMOS field in the UltraDeep layer (1.37 deg$^2$). The five NB datasets allow us to investigate the star-forming galaxies presenting emission-lines at 14 specific redshifts ranging from $z\sim1.6$ down to $z\sim0.05$. We revisit the distribution of large-scale structures and luminosity functions (LFs) for the emission-line galaxies with the large samples of 75,377 emission-line galaxies selected. The redshift revolution of LFs shows that the star formation rate densities (SFRDs) decreases monotonically from $z\sim1.6$, which is consistent with the cosmic SFRD ever known. Our samples of emission-line galaxies covering a sufficiently large survey volume are useful to investigate the evolution of star-forming galaxies since the cosmic noon in a wide range of environments including galaxy clusters, filaments, and voids., Comment: 32 pages, 23 figures, 6 tables, accepted for publication in PASJ. The catalog of emission-line galaxies from HSC-SSP PDR2 data will be available soon at https://hsc-release.mtk.nao.ac.jp/, and the catalog from CHORUS data will be available after the CHORUS data are publicly released

arXiv:2007.08685v2, We present the first discovery of dust-reddened quasars (red quasars) in the high-z universe (z > 5.6). This is a result from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which is based on the sensitive multi-band optical imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We identified four red quasar candidates from the 93 spectroscopically confirmed high-z quasars in the SHELLQs sample, based on detections in the Wide-field Infrared Survey Explorer (WISE) data at 3.4 and 4.6 μm (rest-frame ∼5000–6500 Å). The amount of dust reddening was estimated with spectral energy distribution (SED) fits over optical and mid-infrared wavelengths. Two of the four candidates were found to be red quasars with dust reddening of E(B − V) > 0.1. The remaining SHELLQs quasars without individual WISE detections are significantly fainter in the WISE bands and bluer than the red quasars, although we did detect them in the W1 band in a stacked image. We also conducted the same SED fits for high-z optically-luminous quasars, but no red quasar was found. This demonstrates the power of Subaru HSC to discover high-z red quasars, which are fainter than the limiting magnitudes of past surveys in the rest-frame ultraviolet, due to dust extinction., We are grateful to the referee for his/her useful comments to improve this paper. YM was supported by the Japan Society for the Promotion of Science (JSPS) KAKENHI grant No. JP17H04830 and the Mitsubishi Foundation grant No. 30140. TI acknowledges supports from the JSPS grant No. JP17K14247. The Hyper Suprime-Cam (HSC) collaboration includes the astronomical communities of Japan and Taiwan, and Princeton University. The HSC instrumentation and software were developed by the National Astronomical Observatory of Japan (NAOJ), the Kavli Institute for the Physics and Mathematics of the Universe (Kavli IPMU), the University of Tokyo, the High Energy Accelerator Research Organization (KEK), the Academia Sinica Institute for Astronomy and Astrophysics in Taiwan (ASIAA), and Princeton University. Funding was contributed by the FIRST program from Japanese Cabinet Office, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), the Japan Society for the Promotion of Science (JSPS), Japan Science and Technology Agency (JST), the Toray Science Foundation, NAOJ, Kavli IPMU, KEK, ASIAA, and Princeton University. The Pan-STARRS1 Surveys (PS1) have been made possible through contributions of the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the Max-Planck Society and its participating institutes, the Max Planck Institute for Astronomy, Heidelberg and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hopkins University, Durham University, the University of Edinburgh, Queen’s University Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incorporated, the National Central University of Taiwan, the Space Telescope Science Institute, the National Aeronautics and Space Administration under Grant No. NNX08AR22G issued through the Planetary Science Division of the NASA Science Mission Directorate, the National Science Foundation under Grant No. AST-1238877, the University of Maryland, and Eotvos Lorand University (ELTE). This paper makes use of software developed for the Large Synoptic Survey Telescope. We thank the LSST Project for making their code available as free software at 〈http://dm.lsst.org. This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration.

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., Comment: 10 pages, 5 figures, accepted for publication in AJ

Recent progress in submillimetre surveys by single-dish telescopes allows us to further challenge the consistency between cosmological simulations and observations. In particular, we compare our simulations that include dust formation and destruction with the recent SCUBA-2 surveys (`STUDIES') by putting emphases on basic observational properties of dust emission such as dust temperature, size of infrared (IR)-emitting region, IR luminosity function and IRX--$\beta$ relation. After confirming that our models reproduce the local galaxy properties, we examine the STUDIES sample at $z\approx 1-4$, finding that the simulation reproduces the aforementioned quantities except for the $z\gtrsim 2$ IR luminosity function at the massive end ($\sim 10^{13}$ L$_{\odot}$). This means that the current simulation correctly reproduces the overall scaling between the size and luminosity (or star formation rate) of dusty region, but lacks extreme starburst phenomena at $z\gtrsim 2$. We also discuss extinction curves and possible AGN contribution., Comment: 13 pages, 10 figures, published in MNRAS

We report the largest sample of candidate strong gravitational lenses belonging to the Survey of Gravitationally-lensed Objects in HSC Imaging for group-to-cluster scale (SuGOHI-c) systems. These candidates are compiled from the S18A data release of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) Survey. We visually inspect $\sim39,500$ galaxy clusters, selected from several catalogs, overlapping with the Wide, Deep, and UltraDeep fields, spanning the cluster redshift range $0.05, Comment: Accepted for publication in MNRAS, 19 pages, 11 figures

We present a CO(3-2) survey of selected regions in the M31 disc as part of the JCMT large programme, HARP and SCUBA-2 High-Resolution Terahertz Andromeda Galaxy Survey (HASHTAG). The 12 CO(3-2) fields in this survey cover a total area of 60 square arcminutes, spanning a deprojected radial range of 2 - 14 kpc across the M31 disc. Combining these observations with existing IRAM 30m CO(1-0) observations and JCMT CO(3-2) maps of the nuclear region of M31, as well as dust temperature and star formation rate surface density maps, we are able to explore the radial distribution of the CO(3-2)/CO(1-0) integrated intensity ratio (R31) and its relationship with dust temperature and star formation. We find that the value of R31 between 2 - 9 kpc galactocentric radius is 0.14, significantly lower than what is seen in the nuclear ring at ~1 kpc (R31 ~ 0.8), only to rise again to 0.27 for the fields centred on the 10 kpc star forming ring. We also found that R31 is positively correlated with dust temperature, with Spearman's rank correlation coefficient $\rho$ = 0.55. The correlation between star formation rate surface density and CO(3--2) intensity is much stronger than with CO(1-0), with $\rho$ = 0.54 compared to -0.05, suggesting that the CO(3-2) line traces warmer and denser star forming gas better. We also find that R31 correlates well with star formation rate surface density, with $\rho$ = 0.69., Comment: 21 pages, 10 figures, accepted by MNRAS for publication

We present Keck/LRIS follow-up spectroscopy for 13 photometric candidates of extremely metal-poor galaxies (EMPGs) selected by a machine-learning technique applied to the deep (∼26 AB mag) optical and wide-area (∼500 deg2) Subaru imaging data in the EMPRESS survey. Nine out of the 13 candidates are EMPGs with an oxygen abundance (O/H) less than ∼10% solar value (O/H)⊙, and four sources are contaminants of moderately metal-rich galaxies or no emission-line objects. Notably, two out of the nine EMPGs have extremely low stellar masses and oxygen abundances of 5 × 104–7 × 105 M ⊙ and 2%–3% (O/H)⊙, respectively. With a sample of five EMPGs with (Fe/O) measurements, two (three) of which are taken from this study (the literature), we confirm that two EMPGs with the lowest (O/H) ratios of ∼2% (O/H)⊙ show high (Fe/O) ratios of ∼0.1, close to the solar abundance ratio. Comparing galaxy chemical enrichment models, we find that the two EMPGs cannot be explained by a scenario of metal-poor gas accretion/episodic star formation history due to their low (N/O) ratios. We conclude that the two EMPGs can be reproduced by the inclusion of bright hypernovae and/or hypothetical pair-instability supernovae (SNe) preferentially produced in a metal-poor environment. This conclusion implies that primordial galaxies at z ∼ 10 could have a high abundance of Fe that did not originate from Type Ia SNe with delays and that Fe may not serve as a cosmic clock for primordial galaxies.

We report the serendipitous discovery of HSC J0904$-$0102, a quadruply-lensed Lyman break galaxy (LBG) in the Survey of Gravitationally-lensed Objects in Hyper Suprime-Cam Imaging (SuGOHI). Owing to its point-like appearance, the source was thought to be a lensed active galactic nucleus. We obtained follow-up spectroscopic data with the Gemini Multi-Object Spectrographs on the Gemini South Telescope, which confirmed this to be a lens system. The deflecting foreground galaxy is a typical early-type galaxy at a high redshift of $z_{\ell} = 0.957$ with stellar velocity dispersion $\sigma_v=259\pm56$ km~s$^{-1}$. The lensed source is identified as an LBG at $z_{\rm s} = 3.403$, based on the sharp drop bluewards of Ly$\alpha$ and other absorption features. A simple lens mass model for the system, assuming a singular isothermal ellipsoid, yields an Einstein radius of $\theta_{\rm Ein} = 1. 23^{\prime\prime}$ and a total mass within the Einstein radius of $M_{\rm Ein} = (5.55\pm 0.24) \times 10^{11}M_{\odot}$ corresponding to a velocity dispersion of $\sigma_{\rm SIE}= 283\pm 3$ km~s$^{-1}$, which is in good agreement with the value derived spectroscopically. The most isolated lensed LBG image has a magnification of $\sim 6.5$. In comparison with other lensed LBGs and typical $z\sim4$ LBG populations, HSC J0904$-$0102 is unusually compact, an outlier at $>2\sigma$ confidence. Together with a previously discovered SuGOHI lens, HSC J1152$+$0047, that is similarly compact, we believe that the HSC Survey is extending LBG studies down to smaller galaxy sizes., Comment: Accepted for publication on MNRAS, 10 pages, 6 figures

We report the rest-frame ultraviolet luminosity function of $g$-dropout galaxies in 177 protocluster candidates (PC UVLF) at $z\sim4$ selected in the Hyper Suprime-Cam Subaru Strategic Program. Comparing with the UVLF of field galaxies at the same redshift, we find that the PC UVLF shows a significant excess towards the bright-end. This excess can not be explained by the contribution of only active galactic nuclei, and we also find that this is more significant in higher dense regions. Assuming that all protocluster members are located on the star formation main sequence, the PC UVLF can be converted into a stellar mass function. Consequently, our protocluster members are inferred to have a 2.8 times more massive characteristic stellar mass than that of the field Lyman break galaxies at the same redshift. This study, for the first time, clearly shows that the enhancement in star formation or stellar mass in overdense regions can generally be seen as early as at $z\sim4$. We also estimate the star formation rate density (SFRD) in protocluster regions as $\simeq 6-20\%$ of the cosmic SFRD, based on the measured PC UVLF after correcting for the selection incompleteness in our protocluster sample. This high value suggests that protoclusters make a non-negligible contribution to the cosmic SFRD at $z\sim4$, as previously suggested by simulations. Our results suggest that protoclusters are essential components for the galaxy evolution at $z\sim4$., Comment: 22 pages, 16 figures, Accepted for publication in ApJ

Strong lenses are extremely useful probes of the distribution of matter on galaxy and cluster scales at cosmological distances, but are rare and difficult to find. The number of currently known lenses is on the order of 1,000. We wish to use crowdsourcing to carry out a lens search targeting massive galaxies selected from over 442 square degrees of photometric data from the Hyper Suprime-Cam (HSC) survey. We selected a sample of $\sim300,000$ galaxies with photometric redshifts in the range $0.2 < z_{phot} < 1.2$ and photometrically inferred stellar masses $\log{M_*} > 11.2$. We crowdsourced lens finding on this sample of galaxies on the Zooniverse platform, as part of the Space Warps project. The sample was complemented by a large set of simulated lenses and visually selected non-lenses, for training purposes. Nearly 6,000 citizen volunteers participated in the experiment. In parallel, we used YattaLens, an automated lens finding algorithm, to look for lenses in the same sample of galaxies. Based on a statistical analysis of classification data from the volunteers, we selected a sample of the most promising $\sim1,500$ candidates which we then visually inspected: half of them turned out to be possible (grade C) lenses or better. Including lenses found by YattaLens or serendipitously noticed in the discussion section of the Space Warps website, we were able to find 14 definite lenses, 129 probable lenses and 581 possible lenses. YattaLens found half the number of lenses discovered via crowdsourcing. Crowdsourcing is able to produce samples of lens candidates with high completeness and purity, compared to currently available automated algorithms. A hybrid approach, in which the visual inspection of samples of lens candidates pre-selected by discovery algorithms and/or coupled to machine learning is crowdsourced, will be a viable option for lens finding in the 2020s., Published version

We analyse the physical properties of 121 SNR $\geq$ 5 sub-millimetre galaxies (SMGs) from the STUDIES 450-$\mu$m survey. We model their UV-to-radio spectral energy distributions using MAGPHYS+photo-$z$ and compare the results to similar modelling of 850-$\mu$m-selected SMG sample from AS2UDS, to understand the fundamental physical differences between the two populations at the observed depths. The redshift distribution of the 450-$\mu$m sample has a median of $z$ = 1.85 $\pm$ 0.12 and can be described by strong evolution of the far-infrared luminosity function. The fainter 450-$\mu$m sample has $\sim$14 times higher space density than the brighter 850-$\mu$m sample at $z$ $\lesssim$2, and a comparable space density at $z$ = 2-3, before rapidly declining, suggesting LIRGs are the main obscured population at $z$ $\sim$ 1-2, while ULIRGs dominate at higher redshifts. We construct rest-frame $\sim$ 180-$\mu$m-selected and dust-mass-matched samples at $z$ = 1-2 and $z$ = 3-4 from the 450-$\mu$m and 850-$\mu$m samples, respectively, to probe the evolution of a uniform sample of galaxies spanning the cosmic noon era. Using far-infrared luminosity, dust masses and an optically-thick dust model, we suggest that higher-redshift sources have higher dust densities due to inferred dust continuum sizes which are roughly half of those for the lower-redshift population at a given dust mass, leading to higher dust attenuation. We track the evolution in the cosmic dust mass density and suggest that the dust content of galaxies is governed by a combination of both the variation of gas content and dust destruction timescale., Comment: 20 pages, 7 figures, MNRAS accepted

2I/Borisov is the first-ever observed interstellar comet (and the second detected interstellar object). It was discovered on 30 August 2019 and has a heliocentric orbital eccentricity of ~ 3.35, corresponding to a hyperbolic orbit that is unbound to the Sun. Given that it is an interstellar object, it is of interest to compare its properties -- such as composition and activity -- with the comets in our Solar System. This study reports low-resolution optical spectra of 2I/Borisov. The spectra were obtained by the MDM observatory Hiltner 2.4m telescope/Ohio State Multi-Object Spectrograph (on 1 and 5 November 2019). The wavelength coverage spanned from 3700A to 9200A. The dust continuum reflectance spectra of 2I/Borisov show that the spectral slope is stepper in the blue end of the spectrum (compared to the red). The spectra of 2I/Borisov clearly show CN emission at 3880A, as well as C2 emission at both 4750A and 5150A. Using a Haser model to covert the observed fluxes into estimates for the molecular production rates, we find Q(CN) = 2.4 +/- 0.2 x 10^24 s^-1, and Q(C2) = 5.5 +/- 0.4 x 10^23 s^-1 at the heliocentric distance of 2.145 au. Our Q(CN) estimate is consistent with contemporaneous observations, and the Q(C2) estimate is generally below the upper limits of previous studies. We derived the ratio Q(C2)/Q(CN) = 0.2 +/- 0.1, which indicates that 2I/Borisov is depleted in carbon chain species, but is not empty. This feature is not rare for the comets in our Solar System, especially in the class of Jupiter Family Comets., 8 pages, 3 figures, accepted for publication in Astrophysical Journal Letters

We construct a SCUBA-2 450-$\mu$m map in the COSMOS field that covers an area of 300 arcmin$^{2}$ and reaches a 1$\sigma$ noise level of 0.65 mJy in the deepest region. We extract 256 sources detected at 450 $\mu$m with signal-to-noise ratio $>$ 4.0 and analyze the physical properties of their multi-wavelength counterparts. We find that most of the sources are at $z\lesssim3$, with a median of $z = 1.79^{+0.03}_{-0.15}$. About $35^{+32}_{-25}$% of our sources are classified as starburst galaxies based on their total star-formation rates (SFRs) and stellar masses ($M_{\ast}$). By fitting the far-infrared spectral energy distributions, we find that our 450-$\mu$m-selected sample has a wide range of dust temperatures (20 K $ \lesssim T_{\rm d} \lesssim$ 60 K), with a median of ${T}_{\rm d} = 38.3^{+0.4}_{-0.9}$ K. We do not find a redshift evolution in dust temperature for sources with $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ at $z $10^{12}$ $\rm L_\odot$ population to the SFR density rises dramatically from $z$ = 0 to 2 ($\propto$ ($1+z$)$^{3.9\pm1.1}$) and dominates the total SFR density at $z \gtrsim 2$., Comment: Published in the ApJ

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$��$ 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., 19 pages, 11 figures, 4 tables. Accepted for publication in Publications of the Astronomical Society of Japan (PASJ)

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}$., Comment: 17 pages, 11 figures, Accepted for Publication in ApJ

著者人数: 69名 (所属. 宇宙航空研究開発機構宇宙科学研究所(JAXA)(ISAS): 鈴木, 大介), Number of authors: 69 (Affiliation. Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency(JAXA)(ISAS): Suzuki, Daisuke), Accepted: 2019-09-25, 資料番号: SA1190157000

We present results of the cross-correlation Ly$\alpha$ intensity mapping with Subaru/Hyper Suprime-Cam (HSC) ultra-deep narrowband images and Ly$\alpha$ 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, we further conduct extensive analysis evaluating systematics of large-scale point-spread-function wings, sky subtractions, and unknown errors based on physically uncorrelated signals and sources found in real HSC images and object catalogs, respectively. Removing the systematics, we carefully calculate cross-correlations between Ly$\alpha$ intensity of the narrowband images and the LAEs. We tentatively identify very diffuse Ly$\alpha$ emission with the $\simeq 3\sigma$ ($\simeq 2\sigma$) significance at $\gtrsim$ 100 comoving kpc (ckpc) far from the LAEs at $z=5.7$ ($6.6$), around and probably even beyond a virial radius of star-forming galaxies with $M_\mathrm{h}\sim10^{11}M_\odot$. The diffuse Ly$\alpha$ 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 ($, Comment: 11 pages, 6 figures, accepted for publication in ApJ

We present deep near-infrared spectroscopy of six quasars at 6.15.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., Comment: 20 pages, 10 figures, accepted for publication in ApJ

We present element abundance ratios and ionizing radiation of local young low-mass (∼106 M ⊙) extremely metal-poor galaxies (EMPGs) with a 2% solar oxygen abundance (O/H)⊙ and a high specific star formation rate (sSFR ∼ 300 Gyr−1) and other (extremely) metal-poor galaxies, which are compiled from Extremely Metal-Poor Representatives Explored by the Subaru Survey (EMPRESS) and the literature. Weak emission lines such as [Fe iii] λ4658 and He ii λ4686 are detected in very deep optical spectra of the EMPGs taken with 8 m class telescopes, including Keck and Subaru, enabling us to derive element abundance ratios with photoionization models. We find that neon-to-oxygen and argon-to-oxygen ratios are comparable to those of known local dwarf galaxies and that the nitrogen-to-oxygen abundance ratios (N/O) are lower than 20% (N/O)⊙, consistent with the low oxygen abundance. However, the iron-to-oxygen abundance ratios (Fe/O) of the EMPGs are generally high; the EMPGs with the 2%-solar oxygen abundance show high Fe/O ratios of ∼90%–140% (Fe/O)⊙, which are unlikely to be explained by suggested scenarios of Type Ia supernova iron productions, iron’s dust depletion, and metal-poor gas inflow onto previously metal-riched galaxies with solar abundances. Moreover, the EMPGs with the 2%-solar oxygen abundance have very high He ii λ4686/Hβ ratios of ∼1/40, which are not reproduced by existing models of high-mass X-ray binaries with progenitor stellar masses M ⊙. Comparing stellar-nucleosynthesis and photoionization models with a comprehensive sample of EMPGs identified by this and previous EMPG studies, we propose that both the high Fe/O ratios and the high He ii λ4686/Hβ ratios are explained by the past existence of supermassive (>300 M ⊙) stars, which may evolve into intermediate-mass black holes (≳100 M ⊙).

The DESI Milky Way Survey (MWS) will observe $\ge$8 million stars between $16 < r < 19$ mag, supplemented by observations of brighter targets under poor observing conditions. The survey will permit an accurate determination of stellar kinematics and population gradients; characterize diffuse substructure in the thick disk and stellar halo; enable the discovery of extremely metal-poor stars and other rare stellar types; and improve constraints on the Galaxy's 3D dark matter distribution from halo star kinematics. MWS will also enable a detailed characterization of the stellar populations within 100 pc of the Sun, including a complete census of white dwarfs. The target catalog from the preliminary selection described here is public., 4 pages, 1 figure; published in Res. Notes AAS

We have initiated a new survey for local extremely metal-poor galaxies (EMPGs) with Subaru/Hyper Suprime-Cam (HSC) large-area (~500 deg^2) optical images reaching a 5 sigma limit of ~26 magnitude, about 100 times deeper than the Sloan Digital Sky Survey (SDSS). To select Z/Z_sun6 reported recently. We spectroscopically confirm that 3 out of 10 candidates are truly EMPGs with Z/Z_sun, Comment: 30 pages, 26 figures, and 8 tables; (Revision: metallicities changed due to the use of a more pricise dust correction method, but conclusion does not change, Minor: a table, some figures, and sentences are added for a clear explanation.)

We report our observations of HSC16aayt (SN 2016jiu), which was discovered by the Subaru/Hyper Suprime-Cam (HSC) transient survey conducted as part of Subaru Strategic Program (SSP). It shows very slow photometric evolution and its rise time is more than 100 days. The optical magnitude change in 400 days remains within 0.6 mag. Spectra of HSC16aayt show a strong narrow emission line and we classify it as a Type IIn supernova. The redshift of HSC16aayt is 0.6814 +/- 0.0002 from the spectra. Its host galaxy center is at 5 kpc from the supernova location and HSC16aayt might be another example of isolated Type IIn supernovae, although the possible existence of underlying star forming activity of the host galaxy at the supernova location is not excluded., Comment: 14 pages, 12 figures, 1 table, accepted by The Astrophysical Journal

We present an overview of a deep transient survey of the COSMOS field with the Subaru Hyper Suprime-Cam (HSC). The survey was performed for the 1.77 deg$^2$ ultra-deep layer and 5.78 deg$^2$ deep layer in the Subaru Strategic Program over 6- and 4-month periods from 2016 to 2017, respectively. The ultra-deep layer shows a median depth per epoch of 26.4, 26.3, 26.0, 25.6, and 24.6 mag in $g$, $r$, $i$, $z$, and $y$ bands, respectively; the deep layer is $\sim0.6$ mag shallower. In total, 1,824 supernova candidates were identified. Based on light curve fitting and derived light curve shape parameter, we classified 433 objects as Type Ia supernovae (SNe); among these candidates, 129 objects have spectroscopic or COSMOS2015 photometric redshifts and 58 objects are located at $z > 1$. Our unique dataset doubles the number of Type Ia SNe at $z > 1$ and enables various time-domain analyses of Type II SNe, high redshift superluminous SNe, variable stars, and active galactic nuclei., 17 pages, 17 figures, accepted for publication in PASJ

We present independent optical spectroscopic follow-up of WGD2038-4008/GRAL2038-4008, a background quasar strongly lensed by a foreground elliptical galaxy into four images, recently discovered independently by Agnello et al. [2018, MNRAS, 479, 4345] and Krone-Martins et al. [2018, A&A, 616, L11] thanks to the exquisite spatial resolution of Gaia. The quasar images are bright (i ∼ 19 mag), thus enabling us to reach S/N > 20 for the continuum within 30 min of exposure time with the Andalucia Faint Object Spectrograph and Camera spectrograph mounted on the 2.56-m Nordic Optical Telescope. The flexible scheduling and high sensitivity delivered by Andalucia Faint Object Spectrograph and Camera provide timely redshifts and reveal the nature of the quasar images; both are essential for lensing modelling and cosmography. Our analysis shows a strong emission feature in a data gap in Agnello et al. [2018, MNRAS, 479, 4345], which can be attributed to as an interloper emission line feature from the foreground lensing galaxy, or hinting to a higher redshift of the background quasar. We discuss these two scenarios and outline possible tests to verify these scenarios.

We present the results of a survey of the brightest UV-selected galaxies in protoclusters. These proto-brightest cluster galaxy (proto-BCG) candidates are drawn from 179 overdense regions of $g$-dropout galaxies at $z\sim4$ from the Hyper Suprime-Cam Subaru Strategic Program identified previously as good protocluster candidates. This study is the first to extend the systematic study of the progenitors of BCGs from $z\sim2$ to $z\sim4$. We carefully remove possible contaminants from foreground galaxies and, for each structure, we select the brightest galaxy that is at least 1 mag brighter than the fifth brightest galaxy. We select 63 proto-BCG candidates and compare their properties with those of galaxies in the field and those of other galaxies in overdense structures. The proto-BCG candidates and their surrounding galaxies have different rest-UV color $(i - z)$ distributions to field galaxies and other galaxies in protoclusters that do not host proto-BCGs. In addition, galaxies surrounding proto-BCGs are brighter than those in protoclusters without proto-BCGs. The image stacking analysis reveals that the average effective radius of proto-BCGs is $\sim28\%$ larger than that of field galaxies. The $i-z$ color differences suggest that proto-BCGs and their surrounding galaxies are dustier than other galaxies at $z\sim4$. These results suggest that specific environmental effects or assembly biasses have already emerged in some protoclusters as early as $z \sim 4$, and we suggest that proto-BCGs have different star formation histories than other galaxies in the same epoch., Comment: 18 pages, 11 figures, Accepted for publication in ApJ

With the advent of the Large Synoptic Survey Telescope (LSST), time-domain astronomy will be faced with an unprecedented volume and rate of data. Real-time processing of variables and transients detected by such large-scale surveys is critical to identifying the more unusual events and allocating scarce follow-up resources efficiently. We develop an algorithm to identify these novel events within a given population of variable sources. We determine the distributions of magnitude changes (dm) over time intervals (dt) for a given passband f, pf(dm|dt), and use these distributions to compute the likelihood of a test source being consistent with the population, or an outlier. We demonstrate our algorithm by applying it to the DECam multi-band time-series data of more than 2000 variable stars identified by Saha et al. (2019) in the Galactic Bulge that are largely dominated by long-period variables and pulsating stars. Our algorithm discovers 18 outlier sources in the sample, including a microlensing event, a dwarf nova, and two chromospherically active RS CVn stars, as well as sources in the Blue Horizontal Branch region of the color-magnitude diagram without any known counterparts. We compare the performance of our algorithm for novelty detection with multivariate KDE and Isolation Forest on the simulated PLAsTiCC dataset. We find that our algorithm yields comparable results despite its simplicity. Our method provides an efficient way for flagging the most unusual events in a real-time alert-broker system., 20 pages + Appendix, accepted to ApJ

Strong gravitationally lensed quasars provide powerful means to study galaxy evolution and cosmology. We use Chitah to hunt for new lens systems in the Hyper Suprime$-$Cam Subaru Strategic Program (HSC SSP) S16A. We present 46 lens candidates, of which 3 are previously known. Including 2 additional lenses found by YattaLens, we obtain X-shooter spectra of 6 promising candidates for lens confirmation and redshift measurements. We report new spectroscopic redshift measurements for both the lens and source galaxies in 4 lens systems. We apply the lens modeling software Glee to model our 6 X-shooter lenses uniformly. Through our analysis of the HSC images, we find that HSCJ022622$-$042522, HSCJ115252$+$004733, and HSCJ141136$-$010216 have point-like lensed images, and that the lens light distribution is well aligned with mass distribution within 6 deg. Thanks to the X-shooter spectra, we estimate fluxes on the Baldwin- Phillips-Terlevich (BPT) diagram, and find that HSCJ022622$-$042522 has a probable quasar source, based on the upper limit of the Nii flux intensity. We also measure the FWHM of Ly$\alpha$ emission of HSCJ141136$-$010216 to be $\sim$254 km/s, showing that it is a probable Lyman-$\alpha$ emitter., Comment: 13 pages, 7 figures

We present optical and near-infrared observations of the nearby Type Iax supernova (SN) 2014dt from 14 to 410 days after the maximum light. The velocities of the iron absorption lines in the early phase indicated that SN 2014dt showed slower expansion than the well-observed Type Iax SNe 2002cx, 2005hk and 2012Z. In the late phase, the evolution of the light curve and that of the spectra were considerably slower. The spectral energy distribution kept roughly the same shape after ~100 days, and the bolometric light curve flattened during the same period. These observations suggest the existence of an optically thick component that almost fully trapped the {\gamma}-ray energy from 56 Co decay. These findings are consistent with the predictions of the weak deflagration model, leaving a bound white dwarf remnant after the explosion., Comment: 39 pages, 14 figures, accepted for publication in PASJ

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., Comment: Accepted for publication in ApJ

The extragalactic $\gamma$-ray background (EGB) arises from the accumulation of $\gamma$-ray emissions from resolved and unresolved extragalactic sources as well as diffuse processes. It is important to study the statistical properties of the EGB in the context of cosmological structure formation. Known astrophysical $\gamma$-ray sources such as blazars, star-forming galaxies, and radio galaxies are expected to trace the underlying cosmic matter density distribution. We explore the correlation of the EGB from Fermi-LAT data with the large-scale matter density distribution from the Subaru Hyper Suprime-Cam (HSC) SSP survey. We reconstruct an unbiased surface matter density distribution $\kappa$ at $z 1$ are used for the lensing $\kappa$ reconstruction. The large-scale correlation suggests strong clustering of high-redshift $\gamma$-ray sources such as blazars. However, the inferred bias factor of $4-5$ is larger by about a factor of two than results from other clustering analyses. The final HSC data covering 1,400 squared degrees will play an essential role to determine accurately the blazar bias at $z > 0.5$., Comment: 17 pages, 8 figures, 4 tables, accepted for publication in Physical Review D

We present the largest Cepheid sample in M31 based on the complete Pan-STARRS1 survey of Andromeda (PAndromeda) in the $r_{\mathrm{P1}}$ , $i_{\mathrm{P1}}$ and $g_{\mathrm{P1}}$ bands. We find 2686 Cepheids with 1662 fundamental mode Cepheids, 307 first-overtone Cepheids, 278 type II Cepheids and 439 Cepheids with undetermined Cepheid type. Using the method developed by Kodric et al. (2013) we identify Cepheids by using a three dimensional parameter space of Fourier parameters of the Cepheid light curves combined with a color cut and other selection criteria. This is an unbiased approach to identify Cepheids and results in a homogeneous Cepheid sample. The Period-Luminosity relations obtained for our sample have smaller dispersions than in our previous work. We find a broken slope that we previously observed with HST data in Kodric et al. (2015), albeit with a lower significance., 79 pages, 39 figures, 8 tables, accepted for publication in AJ, K18b is submittted to ApJ, electronic data will be available on CDS

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., 24 pages, 21 figures, accepted in ApJ

We present the Ly$\alpha$ luminosity function (LF) derived from 34 Ly$\alpha$ emitters (LAEs) at $z=7.0$ on the sky of $3.1$ deg$^2$, the largest sample compared to those in the literature obtained at a redshift $z\gtrsim7$. The LAE sample is made by deep large-area Subaru narrowband observations conducted by the Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS) project. The $z=7.0$ Ly$\alpha$ LF of our project is consistent with those of the previous DECam and Subaru studies at the bright and faint ends, respectively, while our $z=7.0$ Ly$\alpha$ LF has uncertainties significantly smaller than those of the previous study results. Exploiting the small errors of our measurements, we investigate the shape of the faint to bright-end Ly$\alpha$ LF. We find that the $z=7.0$ Ly$\alpha$ LF shape can be explained by the steep slope of $\alpha \simeq -2.5$ suggested at $z=6.6$, and that there is no clear signature of a bright-end excess at $z\simeq 7$ claimed by the previous work, which was thought to be made by the ionized bubbles around bright LAEs whose Ly$\alpha$ photons could easily escape from the partly neutral IGM at $z \simeq 7$. We estimate the Ly$\alpha$ luminosity densities (LDs) with Ly$\alpha$ LFs at $z\simeq 6-8$ given by our and the previous studies, and compare the evolution of the UV-continuum LD estimated with dropouts. The Ly$\alpha$ LD monotonically decreases from $z\sim 6$ to $8$, and evolves stronger than the UV-continuum LD, indicative of the Ly$\alpha$ damping wing absorption of the IGM towards the heart of the reionization epoch., Comment: 13 pages, 15 figures, accepted in ApJ

The survey of Lyman $\alpha$ emitters (LAEs) with Subaru Hyper Suprime-Cam, called SILVERRUSH (Ouchi et al.), is producing massive data of LAEs at $z\gtrsim6$. Here we present LAE simulations to compare the SILVERRUSH data. In 162$^3$ comoving Mpc$^3$ boxes, where numerical radiative transfer calculations of reionization were performed, LAEs have been modeled with physically motivated analytic recipes as a function of halo mass. We have examined $2^3$ models depending on the presence or absence of dispersion of halo Ly$\alpha$ emissivity, dispersion of the halo Ly$\alpha$ optical depth, $\tau_\alpha$, and halo mass dependence of $\tau_\alpha$. The unique free parameter in our model, a pivot value of $\tau_\alpha$, is calibrated so as to reproduce the $z=5.7$ Ly$\alpha$ luminosity function (LF). We compare our model predictions with Ly$\alpha$ LFs at $z=6.6$ and $7.3$, LAE angular auto-correlation functions (ACFs) at $z=5.7$ and $6.6$, and LAE fractions in Lyman break galaxies at $5, Comment: PASJ in press, 30 pages, 20 figures

We present 1,201 galaxies at $0.05, Comment: 24 pages, 14 figures. Submitted to ApJ

Our understanding of extra-solar planet systems is highly driven by advances in observations in the past decade. Thanks to high precision spectrograph, we are able to reveal unseen companions to stars with the radial velocity method. High precision photometry from the space, especially with the Kepler mission, enables us to detect planets when they transit their stars and dim the stellar light by merely one percent or smaller. Ultra wide-field, high cadence, continuous monitoring of the Galactic bulge from different sites around the southern hemisphere provides us the opportunity to observe microlensing effects caused by planetary systems from the solar neighborhood, all the way to the Milky Way center. The exquisite AO imaging from ground-based large telescopes, coupled with high-contrast coronagraph, captured the photons directly emitted by planets around other stars. In this article, I present a concise review of the extra-solar planet discoveries, discussing the strengths and weaknesses of the major planetary detection methods, providing an overview of our current understanding of planetary formation and evolution given the tremendous observations delivered by various methods, as well as on-going and planned observation endeavors to provide a clear picture of extra-solar planetary systems., 11 pages, 1 figure, Galaxies accepted

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., 20 pages, 9 figures, 4 tables. Accepted for publication in PASJ

We present an unprecedentedly large catalog consisting of 2,230 > L* Lya emitters (LAEs) at z=5.7 and 6.6 on the 13.8 and 21.2 deg2 sky, respectively, that are identified by the SILVERRUSH program with the first narrowband imaging data of the Hyper Suprime-Cam (HSC) survey. We confirm that the LAE catalog is reliable on the basis of 96 LAEs whose spectroscopic redshifts are already determined by this program and the previous studies. This catalogue is also available on-line. Based on this catalogue, we derive the rest-frame Lya equivalent-width distributions of LAEs at z~5.7-6.6 that are reasonably explained by the exponential profiles with the scale lengths of ~120-170A, showing no significant evolution from z~5.7 to z~6.6. We find that 275 LAEs with a large equivalent width (LEW) of >240A are candidates of young-metal poor galaxies and AGNs. We also find that the fraction of LEW LAEs to all ones is 4% and 21% at z~5.7 and z~6.6, respectively. Our LAE catalog includes 11 Lya blobs (LABs) that are LAEs with spatially extended Lya emission whose profile is clearly distinguished from those of stellar objects at the >~ 3sigma level. The number density of the LABs at z=6-7 is ~ 10^-7-10^-6 Mpc^-3, being ~ 10-100 times lower than those claimed for LABs at z~ 2-3, suggestive of disappearing LABs at z>~6, albeit with the different selection methods and criteria for the low and high-z LABs., 19 pages, 11 figures. Accepted for publication in PASJ special issue

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, Comment: 20 pages, 10 figures, submitted to PASJ