Source author record

Toshihiro Kawaguchi

Toshihiro Kawaguchi 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

33works
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

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.

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

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.

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.

preprint2020arXiv

Application of X-Ray Clumpy Torus Model (XCLUMPY) to 10 Obscured Active Galactic Nuclei Observed with Suzaku and NuSTAR

We apply XCLUMPY, an X-ray spectral model from a clumpy torus in an active galactic nucleus (AGN), to the broadband X-ray spectra of 10 obscured AGNs observed with both Suzaku and NuSTAR. The infrared spectra of these AGNs were analyzed with the CLUMPY code. Since XCLUMPY adopts the same clump distribution as that in the CLUMPY, we can directly compare the torus parameters obtained from the X-ray spectra and those from the infrared ones. The torus angular widths determined from the infrared spectra ($σ_{\mathrm{IR}}$) are systematically larger than those from the X-ray data ($σ_{\mathrm{X}}$); the difference ($σ_{\mathrm{IR}}-σ_{\mathrm{X}}$) correlates with the inclination angle determined from the X-ray spectrum. These results can be explained by the contribution from dusty polar outflows to the observed infrared flux, which becomes more significant at higher inclinations (more edge-on views). The ratio of the hydrogen column density and V-band extinction in the line of sight absorber shows large scatter ($\simeq$1 dex) around the Galactic value, suggesting that a significant fraction of AGNs have dust-rich circumnuclear environments.

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

Revisiting Soltan's argument based on a semi-analytical model for galaxy and black hole evolution

We show the significance of the super-Eddington accretion for the cosmic growth of supermassive black holes (SMBHs) with a semi-analytical model for galaxy and black hole evolution. The model explains various observed properties of galaxies and active galactic nuclei at a wide redshift range. By tracing the growth history of individual SMBHs, we find that the fraction of the SMBH mass acquired during the super-Eddington accretion phases to the total SMBH mass becomes larger for less massive black holes and at higher redshift. Even at z = 0, SMBHs with > 1e+9 Msun have acquired more than 50% of their mass by super-Eddington accretions, which is apparently inconsistent with classical Soltan's argument. However, the mass-weighted radiation efficiency of SMBHs with > 1e+8 Msun obtained with our model, is about 0.08 at z = 0, which is consistent with Soltan's argument within the observational uncertainties. We, therefore, conclude that Soltan's argument cannot reject the possibility that SMBHs are grown mainly by super-Eddington accretions.

preprint2020arXiv

Semi-analytic modelling of AGNs: auto-correlation function and halo occupation

The spatial clustering of active galactic nuclei (AGNs) is considered to be one of the important diagnostics for the understanding of the underlying processes behind their activities complementary to measurements of the luminosity function (LF). We analyse the AGN clustering from a recent semi-analytic model performed on a large cosmological $N$-body simulation covering a cubic gigaparsec comoving volume. We have introduced a new time-scale of gas accretion on to the supermassive black holes to account for the loss of the angular momentum on small scales, which is required to match the faint end of the observed X-ray LF. The large simulation box allows us accurate determination of the auto-correlation function of the AGNs. The model prediction indicates that this time-scale plays a significant role in allowing massive haloes to host relatively faint population of AGNs, leading to a higher bias factor for those AGNs. The model predictions are in agreement with observations of X-ray selected AGNs in the luminosity range $10^{41.5}~\mathrm{erg} \ \mathrm{s}^{-1} \leq L_{2-10\mathrm{keV}} \leq 10^{44.5}~\mathrm{erg} \ \mathrm{s}^{-1}$, with the typical host halo mass of $10^{12.5-13.5} h^{-1}\,{\rm M}_{\odot}$ at $z \lesssim 1$. This result shows that the observational clustering measurements impose an independent constraint on the accretion time-scale complementary to the LF measurements. Moreover, we find that not only the effective halo mass corresponding to the overall bias factor, but the extended shape of the predicted AGN correlation function shows remarkable agreement with those from observations. Further observational efforts towards the low luminosity end at $z \sim 1$ would give us stronger constraints on the triggering mechanisms of AGN activities through their clustering.

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

X-Ray Spectral Model from Clumpy Torus and Its Application to Circinus Galaxy

