Search

We present the study on the relationship between SMBHs and their host galaxies using our variability-selected AGN sample ($i_\mathrm{AB} \leq 25.9,\ z \leq 4.5$) constructed from the HSC-SSP Ultra-deep survey in the COSMOS field. We estimated the BH mass ($M_\mathrm{BH}=10^{5.5-10}\ M_{\odot}$) based on the single-epoch virial method and the total stellar mass ($M_\mathrm{star}=10^{10-12}\ M_{\odot}$) by separating the AGN component with SED fitting. We found that the redshift evolution of the BH-stellar mass ratio ($M_\mathrm{BH}/M_\mathrm{star}$) depends on the $M_\mathrm{BH}$ which is caused by the no significant correlation between $M_\mathrm{BH}$ and $M_\mathrm{star}$. Variable AGNs with massive SMBHs ($M_\mathrm{BH}>10^{9}\ M_{\odot}$) at $1.5\sim1\%$) than the BH-bulge ratios ($M_\mathrm{BH}/M_\mathrm{bulge}$) observed in the local universe for the same BH range. This implies that there is a typical growth path of massive SMBHs which is faster than the formation of the bulge component as final products seen in the present day. For the low-mass SMBHs ($M_\mathrm{BH}<10^{8}\ M_{\odot}$) at $0.5 4$. We interpret that host galaxies harboring less massive SMBHs at intermediate redshift have already acquired sufficient stellar mass, although high-z galaxies are still in the early stage of galaxy formation relative to those at the intermediate/local universe., Comment: Accepted for publication in ApJ, 14 pages, 5 figures

SDSS J1430+2303 has been argued to possess a supermassive black hole binary which is predicted to merge within a few months or three years from January 2022. We conducted follow-up optical spectroscopic observations of SDSS J1430+2303 with KOOLS-IFU on Seimei Telescope in May, June, and July 2022, and April 2023. The observed spectrum around $\mathrm{H}\mathrm{\alpha}$ shows a central broad component $\sim 10^3\ \mathrm{km\ s^{-1}}$ blueshifted from the narrow H$\mathrm{\alpha}$ line as well as the broader double-peaked component with a separation of $\sim\pm 5\times10^3\ \mathrm{km\ s^{-1}}$, similar to the spectrum reported in January 2022. We investigate the variability of the complex broad $\mathrm{H}\mathrm{\alpha}$ emission line relative to the continuum over the observation period. The continuum-normalized relative flux of the central broad component shows the increasing trend from May to July 2022 which is interpreted to be caused by the decrease of the continuum as also supported by damping of the X-ray, UV, and optical light curves observed for the same period. From July 2022 to April 2023, however, the central broad component decreased significantly. For the relative flux of the broader double-peaked component, on the other hand, no significant change appears at any epoch. These results suggest that the complicated broad line profile of SDSS J1430+2303 is generated from at least two distinct regions. While the central broad component originates from a broad line region, the broader double-peaked component arises in the vicinity of the continuum source., Comment: Published in Publications of the Astronomical Society of Japan. 6 pages, 3 figures

We performed an experimental verification of a coronagraph. As a result, we confirmed that, at the focal region where the planetary point spread function exists, the coronagraph system mitigates the raw contrast of a star-planet system by at least $1\times10^{-5}$ even for the 1-$\lambda/D$ star-planet separation. In addition, the verified coronagraph keeps the shapes of the off-axis point spread functions when the setup has the source angular separation of 1$\lambda/D$. The low-order wavefront error and the non-zero extinction ratio of the linear polarizer may affect the currently confirmed contrast. The sharpness of the off-axis point spread function generated by the sub-$\lambda/D$ separated sources is promising for the fiber-based observation of exoplanets. The coupling efficiency with a single mode fiber exceeds 50% when the angular separation is greater than 3--4$\times 10^{-1}\lambda/D$. For sub-$\lambda/D$ separated sources, the peak positions (obtained with Gaussian fitting) of the output point spread functions are different from the angular positions of sources; the peak position moved from about $0.8\lambda/D$ to $1.0\lambda/D$ as the angular separation of the light source varies from $0.1\lambda/D$ to $1.0\lambda/D$. The off-axis throughput including the fiber-coupling efficiency (with respect to no focal plane mask) is about 40% for 1-$\lambda/D$ separated sources and 10% for 0.5-$\lambda/D$ separated ones (excluding the factor of the ratio of pupil aperture width and Lyot stop width), where we assumed a linear-polarized-light injection. In addition, because this coronagraph can remove point sources on a line in the sky, it has another promising application for high-contrast imaging of exoplanets in binary systems., Comment: 18 pages, 10 figures, accepted for the Publications of the Astronomical Society of the Pacific

We report the detection of an ionized gas outflow from an $X$-ray active galactic nucleus (AGN) hosted in a massive quiescent galaxy in a protocluster at $z=3.09$ (J221737.29+001823.4). It is a type-2 QSO with broad ($W_{80}>1000$ km s$^{-1}$) and strong ($\log (L_{\rm [OIII]}$ / erg s$^{-1})\approx43.4$) [O {\footnotesize III}]$\lambda\lambda$4959,5007 emission lines detected by slit spectroscopy in three-position angles using Multi-Object Infra-Red Camera and Spectrograph (MOIRCS) on the Subaru telescope and the Multi-Object Spectrometer For Infra-Red Exploration (MOSFIRE) on the Keck-I telescope. In the all slit directions, [O {\footnotesize III}] emission is extended to $\sim15$ physical kpc and indicates a powerful outflow spreading over the host galaxy. The inferred ionized gas mass outflow rate is $\rm 22\pm3~M_{\odot}~yr^{-1}$. Although it is a radio source, according to the line diagnostics using H$\beta$, [O {\footnotesize II}], and [O {\footnotesize III}], photoionization by the central QSO is likely the dominant ionization mechanism rather than shocks caused by radio jets. On the other hand, the spectral energy distribution of the host galaxy is well characterized as a quiescent galaxy that has shut down star formation by several hundred Myr ago. Our results suggest a scenario that QSOs are powered after the shut-down of the star formation and help to complete the quenching of massive quiescent galaxies at high redshift., Comment: 18 pages 6 figures, accepted for publication in ApJ

