Source author record

Yoshiki Toba

Yoshiki Toba appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

44works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

44 published item(s)

preprint2026arXiv

A z ~ 0.4 Galaxy Reflecting the high-redshift Little Red Dots: An Extended Starburst with an Overmassive Black Hole

One of the most remarkable discoveries of JWST is a population of compact, red sources at z > 4, commonly referred to as Little Red Dots (LRDs). Spectroscopic identifications reported that most LRDs are active galactic nuclei (AGNs), which are preferentially found around z ~ 6 and could imply a key phase in the formation and growth of black holes (BHs) in the early universe. Photometric surveys at lower redshift have recently been carried out to trace their evolution across cosmic time, and a small number of LRD-like galaxies have been spectroscopically identified at both Cosmic Noon and in the local universe. Here we report the discovery of one of the lowest-z LRD-like galaxies, J204837.26-002437.2 (hereafter J2048) at z = 0.4330, using new Gemini-N/GMOS IFU observations combined with archival multi-band photometric SED data. The GMOS data reveal extended blue emission from starburst with a star formation rate of 400 solar mass per year, together with an extended, highly fast ionized outflow. This is the first spectroscopic confirmation of extended host emission and outflow in an LRD-like galaxy, providing a unique laboratory for understanding the nature of their high-redshift counterparts. Moreover, J2048 would host an extremely overmassive BH with a BH-to-stellar mass ratio of ~ 60%, with the BH mass and host stellar mass estimated to be 10^10.2 and 10^10.4 solar masses, respectively. We discuss the origin and evolutionary fate of J2048, and the implications that such low-z analogs have for interpreting the properties of high-z LRDs.

preprint2024arXiv

A population of Optically Quiescent Quasars from WISE and SDSS

The growth of active galactic nuclei (AGN) occurs under some form of obscuration in a large fraction of the population. The difficulty in constraining this population leads to high uncertainties in cosmic X-ray background and galaxy evolution models. Using an SDSS-WISE cross-match, we target infrared luminous AGN ($W1-W2$ > 0.8, and monochromatic rest-frame luminosity above $λL_λ$(12$μ$m) $\approx$ 3 $\times$ 10$^{44}$ erg s$^{-1}$), but with passive galaxy-like optical spectra (Optically Quiescent Quasars; OQQs). We find 47 objects that show no significant [O III]$λ$5007 emission, a typically strong AGN optical emission line. As a comparison sample, we examine SDSS-selected Type 2 quasars (QSO2s), which show a significant [O III]$λ$5007 line by definition. We find a 1:16 ratio of OQQs compared to QSO2s, suggesting that the OQQ duty cycle is likely much shorter than that of QSO2s (though selection biases are not fully quantified). We consider observed properties in comparison with other galaxy types, and examine them for consistency with theories on their intrinsic nature: chiefly (a) a high covering factor for surrounding obscuring matter, preventing the detection of high-ionisation emission lines - `cocooned AGN'; or (b) ionised gas being absent on the kpc scales of the Narrow Line Region (NLR), perhaps due to a `switching on' or `young' AGN. OQQs do not obviously fit the standard paradigm for merger-driven AGN and host galaxy evolution, implying we may be missing part of the flow of AGN evolution.

preprint2022arXiv

A Wide and Deep Exploration of Radio Galaxies with Subaru HSC (WERGS). VII. Redshift Evolution of Radio Galaxy Environments at z=0.3-1.4

We examine the redshift evolution of density environments around 2,163 radio galaxies with the stellar masses of $\sim10^{9}-10^{12} M_\odot$ between redshifts of $z=0.3-1.4$, based on the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) and Faint Images of the Radio Sky at Twenty-cm (FIRST). We use the $k$-nearest neighbor method to measure the local galaxy number density around our radio galaxy sample. We find that the overdensities of the radio galaxies are weakly but significantly anti-correlated with redshift. This is consistent with the known result that the relative abundance of less-massive radio galaxies increases with redshift, because less-massive radio galaxies reside in relatively low density regions. Massive radio galaxies with stellar mass of $M_* >10^{11}M_\odot$ are found in high density environments compared with the control sample galaxies with radio-non-detection and matched-stellar-mass. Less-massive radio galaxies with $M_* <10^{11}M_\odot$ reside in average density environments. The fraction of the radio galaxies associated with the neighbors within a typical major merger scale, $<70$ kpc, is higher than (comparable to) that of the control galaxies at $M_* >10^{11}M_\odot$ ($M_* <10^{11}M_\odot$). We also find that the local densities around the radio galaxies are anti-correlated with the radio luminosities and black hole mass accretion rates at fixed stellar mass. These findings suggest that massive radio galaxies have matured through galaxy mergers in the past, and have supermassive black holes whose mass accretion almost ceased at $z>1.4$, while less-massive radio galaxies undergo active accretion just at this epoch, as they have avoided such merger events.

preprint2022arXiv

Detection of extended millimeter emission in the host galaxy of 3C273 and its implications for QSO feedback via high dynamic range ALMA imaging

We estimate the amount of negative feedback energy injected into the ISM of the host galaxy of 3C273, a prototypical radio loud quasar. We obtained 93, 233 and 343 GHz continuum images with the Atacama Large Millimeter/Sub-millimeter Array (ALMA). After self calibration and point source subtraction, we reach an image dynamic range of $\sim 85000$ at 93\ GHz, $\sim 39000$ at 233\ GHz and $\sim 2500$ at 343\ GHz. These are currently the highest image dynamic range obtained using ALMA. We detect spatially extended millimeter emission associated with the host galaxy, cospatial with the Extended Emission Line Region (EELR) observed in the optical. The millimeter spectral energy distribution and comparison with centimeter data show that the extended emission cannot be explained by dust thermal emission, synchrotron or thermal bremsstrahlung arising from massive star formation. We interpret the extended millimeter emission as thermal bremsstrahlung from gas directly ionized by the central source. The extended flux indicates that at least $\sim 7\%$ of the bolometric flux of the nuclear source was used to ionize atomic hydrogen in the host galaxy. The ionized gas is estimated to be as massive as $10^{10}$ to $10^{11}\ \mathrm{M_\odot}$, but the molecular gas fraction with respect to the stellar mass is consistent with other ellipticals, suggesting that direct ionization ISM by the QSO may not be sufficient to suppress star formation, or we are witnessing a short timescale before negative feedback becomes observable. The discovery of a radio counterpart to EELRs provides a new pathway to studying the QSO-host ISM interaction.

preprint2022arXiv

Do blue galaxy-clusters have hot intracluster gas?

We present herein a systematic X-ray analysis of blue galaxy-clusters at $z=0.84$ discovered by the Subaru telescope. The sample consisted of 43 clusters identified by combining red-sequence and blue-cloud surveys, covering a wide range of emitter fractions (i.e., 0.3--0.8). The spatial extent of the over-density region of emitter galaxies was approximately 1~Mpc in radius. The average cluster mass was estimated as $0.6 (<1.5)\times10^{14}~{\rm M_\odot}$ from the stacked weak-lensing measurement. We analyzed the XMM-Newton archival data, and measured the X-ray luminosity of the hot intracluster medium. As a result, diffuse X-ray emission was marginally detected in 14 clusters, yielding an average luminosity of $5\times 10^{42}~{\rm erg\,s^{-1}}$. On the contrary, it was not significant in 29 clusters. The blue clusters were significantly fainter than the red-dominated clusters, and the X-ray luminosity did not show any meaningful correlation either with emitter fraction or richness. The X-ray surface brightness was low, but the amount of gas mass was estimated to be comparable to that observed in the $10^{13-14}~{\rm M_{\odot}}$ cluster. Based on the results, we suggest that the blue clusters are at the early formation stage, and the gas is yet to be compressed and heated up to produce appreciable X-rays. Follow-up spectroscopic measurements are essential to clarify the dynamical status and co-evolution of galaxies and hot gas in the blue clusters.