We construct an X-ray spectral model from the clumpy torus in an active galactic nucleus (AGN), designated as "XCLUMPY", utilizing the Monte Carlo simulation for Astrophysics and Cosmology framework (MONACO: Odaka et al. 2011, 2016). The adopted geometry of the torus is the same as that in Nenkova et al. (2008), who assume a power law distribution of clumps in the radial direction and a normal distribution in the elevation direction. We investigate the dependence of the X-ray continuum and Fe K$α$ fluorescence line profile on the torus parameters. Our model is compared with other torus models: MYTorus model (Murphy & Yaqoob 2009), Ikeda model (Ikeda et al. 2009), and CTorus model (Liu & Li 2014). As an example, we also present the results applied to the broadband X-ray spectra of the Circinus galaxy observed with XMM-Newton, Suzaku, and NuSTAR. Our model can well reproduce the data, yielding a hydrogen column density along the equatorial plane $N_{\mathrm{H}}^{\mathrm{Equ}} = 9.08_{-0.08}^{+0.14} \times 10^{24}$ cm$^{-2}$, a torus angular width $σ= 14.7_{-0.39}^{+0.44}$ degree, and a 2--10 keV luminosity $\log L_{2-10}/\mathrm{erg \ s^{-1}} = 42.8$. These results are discussed in comparison with the observations in other wavelengths.

preprint2016arXiv

Formation of the Andromeda Giant Stream: Asymmetric Structure and Disc Progenitor

We focus on the evidence of a past minor merger discovered in the halo of the Andromeda galaxy (M31). Previous N-body studies have enjoyed moderate success in producing the observed giant stellar stream (GSS) and stellar shells in M31's halo. The observed distribution of stars in the halo of M31 shows an asymmetric surface brightness profile across the GSS; however, the effect of the morphology of the progenitor galaxy on the internal structure of the GSS requires further investigation in theoretical studies. To investigate the physical connection between the characteristic surface brightness in the GSS and the morphology of the progenitor dwarf galaxy, we systematically vary the thickness, rotation velocity and initial inclination of the disc dwarf galaxy in N-body simulations. The formation of the observed structures appears to be dominated by the progenitor's rotation. Besides reproducing the observed GSS and two shells in detail, we predict additional structures for further observations. We predict the detectability of the progenitor's stellar core in the phase-space density distribution, azimuthal metallicity gradient of the western shell-like structure and an additional extended shell in the north-western direction that may constrain the properties of the progenitor galaxy.

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

Theoretical reevaluations of the black hole mass - bulge mass relation - I. Effect of the seed black hole mass

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

preprint2016arXiv

X-ray Spectral Model of Reprocessing by Smooth and Clumpy Molecular Tori in Active Galactic Nuclei with the MONACO framework

We construct an X-ray spectral model of reprocessing by a torus in an active galactic nucleus (AGN) with a Monte Carlo simulation framework MONACO. Two torus geometries of smooth and clumpy cases are considered and compared. In order to reproduce a Compton shoulder accurately, MONACO includes not only free electron scattering but also bound electron scattering. Raman and Reyleigh scattering are also treated, and scattering cross sections dependent on chemical states of hydrogen and helium are included. Doppler broadening by turbulence velocity can be implemented. Our model gives consistent results with other available models, such as MYTorus, except for differences due to different physical parameters and assumptions. We studied the dependence on torus parameters for Compton shoulder, and found that a intensity ratio of Compton shoulder to line core mainly depends on the column density, inclination angle, and metal abundance. For instance, an increase of metal abundance makes the Compton shoulder relatively weak. Also, shape of Compton shoulder depends on the column density. Furthermore, these dependences become different between smooth and clumpy cases. Then, we discuss the possibility of ASTRO-H SXS spectroscopy of Compton shoulder in AGN reflection spectra.

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.

preprint2015arXiv

Near-Infrared Spectroscopy of Quasars at z~3 and Estimates of Their Supermassive Black Hole Masses

We present the results of new infrared spectroscopic observations of 37 quasars at z~3, selected based on the optical r'-band magnitude and the availability of nearby bright stars for future imaging follow-up with Adaptive Optics system. The supermassive black hole (SMBH) masses (M_BH) were successfully estimated in 28 out of 37 observed objects from the combination of the H_beta emission linewidth and continuum luminosity at rest-frame 5100A. Comparing these results with those from previous studies of quasars with similar redshift, our sample exhibited slightly lower (~ -0.11 dex in median) Eddington ratios; and, the SMBH masses are slightly (~ 0.38 dex in median) higher. The SMBH growth time, t_grow, was calculated by dividing the estimated SMBH mass by the mass accretion rate measured using optical luminosity. We found, given reasonable assumptions, that t_grow was smaller than the age of the universe at the redshift of individual quasars for a large fraction of observed sources, suggesting that the SMBHs in many of our observed quasars are in growing phase with high accretion rates. A comparison of the SMBH masses derived from our H_beta data and archived CIV data indicated considerable scattering, as indicated in previous studies. All quasars with measured SMBH masses have at least one nearby bright star, such that they are suitable targets for adaptive optics observations to study the mass relationship between SMBHs and host galaxies' stellar component at high redshift.