The morphological study is crucial to investigate the connections between active galactic nuclei (AGN) activities and the evolution of galaxies. Substantial studies have found that radiative-mode AGNs primarily reside in disk galaxies, questioning the merger-driven mechanism of AGN activities. In this study, through S{\'e}rsic profile fitting and non-parametric morphological parameter measurements, we investigated the morphology of host galaxies of 485 optical variability-selected low luminosity AGNs at $z\lesssim4.26$ in the COSMOS field. We analyzed high-resolution images of the Hubble Space Telescope to measure these morphological parameters. We only successfully measured the morphological parameters for 76 objects and most AGN hosts ($\sim70\%$) were visually compact point-like sources. We examined the obtained morphological information as a function of redshift and compared them with literature data. We found that these AGN host galaxies showed no clear morphological preference. However, the merger rate increased with the higher hosts' SFRs and AGN luminosity. Interestingly, we found ongoing star formation consistent with the typical star forming populations in both elliptical and spiral galaxies while these two types of galaxies were more symmetric than normal star forming galaxies. These results suggested that optical variability-selected AGNs have} higher probabilities to reside in elliptical galaxies than infrared-selected AGNs (IR-AGNs), whose host galaxies had a strong disk-dominance, and supported recent studies that the AGN feedback could enhance star forming activities in host galaxies., Comment: Accepted for publication in ApJ

We conducted a deep spectroscopic survey, named SSA22-HIT, in the SSA22 field with the DEep Imaging MultiObject Spectrograph (DEIMOS) on the Keck telescope, designed to tomographically map high-z HI gas through analysis of Lya absorption in background galaxies' spectra. In total, 198 galaxies were spectroscopically confirmed at 2.5 < z < 6 with a few low-z exceptions in the 26 x 15 arcmin^2 area, of which 148 were newly determined in this study. Our redshift measurements were merged with previously confirmed redshifts available in the 34 x 27 arcmin^2 area of the SSA22 field. This compiled catalog containing 730 galaxies of various types at z > 2 is useful for various applications, and it is made publicly available. Our SSA22-HIT survey has increased by approximately twice the number of spectroscopic redshifts of sources at z > 3.2 in the observed field. From a comparison with publicly available redshift catalogs, we show that our compiled redshift catalog in the SSA22 field is comparable to those among major extragalactic survey fields in terms of a combination of wide area and high surface number density of objects at z > 2. About 40 % of the spectroscopically confirmed objects in SSA22-HIT show reasonable quality of spectra in the wavelengths shorter than Lya when a sufficient amount of smoothing is adopted. Our data set enables us to make the HI tomographic map at z > 3, which we present in a parallel study., Comment: Accepted for publication in Astronomical Journal, 29 pages, 17 figures, 7 tables

We searched for massive galaxy population in the known large-scale high-density structure of Lyman$~\alpha$ emitters (LAEs) at $z=2.39$ near the radio galaxy 53W002 by using $B,~V,~i^\prime,~J,~H,~{\rm{and}}~K_s$-bands imaging data taken with Suprime-Cam and MOIRCS on the Subaru telescope. We selected 62 protocluster member candidates by their $JHK_s$-band colors and Spectral Energy Distribution (SED) fitting analysis ($JHK_s$-selected galaxies) in our survey field of $70.2~{\rm{arcmin}}^{2}$, and compared their physical properties estimated from the SED fitting with a comparison sample in the COSMOS field. We found significant number density excesses for the $JHK_s$-selected galaxies in the 53W002 field at $K_s<22.25,~J-K_s>2,~{\rm{or}}~V-K_s>4$. In particular the number density of the $JHK_s$-selected galaxies with $K_s<22.25~{\rm{and}}~J-K_s>2$ in the 53W002 field is nine times higher than the comparison sample. Most of those with $K_s<22.25~{\rm{and}}~J-K_s>2$ are massive galaxies with $M_s>10^{11}~M_\odot$, and their sSFRs of $10^{-11}$--$10^{-10}~\rm{yr^{-1}}$ suggest that the star formation has not yet stopped completely. We also found a density excess of quiescent galaxies with $M_s=5\times10^{10}$--$10^{11}~M_\odot~{\rm{and}}~{\rm{sSFR}}<10^{-11}~\rm{yr^{-1}}$ as well as that of low-mass galaxies with $M_s=10^{9.75}$--$10^{10}~M_\odot$ and various sSFRs. The massive galaxies with $M_s>10^{11}~M_\odot$ are not located at the density peaks of LAEs, but widely distributed along the similar direction with the structure of LAEs over $\sim15$--$20$ comoving Mpc. On the other hand, the quiescent galaxies with ${\rm{sSFR}}<10^{-11}~\rm{yr^{-1}}$ clearly avoid the structure of LAEs. Our results suggest that massive galaxies also exist in this protocluster discovered by the moderate overdensity of LAEs and their star formation activities depend on location in the protocluster., Comment: 33 pages, 26 figures, accepted for publication in ApJ