preprint2022arXiv

Extreme nature of four blue-excess dust-obscured galaxies revealed by optical spectroscopy

We report optical spectroscopic observations of four blue-excess dust-obscured galaxies (BluDOGs) identified by Subaru Hyper Suprime-Cam. BluDOGs are a sub-class of dust-obscured galaxies (DOGs, defined with the extremely red color $(i-[22])_{\rm AB} \geq 7.0$; Toba et al. 2015), showing a significant flux excess in the optical $g$- and $r$-bands over the power-law fits to the fluxes at the longer wavelengths. Noboriguchi et al. (2019) has suggested that BluDOGs may correspond to the blowing-out phase involved in a gas-rich major merger scenario. However the detailed properties of BluDOGs are not understood because of the lack of spectroscopic information. In this work, we carry out deep optical spectroscopic observations of four BluDOGs using Subaru/FOCAS and VLT/FORS2. The obtained spectra show broad emission lines with extremely large equivalent widths, and a blue wing in the CIV line profile. The redshifts are between 2.2 and 3.3. The averaged rest-frame equivalent widths of the CIV lines are $160\pm33$ $\mathrm{\mathring{A}}$, $\sim$7 times higher than the average of a typical type-1 quasar. The FWHMs of their velocity profiles are between 1990 and 4470 ${\rm km\ s^{-1}}$, and their asymmetric parameters are 0.05 and 0.25. Such strong CIV lines significantly affect the broad-band magnitudes, which is partly the origin of the blue excess seen in the spectral energy distribution of BluDOGs. Their estimated supermassive black hole masses are $1.1\times10^8 < M_{\rm BH}/M_\odot < 5.5 \times 10^8$. The inferred Eddington ratios of the BluDOGs are higher than 1 ($1.1< λ_{\rm Edd} < 3.8$), suggesting that the BluDOGs are in a rapidly evolving phase of supermassive black holes.

preprint2022arXiv

Metallicity-PAH Relation of MIR-selected Star-forming Galaxies in AKARI North Ecliptic Pole-wide Survey

We investigate the variation in the mid-infrared spectral energy distributions of 373 low-redshift ($z<0.4$) star-forming galaxies, which reflects a variety of polycyclic aromatic hydrocarbon (PAH) emission features. The relative strength of PAH emission is parameterized as $q_\mathrm{PAH}$, which is defined as the mass fraction of PAH particles in the total dust mass. With the aid of continuous mid-infrared photometric data points covering 7-24$μ$m and far-infrared flux densities, $q_\mathrm{PAH}$ values are derived through spectral energy distribution fitting. The correlation between $q_\mathrm{PAH}$ and other physical properties of galaxies, i.e., gas-phase metallicity ($12+\mathrm{log(O/H)}$), stellar mass, and specific star-formation rate (sSFR) are explored. As in previous studies, $q_\mathrm{PAH}$ values of galaxies with high metallicity are found to be higher than those with low metallicity. The strength of PAH emission is also positively correlated with the stellar mass and negatively correlated with the sSFR. The correlation between $q_\mathrm{PAH}$ and each parameter still exists even after the other two parameters are fixed. In addition to the PAH strength, the application of metallicity-dependent gas-to-dust mass ratio appears to work well to estimate gas mass that matches the observed relationship between molecular gas and physical parameters. The result obtained will be used to calibrate the observed PAH luminosity-total infrared luminosity relation, based on the variation of MIR-FIR SED, which is used in the estimation of hidden star formation.

preprint2022arXiv

On the Connection between Supermassive Black Hole and Galaxy Growth in the Reionization Epoch

The correlation between the mass of supermassive black holes (SMBHs; $\mathcal{M}_{\rm BH}$) and their host galaxies ($\mathcal{M}_\star$) in the reionization epoch provides valuable constraints on their early growth. High-redshift quasars typically have a $\mathcal{M}_{\rm BH}$/$\mathcal{M}_\star$ ratio significantly elevated in comparison to the local value. However, the degree to which this apparent offset is driven by observational biases is unclear for the most distant quasars. To address this issue, we model the sample selection and measurement biases for a compilation of 20 quasars at $z\sim6$ with host properties based on ALMA observations. We find that the observed distribution of quasars in the $\mathcal{M}_{\rm BH} - \mathcal{M}_\star$ plane can be reproduced by assuming that the underlying SMBH population at $z\sim6$ follows the relationship in the local universe. However, a positive or even a negative evolution in $\mathcal{M}_{\rm BH}$/$\mathcal{M}_\star$ can also explain the data, depending on whether the intrinsic scatter evolves and the strength of various systematic uncertainties. To break these degeneracies, an improvement in the accuracy of mass measurements and an expansion of the current sample to lower $\mathcal{M}_{\rm BH}$ limits are needed. Furthermore, assuming a radiative efficiency of 0.1 and quasar duty cycles estimated from the active SMBH fraction, significant outliers in $\mathcal{M}_{\rm BH}$/$\mathcal{M}_\star$ tend to move toward the local relation given their instantaneous BH growth rate and star formation rate. This may provide evidence for a self-regulated SMBH-galaxy coevolution scenario that is in place at $z\sim6$, with AGN feedback being a possible driver.

preprint2022arXiv

Origin and evolution of dust-obscured galaxies in galaxy mergers

Dust Obscured Galaxies (DOGs), which are observationally characterized as faint in the optical and bright in the infrared, are the final stage of galaxy mergers and are essential objects in the evolution of galaxies and active galactic nuclei (AGNs). However, the relationship between torus-scale gas dynamics around AGNs and DOGs lifetime remain unclear. We obtained evolution of the spectral energy distributions (SEDs) of a galaxy merger system with AGN feedback, from post-processed pseudo-observations based on an N-body/Smoothed Particle Hydrodynamics (SPH) simulation. We focused on a late stage merger of two identical galaxies with a supermassive black hole (SMBH) of 10$^8$ M$_\odot$. We found that the infrared luminosity of the system reaches ultra- and hyper-luminous infrared galaxy classes (10$^{12}$ and 10$^{13}$ L$_\odot$, respectively). The DOGs phase corresponds to a state in which the AGNs are buried in dense gas and dust, with the infrared luminosity exceeding 3.3 $\times$ 10$^{12}$ L$_\odot$. We also identified the sub-categories of DOGs, namely bump and power-law DOGs from the SEDs and their evolution. The bump DOGs tend to evolve to power-law DOGs on several Myrs. We found that contribution from the hot dust around the nucleus in the infrared radiation is essential for identifying the system as a power-law DOG; the gas and dust distribute non-spherically around the nucleus, therefore, the observed properties of DOGs depend on the viewing angle. In our model, the lifetime of merger-driven DOGs is less than 4 Myrs, suggesting that the observed DOGs phase is a brief aspect of galaxy mergers.

preprint2022arXiv

Passive spiral galaxies deeply captured by Subaru Hyper Suprime-Cam