preprint2014arXiv

Hunting A Wandering Supermassive Black Hole in M31 Halo -- Hermitage of Black Hole

In the hierarchical structure formation scenario, galaxies enlarge through multiple merging events with less massive galaxies. In addition, the Magorrian relation indicates that almost all galaxies are occupied by a central supermassive black hole (SMBH) of mass $10^{-3}$ of its spheroidal component. Consequently, SMBHs are expected to wander in the halos of their host galaxies following a galaxy collision, although evidence of this activity is currently lacking. We investigate a current plausible location of an SMBH wandering in the halo of the Andromeda galaxy (M31). According to theoretical studies of $N$-body simulations, some of the many substructures in the M31 halo are remnants of a minor merger occurring about 1 Gyr ago. First, to evaluate the possible parameter space of the infalling orbit of the progenitor, we perform numerous parameter studies using a Graphics Processing Unit (GPU) cluster. To reduce uncertainties in the predicted position of the expected SMBH, we then calculate the time evolution of the SMBH in the progenitor dwarf galaxy from $N$-body simulations using the plausible parameter sets. Our results show that the SMBH lies within the halo ($\sim$20--50 kpc from the M31 center), closer to the Milky Way than the M31 disk. Furthermore, the predicted current positions of the SMBH were restricted to an observational field of $0\degr.6 \times 0\degr.7$ in the northeast region of the M31 halo. We also discuss the origin of the infalling orbit of the satellite galaxy and its relationships with the recently discovered vast thin disk plane of satellite galaxies around M31.

preprint2014arXiv

Relics of Galaxy Merging: Observational Predictions for a Wandering Massive Black Hole and Accompanying Star Cluster in the Halo of M31

Galaxies and massive black holes (BHs) presumably grow via galactic merging events and subsequent BH coalescence. As a case study, we investigate the merging event between the Andromeda galaxy (M31) and a satellite galaxy. We compute the expected observational appearance of the massive BH that was at the center of the satellite galaxy prior to the merger, and is currently wandering in the M31 halo. We demonstrate that a radiatively inefficient accretion flow with a bolometric luminosity of a few tens of solar luminosities develops when Hoyle-Lyttleton accretion onto the BH is assumed. We compute the associated broadband spectrum and show that the radio band (observable with EVLA, ALMA and SKA) is the best frequency range to detect the emission. We also evaluate the mass and the luminosity of the stars bound by the wandering BH and find that such a star cluster is sufficiently luminous that it could correspond to one of the star clusters found by the PAndAS survey. The discovery of a relic massive BH wandering in a galactic halo will provide a direct means to investigate in detail the coevolution of galaxies and BHs. It also means a new population of BHs (off-center massive BHs), and offers targets for clean BH imaging that avoids strong interstellar scattering in the center of galaxies.

preprint2014arXiv

Statistical Properties of Multi-epoch Spectral Variability of SDSS Stripe82 Quasars

We investigate the UV-optical (longward of Ly$α$ 1216Å) spectral variability of nearly 9000 quasars ($0<z<4$) using multi-epoch photometric data within the SDSS Stripe 82 region. The regression slope in the flux-flux space of a quasar light curve directly measures the color of the flux difference spectrum, then the spectral shape of the flux difference spectra can be derived by taking a careful look at the redshift dependence of the regression slopes. First, we confirm that the observed quasar spectrum becomes bluer when the quasar becomes brighter. We infer the spectral index of the composite difference spectrum as $α_ν^{\text{dif}}\sim +1/3$ (in the form of $f_ν\propto ν^{α_ν}$), which is significantly bluer than that of the composite spectrum $α_ν^{\text{com}}\sim -0.5$. We also show that the continuum variability cannot be explained by the accretion disk models with varying mass accretion rate. Second, we examine the effects of broad emission line variability on the color-redshift space. The variability of the "Small Blue Bump" is extensively discussed. We show that the low-ionization lines of MgII and FeII are less variable compared to Balmer emission lines and high-ionization lines, and the Balmer continuum is the dominant variable source around $\sim 3000$Å. These results are compared with previous studies, and the physical mechanisms of the variability of the continuum and emission lines are discussed.

