Researcher profile

Yoshiki Toba

Yoshiki Toba contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
33works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

33 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&#39;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 &#34;outliers&#34; 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&#34; to <~0.&#34;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&#39;&#39;.8 \times 0&#39;&#39;.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&#34; 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 &#34;Non-merging&#39;&#39; 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 &#34;non-merging&#34; 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 &#34;non-merger&#34; 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&#39;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 &#34;AGN feedback&#34;, 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.