This paper presents a thousand passive spiral galaxy samples at $z=$ 0.01-0.3 based on a combined analysis of the Third Public Data Release of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP PDR3) and the GALEX-SDSS-WISE Legacy Catalog (GSWLC-2). Among 54871 $gri$ galaxy cutouts taken from the HSC-SSP PDR3 over 1072 deg$^2$, we conducted a search with deep-learning morphological classification for candidates of passive spirals below the star-forming main sequence derived by UV to mid-IR SED fitting in the GSWLC-2. We then classified the candidates into 1100 passive spirals and 1141 secondary samples based on visual inspections. Most of the latter cases are considered to be passive ringed S0 or pseudo-ringed galaxies. The remainder of these secondary samples has ambiguous morphologies, including two peculiar objects with diamond-shaped stellar wings. The selected passive spirals have a similar distribution to the general quiescent galaxies on the EW$_\mathrm{Hδ}$-D$_n$4000 diagram and concentration indices. Moreover, we detected an enhanced passive fraction of spiral galaxies in X-ray clusters. Passive spirals in galaxy clusters are preferentially located in the midterm or late infall phase on the phase-space diagram, supporting the ram pressure scenario, which has been widely advocated in previous studies. The source catalog and $gri$-composite images are available on the HSC-SSP PDR3 website (https://hsc.mtk.nao.ac.jp/ssp/data-release/). Future updates, including integration with a citizen science project dedicated to the HSC data, will achieve more effective and comprehensive classifications.

preprint2022arXiv

Where's Swimmy?: Mining unique color features buried in galaxies by deep anomaly detection using Subaru Hyper Suprime-Cam data

We present the Swimmy (Subaru WIde-field Machine-learning anoMalY) survey program, a deep-learning-based search for unique sources using multicolored ($grizy$) imaging data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). This program aims to detect unexpected, novel, and rare populations and phenomena, by utilizing the deep imaging data acquired from the wide-field coverage of the HSC-SSP. This article, as the first paper in the Swimmy series, describes an anomaly detection technique to select unique populations as "outliers" from the data-set. The model was tested with known extreme emission-line galaxies (XELGs) and quasars, which consequently confirmed that the proposed method successfully selected 60-70% of the quasars and 60% of the XELGs without labeled training data. In reference to the spectral information of local galaxies at $z=$0.05-0.2 obtained from the Sloan Digital Sky Survey, we investigated the physical properties of the selected anomalies and compared them based on the significance of their outlier values. The results revealed that XELGs constitute notable fractions of the most anomalous galaxies, and certain galaxies manifest unique morphological features. In summary, a deep anomaly detection is an effective tool that can search rare objects, and ultimately, unknown unknowns with large data-sets. Further development of the proposed model and selection process can promote the practical applications required to achieve specific scientific goals.

preprint2021arXiv

A Wide and Deep Exploration of Radio Galaxies with Subaru HSC (WERGS). VI. Distant Filamentary Structures Pointed by High-z Radio Galaxies at z~4

We present the environmental properties around high-$z$ radio galaxies (HzRGs) at $z\sim4$, which have been poorly investigated because of their rarity. We use the largest samples of HzRGs and $g$-dropout galaxy overdense regions at $z\sim4$, which were constructed from Hyper Suprime-Cam Subaru Strategic Program, to characterize the HzRG environments statistically. We measure the $g$-dropout galaxy overdensities around 21 HzRGs whose rest-frame 1.4 GHz radio luminosities ($L_{1.4\mathrm{GHz}}$) are $10^{26-27}$ W Hz$^{-1}$. We find that the overdensities around the faint HzRGs with $L_{1.4\mathrm{GHz}}\sim10^{26.0-26.5}$ W Hz$^{-1}$ tend to be higher than that of the $g$-dropout galaxies. On the other hand, no significant difference of density environments is found between the luminous HzRGs with $ L_{1.4\mathrm{GHz}}\sim10^{26.5-27.0} $ W Hz$^{-1}$ and the $g$-dropout galaxies. The HzRGs are found to occupy more massive halos than $g$-dropout galaxies through a cross-correlation between the HzRGs and $g$-dropout galaxies. This trend is more pronounced in the faint HzRGs. These results are consistent with a scenario where HzRGs get older and more massive as the radio-luminosity decreases. The HzRGs are expected to trace the progenitors of local cluster halos from their calculated halo mass. In addition, we find that surrounding galaxies tend to distribute along the radio-jet major axis of the HzRGs at angular distances less than $\lesssim500$ physical kpc. Our findings imply the onset of the filamentary structures around the HzRGs at $z\sim4$.

preprint2021arXiv

An Active Galactic Nucleus Recognition Model based on Deep Neural Network

To understand the cosmic accretion history of supermassive black holes, separating the radiation from active galactic nuclei (AGNs) and star-forming galaxies (SFGs) is critical. However, a reliable solution on photometrically recognising AGNs still remains unsolved. In this work, we present a novel AGN recognition method based on Deep Neural Network (Neural Net; NN). The main goals of this work are (i) to test if the AGN recognition problem in the North Ecliptic Pole Wide (NEPW) field could be solved by NN; (ii) to shows that NN exhibits an improvement in the performance compared with the traditional, standard spectral energy distribution (SED) fitting method in our testing samples; and (iii) to publicly release a reliable AGN/SFG catalogue to the astronomical community using the best available NEPW data, and propose a better method that helps future researchers plan an advanced NEPW database. Finally, according to our experimental result, the NN recognition accuracy is around 80.29% - 85.15%, with AGN completeness around 85.42% - 88.53% and SFG completeness around 81.17% - 85.09%.

preprint2021arXiv

Galaxy Morphologies Revealed with Subaru HSC and Super-Resolution Techniques I: Major Merger Fractions of L_UV~3-15 L_UV* Dropout Galaxies at z~4-7

We perform a super-resolution analysis of the Subaru Hyper Suprime-Cam (HSC) images to estimate the major merger fractions of z~4-7 dropout galaxies at the bright end of galaxy UV luminosity functions (LFs). Our super-resolution technique improves the spatial resolution of the ground-based HSC images, from ~1" to <~0."1, which is comparable to that of the Hubble Space Telescope, allowing us to identify z~4-7 bright major mergers at a high completeness value of >~90%. We apply the super-resolution technique to 6412, 16, 94, and 13 very bright dropout galaxies at z~4, 5, 6, and 7, respectively, in a UV luminosity range of L_UV~3-15 L_UV* corresponding to -24<~M_UV<~-22. The major merger fractions are estimated to be f_merger~10-20% at z~4 and ~50-70% at z~5-7, which shows no f_merger difference compared to those of a control faint galaxy sample. Based on the f_merger estimates, we verify contributions of source blending effects and major mergers to the bright-end of double power-law (DPL) shape of z~4-7 galaxy UV LFs. While these two effects partly explain the DPL shape at L_UV~3-10 L_UV*, the DPL shape cannot be explained at the very bright end of L_UV>~10 L_UV*, even after the AGN contribution is subtracted. The results support scenarios in which other additional mechanisms such as insignificant mass quenching and low dust obscuration contribute to the DPL shape of galaxy UV LFs.

preprint2021arXiv

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) XII. Extended [C II] Structure (Merger or Outflow) in a z = 6.72 Red Quasar