We report a massive quiescent galaxy at $z_{\rm spec}=3.0922^{+0.008}_{-0.004}$ spectroscopically confirmed at a protocluster in the SSA22 field by detecting the Balmer and Ca {\footnotesize II} absorption features with multi-object spectrometer for infrared exploration (MOSFIRE) on the Keck I telescope. This is the most distant quiescent galaxy confirmed in a protocluster to date. We fit the optical to mid-infrared photometry and spectrum simultaneously with spectral energy distribution (SED) models of parametric and nonparametric star formation histories (SFH). Both models fit the observed SED well and confirm that this object is a massive quiescent galaxy with the stellar mass of $\log(\rm M_{\star}/M_{\odot}) = 11.26^{+0.03}_{-0.04}$ and $11.54^{+0.03}_{-0.00}$, and star formation rate of $\rm SFR/M_{\odot}~yr^{-1} <0.3$ and $=0.01^{+0.03}_{-0.01}$ for parametric and nonparametric models, respectively. The SFH from the former modeling is described as an instantaneous starburst while that of the latter modeling is longer-lived but both models agree with a sudden quenching of the star formation at $\sim0.6$ Gyr ago. This massive quiescent galaxy is confirmed in an extremely dense group of galaxies predicted as a progenitor of a brightest cluster galaxy formed via multiple mergers in cosmological numerical simulations. We newly find three plausible [O III]$\lambda$5007 emitters at $3.0791\leq z_{\rm spec}\leq3.0833$ happened to be detected around the target. Two of them just between the target and its nearest massive galaxy are possible evidence of their interactions. They suggest the future strong size and stellar mass evolution of this massive quiescent galaxy via mergers., Comment: 13 pages, 7 figures (in main text. 21 pages including appendices), accepted for publication in ApJ

We use the SPace Infrared telescope for Cosmology and Astrophysics (SPICA) project as a template to demonstrate how deep spectrophotometric surveys covering large cosmological volumes over extended fields (1-15 square degrees) with a mid-IR imaging spectrometer (17-36 micron) in conjunction with deep 70 micron photometry with a far-IR camera, at wavelengths which are not affected by dust extinction can answer the most crucial questions in current galaxy evolution studies. A SPICA-like mission will be able for the first time to provide an unobscured three dimensional (3-D, i.e. x, y and redshift z) view of galaxy evolution back to an age of the Universe of less than ~2 Gyrs, in the mid-IR rest-frame. This survey strategy will produce a full census of the Star formation Rate (SFR) in the Universe, using Polycyclic Aromatic Hydrocarbons (PAH) bands and fine-structure ionic lines, reaching the characteristic knee of the galaxy luminosity function, where the bulk of the population is distributed, at any redshift up to z ~3.5. Deep follow-up pointed spectroscopic observations with grating spectrometers { onboard the satellite}, across the full IR spectral range (17-210 micron), would simultaneously measure Black Hole Accretion Rate (BHAR), from high-ionization fine-structure lines, and SFR, from PAH and low- to mid-ionization lines in thousands of galaxies from solar to low metallicities, down to the knee of their luminosity functions. The analysis of the resulting atlas of IR spectra will reveal the physical processes at play in evolving galaxies across cosmic time, especially its heavily dust-embedded phase during the activity peak at the cosmic noon (z ~1-3), through IR emission lines and features that are insensitive to the dust obscuration., Comment: Paper accepted for publication on PASA on 5th March 2021, as part of the SPICA Special Issue

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

Using a sample of galaxies at $z\approx 3$ with detected Lyman Continuum (LyC) leakage in the SSA22 field, we attempt to verify a proposed indirect method for identifying cases with high LyC escape fraction $f_\mathrm{esc}$ based on measurements of the H$\beta$ equivalent width (EW) and the $\beta$ slope of the UV continuum. To this end, we present Keck/MOSFIRE H$\beta$ emission line flux measurements of LyC galaxies (LCGs) at spectroscopic redshifts $z_{\mathrm{spec}} \sim 3.3$, Lyman break galaxies (LBGs) at photometric redshifts $z_{\mathrm{phot}} = 2.7$--$3.7$, and Ly$\alpha$ emitters at $z_{\mathrm{phot}} = 3.1$. We also reconfirm the spectroscopic redshifts and measure the H$\beta$ emission line fluxes from 2 LCGs and 6 LBGs. For the LCG in our sample with the most extreme $f_\mathrm{esc}$, as revealed by the direct detection of LyC photons, we find that the EW(H$\beta$)--$\beta$ method gives a broadly consistent estimate for $f_\mathrm{esc}$, although the error bars remain very large. We also discuss how a combination of $f_\mathrm{esc}$ measurements based on direct and indirect methods can shed light on the LyC escape mechanism and the anisotropy of the leakage., Comment: 21 pages, 5figures, published in MNRAS

We analyzed the young (2.8-Myr-old) binary system FS Tau A using near-infrared (H-band) high-contrast polarimetry data from Subaru/HiCIAO and sub-millimeter CO (J=2-1) line emission data from ALMA. Both the near-infrared and sub-millimeter observations reveal several clear structures extending to $\sim$240 AU from the stars. Based on these observations at different wavelengths, we report the following discoveries. One arm-like structure detected in the near-infrared band initially extends from the south of the binary with a subsequent turn to the northeast, corresponding to two bar-like structures detected in ALMA observations with an LSRK velocity of 1.19-5.64 km/s. Another feature detected in the near-infrared band extends initially from the north of the binary, relating to an arm-like structure detected in ALMA observations with an LSRK velocity of 8.17-16.43 km/s. From their shapes and velocities, we suggest that these structures can mostly be explained by two streamers that connect the outer circumbinary disk and the central binary components. These discoveries will be helpful for understanding the evolution of streamers and circumstellar disks in young binary systems., Comment: 20 pages, 11 figures, accepted for publication in The Astrophysical Journal