preprint2013arXiv

Innermost structure and near-infrared radiation of dusty clumpy tori in active galactic nuclei

The dusty clumpy torus surrounds the central black hole (BH) and the accretion disk in active galactic nuclei, and governs the growth of super-massive BHs via gas fueling towards the central engine. Near-infrared (NIR) monitoring observations have revealed that the torus inner radius is determined by the dust sublimation process. However, the observed radii are systematically smaller than the theoretical predictions by a factor of three. We take into account the anisotropic illumination by the central accretion disk to the torus, and calculate the innermost structure of the torus and the NIR time variablity. We then show that the anisotropy naturally solves the systematic descrepancy and that the viewing angle is the primary source to produce an object-to-object scatter of the NIR time delay. Dynamics of clumps at the innermost region of the torus will be unveiled via future high-resolution X-ray spectroscopy (e.g., Astro-H).

preprint2011arXiv

Combined Spectral and Timing Analysis of the Black Hole Candidate MAXI J1659-152 Discovered by MAXI and Swift

We report on X-ray spectral and timing results of the new black hole candidate (BHC) MAXI J1659-152 with the orbital period of 2.41 hours (shortest among BHCs) in the 2010 outburst from 65 Rossi X-ray Timing Explorer (RXTE) observations and 8 simultaneous Swift and RXTE observations. According to the definitions of the spectral states in Remillard & McClintock (2006), most of the observations have been classified into the intermediate state. All the X-ray broadband spectra can be modeled by a multi-color disk plus a power-law with an exponential cutoff or a multi-color disk plus a Comptonization component. During the initial phase of the outburst, a high energy cutoff was visible at 30-40 keV. The innermost radius of the disk gradually decreased by a factor of more than 3 from the onset of the outburst and reached a constant value of 35 d_10 cos i^-1/2 km, where d_10 is the distance in units of 10 kpc and $i$ is the inclination. The type-C quasi-periodic oscillation (QPO) frequency varied from 1.6 Hz to 7.3 Hz in association with a change of the innermost radius, while the innermost radius remained constant during the type-B QPO detections at 1.6-4.1 Hz. Hence, we suggest that the origin of the type-B QPOs is different from that of type-C QPOs, the latter of which would originate from the disk truncation radius. Assuming the constant innermost radius in the latter phase of the outburst as the innermost stable circular orbit, the black hole mass in MAXI J1659-152 is estimated to be 3.6-8.0 M_solar for a distance of 5.3-8.6 kpc and an inclination angle of 60-75 degrees.

preprint2011arXiv

Near-Infrared Reverberation by Dusty Clumpy Tori in Active Galactic Nuclei

According to recent models, the accretion disk and black hole in active galactic nuclei are surrounded by a clumpy torus. We investigate the NIR flux variation of the torus in response to a UV flash for various geometries. Anisotropic illumination by the disk and the torus self-occultation contrast our study with earlier works. Both the waning effect of each clump and the torus self-occultation selectively reduce the emission from the region with a short delay. Therefore, the NIR delay depends on the viewing angle (where a more inclined angle leads to a longer delay) and the time response shows an asymmetric profile with a negative skewness, opposing to the results for optically thin tori. The range of the computed delay coincides with the observed one, suggesting that the viewing angle is primarily responsible for the scatter of the observed delay. We also propose that the red NIR-to-optical color of type-1.8/1.9 objects is caused by not only the dust extinction but also the intrinsically red color. Compared with the modest torus thickness, both a thick and a thin tori display the weaker NIR emission. A selection bias is thus expected such that NIR-selected AGNs tend to possess moderately thick tori. A thicker torus shows a narrower and more heavily skewed time profile, while a thin torus produces a rapid response. A super-Eddington accretion rate leads to a much weaker NIR emission due to the disk self-occultation and the disk truncation by the self-gravity. A long delay is expected from an optically thin and/or a largely misaligned torus. A very weak NIR emission, such as in hot-dust-poor active nuclei, can arise from a geometrically thin torus, a super-Eddington accretion rate or a slightly misaligned torus.

preprint2010arXiv