We present ALMA [C II] 158 $μ$m line and far-infrared (FIR) continuum emission observations toward HSC J120505.09$-$000027.9 (J1205$-$0000) at $z = 6.72$ with the beam size of $\sim 0''.8 \times 0''.5$ (or 4.1 kpc $\times$ 2.6 kpc), the most distant red quasar known to date. Red quasars are modestly reddened by dust, and are thought to be in rapid transition from an obscured starburst to an unobscured normal quasar, driven by powerful active galactic nucleus (AGN) feedback which blows out a cocoon of interstellar medium (ISM). The FIR continuum of J1205$-$0000 is bright, with an estimated luminosity of $L_{\rm FIR} \sim 3 \times 10^{12}~L_\odot$. The [C II] line emission is extended on scales of $r \sim 5$ kpc, greater than the FIR continuum. The line profiles at the extended regions are complex and broad (FWHM $\sim 630-780$ km s$^{-1}$). Although it is not practical to identify the nature of this extended structure, possible explanations include (i) companion/merging galaxies and (ii) massive AGN-driven outflows. For the case of (i), the companions are modestly star-forming ($\sim 10~M_\odot$ yr$^{-1}$), but are not detected by our Subaru optical observations ($y_{\rm AB,5σ} = 24.4$ mag). For the case of (ii), our lower-limit to the cold neutral outflow rate is $\sim 100~M_\odot$ yr$^{-1}$. The outflow kinetic energy and momentum are both much smaller than what predicted in energy-conserving wind models, suggesting that the AGN feedback in this quasar is not capable of completely suppressing its star formation.

preprint2021arXiv

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). XVI. 69 New Quasars at 5.8 < z < 7.0

We present the spectroscopic discovery of 69 quasars at 5.8 < z < 7.0, drawn from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) imaging survey data. This is the 16th publication from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, and completes identification of all but the faintest candidates (i.e., i-band dropouts with zAB < 24 and y-band detections, and z-band dropouts with yAB < 24) with Bayesian quasar probability Pq > 0.1 in the HSC-SSP third public data release (PDR3). The sample reported here also includes three quasars with Pq < 0.1 at z ~ 6.6, which we selected in an effort to completely cover the reddest point sources with simple color cuts. The number of high-z quasars discovered in SHELLQs has now grown to 162, including 23 type-II quasar candidates. This paper also presents identification of seven galaxies at 5.6 < z < 6.7, an [O III] emitter at z = 0.954, and 31 Galactic cool stars and brown dwarfs. High-z quasars and galaxies comprise 75 % and 16 % respectively of all the spectroscopic SHELLQs objects that pass our latest selection algorithm with the PDR3 photometry. That is, a total of 91 % of the objects lie at z > 5.6. This demonstrates that the algorithm has very high efficiency, even though we are probing an unprecedentedly low-luminosity population down to M1450 ~ -21 mag.

preprint2021arXiv

The eROSITA Final Equatorial-Depth Survey (eFEDS): A multiwavelength view of WISE mid-infrared galaxies/active galactic nuclei

We investigate the physical properties--such as the stellar mass, SFR, IR luminosity, X-ray luminosity, and hydrogen column density--of MIR galaxies and AGN at $z < 4$ in the 140 deg$^2$ field observed by SRG/eROSITA through the eFEDS survey. By cross-matching the WISE 22 $μ$m (W4)-detected sample and the eFEDS X-ray point-source catalog, we find that 692 extragalactic objects are detected by eROSITA. We have compiled a multiwavelength dataset. We have also performed (i) an X-ray spectral analysis, (ii) SED fitting using X-CIGALE, (iii) 2D image-decomposition analysis using Subaru HSC images, and (iv) optical spectral fitting with QSFit to investigate the AGN and host-galaxy properties. For 7,088 WISE W4 objects that are undetected by eROSITA, we have performed an X-ray stacking analysis to examine the typical physical properties of these X-ray faint and/or probably obscured objects. We find that (i) 82% of the eFEDS-W4 sources are classified as X-ray AGN with $\log\,L_{\rm X} >$ 42 erg s$^{-1}$; (ii) 67% and 24% of the objects have $\log\,(L_{\rm IR}/L_{\odot}) > 12$ and 13, respectively; (iii) the relationship between $L_{\rm X}$ and the 6 $μ$m luminosity is consistent with that reported in previous works; and (iv) the relationship between the Eddington ratio and $N_{\rm H}$ for the eFEDS-W4 sample and a comparison with a model prediction from a galaxy-merger simulation indicates that approximately 5% of the eFEDS-W4 sources in our sample are likely to be in an AGN-feedback phase, in which strong radiation pressure from the AGN blows out the surrounding material from the nuclear region. Thanks to the wide area coverage of eFEDS, we have been able to constrain the ranges of the physical properties of the WISE W4 sample of AGNs at $z < 4$, providing a benchmark for forthcoming studies on a complete census of MIR galaxies selected from the full-depth eROSITA all-sky survey.

preprint2021arXiv

The eROSITA Final Equatorial-Depth Survey (eFEDS): The AGN Catalogue and its X-ray Spectral Properties

Context: After the successful launch of the Spectrum-Roentgen-Gamma (SRG) mission in July 2019, eROSITA, the soft X-ray instrument aboard SRG, performed scanning observations of a large contiguous field, namely the eROSITA Final Equatorial Depth Survey (eFEDS), ahead of the planned four-year all-sky survey. eFEDS yielded a large sample of X-ray sources with very rich multi-band photometric and spectroscopic coverage. Aims: We present here the eFEDS Active Galactic Nuclei (AGN) catalog and the eROSITA X-ray spectral properties of the eFEDS sources. Methods: Using a Bayesian method, we perform a systematic X-ray spectral analysis for all eFEDS sources. The appropriate model is chosen based on the source classification and the spectral quality, and, in the case of AGN, including the possibility of intrinsic (rest-frame) absorption and/or soft excess emission. Hierarchical Bayesian modeling (HBM) is used to estimate the spectral parameter distribution of the sample. Results: X-ray spectral properties are presented for all eFEDS X-ray sources. There are 21952 candidate AGN, which comprise 79% of the eFEDS sample. Despite a large number of faint sources with low photon counts, our spectral fitting provides meaningful measurements of fluxes, luminosities, and spectral shapes for a majority of the sources. This AGN catalog is dominated by X-ray unobscured sources, with an obscured (logNH>21.5) fraction of 10% derived by HBM. The power-law slope of the catalog can be described by a Gaussian distribution of 1.94+-0.22. Above a photon counts threshold of 500, nine out of 50 AGN have soft excess detected. For the sources with blue UV to optical color (type-I AGN), the X-ray emission is well correlated with the UV emission with the usual anti-correlation between the X-ray to UV spectral slope α_{OX} and the UV luminosity.

preprint2021arXiv

Third Data Release of the Hyper Suprime-Cam Subaru Strategic Program

The paper presents the third data release of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), a wide-field multi-band imaging survey with the Subaru 8.2m telescope. HSC-SSP has three survey layers (Wide, Deep, and UltraDeep) with different area coverages and depths, designed to address a wide array of astrophysical questions. This third release from HSC-SSP includes data from 278 nights of observing time and covers about 670 square degrees in all five broad-band filters at the full depth ($\sim26$~mag at $5σ$) in the Wide layer. If we include partially observed area, the release covers 1,470 square degrees. The Deep and UltraDeep layers have $\sim80\%$ of the originally planned integration times, and are considered done, as we have slightly changed the observing strategy in order to compensate for various time losses. There are a number of updates in the image processing pipeline. Of particular importance is the change in the sky subtraction algorithm; we subtract the sky on small scales before the detection and measurement stages, which has significantly reduced false detections. Thanks to this and other updates, the overall quality of the processed data has improved since the previous release. However, there are limitations in the data (for example, the pipeline is not optimized for crowded fields), and we encourage the user to check the quality assurance plots as well as a list of known issues before exploiting the data. The data release website is https://hsc-release.mtk.nao.ac.jp/.

preprint2020arXiv

Dual supermassive black holes at close separation revealed by the Hyper Suprime-Cam Subaru Strategic Program

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

preprint2020arXiv