We search for galaxies with a strong Balmer break (Balmer Break Galaxies; BBGs) at $z \sim 6$ over a 0.41 deg$^2$ effective area in the COSMOS field. Based on rich imaging data, including data obtained with the Atacama Large Millimeter/submillimeter Array (ALMA), three candidates are identified by their extremely red $K - [3.6]$ colors as well as by non-detection in X-ray, optical, far-infrared (FIR), and radio bands. The non-detection in the deep ALMA observations suggests that they are not dusty galaxies but BBGs at $z \sim 6$, although contamination from Active Galactic Nuclei (AGNs) at $z \sim 0$ cannot be completely ruled out for the moment. Our spectral energy distribution (SED) analyses reveal that the BBG candidates at $z \sim 6$ have stellar masses of $\approx 5 \times 10^{10} M_{\odot}$ dominated by old stellar populations with ages of $\gtrsim 700$ Myr. Assuming that all the three candidates are real BBGs at $z \sim 6$, we estimate the stellar mass density (SMD) to be $2.4^{+2.3}_{-1.3} \times 10^{4} M_{\odot}$ Mpc$^{-3}$. This is consistent with an extrapolation from the lower redshift measurements. The onset of star formation in the three BBG candidates is expected to be several hundred million years before the observed epoch of $z \sim 6$. We estimate the star-formation rate density (SFRD) contributed by progenitors of the BBGs to be 2.4 -- 12 $\times 10^{-5} M_{\odot}$ yr$^{-1} $Mpc$^{-3}$ at $z > 14$ (99.7\% confidence range). Our result suggests a smooth evolution of the SFRD beyond $z = 8$., Comment: 29 pages, 16 figures, 4 tables, accepted for publication in ApJ

The SR24 multi-star system hosts both circumprimary and circumsecondary disks, which are strongly misaligned from each other. The circumsecondary disk is circumbinary in nature. Interestingly, both disks are interacting, and they possibly rotate in opposite directions. To investigate the nature of this unique twin disk system, we present 0.''1 resolution near-infrared polarized intensity images of the circumstellar structures around SR24, obtained with HiCIAO mounted on the Subaru 8.2 m telescope. Both the circumprimary disk and the circumsecondary disk are resolved and have elongated features. While the position angle of the major axis and radius of the NIR polarization disk around SR24S are 55$^{\circ}$ and 137 au, respectively, those around SR24N are 110$^{\circ}$ and 34 au, respectively. With regard to overall morphology, the circumprimary disk around SR24S shows strong asymmetry, whereas the circumsecondary disk around SR24N shows relatively strong symmetry. Our NIR observations confirm the previous claim that the circumprimary and circumsecondary disks are misaligned from each other. Both the circumprimary and circumsecondary disks show similar structures in $^{12}$CO observations in terms of its size and elongation direction. This consistency is because both NIR and $^{12}$CO are tracing surface layers of the flared disks. As the radius of the polarization disk around SR24N is roughly consistent with the size of the outer Roche lobe, it is natural to interpret the polarization disk around SR24N as a circumbinary disk surrounding the SR24Nb-Nc system., Comment: 14 pages, 5 figures, accepted for publication in AJ

The CFHT Large Area U-band Deep Survey (CLAUDS) uses data taken with the MegaCam mosaic imager on CFHT to produce images of 18.60 deg2 with median seeing of FWHM=0.92 arcsec and to a median depth of U = 27.1 AB (5 sigma in 2 arcsec apertures), with selected areas that total 1.36 deg2 reaching a median depth of U=27.7 AB. These are the deepest U-band images assembled to date over this large an area. These data are located in four fields also imaged to comparably faint levels in grizy and several narrowband filters as part of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP). These CFHT and Subaru datasets will remain unmatched in their combination of area and depth until the advent of the Large Synoptic Survey Telescope (LSST). This paper provides an overview of the scientific motivation for CLAUDS and gives details of the observing strategy, observations, data reduction, and data merging with the HSC-SSP. Three early applications of these deep data are used to illustrate the potential of the dataset: deep U-band galaxy number counts, z~3 Lyman break galaxy (LBG) selection, and photometric redshifts improved by adding CLAUDS U to the Subaru HSC grizy photometry., Comment: MNRAS, in press

We report results of a search for galaxies at z>3 with Lyman continuum (LyC) emission using a narrow-band filter NB359 with Subaru / Suprime-Cam in a ~800 arcmin$^2$ blank field around the GOODS-N. We use 103 star-forming galaxies (SFGs) and 8 AGNs with spectroscopic redshifts in a range between 3.06 and 3.5, and 157 photometrically selected z=3.1 Lyman $\alpha$ emitter (LAE) candidates as the targets. After removing galaxies spectroscopically confirmed to be contaminated by foreground sources, we found two SFGs and one AGN as candidate LyC emitting sources among the targets with spectroscopic redshifts. Among LAE candidates, five sources are detected in the NB359 image, and three among them may be contaminated by foreground sources. We compare the sample galaxies in the GOODS-N with those in the SSA22, where a prominent protocluster at z=3.1 is known and the LyC search using the same NB359 filter has been made. Frequency of galaxies with LyC leakage in the SSA22 field may be about two times higher than the galaxies in the GOODS-N with the sample UV magnitude range, although the numbers of LyC detections in these fields are too small to make a statistically significant conclusion. By combining the sample galaxies in these fields, we place the 3$\sigma$ upper limits of the observed LyC-to-UV flux density ratio and LyC escape fraction for galaxies at z=3.1 with absolute UV magnitude $M_{UV} < -18.8$ as ($f_{LyC}/f_{UV})_{obs} < 0.036$ and $f_{esc}^{abs}<8$\%, respectively., Comment: 23 pages, 19 figures. Accepted for publication in MNRAS