Long-Term Spectral Variations of Ultraluminous X-ray Sources in the Interacting Galaxy Systems M51 and NGC4490/85

Variable ultraluminous X-ray sources (ULXs), which are considered to be black hole binaries (BHBs), are known to show state transitions similarly to Galactic BHBs. However, the relation between the ULX states and the Galactic BHB states is still unclear primarily due to less well-understood behaviors of ULXs in contrast to the Galactic BHBs. Here, we report a statistical X-ray spectral study of 34 energy spectra from seven bright ULXs in the interacting galaxy systems M51 and NGC4490/85, using archive data from multiple Chandra and XMM-Newton observations spanning for a few years. In order to compare with Galactic BHB states, we applied representative spectral models of BHBs; a power-law (PL), a multi-color disk black body (MCD), and a slim disk model to all the ULX spectra. We found a hint of a bimodal structure in the luminosity distribution of the samples, suggesting that ULXs have two states with typical luminosities of 3-6*10^{39} and 1.5-3*10^{39} ergs/s. Most spectra in the brighter state are explained by the MCD or the slim disk model, whereas those in the fainter state are explained by the PL model. In particular, the slim disk model successfully explains the observed spectral variations of NGC4490/85 ULX-6 and ULX-8 by changes of the mass accretion rate to a black hole of an estimated mass of <40 Msun. From the best-fit model parameters of each state, we speculate that the brighter state in these two ULXs corresponds to the brightest state of Galactic BHBs, which is often called the ``apparently standard state''. The fainter state of the ULXs has a PL shaped spectrum, but the photon index range is much wider than that seen in any single state of Galactic BHBs. We thus speculate that it is a state unique to ULXs. Some sources show much fainter and steeper spectra than the faint state, which we identified as another state.

preprint2010arXiv

Orientation Effects on the Inner Region of Dusty Torus of Active Galactic Nuclei

A sublimation process governs the innermost region of the dusty torus of active galactic nuclei. However, the observed inner radius of the torus is systematically smaller than the expected radius by a factor of ~ 1/3. We show that the anisotropy of the emission from accretion disks resolves this conflict naturally and quantitatively. An accretion disk emits lesser radiation in the direction closer to its equatorial plane (i.e., to the torus). We find that the anisotropy makes the torus inner region closer to the central black hole and concave. Moreover, the innermost edge of the torus may connect with the outermost edge of the disk continuously. Considering the anisotropic emission of each clump in the torus, we calculate the near-infrared flux variation in response to a UV flash. For an observer at the polar angle theta_obs = 25 deg, the centroid of the time delay is found to be 37% of the delay expected in the case of isotropic illumination, which explains the observed systematic deviation.

preprint2006arXiv

A Bar Fuels a Super-Massive Black Hole?: Host Galaxies of Narrow-Line Seyfert 1 Galaxies

We present optical images of nearby 50 narrow-line Seyfert 1 galaxies (NLS1s) which cover all the NLS1s at z<0.0666 and $δ\ge -25^{\circ}$ known at the time of 2001. Among the 50 NLS1s, 40 images are newly obtained by our observations and 10 images are taken from archive data. Motivated by the hypothesis that NLS1s are in an early phase of a super-massive black hole (BH) evolution, we present a study of NLS1 host galaxy morphology to examine trigger mechanism(s) of active galactic nuclei (AGNs) by seeing the early phase of AGN. With these images, we made morphological classification by eye inspection and by quantitative method, and found a high bar frequency of the NLS1s in the optical band; the bar frequency is $85 \pm 7% among disk galaxies ($64- 71$% in total sample) which is more frequent than that (40-70%) of broad-line Seyfert 1 galaxies (BLS1s) and normal disk galaxies, though the significance is marginal. Our results confirm the claim by Crenshaw et al. (2003) with a similar analysis for 19 NLS1s. The frequency is comparable to that of HII/starburst galaxies. We also examined the bar frequency against width of the broad H$β$ emission line, Eddington ratio, and black hole mass, but no clear trend is seen. Possible implications such as an evolutionary sequence from NLS1s to BLS1 are discussed briefly.

preprint2006arXiv

Do Ultraluminous X-ray Sources Really Contain Intermediate-mass Black Holes?