Nature of Compton-thick Active Galactic Nuclei in "Non-merging'' Luminous Infrared Galaxies UGC 2608 and NGC 5135 Revealed with Broadband X-ray Spectroscopy

We have analyzed the broadband X-ray spectra of active galactic nuclei (AGNs) in two "non-merging" luminous infrared galaxies (LIRGs) UGC 2608 and NGC 5135, utilizing the data of NuSTAR, Suzaku, XMM-Newton, and Chandra. Applying the X-ray clumpy-torus model (XCLUMPY: Tanimoto et al. 2019), we find that both sources have similar spectra characterized by Compton-thick (CT) absorption ($N_{\rm H} \sim$ 5-7 $\times$ $10^{24}$ cm$^{-2}$) and small torus angular width ($σ$ $<$ 20$^{\circ}$). The intrinsic 2-10 keV luminosities are $3.9^{+2.2}_{-1.7}$ $\times$ $10^{43}$ erg s$^{-1}$ (UGC 2608) and $2.0^{+3.3}_{-1.0}$ $\times$ $10^{43}$ erg s$^{-1}$ (NGC 5135). The [O IV]-to-nuclear-12 $μ$m luminosity ratios are larger than those of typical Seyferts, which are consistent with the torus covering factors ($C_{\rm T} \lesssim$ 0.7) estimated from the torus angular widths and column densities by X-ray spectroscopy. The torus covering factors and Eddington ratios ($λ_{\rm Edd} \sim$ 0.1) follow the relation found by Ricci et al. (2017c) for local AGNs, implying that their tori become geometrically thin due to significant radiation pressure of the AGN that blows out some part of the tori. These results indicate that the CT AGNs in these "non-merger" LIRGs are just a normal AGN population seen edge-on through a large line-of-sight column density. They are in contrast to the buried CT AGNs in late-stage mergers that have large torus covering factors even at large Eddington ratios.

preprint2020arXiv

Photometric Redshifts in the North Ecliptic Pole Wide Field based on a Deep Optical Survey with Hyper Suprime-Cam

The $AKARI$ space infrared telescope has performed near- to mid-infrared (MIR) observations on the North Ecliptic Pole Wide (NEPW) field (5.4 deg$^2$) for about one year. $AKARI$ took advantage of its continuous nine photometric bands, compared with NASA's $Spitzer$ and WISE space telescopes, which had only four filters with a wide gap in the MIR. The $AKARI$ NEPW field lacked deep and homogeneous optical data, limiting the use of nearly half of the IR sources for extra-galactic studies owing to the absence of photometric redshifts (photo-zs). To remedy this, we have recently obtained deep optical imaging over the NEPW field with 5 bands ($g$, $r$, $i$, $z$, and $Y$) of the Hyper Suprime-Camera (HSC) on the Subaru 8m telescope. We optically identify AKARI-IR sources along with supplementary $Spitzer$ and WISE data as well as pre-existing optical data. In this work, we derive new photo-zs using a $χ^2$ template-fitting method code ($Le$ $Phare$) and reliable photometry from 26 selected filters including HSC, $AKARI$, CFHT, Maidanak, KPNO, $Spitzer$ and WISE data. We take 2026 spectroscopic redshifts (spec-z) from all available spectroscopic surveys over the NEPW to calibrate and assess the accuracy of the photo-zs. At z < 1.5, we achieve a weighted photo-z dispersion of $σ_{Δ{z/(1+z)}}$ = 0.053 with $η$ = 11.3% catastrophic errors.

preprint2020arXiv

Properties of environment around AGN and luminous galaxy pairs through HSC wide survey

We investigated the properties of AGN environments, particularly environments where the association of luminous galaxies (LGs) is found within 4~Mpc from AGNs with redshifts of 0.8 -- 1.1. For comparison, three additional AGN environments, (namely, AGNs of all types, type~1 AGNs with X-ray and/or radio detection, and type~2 AGNs) and an environment of blue $M_{*}$, characteristic luminosity of the Schechter function, galaxies were investigated. The cross-correlation function with the surrounding galaxies was measured and compared between the AGN and blue galaxy samples. We also compared the distributions of color, absolute magnitude, and stellar mass of the galaxies around such target objects. The properties of clusters detected using surrounding galaxies selected based on a photometric redshift were examined and compared for different samples. The target AGNs were drawn from the Million Quasars (MILLIQUAS) catalog, and the blue galaxies were drawn from six redshift survey catalogs (SDSS, WiggleZ, DEEP2, VVDS, VIPERS, and PRIMUS). The galaxies used as a measure of the environment around the targets are drawn from S18a internal data released by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We found that, among the five AGN and blue galaxy samples considered, the environment of AGN-LG pairs is the most enriched with luminous galaxies. We also found an enhancement in the number of mass-selected clusters in the AGN-LG pair sample against those in the other samples. The results obtained in this study indicate that existence of multiple clusters is the major driver in the association of AGNs and LGs, rather than a single large-mass dark matter halo hosting the AGN.

preprint2020arXiv

SCUBA-2 Ultra Deep Imaging Eao Survey (Studies) III: Multi-wavelength properties, luminosity functions and preliminary source catalog of 450-$μ$m-selected galaxies

We construct a SCUBA-2 450-$μ$m map in the COSMOS field that covers an area of 300 arcmin$^{2}$ and reaches a 1$σ$ noise level of 0.65 mJy in the deepest region. We extract 256 sources detected at 450 $μ$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-$μ$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<3$. However, we find a moderate correlation where dust temperature increases with the deviation from the SFR-$M_{\ast}$ relation. The increase in dust temperature also correlates with optical morphology, which is consistent with merger-triggered starbursts in sub-millimeter galaxies. Our galaxies do not show the tight IRX-$β_{\rm UV}$ correlation that has been observed in the local Universe. We construct the infrared luminosity functions of our 450-$μ$m sources and measure their comoving SFR densities. The contribution of the $L_{\rm IR}$ > $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$.

preprint2020arXiv

SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) IV: Spatial clustering and halo masses of 450-$μ$m-selected sub-millimeter galaxies

We analyze an extremely deep 450-$μ$m image ($1σ=0.56$\,mJy\,beam$^{-1}$) of a $\simeq 300$\,arcmin$^{2}$ area in the CANDELS/COSMOS field as part of the SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). We select a robust (signal-to-noise ratio $\geqslant 4$) and flux-limited ($\geqslant 4$\,mJy) sample of 164 sub-millimeter galaxies (SMGs) at 450-$μ$m that have $K$-band counterparts in the COSMOS2015 catalog identified from radio or mid-infrared imaging. Utilizing this SMG sample and the 4705 $K$-band-selected non-SMGs that reside within the noise level $\leqslant 1$\,mJy\,beam$^{-1}$ region of the 450-$μ$m image as a training set, we develop a machine-learning classifier using $K$-band magnitude and color-color pairs based on the thirteen-band photometry available in this field. We apply the trained machine-learning classifier to the wider COSMOS field (1.6\,deg$^{2}$) using the same COSMOS2015 catalog and identify a sample of 6182 450-$μ$m SMG candidates with similar colors. The number density, radio and/or mid-infrared detection rates, redshift and stellar mass distributions, and the stacked 450-$μ$m fluxes of these SMG candidates, from the S2COSMOS observations of the wide field, agree with the measurements made in the much smaller CANDELS field, supporting the effectiveness of the classifier. Using this 450-$μ$m SMG candidate sample, we measure the two-point autocorrelation functions from $z=3$ down to $z=0.5$. We find that the 450-$μ$m SMG candidates reside in halos with masses of $\simeq (2.0\pm0.5) \times10^{13}\,h^{-1}\,\rm M_{\odot}$ across this redshift range. We do not find evidence of downsizing that has been suggested by other recent observational studies.