We investigate the stellar population of star-forming galaxies at $z \sim 4$ by focusing on their slope of rest-frame ultraviolet (UV) continuum, $\beta$ where $f_{\lambda} \propto \lambda^{\beta}$. We investigate the sample of bright Lyman Break Galaxies (LBGs) with $i' \leq 26.0$ in the Subaru/XMM-Newton Deep Survey field by using the SED fitting analysis. We find that the apparently redder ($\beta_{\mathrm{obs}} > -1.73$) LBGs tend to be dusty ($\mathrm{Av} > 1.0$), have the young stellar population ($\beta_{\mathrm{int}} < -2.42$), and the intrinsically active star-forming galaxies (SFR $\gtrsim$ a few $\times 10^{2}$ $\mathrm{M_{\odot}} \mathrm{yr^{-1}}$). It means that a significant fraction of the UV-selected LBGs at $z \sim 4$ is on-going active and dust obscured star-forming galaxies. We compare the IR to UV luminosity ratio assuming the dust attenuation laws with the sub-millimeter observations from previous works. The result suggests that the Calzetti-like dust attenuation law is preferable for the active and dusty star-forming LBGs at $z = 4$. We also find that an extrapolation of the $\beta_{\mathrm{int}}$-$M_{\mathrm{UV,int}}$ relation toward the fainter magnitude range below our sample magnitude limit intersects the $\beta_{\mathrm{obs}}$-$M_{\mathrm{UV,obs}}$ relation previously obtained in the deeper narrow-area observations at $M_{\mathrm{UV}} = -18.9$ and $\beta = -1.94$, which coincides with the break point of $\beta_{\mathrm{obs}}$-$M_{\mathrm{UV,obs}}$ relation observed so far. The coincidence suggest that we see the almost dust-free population at $M_{\mathrm{UV,obs}} \gtrsim -18.9$., Comment: 20 pages, 20 figures

We present H-band polarized scattered light imagery and JHK high-contrast spectroscopy of the protoplanetary disk around HD 163296 observed with the HiCIAO and SCExAO/CHARIS instruments at Subaru Observatory. The polarimetric imagery resolve a broken ring structure surrounding HD 163296 that peaks at a distance along the major axis of 0.65 (66 AU) and extends out to 0.98 (100 AU) along the major axis. Our 2011 H-band data exhibit clear axisymmetry, with the NW- and SE- side of the disk exhibiting similar intensities. Our data are clearly different than 2016 epoch H-band observations from VLT/SPHERE that found a strong 2.7x asymmetry between the NW- and SE-side of the disk. Collectively, these results indicate the presence of time variable, non-azimuthally symmetric illumination of the outer disk. Based on our 3D-MCRT modeling of contemporaneous IR spectroscopic and H-band polarized intensity imagery of the system, we suggest that while the system could plausibly host an inclined inner disk component, such a component is unlikely to be responsible for producing the observed time-dependent azimuthal variations in the outer scattered light disk of the system. While our SCExAO/CHARIS data are sensitive enough to recover the planet candidate identified from NIRC2 in the thermal IR, we fail to detect an object with a corresponding JHK brightness estimated from the atmospheric models of Baraffe et al. 2003. This suggests that the candidate is either fainter in JHK bands than model predictions, possibly due to extinction from the disk or atmospheric dust/clouds, or that it is an artifact of the dataset/data processing. Our SCExAO/CHARIS data lower the IR mass limits for planets inferred at larger stellocentric separations; however, these ALMA-predicted protoplanet candidates are currently still consistent with direct imaging constraints., Comment: 19 pages, 13 figures, 2 tables

We conducted high-contrast polarimetry observations of T Tau in the H-band, using the HiCIAO instrument mounted on the Subaru Telescope, revealing structures as near as 0.$\arcsec$1 from the stars T Tau N and T Tau S. The whole T Tau system is found to be surrounded by nebula-like envelopes, and several outflow-related structures are detected in these envelopes. We analyzed the detailed polarization patterns of the circumstellar structures near each component of this triple young star system and determined constraints on the circumstellar disks and outflow structures. We suggest that the nearly face-on circumstellar disk of T Tau N is no larger than 0.$\arcsec$8, or 117 AU, in the northwest, based on the existence of a hole in this direction, and no larger than 0.$\arcsec$27, or 40 AU, in the south. A new structure "N5" extends to about 0.$\arcsec$42, or 59 AU, on the southwest of the star, believed to be part of the disk. We suggest that T Tau S is surrounded by a highly inclined circumbinary disk with a radius of about 0.$\arcsec$3, or 44 AU, with a position angle of about 30$^\circ$, that is misaligned with the orbit of the T Tau S binary. After analyzing the positions and polarization vector patterns of the outflow-related structures, we suggest that T Tau S should trigger the well-known E-W outflow, and is also likely to be responsible for a southwest precessing outflow "coil" and a possible south outflow., Comment: 12 pages, 6 figures, accepted by The Astrophysical Journal