An open question remains whether Ultraluminous X-ray Sources (ULXs) really contain intermediate-mass black holes (IMBHs). We carefully investigated the XMM-Newton EPIC spectra of the four ULXs that were claimed to be strong candidates of IMBHs by several authors. We first tried fitting by the standard spectral model of disk blackbody (DBB) + power-law (PL), finding good fits to all of the data, in agreement with others. We, however, found that the PL component dominates the DBB component at $\sim$ 0.3 to 10 keV. Thus, the black hole parameters derived solely from the minor DBB component are questionable. Next, we tried to fit the same data by the ``$p$-free disk model'' without the PL component, assuming the effective temperature profile of $T_{\rm eff} \propto r^{-p}$ where $r$ is the disk radius. Interestingly, in spite of one less free model parameters, we obtained similarly good fits with much higher innermost disk temperatures, $1.8 < kT_{\rm in} < 3.2$ keV. More importantly, we obtained $p \sim 0.5$, just the value predicted by the slim (super-critical) disk theory, rather than $p = 0.75$ that is expected from the standard disk model. The estimated black hole masses from the $p$-free disk model are much smaller; $M\ltsim 40 M_\odot$. Furthermore, we applied a more sophisticated slim disk model by Kawaguchi (2003, ApJ, 593,69), and obtained good fits with roughly consistent black hole masses. We thus conclude that the central engines of these ULXs are super-critical accretion flows to stellar-mass black holes.

preprint2006arXiv

Systematic Effects in Measurement of Black Hole Masses by Emission-Line Reverberation of Active Galactic Nuclei: Eddington Ratio and Inclination

We try to identify the sources and systematics of the uncertainty of the masses of the central black holes in active galactic nuclei (AGNs) which have been obtained by reverberation-mapping methods. We characterize the broad H$β$ emission-line profiles by the ratio of their full-width at half maximum (FWHM) to their line dispersion, i.e., the second moment of the line profile. We use this parameter to separate the reverberation-mapped AGNs into two populations, the first with narrower H$β$ lines with relatively extended wings, and the second with broader lines that are relatively flat-topped. The first population is characterized by higher Eddington ratios than the second. Within each population, we calibrate the black-hole mass scale by comparison of the reverberation-based mass with that predicted by the bulge velocity dispersion. We also use the distribution of ratios of the reverberation-based mass to the velocity-dispersion mass prediction in a comparison with a ``generalized thick disk'' model in order to see if inclination can plausibly account for the observed distribution. We find that the line dispersion is a less biased parameter in general than FWHM for black hole mass estimation, although we show that it is possible to empirically correct for the bias introduced by using FWHM to characterize the emission-line width. We also argue that inclination effects are apparent in a small subset of objects with the narrowest emission lines. Our principal conclusion is that the H$β$ profile is sensitive primarily to Eddington ratio, but that inclination effects play a role in some cases. (abridged)

preprint2000arXiv

Slim Disk Model for Soft X-Ray Excess and Variability of Narrow-Line Seyfert 1 Galaxies

Narrow-line Seyfert 1 galaxies (NLS1s) exhibit extreme soft X-ray excess and large variability. We argue that both features can be basically accounted for by the slim disk model. We assume that a central black-hole mass in NLS1 is relatively small, $M \sim 10^{5-7}M_\odot$, and that a disk shines nearly at the Eddington luminosity, $L_{\rm E}$. Then, the disk becomes a slim disk and exhibits the following distinctive signatures: (1) The disk luminosity (particularly of X-rays) is insensitive to mass-flow rates, $\dot M$, since the generated energy is partly carried away to the black hole by trapped photons in accretion flow. (2) The spectra are multi-color blackbody. The maximum blackbody temperature is $T_{\rm bb} \simeq 0.2(M/10^5 M_\odot)^{-1/4}$ keV, and the size of the blackbody emitting region is small, $r_{\rm bb} \lsim 3 r_{\rm S}$ (with $r_{\rm S}$ being Schwarzschild radius) even for a Schwarzschild black hole. (3) All the ASCA observation data of NLS1s fall onto the region of $\dot M/(L_{\rm E}/c^2)>10$ (with $L_{\rm E}$ being the Eddington luminosity) on the ($r_{\rm bb},T_{\rm bb}$) plane, supporting our view that a slim disk emits soft X-rays at $\sim L_{\rm E}$ in NLS1s. (4) Magnetic energy can be amplified, at most, up to the equipartition value with the trapped radiation energy which greatly exceeds radiation energy emitted from the disk. Hence, energy release by consecutive magnetic reconnection will give rise to substantial variability in soft X-ray emission.