preprint2020arXiv

Search for Optically Dark Infrared Galaxies without Counterparts of Subaru Hyper Suprime-Cam in the AKARI North Ecliptic Pole Wide Survey Field

We present the physical properties of AKARI sources without optical counterparts in optical images from the Hyper Suprime-Cam (HSC) on the Subaru telescope. Using the AKARI infrared (IR) source catalog and HSC optical catalog, we select 583 objects that do not have HSC counterparts in the AKARI North Ecliptic Pole (NEP) wide survey field ($\sim 5$ deg$^{2}$). Because the HSC limiting magnitude is deep ($g_{\rm AB}$ $\sim 28.6$), these are good candidates for extremely red star-forming galaxies (SFGs) and/or active galactic nuclei (AGNs), possibly at high redshifts. We compile multi-wavelength data out to 500 $μ$m and use it for Spectral Energy Distribution (SED) fitting with CIGALE to investigate the physical properties of AKARI galaxies without optical counterparts. We also compare their physical quantities with AKARI mid-IR selected galaxies with HSC counterparts. The estimated redshifts of AKARI objects without HSC counterparts range up to $z\sim 4$, significantly higher than that of AKARI objects with HSC counterparts. We find that: (i) 3.6 $-$ 4.5 $μ$m color, (ii) AGN luminosity, (iii) stellar mass, (iv) star formation rate, and (v) $V$-band dust attenuation in the interstellar medium of AKARI objects without HSC counterparts are systematically larger than those of AKARI objects with counterparts. These results suggest that our sample includes luminous, heavily dust-obscured SFGs/AGNs at $z\sim 1-4$ that are missed by previous optical surveys, providing very interesting targets for the coming James Webb Space Telescope era.

preprint2020arXiv

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). IX. Identification of Two Red Quasars at z > 5.6

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 spectroscopically confirmed 93 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 um (rest-frame ~5000-6500 A). 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.

preprint2020arXiv

Subaru Hyper Suprime-Cam View of Quasar Host Galaxies at z < 1

Active galactic nuclei (AGNs) are key for understanding the coevolution of galaxies and supermassive black holes (SMBHs). AGN activity is thought to affect the properties of their host galaxies, via a process called "AGN feedback", which drives the co-evolution. From a parent sample of 1151 z < 1 type-1 quasars from the Sloan Digital Sky Survey quasar catalog, we detected host galaxies of 862 of them in the high-quality grizy images of the Subaru Hyper Suprime-Cam (HSC) survey. The unprecedented combination of the survey area and depth allows us to perform a statistical analysis of the quasar host galaxies, with small sample variance. We fit the radial image profile of each quasar as a linear combination of the point spread function and the Sersic function, decomposing the images into the quasar nucleus and the host galaxy components. We found that the host galaxies are massive, with stellar mass Mstar > 10^(10) Msun, and are mainly located on the green valley. This trend is consistent with a scenario in which star formation of the host galaxies is suppressed by AGN feedback, that is, AGN activity may be responsible for the transition of these galaxies from the blue cloud to the red sequence. We also investigated the SMBH mass to stellar mass relation of the z < 1 quasars, and found a consistent slope with the local relation, while the SMBHs may be slightly undermassive. However, the above results are subject to our sample selection, which biases against host galaxies with low masses and/or large quasar-to-host flux ratios.

preprint2020arXiv

Tracing the co-evolution path of super massive black holes and spheroids with AKARI-selected ultra-luminous IR galaxies at intermediate redshifts

We present the stellar population and ionized-gas outflow properties of ultra-luminous IR galaxies (ULIRGs) at $z=$ 0.1-1.0, which are selected from AKARI FIR all-sky survey. We construct a catalog of 1077 ULIRGs to examine feedback effect after major mergers. 202 out of the 1077 ULIRGs are spectroscopically identified by SDSS and Subaru/FOCAS observations. Thanks to deeper depth and higher resolution of AKARI compared to the previous IRAS survey, and reliable identification from WISE MIR pointing, the sample is unique in identifying optically-faint (i$\sim$20) IR-bright galaxies, which could be missed in previous surveys. A self-consistent spectrum-SED decomposition method, which constrains stellar population properties in SED modeling based on spectral fitting results, has been employed for 149 ULIRGs whose optical continua are dominated by host galaxies. They are massive galaxies ($M_{\rm star}\sim10^{11}$-$10^{12}$ M$_{\odot}$), associated with intense star formation activities (SFR $\sim$ 200-2000 M$_{\odot}$ yr$^{-1}$). The sample covers a range of AGN bolometric luminosity of $10^{10}$-$10^{13}$ L$_{\odot}$, and the outflow velocity measured from [OIII] 5007A line shows a correlation with AGN luminosity. Eight galaxies show extremely fast outflows with velocity up to 1500-2000 km s$^{-1}$. However, the co-existence of vigorous starbursts and strong outflows suggests the star formation has not been quenched during the ULIRG phase. By deriving stellar mass and mass fraction of young stellar population, we find no significant discrepancies between stellar properties of ULIRGs with weak and powerful AGNs. The results are not consistent with the merger-induced evolutionary scenario, which predicts that SF-dominated ULIRGs show smaller stellar mass and younger stellar populations compared to AGN-dominated ULIRGs.

preprint2020arXiv

Wide and Deep Exploration of Radio Galaxies with Subaru HSC (WERGS). III. Discovery of a z = 4.72 Radio Galaxy with Lyman Break Technique

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$α$ 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$α$ 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: $α^{325}_{1400}$ of $-1.1$ and $α^{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.

preprint2019arXiv

NuSTAR Discovery of a Compton-thick Dust-obscured Galaxy WISE J0825+3002

We report the discovery of a Compton-thick (CT) dust-obscured galaxy (DOG) at $z$ = 0.89, WISE J082501.48+300257.2 (WISE0825+3002), observed by Nuclear Spectroscopic Telescope Array (NuSTAR). X-ray analysis with the XCLUMPY model revealed that hard X-ray luminosity in the rest-frame 2-10 keV band of WISE0825+3002 is $L_{\rm X}$ (2-10 keV) = $4.2^{+2.8}_{-1.6} \times 10^{44}$ erg s$^{-1}$ while its hydrogen column density is $N_{\rm H}$ = $1.0^{+0.8}_{-0.4} \times 10^{24}$ cm$^{-2}$, indicating that WISE0825+3002 is a mildly CT active galactic nucleus (AGN). We performed the spectral energy distribution (SED) fitting with CIGALE to derive its stellar mass, star formation rate, and infrared luminosity. The estimated Eddington ratio based on stellar mass and integration of the best-fit SED of AGN component is $λ_{\rm Edd}$ = 0.70, which suggests that WISE0825+3002 harbors an actively growing black hole behind a large amount of gas and dust. We found that the relationship between luminosity ratio of X-ray and 6 $μ$m, and Eddington ratio follows an empirical relation for AGNs reported by Toba et al. (2019a).

preprint2019arXiv

SOFIA/HAWC+ View of an Extremely Luminous Infrared Galaxy, WISE1013+6112

We present far-infrared (FIR) properties of an extremely luminous infrared galaxy (ELIRG) at $z_{\rm spec}$ = 3.703, WISE J101326.25+611220.1 (WISE1013+6112). This ELIRG is selected as an IR-bright dust-obscured galaxy (DOG) based on the photometry from the Sloan digital sky survey (SDSS) and wide-field infrared survey explorer (WISE). In order to derive its accurate IR luminosity, we perform follow-up observations at 89 and 154 $μ$m using the high-resolution airborne wideband camera-plus (HAWC+) on board the 2.7-m stratospheric observatory for infrared astronomy (SOFIA) telescope. We conduct spectral energy distribution (SED) fitting with CIGALE using 15 photometric data (0.4-1300 $μ$m). We successfully pin down FIR SED of WISE1013+6112 and its IR luminosity is estimated to be $L_{\rm IR}$ = (1.62 $\pm$ 0.08) $\times 10^{14}$ $L_{\odot}$, making it one of the most luminous IR galaxies in the universe. We determine the dust temperature of WISE1013+6112 is $T_{\rm dust}$ = 89 $\pm$ 3 K, which is significantly higher than that of other populations such as SMGs and FIR-selected galaxies at similar IR luminosities. The resultant dust mass is $M_{\rm dust}$ = (2.2 $\pm$ 0.1) $\times 10^{8}$ $M_{\odot}$. This indicates that WISE1013+6112 has a significant active galactic nucleus (AGN) and star-forming activity behind a large amount of dust.

preprint2019arXiv

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) VIII. A less biased view of the early co-evolution of black holes and host galaxies

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.