We present $H$- and $K_{\rm s}$-bands observations of the LkH$\alpha$ 330 disk with a multi-band detection of the large gap and spiral-like structures. The morphology of the outer disk ($r\sim$$0\farcs3$) at PA=0--45$^\circ$ and PA=180--290$^\circ$ are likely density wave-induced spirals and comparison between our observational results and simulations suggests a planet formation. We have also investigated the azimuthal profiles at the ring and the outer-disk regions as well as radial profiles in the directions of the spiral-like structures and semi-major axis. Azimuthal analysis shows a large variety in wavelength and implies that the disk has non-axisymmetric dust distributions. The radial profiles in the major-axis direction (PA=$271^\circ$) suggest that the outer region ($r\geq0\farcs25$) may be influenced by shadows of the inner region of the disk. The spiral-like directions (PA=10$^\circ$ and 230$^\circ$) show different radial profiles, which suggests that the surfaces of the spiral-like structures are highly flared and/or have different dust properties. Finally, a color-map of the disk shows a lack of an outer eastern region in the $H$-band disk, which may hint the presence of an inner object that casts a directional shadow onto the disk., Comment: 12pages, 16 figures, 2 tables, accepted for publication in AJ

Large exoplanet surveys have successfully detected thousands of exoplanets to-date. Utilizing these detections and non-detections to constrain our understanding of the formation and evolution of planetary systems also requires a detailed understanding of the basic properties of their host stars. We have determined the basic stellar properties of F, K, and G stars in the Strategic Exploration of Exoplanets and Disks with Subaru (SEEDS) survey from echelle spectra taken at the Apache Point Observatory's 3.5m telescope. Using ROBOSPECT to extract line equivalent widths and TGVIT to calculate the fundamental parameters, we have computed Teff, log(g), vt, [Fe/H], chromospheric activity, and the age for our sample. Our methodology was calibrated against previously published results for a portion of our sample. The distribution of [Fe/H] in our sample is consistent with that typical of the Solar neighborhood. Additionally, we find the ages of most of our sample are $< 500 Myrs$, but note that we cannot determine robust ages from significantly older stars via chromospheric activity age indicators. The future meta-analysis of the frequency of wide stellar and sub-stellar companions imaged via the SEEDS survey will utilize our results to constrain the occurrence of detected co-moving companions with the properties of their host stars., Comment: 22 Pages, 10 Figures, 5 Tables. Published in MNRAS

Using galaxies as background light sources to map intervening Lya absorption is a novel approach to study the interplay among galaxies, the circum-galactic medium (CGM), and the intergalactic medium (IGM). Introducing a new measure of z = 3.1$ HI Lya absorption relative to the cosmic mean, Delta_NB497, estimated from photometric data of star-forming galaxies at 3.3 < z < 3.5, we have made two-dimensional Delta_NB497 maps in the z = 3.1 SSA22 proto-cluster region and two control fields (SXDS and GOODS-N fields) with a spatial resolution of ~ 5 comoving Mpc. The Delta_NB497 measurements in the SSA22 field are systematically larger than those in the control fields, and this HI absorption enhancement extends more than 50 comoving Mpc. The field-averaged (i.e., ~50 comoving Mpc scale) Delta_NB497 and the overdensity of Lya emitters (LAEs) seem to be correlated, while there is no clear dependency of the Delta_NB497 on the local LAE overdensity in a few comoving Mpc scale. These results suggest that diffuse HI gas spreads out in/around the SSA22 proto-cluster. We have also found an enhancement of Delta_NB497 at a projected distance < 100 physical kpc from the nearest z = 3.1 galaxies at least in the SSA22 field, which is probably due to HI gas associated with the CGM of individual galaxies. The HI absorption enhancement in the CGM-scale tends to be weaker around galaxies with stronger Lya emission, which suggests that the Lya escape fraction from galaxies depends on hydrogen neutrality in the CGM., Comment: 13 pages, 11 figures, 1 table, accepted to MNRAS

By performing non-masked polarization imaging with Subaru/HiCIAO, polarized scattered light from the inner region of the disk around the GG Tau A system was successfully detected in the $H$ band with a spatial resolution of approximately 0.07$\arcsec$, revealing the complicated inner disk structures around this young binary. This paper reports the observation of an arc-like structure to the north of GG Tau Ab and part of a circumstellar structure that is noticeable around GG Tau Aa extending to a distance of approximately 28 AU from the primary star. The speckle noise around GG Tau Ab constrains its disk radius to <13 AU. Based on the size of the circumbinary ring and the circumstellar disk around GG Tau Aa, the semi-major axis of the binary's orbit is likely to be 62 AU. A comparison of the present observations with previous ALMA and near-infrared (NIR) H$_2$ emission observations suggests that the north arc could be part of a large streamer flowing from the circumbinary ring to sustain the circumstellar disks. According to the previous studies, the circumstellar disk around GG Tau Aa has enough mass and can sustain itself for a duration sufficient for planet formation; thus, our study indicates that planets can form within close (separation $\lesssim$ 100 AU) young binary systems., Comment: Accepted for publication in AJ, 12 pages, 5 figures

We present H-band (1.6 {\mu}m) scattered light observations of the transitional disk RX J1615.3-3255, located in the ~1 Myr old Lupus association. From a polarized intensity image, taken with the HiCIAO instrument of the Subaru Telescope, we deduce the position angle and the inclination angle of the disk. The disk is found to extend out to 68 $\pm$ 12 AU in scattered light and no clear structure is observed. Our inner working angle of 24 AU does not allow us to detect a central decrease in intensity similar to that seen at 30 AU in the 880 {\mu}m continuum observations. We compare the observations with multiple disk models based on the Spectral Energy Distribution (SED) and submm interferometry and find that an inner rim of the outer disk at 30 AU containing small silicate grains produces a polarized intensity signal which is an order of magnitude larger than observed. We show that a model in which the small dust grains extend smoothly into the cavity found for large grains is closer to the actual H-band observations. A comparison of models with different dust size distributions suggests that the dust in the disk might have undergone significant processing compared to the interstellar medium., Comment: 8 pages, 7 figures, 4 tables. Accepted for publication in A&A

We find a new substellar companion to the Pleiades member star, Pleiades HII 3441, using the Subaru telescope with adaptive optics. The discovery is made as part of the high-contrast imaging survey to search for planetary-mass and substellar companions in the Pleiades and young moving groups. The companion has a projected separation of 0".49 +/- 0".02 (66 +/- 2 AU) and a mass of 68 +/- 5 M_J based on three observations in the J-, H-, and K_S-band. The spectral type is estimated to be M7 (~2700 K), and thus no methane absorption is detected in the H band. Our Pleiades observations result in the detection of two substellar companions including one previously reported among 20 observed Pleiades stars, and indicate that the fraction of substellar companions in the Pleiades is about 10.0 +26.1/-8.8 %. This is consistent with multiplicity studies of both the Pleiades stars and other open clusters., Comment: Main text (14 pages, 4 figures, 4 tables), and Supplementary data (8 pages, 3 tables). Accepted for Publications of Astronomical Society of Japan

We present new analysis of multi-epoch, H-band, scattered light images of the AB Aur system. We used a Monte Carlo, radiative transfer code to simultaneously model the system's SED and H-band polarized intensity imagery. We find that a disk-dominated model, as opposed to one that is envelope dominated, can plausibly reproduce AB Aur's SED and near-IR imagery. This is consistent with previous modeling attempts presented in the literature and supports the idea that at least a subset of AB Aur's spirals originate within the disk. In light of this, we also analyzed the movement of spiral structures in multi-epoch H-band total light and polarized intensity imagery of the disk. We detect no significant rotation or change in spatial location of the spiral structures in these data, which span a 5.8 year baseline. If such structures are caused by disk-planet interactions, the lack of observed rotation constrains the location of the orbit of planetary perturbers to be >47 AU., Comment: 8 pages, 3 figures, 1 table, Accepted to ApJ

We present the first near-infrared scattered-light detection of the transitional disk around V1247 Ori, which was obtained using high-resolution polarimetric differential imaging observations with Subaru/HiCIAO. Our imaging in the H band reveals the disk morphology at separations of ~0.14"-0.86" (54-330 au) from the central star. The polarized intensity (PI) image shows a remarkable arc-like structure toward the southeast of the star, whereas the fainter northwest region does not exhibit any notable features. The shape of the arm is consistent with an arc of 0.28" $\pm$ 0.09" in radius (108 au from the star), although the possibility of a spiral arm with a small pitch angle cannot be excluded. V1247 Ori features an exceptionally large azimuthal contrast in scattered, polarized light; the radial peak of the southeastern arc is about three times brighter than the northwestern disk measured at the same distance from the star. Combined with the previous indication of an inhomogeneous density distribution in the gap at $\lesssim$46 au, the notable asymmetry in the outer disk suggests the presence of unseen companions and/or planet-forming processes ongoing in the arc., Comment: 21 pages, 5 figures, accepted for publication in PASJ

We present high-contrast observations of 68 young stellar objects (YSOs) explored as part of the SEEDS survey on the Subaru telescope. Our targets are very young ($<$10 Myr) stars, which often harbor protoplanetary disks where planets may be forming. We achieve a typical contrast of $\sim$$10^{-4}$--$10^{-5.5}$ at an angular distance of 1\arcsec\ from the central star, corresponding to typical mass sensitivities (assuming hot-start evolutionary models) of $\sim$10 ${\rm M_J}$ at 70 AU and $\sim$6 ${\rm M_J}$ at 140 AU. We detected a new stellar companion to HIP 79462 and confirmed the substellar objects GQ Lup b and ROXs 42B b. An additional six companion candidates await follow-up observations to check for common proper motion. Our SEEDS YSO observations probe the population of planets and brown dwarfs at the very youngest ages; these may be compared to the results of surveys targeting somewhat older stars. Our sample and the associated observational results will help enable detailed statistical analyses of giant planet formation., Comment: 31pages, 13 figures, accepted to AJ

We study the extent to which very bright (-23.0 < MUV < -21.75) Lyman-break selected galaxies at redshifts z~7 display detectable Lya emission. To explore this issue, we have obtained follow-up optical spectroscopy of 9 z~7 galaxies from a parent sample of 24 z~7 galaxy candidates selected from the 1.65 sq.deg COSMOS-UltraVISTA and SXDS-UDS survey fields using the latest near-infrared public survey data, and new ultra-deep Subaru z'-band imaging (which we also present and describe in this paper). Our spectroscopy has yielded only one possible detection of Lya at z=7.168 with a rest-frame equivalent width EW_0 = 3.7 (+1.7/-1.1) Angstrom. The relative weakness of this line, combined with our failure to detect Lya emission from the other spectroscopic targets allows us to place a new upper limit on the prevalence of strong Lya emission at these redshifts. For conservative calculation and to facilitate comparison with previous studies at lower redshifts, we derive a 1-sigma upper limit on the fraction of UV bright galaxies at z~7 that display EW_0 > 50 Angstrom, which we estimate to be < 0.23. This result may indicate a weak trend where the fraction of strong Lya emitters ceases to rise, and possibly falls between z~6 and z~7. Our results also leave open the possibility that strong Lya may still be more prevalent in the brightest galaxies in the reionization era than their fainter counterparts. A larger spectroscopic sample of galaxies is required to derive a more reliable constraint on the neutral hydrogen fraction at z~7 based on the Lya fraction in the bright galaxies., Comment: 20 pages, 7 figures, accepted for publication in ApJ