preprint2016arXiv

Near-infrared spectroscopy of nearby Seyfert galaxies: Is there evidence for the shock excitation in narrow-line regions?

One of the unsettled, important problems in active galactic nuclei (AGNs) is the major ionization mechanism of gas clouds in their narrow-line regions (NLRs). In order to investigate this issue, we present our $J$-band spectroscopic observations for a sample of 26 nearby Seyfert galaxies. In our study, we use the flux ratio of the following two forbidden emission lines, [Fe II]1.257$μ$m and [P II]1.188$μ$m because it is known that this ratio is sensitive to the ionization mechanism. We obtain the [Fe II]/[P II] flux ratio or its lower limit for 19 objects. In addition to our data, we compile this flux ratio (or its lower limit) of 23 nearby Seyfert galaxies from the literature. Based on the collected data, we find that three Seyfert galaxies show very large lower limits of [Fe II]/[P II] flux ratios ($\gtrsim$10); NGC 2782, 5005, and Mrk 463. It is thus suggested that the contribution of the fast shock in the gas excitation is significantly large for them. However, more than half of Seyfert galaxies in our sample show moderate [Fe II]/[P II] flux ratios ($\sim$2), that is consistent to the pure photoionization by power-law ionizing continuum emission. We also find that the [Fe II]/[P II] flux ratio shows no clear correlation with the radio loudness, suggesting that the radio jet is not the primary origin of shocks in NLRs of Seyfert galaxies.

preprint2016arXiv

Revisiting the Completeness and the Luminosity Function in High-Redshift Low-Luminosity Quasar Surveys

Recent studies have derived quasar luminosity functions (QLFs) at various redshifts. However, the faint side of the QLF at high redshifts is still too uncertain. An accurate estimate of the survey completeness is essential to derive an accurate QLF for use in studying the luminosity-dependent density evolution of the quasar population. Here we investigate how the luminosity dependence of quasar spectra (the Baldwin effect) and the attenuation model for the inter-galactic medium (IGM) affect the completeness estimates. For this purpose, we revisit the completeness of quasar surveys specifically at $z\sim4-5$, using the COSMOS images observed with Subaru/Suprime-Cam. As the result, we find that the completeness estimates are sensitive to the luminosity dependence of the quasar spectrum and difference in the IGM attenuation models. At $z\sim4$, the number density of quasars when we adopt the latest IGM model and take the luminosity dependence of spectra into account are $(3.49\pm1.62)\times10^{-7}$ Mpc$^{-3}$ mag$^{-1}$ for $-24.09<M_{1450}<-23.09$ and $(5.24\pm2.13)\times10^{-7}$ Mpc$^{-3}$ mag$^{-1}$ for $-23.09<M_{1450}<-22.09$, respectively, which are $\sim$$24$$\%$ lower than that estimated by the conventional method. On the other hand, at $z\sim5$, a $1σ$ confidence upper limit of the number density at $-24.5<M_{1450}<-22.5$ in our new estimates is $\sim$$43$$\%$ higher than that estimated previously. The results suggest that the luminosity dependence of the quasar spectrum and the IGM model are important for deriving accurate number density of high-$z$ quasars. Even taking these effects into account, the inferred luminosity-dependent density evolution of quasars is consistent with the AGN down-sizing evolutionary picture.

preprint2016arXiv

Search for Hyper Infrared-Luminous Dust Obscured Galaxies selected with WISE and SDSS

We aim to search for hyperliminous infrared (IR) galaxies (HyLIRGs) with IR luminosity $L_{\rm IR}$ $>$ 10$^{13}$ $L_{\odot}$ by applying the selection method of Dust Obscured Galaxies (DOGs). They are spatially rare but could correspond to a maximum phase of cosmic star formation and/or active galactic nucleus (AGN) activity, hence they are a crucial population for understanding the star formation and mass assembly history of galaxies. Combining the optical and IR catalogs obtained from Sloan Digital Sky Survey (SDSS) and Wide-field Infrared Survey Explorer (WISE), we performed the extensive HyLIRGs survey; we selected 5,311 IR-bright DOGs with $i$ -- [22] $>$ 7.0 and flux at 22 $μ$m $>$ 3.8 mJy in 14,555 deg$^2$, where $i$ and [22] are $i$-band and 22 $μ$m AB magnitudes, respectively. Among them, 67 DOGs have reliable spectroscopic redshifts that enable us to estimate their total IR luminosity based on the SED fitting. Consequently, we successfully discovered 24 HyLIRGs among the 67 spectroscopically-confirmed DOGs. We found that (i) $i$ - [22] color of IR-bright DOGs correlates with the total IR luminosity and (ii) surface number density of HyLIRGs is $>$ 0.17 deg$^{-2}$. A high fraction ($\sim$ 73%) of IR-bright DOGs with $i$ - [22] $>$ 7.5 shows $L_{\rm IR}$ $>$ 10$^{13}$ $L_{\odot}$, and the DOGs criterion we adopted could be independently-effective against "W1W2-dropout method" based on four WISE bands, for searching hyper IR luminous populations of galaxies.

preprint2016arXiv

Subaru high-z exploration of low-luminosity quasars (SHELLQs). I. Discovery of 15 quasars and bright galaxies at 5.7 < z < 6.9

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

preprint2016arXiv

X-ray and infrared diagnostics of nearby active galactic nuclei with MAXI and AKARI

Nearby active galactic nuclei were diagnosed in the X-ray and mid-to-far infrared wavelengths, with Monitor of All-sky X-ray Image (MAXI) and the Japanese infrared observatory AKARI, respectively. Among the X-ray sources listed in the second release of the MAXI all-sky X-ray source catalog, 100 ones are currently identified as a non-blazar-type active galactic nucleus. These include 95 Seyfert galaxies and 5 quasars, and they are composed of 73 type-1 and 27 type-2 objects. The AKARI all-sky survey point source catalog was searched for their mid- and far-infrared counterparts at 9, 18, and 90 $μ$m. As a result, 69 Seyfert galaxies in the MAXI catalog (48 type-1 and 21 type-2 ones) were found to be detected with AKARI. The X-ray (3-4 keV and 4-10 keV) and infrared luminosities of these objects were investigated, together with their color information. Adopting the canonical photon index, $Γ= 1.9$, of the intrinsic X-ray spectrum of the Seyfert galaxies, the X-ray hardness ratio between the 3-4 and 4-10 keV ranges derived with MAXI was roughly converted into the absorption column density. After the X-ray luminosity was corrected for absorption from the estimated column density, the well-known X-ray-to-infrared luminosity correlation was confirmed at least in the Compton-thin regime. In contrast, NGC 1365, only one Compton-thick object in the MAXI catalog, was found to deviate from the correlation toward a significantly lower X-ray luminosity by nearly an order of magnitude. It was verified that the relation between the X-ray hardness below 10 keV and X-ray-to-infrared color acts as an effective tool to pick up Compton-thick objects. The difference in the infrared colors between the type-1 and type-2 Seyfert galaxies and its physical implication on the classification and unification of active galactic nuclei were briefly discussed.

preprint2015arXiv

A New Sample of Obscured AGNs Selected from the XMM-Newton and AKARI Surveys

We report a new sample of obscured active galactic nuclei (AGNs) selected from the XMM serendipitous source and AKARI point-source catalogs. We match X-ray sources with infrared (18 and 90 micron) sources located at |b|>10 deg to create a sample consisting of 173 objects. Their optical classifications and absorption column densities measured by X-ray spectra are compiled and study efficient selection criteria to find obscured AGNs. We apply the criteria (1) X-ray hardness ratio defined by using the 2-4.5 keV and 4.5-12 keV bands >-0.1 and (2) EPIC-PN count rate (CR) in the 0.2-12 keV to infrared flux ratio CR/F90<0.1 or CR/F18<1, where F18 and F90 are infrared fluxes at 18 and 90 micron in Jy, respectively, to search for obscured AGNs. X-ray spectra of 48 candidates, for which no X-ray results have been published, are analyzed and X-ray evidence for the presence of obscured AGNs such as a convex shape X-ray spectrum indicative of absorption of NH~10^{22-24} cm^-2, a very flat continuum, or a strong Fe-K emission line with an equivalent width of >700 eV is found in 26 objects. Six among them are classified as Compton-thick AGNs, and four are represented by either Compton-thin or Compton-thick spectral models. The success rate of finding obscured AGNs combining our analysis and the literature is 92% if the 18 micron condition is used. Of the 26 objects, 4 are optically classified as an HII nucleus and are new "elusive AGNs" in which star formation activity likely overwhelms AGN emission in the optical and infrared bands.

preprint2015arXiv

Evolution of mid-infrared galaxy luminosity functions from the entire AKARI NEP-Deep field with new CFHT photometry

We present infrared galaxy luminosity functions (LFs) in the AKARI North Ecliptic Pole (NEP) deep field using recently-obtained, wider CFHT optical/near-IR images. AKARI has obtained deep images in the mid-infrared (IR), covering 0.6 deg$^2$ of the NEP deep field. However, our previous work was limited to the central area of 0.25 deg$^2$ due to the lack of optical coverage of the full AKARI NEP survey. To rectify the situation, we recently obtained CFHT optical and near-IR images over the entire AKARI NEP deep field. These new CFHT images are used to derive accurate photometric redshifts, allowing us to fully exploit the whole AKARI NEP deep field. AKARI's deep, continuous filter coverage in the mid-IR wavelengths (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, and 24$μ$m) exists nowhere else, due to filter gaps of other space telescopes. It allows us to estimate restframe 8$μ$m and 12$μ$m luminosities without using a large extrapolation based on spectral energy distribution (SED) fitting, which was the largest uncertainty in previous studies. Total infrared luminosity (TIR) is also obtained more reliably due to the superior filter coverage. The resulting restframe 8$μ$m, 12$μ$m, and TIR LFs at $0.15<z<2.2$ are consistent with previous works, but with reduced uncertainties, especially at the high luminosity-end, thanks to the wide field coverage. In terms of cosmic infrared luminosity density ($Ω_{\mathrm{IR}}$), we found that the $Ω_{\mathrm{IR}}$ evolves as $\propto (1+z)^{4.2\pm 0.4}$.

preprint2015arXiv

Hyper-luminous Dust Obscured Galaxies discovered by the Hyper Suprime-Cam on Subaru and WISE

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

preprint2014arXiv

Luminosity and Redshift Dependence of the Covering Factor of AGNs viewed with WISE and SDSS

In this work, we investigate the dependence of the covering factor (CF) of active galactic nuclei (AGNs) (i) on the mid-infrared (MIR) luminosity and (ii) on the redshift. We constructed 12- and 22-micron luminosity functions (LFs) at 0.006 < z < 0.3 using the Wide-field Infrared Survey Explorer} (WISE) data. Combining the WISE catalog with the Sloan Digital Sky Survey (SDSS) spectroscopic data, we selected 223,982 galaxies at 12 micron and 25,721 galaxies at 22 micron for spectroscopic classification. We then identified 16,355 AGNs at 12 micron and 4,683 AGNs at 22 micron by their optical emission lines and cataloged classifications in the SDSS. Following that, we estimated the CF as the fraction of type 2 AGN in all AGNs whose MIR emissions are dominated by the active nucleus (not their host galaxies) based on their MIR colors. We found that (i) the CF decreased with increasing MIR luminosity, regardless of the choice of type 2 AGN classification criteria, and (ii) the CF did not change significantly with the redshift for z < 0.2. Furthermore, we carried out various tests to determine the influence of selection bias and confirmed similar dependences exist even when taking these uncertainties into account. The luminosity dependence of the CF can be explained by the receding torus model, but the "modified" receding torus model gives a slightly better fit, as suggested by Simpson.

preprint2013arXiv

The 9 and 18 Micron Luminosity Function of Various Types of Galaxies with AKARI: Implication for the Dust Torus Structure of AGN

We present the 9 and 18 micron luminosity functions (LFs) of galaxies at 0.006 < z < 0.8 (with an average redshift of ~ 0.04) using the AKARI mid-infrared all-sky survey catalog. We selected 243 galaxies at 9 micron and 255 galaxies at 18 micron from the Sloan Digital Sky Survey (SDSS) spectroscopy region. These galaxies were then classified by their optical emission lines, such as the line width of H_alpha or by their emission line ratios of [OIII]/H_beta and [NII]/H_alpha into five types: Type 1 active galactic nuclei (AGN) (Type 1); Type 2 AGN (Type 2); low-ionization narrow emission line galaxies (LINER); galaxies with both star formation and narrow-line AGN activity (composite galaxies); and star-forming galaxies (SF). We found that (i) the number density ratio of Type 2 to Type 1 AGNs is 1.73 +/- 0.36, which is larger than a result obtained from the optical LF and (ii) this ratio decreases with increasing 18 micron luminosity.

preprint2011arXiv

Outflow in Overlooked Luminous Quasar: Subaru Observations of AKARI J1757+5907

We present Subaru observations of the newly discovered luminous quasar AKARI J1757+5907, which shows an absorption outflow in its spectrum. The absorption consists of 9 distinct troughs, and our analysis focuses on the troughs at ~ -1000$ km s^{-1} for which we can measure accurate column densities of He I*, Fe II and Mg II. We use photoionization models to constrain the ionization parameter, total hydrogen column density, and the number density of the outflowing gas. These constraints yield lower limits for the distance, mass flow rate and kinetic luminosity for the outflow of 3.7 kpc, 70 M_{sun} yr^{-1}, and 2.0 x 10^{43} ergs s^{-1}, respectively. Such mass flow rate value can contribute significantly to the metal enrichment of the intra-cluster medium. We find that this moderate velocity outflow is similar to those recently discovered in massive post-starburst galaxies. Finally, we describe the scientific potential of future observations targeting this object.