We report on the identification of the old stellar population galaxy candidates at z > 5. We developed a new infrared color selection scheme to isolate galaxies with the strong Balmer breaks at z > 5, and applied it to the ultra-deep and wide infrared survey data from the Spitzer Extended Deep Survey (SEDS) and the UKIRT Infrared Deep Sky Survey. The eight objects satisfying K - [3.6] > 1.3 and K - [3.6] > 2.4 ([3.6] - [4.5]) + 0.6 are selected in the 0.34 deg^2 SEDS Ultra Deep Survey field. Rich multi-wavelength imaging data from optical to far-infrared are also used to reject blending sources and strong nebular line emitters, and we finally obtained the three most likely evolved galaxies at z > 5. Their stacked spectral energy distribution is fitted well with the old stellar population template with M_{*} = (7.5+-1.5) x 10^{10} Msun, star formation rate = 0.9 +- 0.2 Msun yr^{-1}, dust A_V < 1, and age = 0.7+-0.4 Gyr at z = 5.7+-0.6, where the dusty star-forming galaxies at z ~ 2.8 are disfavored because of the faintness in the 24um. The stellar mass density of these evolved galaxy candidates, (6+-4) x 10^4 Msun Mpc^{-3}, is much lower than that of star-forming galaxies, but the non-zero fraction suggests that initial star-formation and quenching have been completed by z ~ 6., Comment: 13 pages, 9 figures, 4 tables, and accepted for publication in PASJ

A radial velocity (RV) survey for intermediate-mass giants has been operated for over a decade at Okayama Astrophysical Observatory (OAO). The OAO survey has revealed that some giants show long-term linear RV accelerations (RV trends), indicating the presence of outer companions. Direct imaging observations can help clarify what objects generate these RV trends. We present the results of high-contrast imaging observations or six intermediate-mass giants with long-term RV trends using the Subaru Telescope and HiCIAO camera. We detected co-moving companions to $\gamma$ Hya B ($0.61^{+0.12}_{-0.14} M_\odot$), HD 5608 B ($0.10 \pm 0.01 M_\odot$), and HD 109272 B ($0.28 \pm 0.06 M_\odot$). For the remaining targets($\iota$ Dra, 18 Del, and HD 14067) we exclude companions more massive than 30-60 $M_\mathrm{Jup}$ at projected separations of 1arcsec-7arcsec. We examine whether these directly imaged companions or unidentified long-period companions can account for the RV trends observed around the six giants. We find that the Kozai mechanism can explain the high eccentricity of the inner planets $\iota$ Dra b, HD 5608 b, and HD 14067 b., Comment: 33 pases, 9 figures, accepted to ApJ

We present high-contrast H-band polarized intensity images of the transitional disk around the young solar-like star LkCa 15. By utilizing Subaru/HiCIAO for polarimetric differential imaging, both the angular resolution and the inner working angle reach 0.07" and r=0.1", respectively. We obtained a clearly resolved gap (width <~ 27 AU) at ~ 48 AU from the central star. This gap is consistent with images reported in previous studies. We also confirmed the existence of a bright inner disk with a misaligned position angle of 13+/-4 degree with respect to that of the outer disk, i.e., the inner disk is possibly warped. The large gap and the warped inner disk both point to the existence of a multiple planetary system with a mass of <~1Mjup., Comment: 6 pages, 3 figures, accepted in PASJ on December 14, 2015

We present the largest to date sample of hydrogen Lyman continuum (LyC) emitting galaxy candidates at any redshift, with 18 Ly$\alpha$ Emitters (LAEs) and 7 Lyman Break Galaxies (LBGs), obtained from the SSA22 field with Subaru/Suprime-Cam. The sample is based on the 159 LAEs and 136 LBGs observed in the field, all with spectroscopically confirmed redshifts, and these LyC candidates are selected as galaxies with counterpart in a narrow-band filter image which traces LyC at $z\geq 3.06$. Many LyC candidates show a spatial offset between the rest-frame non-ionizing ultraviolet (UV) detection and the LyC-emitting substructure or between the Ly$\alpha$ emission and LyC. The possibility of foreground contamination complicates the analysis of the nature of LyC emitters, although statistically it is highly unlikely that all candidates in our sample are contaminated by foreground sources. Many viable LyC LAE candidates have flux density ratios inconsistent with standard models, while also having too blue UV slopes to be foreground contaminants. Stacking reveals no significant LyC detection, suggesting that there is a dearth of objects with marginal LyC signal strength, perhaps due to a bimodality in the LyC emission. The foreground contamination-corrected $3\sigma$ upper limits of the observed average flux density ratios are $f_{LyC}/f_{UV}<0.08$ from stacking LAEs and $f_{LyC}/f_{UV}<0.02$ from stacking LBGs. There is a sign of a positive correlation between LyC and Ly$\alpha$, suggesting that both types of photons escape via a similar mechanism. The LyC detection rate among proto-cluster LBGs is seemingly lower compared to the field., Comment: 24 pages, 23 figures, 4 tables. Accepted