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.

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.

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.

We report a discovery of $z = 4.72$ radio galaxy, HSC J083913.17+011308.1, by using the Lyman break technique with the Hyper Suprime-Cam Subaru Strategic Survey (HSC-SSP) catalog for VLA FIRST radio sources. The number of known high-$z$ radio galaxies (HzRGs) at $z > 3$ is quite small to constrain the evolution of HzRGs so far. The deep and wide-area optical survey by HSC-SSP enables us to apply the Lyman break technique to a large search for HzRGs. For an HzRG candidate among pre-selected $r$-band dropouts with a radio detection, a follow-up optical spectroscopy with GMOS/Gemini has been performed. The obtained spectrum presents a clear Ly$\alpha$ emission line redshifted to $z=4.72$. The SED fitting analysis with the rest-frame UV and optical photometries suggests the massive nature of this HzRG with $\log{M_*/M_{\odot}} = 11.4$. The small equivalent width of Ly$\alpha$ and the moderately red UV colors indicate its dusty host galaxy, implying a chemically evolved and dusty system. The radio spectral index does not meet a criterion for an ultra-steep spectrum: $\alpha^{325}_{1400}$ of $-1.1$ and $\alpha^{150}_{1400}$ of $-0.9$, demonstrating that the HSC-SSP survey compensates for a sub-population of HzRGs which are missed in surveys focusing on an ultra-steep spectral index.

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.

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 ($\sigma_{\mathrm{IR}}$) are systematically larger than those from the X-ray data ($\sigma_{\mathrm{X}}$); the difference ($\sigma_{\mathrm{IR}}-\sigma_{\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.

We investigate the correlation between rest-frame UV flux variability of broad absorption line (BAL) quasars and their variability in BAL equivalent widths (EWs) in a various timescale from $<10$~days to a few years in the quasar rest-frame. We use the data sets of BAL EWs taken by the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project and photometric data taken by the intermediate Palomar Transient Factory (iPTF) in $g$ and $R$-bands and the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) in $grizy$ bands. Our results are summarized as below; (1) the distributions of flux variability versus BAL variability show weak, moderate, or a strong positive correlation, (2) there is no significant difference in flux variability amplitudes between BAL quasar with significant short timescale EW variability (called class S1) and without (class S2), (3) in all time scales considered in this paper, the class S1 quasars show systematically larger BAL variability amplitudes than those of the class S2 quasars, and (4) there are possible correlations between BAL variability and physical parameters of the quasars such as black hole masses (moderate positive), Eddington ratios, and accretion disk temperature (strong negative) in the class S2 quasars. These results indicate that the BAL variability requires changing in the ionizing continuum and an ancillary mechanism such as variability in X-ray shielding gas located at the innermost region of an accretion disk.

We studied the final phases of galactic mergers, focusing on interactions between supermassive black holes (SMBHs) and the interstellar medium in a central sub-kpc region, using an N-body/hydrodynamics code. This numerical experiment aims to understand the fate of the gas supplied by mergers of two or more galaxies with SMBHs, whose masses are $10^7 M_\odot$. We observed that the mass accretion rate to one SMBH exceeds the Eddington accretion rate when the distance between two black holes (BHs) rapidly decreases. However, this rapid accretion phase does not last for more than $10^7$ yrs, and it drops to $\sim$ 10% of the Eddington rate in the quasi-steady accretion phase. The rapid accretion is caused by the angular momentum transfer from the gas to the stellar component, and the moderate accretion in the quasi-steady phase is caused {by turbulent viscosity and gravitational torque in the disk. The second merger event enhances the mass accretion to the BHs; however, this phase takes place on a similar timescale to the first merger event. We also found that the AGN feedback and the mass accretion to BHs can coexist in the central region of merged galaxies, if the amount of feedback energy is given as $(2 \times 10^{-4} - 2 \times 10^{-3} )\dot{M} c^2$, where $\dot{M}$ is the accretion rate to $r= 1$ pc. The accretion rate is suppressed by $\sim$ 1/50 in the quasi-steady accretion phase for $0.02 \dot{M} c^2$. The fraction of the gas that finally falls to each BH is approximately 5-7% of the supplied total gas mass ($10^8 M_\odot$), and 15-20% of the gas forms a circumnuclear gas inside 100 pc. This remnant gas heavily obscures the luminous phase of the active galactic nuclei (AGN) during merger events, and the moderate AGN feedback does not alter this property.

We report the discovery of 28 quasars and 7 luminous galaxies at 5.7 $\le$ z $\le$ 7.0. This is the tenth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The total number of spectroscopically identified objects in SHELLQs has now grown to 93 high-z quasars, 31 high-z luminous galaxies, 16 [O III] emitters at z ~ 0.8, and 65 Galactic cool dwarfs (low-mass stars and brown dwarfs). These objects were found over 900 deg2, surveyed by HSC between 2014 March and 2018 January. The full quasar sample includes 18 objects with very strong and narrow Ly alpha emission, whose stacked spectrum is clearly different from that of other quasars or galaxies. While the stacked spectrum shows N V 1240 emission and resembles that of lower-z narrow-line quasars, the small Ly alpha width may suggest a significant contribution from the host galaxies. Thus these objects may be composites of quasars and star-forming galaxies.

We present deep near-infrared spectroscopy of six quasars at 6.1<z<6.7 with VLT/X-Shooter and Gemini-N/GNIRS. Our objects, originally discovered through a wide-field optical survey with the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP), have the lowest luminosities (-25.5< M1450<-23.1 mag) of the z>5.8 quasars with measured black hole masses. From single-epoch mass measurements based on MgII2798, we find a wide range in black hole masses, from M_BH=10^7.6 to 10^9.3 Msun. The Eddington ratios L_bol/L_Edd range from 0.16 to 1.1, but the majority of the HSC quasars are powered by M_BH=10^9 Msun supermassive black holes (SMBHs) accreting at sub-Eddington rates. The Eddington ratio distribution of the HSC quasars is inclined to lower accretion rates than those of Willott et al. (2010a), who measured the black hole masses for similarly faint z=6 quasars. This suggests that the global Eddington ratio distribution is wider than has previously been thought. The presence of M_BH=10^9 Msun SMBHs at z=6 cannot be explained with constant sub-Eddington accretion from stellar remnant seed black holes. Therefore, we may be witnessing the first buildup of the most massive black holes in the first billion years of the universe, the accretion activity of which is transforming from active growth to a quiescent phase. Measurements of a larger complete sample of z>6 low-luminosity quasars, as well as deeper observations with future facilities will enable us to better understand the early SMBH growth in the reionization epoch.

We have obtained three-dimensional maps of the universe in $\sim200\times200\times80$ comoving Mpc$^3$ (cMpc$^3$) volumes each at $z=5.7$ and $6.6$ based on a spectroscopic sample of 179 galaxies that achieves $\gtrsim80$\% completeness down to the Ly$\alpha$ luminosity of $\log(L_{\rm Ly\alpha}/[\mathrm{erg\ s^{-1}}])=43.0$, based on our Keck and Gemini observations and the literature. The maps reveal filamentary large-scale structures and two remarkable overdensities made out of at least 44 and 12 galaxies at $z=5.692$ (z57OD) and $z=6.585$ (z66OD), respectively, making z66OD the most distant overdensity spectroscopically confirmed to date with $>10$ spectroscopically confirmed galaxies. We compare spatial distributions of submillimeter galaxies at $z\simeq 4-6$ with our $z=5.7$ galaxies forming the large-scale structures, and detect a $99.97\%$ signal of cross correlation, indicative of a clear coincidence of dusty star-forming galaxy and dust unobscured galaxy formation at this early epoch. The galaxies in z57OD and z66OD are actively forming stars with star formation rates (SFRs) $\gtrsim5$ times higher than the main sequence, and particularly the SFR density in z57OD is 10 times higher than the cosmic average at the redshift (a.k.a. the Madau-Lilly plot). Comparisons with numerical simulations suggest that z57OD and z66OD are protoclusters that are progenitors of the present-day clusters with halo masses of $\sim10^{14}\ \mathrm{M_\odot}$.

We show the Eddington ratio distributions of supermassive black holes at a wide redshift range (0 < z < 8) obtained with a semi-analytic model of galaxy formation. The distribution is broadly consistent with observational estimates at low redshift. We find that the growth rate of black holes at higher redshift is more likely to exceed the Eddington limit because the typical gas fraction of the host galaxies is higher at higher redshift. We also find that the super- Eddington growth is more common for less massive supermassive black holes, supporting an idea that supermassive black holes have been formed via super-Eddington accretion. These results indicate the "slowing down" of cosmic growth of supermassive black holes: the growth of supermassive black holes with a higher Eddington ratio peaks at higher redshift. We also show the effect of the sample selection on the shape of the Eddington ratio distribution functions and find that shallower observations will miss active galactic nuclei with not only the smaller but also higher Eddington ratios.

We present physical properties of radio galaxies (RGs) with $f_{\rm 1.4 GHz} >$ 1 mJy discovered by Subaru Hyper Supreme-Cam (HSC) and VLA Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) survey. For 1056 FIRST RGs at $0 < z \leq 1.7$ with HSC counterparts in about 100 deg$^2$, we compiled multi-wavelength data of optical, near-infrared (IR), mid-IR, far-IR, and radio (150 MHz). We derived their color excess ($E (B-V)_{*}$), stellar mass, star formation rate (SFR), IR luminosity, the ratio of IR and radio luminosity ($q_{\rm IR}$), and radio spectral index ($\alpha_{\rm radio}$) that are derived from the SED fitting with CIGALE. We also estimated Eddington ratio based on stellar mass and integration of the best-fit SEDs of AGN component. We found that $E (B-V)_{*}$, SFR, and IR luminosity clearly depend on redshift while stellar mass, $q_{\rm IR}$, and $\alpha_{\rm radio}$ do not significantly depend on redshift. Since optically-faint ($i_{\rm AB} \geq 21.3$) RGs that are newly discovered by our RG survey tend to be high redshift, they tend to not only have a large dust extinction and low stellar mass but also have high SFR and AGN luminosity, high IR luminosity, and high Eddington ratio compared to optically-bright ones. The physical properties of a fraction of RGs in our sample seem to differ from a classical view of RGs with massive stellar mass, low SFR, and low Eddington ratio, demonstrating that our RG survey with HSC and FIRST provides us curious RGs among entire RG population.

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$\alpha$ 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 $\sigma = 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.

We report on the optical properties of infrared (IR)-bright dust-obscured galaxies (DOGs) that are defined as (i-[22])_AB > 7.0. Because supermassive black holes (SMBHs) in IR-bright DOGs are expected to be rapidly growing in the major merger scenario, they provide useful clues for understanding the co-evolution of SMBHs and their host galaxies. However, the optical properties of IR-bright DOGs remain unclear because the optical emission of a DOG is very faint. By combining ~ 105 deg^2 images of the optical, near-IR, and mid-IR data obtained from the Subaru Hyper Suprime-Cam (HSC) survey, the VISTA VIKING survey, and the WISE all-sky survey, respectively, 571 IR-bright DOGs were selected. We found that IR-bright DOGs show a redder (g-z)_AB color than other populations of dusty galaxies, such as ultra-luminous IR galaxies (ULIRGs) at a similar redshift, with a significantly large dispersion. Among the selected DOGs, star-formation (SF) dominated DOGs show a relatively red color, while active galactic nucleus (AGN) dominated DOGs show a rather blue color in optical. This result is consistent with the idea that the relative AGN contribution in the optical emission becomes more significant at a later stage in the major merger scenario. We discovered eight IR-bright DOGs showing a significant blue excess in blue HSC bands (BluDOGs). This blue excess can be interpreted as a leaked AGN emission that is either a directly leaking or a scattered AGN emission, as proposed for some blue-excess Hot DOGs in earlier studies.

We apply a new X-ray clumpy torus model called XCLUMPY (Tanimoto et al. 2019), where the clump distribution is assumed to be the same as in the infrared clumpy torus model (CLUMPY) by Nenkova et al. 2008, to the broadband X-ray spectra of type-1 active galactic nuclei (AGNs) for the first time. We analyze the archival data of IC 4329A and NGC 7469 observed with NuSTAR/Suzaku and NuSTAR/XMM-Newton, respectively, whose infrared spectra were studied with CLUMPY by Ichikawa et al. 2015 and optical extinctions ($A_{\rm V}$) of the tori were estimated. We consider two models, invoking (Model 1) a relativistic reflection component from the accretion disk and (Model 2) a partial absorber. Assuming that the narrow Fe K$\alpha$ emission line at 6.4 keV originates from the torus, we separate the contribution of the torus reflection components in the total spectra. Our models yield equatorial hydrogen column densities of the tori to be $N^{\rm Equ}_{\rm H} =$ (0.53--1.43) $\times 10^{23}~\rm cm^{-2}$ and $N^{\rm Equ}_{\rm H} =$ (0.84--1.43) $\times 10^{24}~\rm cm^{-2}$, for IC 4329A and NGC 7469, respectively. We find that the $N_{\rm H}/A_{\rm V}$ ratios in the tori are by factors of 25--68 (IC 4329A) and 2.4--3.9 (NGC 7469) smaller than that in the Galactic interstellar medium (ISM). These results suggest that a non-negligible fraction of AGNs are 'dust-rich' compared with the Galactic ISM, as opposite to the general trend previously reported in many obscured AGNs.

We report the discovery of a quasar at z = 7.07, which was selected from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. This quasar, HSC J124353.93+010038.5, has an order of magnitude lower luminosity than do the other known quasars at z > 7. The rest-frame ultraviolet absolute magnitude is M1450 = -24.13 +/- 0.08 mag and the bolometric luminosity is Lbol = (1.4 +/- 0.1) x 10^{46} erg/s. Its spectrum in the optical to near-infrared shows strong emission lines, and shows evidence for a fast gas outflow, as the C IV line is blueshifted and there is indication of broad absorption lines. The Mg II-based black hole mass is Mbh = (3.3 +/- 2.0) x 10^8 Msun, thus indicating a moderate mass accretion rate with an Eddington ratio 0.34 +/- 0.20. It is the first z > 7 quasar with sub-Eddington accretion, besides being the third most distant quasar, known to date. The luminosity and black hole mass are comparable to, or even lower than, those measured for the majority of low-z quasars discovered by the Sloan Digital Sky Survey, and thus this quasar likely represents a z > 7 counterpart to quasars commonly observed in the low-z universe.

In order to construct a sample of ultra-luminous infrared galaxies (ULIRGs, with infrared luminosity, $L_{\rm IR} > 10^{12}$ L$_{\odot}$) at 0.5 < z < 1, we are conducting an optical follow-up program for bright 90-$\mu$m FIR sources with a faint optical (i < 20) counterpart selected in the AKARI Far-Infrared Surveyor (FIS) Bright Source catalog (Ver.2). AKARI-FIS-V2 J0916248+073034, identified as a ULIRG at z = 0.49 in the spectroscopic follow-up observation, indicates signatures of an extremely strong outflow in its emission line profiles. Its [OIII] 5007\AA\ emission line shows FWHM of 1830 km s$^{-1}$ and velocity shift of -770 km s$^{-1}$ in relative to the stellar absorption lines. Furthermore, low-ionization [OII] 3726\AA\ 3729\AA\ doublet also shows large FWHM of 910 km s$^{-1}$ and velocity shift of -380 km s$^{-1}$. After the removal of an unresolved nuclear component, the long-slit spectroscopy 2D image possibly shows that the outflow extends to radius of 4 kpc. The mass outflow and energy ejection rates are estimated to be 500 M$_{\odot}$ yr$^{-1}$ and $4\times10^{44}$ erg s$^{-1}$, respectively, which imply that the outflow is among the most powerful ones observed in ULIRGs and QSOs at 0.3 < z < 1.6. The co-existence of the strong outflow and intense star formation (star formation rate of 990 M$_{\odot}$ yr$^{-1}$) indicates that the feedback of the strong outflow has not severely affect the star-forming region of the galaxy.

We present new measurements of the quasar luminosity function (LF) at $z \sim 6$, over an unprecedentedly wide range of the rest-frame ultraviolet luminosity $M_{1450}$ from $-30$ to $-22$ mag. This is the fifth in a series of publications from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The LF was calculated with a complete sample of 110 quasars at $5.7 \le z \le 6.5$, which includes 48 SHELLQs quasars discovered over 650 deg$^2$, and 63 brighter quasars discovered by the Sloan Digital Sky Survey and the Canada-France-Hawaii Quasar Survey (including one overlapping object). This is the largest sample of $z \sim 6$ quasars with a well-defined selection function constructed to date, and has allowed us to detect significant flattening of the LF at its faint end. A double power-law function fit to the sample yields a faint-end slope $\alpha = -1.23^{+0.44}_{-0.34}$, a bright-end slope $\beta = -2.73^{+0.23}_{-0.31}$, a break magnitude $M_{1450}^* = -24.90^{+0.75}_{-0.90}$, and a characteristic space density $\Phi^* = 10.9^{+10.0}_{-6.8}$ Gpc$^{-3}$ mag$^{-1}$. Integrating this best-fit model over the range $-18 < M_{1450} < -30$ mag, quasars emit ionizing photons at the rate of $\dot{n}_{\rm ion} = 10^{48.8 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$ at $z = 6.0$. This is less than 10 % of the critical rate necessary to keep the intergalactic medium ionized, which indicates that quasars are not a major contributor to cosmic reionization.

We present the latest results of a semi-analytic model of galaxy formation, "New Numerical Galaxy Catalogue", which is combined with large cosmological N-body simulations. This model can reproduce statistical properties of galaxies at z < 6.0. We focus on the properties of active galactic nuclei (AGNs) and supermassive black holes, especially on the accretion timescale onto black holes. We find that the number density of AGNs at z < 1.5 and at hard X-ray luminosity 10^{ 44 }< erg/s is underestimated compared with recent observational estimates when we assume the exponentially decreasing accretion rate and the accretion timescale which is proportional to the dynamical time of the host halo or the bulge, as is often assumed in semi-analytic models. We show that to solve this discrepancy, the accretion timescale of such less luminous AGNs instead should be a function of the black hole mass and the accreted gas mass. This timescale can be obtained from a phenomenological modelling of the gas angular momentum loss in the circumnuclear torus and/or the accretion disc. Such models predict a longer accretion timescale for less luminous AGNs at z < 1.0 than bright QSOs whose accretion timescale would be 10^{ 7-8 } yr. With this newly introduced accretion timescale, our model can explain the observed luminosity functions of AGNs at z < 6.0.

We report the result of optical identifications of FIRST radio sources with the Hyper Suprime-Cam Subaru Strategic Program survey (HSC-SSP). The positional cross-match within 1" between the FIRST and HSC-SSP catalogs (i ~< 26) produced more than 3600 optical counterparts in the 156 deg^2 of the HSC-SSP field. The matched counterparts account for more than 50% of the FIRST sources in the search field, which substantially exceed previously reported fractions of SDSS counterparts (i ~< 22) of ~30%. Among the matched sample, 9% are optically unresolved sources such as radio-loud quasars. The optically faint (i > 21) radio galaxies (RGs) show that the fitting linear function of the 1.4 GHz source counts has a slope that is flatter than that of the bright RGs, while optically faint radio quasars show a slope steeper than that of bright radio quasars. The optically faint RGs show a flat slope in the i-band number counts down to 24 mag, implying either less massive or distant radio-active galactic nuclei (AGNs) beyond 24 mag. The photometric redshift and the comparison of colors with the galaxy models show that most of the matched RGs are distributed at redshifts from 0 to 1.5. The optically faint sample includes the high radio-loudness sources that are not seen in the optically bright sample. Such sources are located at redshift z > 1. This study gives ~1500 radio AGNs lying at the optically faint end and high-redshift regime not probed by previous searches.

We analyzed the spectral shape of the Compton shoulder around the neutral Fe-K$_\alpha$ line of the Compton-thick type II Seyfert nucleus of the Circinus galaxy. The characteristics of this Compton shoulder with respect to the reflected continuum and Fe-K$_\alpha$ line core intensity are a powerful diagnostics tool for analyzing the structure of the molecular tori, which obscure the central engine. We applied our Monte-Carlo-based X-ray reflection spectral model to the Chandra High Energy Transmission Grating data and successfully constrained the various spectral parameters independently, using only the spectral data only around the Fe-K$_\alpha$ emission line. The obtained column density and inclination angle are consistent with the previous observations and the Compton-thick type II Seyfert picture. In addition, we determined the metal abundance of the molecular torus for the case of the smooth and clumpy torus to be 1.75$^{+0.19}_{-0.17}$ and 1.74$\pm$0.16 solar abundance, respectively. Such slightly over-solar abundance can be useful information for discussing the star formation rate in the molecular tori of active galactic nuclei.

We report the initial result of an adaptive-optics assisted, optical integral-field-unit observation on IRAS04576+0912, the nearest (z=0.039) active galactic nucleus with a prominent blueshift/tail in [O III] emission from a sample of such objects that we have collected from the literature. We aim at addressing the putative quasar-mode feedback process with Subaru/Kyoto 3D II+AO188. The optical waveband (6400--7500 AA) enables us to measure the gas density via the [S II] doublets, in contrast to earlier Near-IR studies. Since the fast [O III] outflow happens only around rapidly growing central black holes, this object is suitable for investigating the black hole-galaxy coevolution. The obtained data cube exhibits blue tail in the [S II] emission at many lenslets. By fitting the spectrum with the high excess flux at the [S II] blue tail, we find the fast (~ 860 km/s), dense (>3000/cc), wide-angle and offset outflow in central 100-pc scales. Although the large opening angle and the high gas outflow-to-accretion ratio may favour the feedback hypothesis, the inferred kinetic power injection rate of this ionized gas outflow seems insufficient to influence the whole host galaxy. A conventional assumption of a low density must have overestimated the feedback process.

We report discovery of 41 new high-z quasars and luminous galaxies, which were spectroscopically identified at 5.7 < z < 6.9. This is the fourth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, based on the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We selected the photometric candidates by a Bayesian probabilistic algorithm, and then carried out follow-up spectroscopy with the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous papers, we have now spectroscopically identified 137 extremely-red HSC sources over about 650 deg2, which include 64 high-z quasars, 24 high-z luminous galaxies, 6 [O III] emitters at z ~ 0.8, and 43 Galactic cool dwarfs (low-mass stars and brown dwarfs). The new quasars span the luminosity range from M1450 ~ -26 to -22 mag, and continue to populate a few magnitude lower luminosities than have been probed by previous wide-field surveys. In a companion paper, we derive the quasar luminosity function at z ~ 6 over an unprecedentedly wide range of M1450 ~ -28 to -21 mag, exploiting the SHELLQs and other survey outcomes.

Narrow-line regions excited by active galactic nuclei (AGN) are important for studying AGN photoionization and feedback. Their strong [O III] lines can be detected with broadband images, allowing morphological studies of these systems with large-area imaging surveys. We develop a new technique to reconstruct the [O III] images using the Subaru Hyper Suprime-Cam (HSC) Survey aided with spectra from the Sloan Digital Sky Survey (SDSS). The technique involves a careful subtraction of the galactic continuum to isolate emission from the [O III]$\lambda$5007 and [O III]$\lambda$4959 lines. Compared to traditional targeted observations, this technique is more efficient at covering larger samples with less dedicated observational resources. We apply this technique to an SDSS spectroscopically selected sample of 300 obscured AGN at redshifts 0.1 - 0.7, uncovering extended emission-line region candidates with sizes up to tens of kpc. With the largest sample of uniformly derived narrow-line region sizes, we revisit the narrow-line region size-luminosity relation. The area and radii of the [O III] emission-line regions are strongly correlated with the AGN luminosity inferred from the mid-infrared (15 $\mu$m rest-frame) with a power-law slope of $0.62^{+0.05}_{-0.06} \pm 0.10$ (statistical and systemic errors), consistent with previous spectroscopic findings. We discuss the implications for the physics of AGN emission-line region and future applications of this technique, which should be useful for current and next-generation imaging surveys to study AGN photoionization and feedback with large statistical samples.

We construct a sample of X-ray bright optically faint active galactic nuclei by combining Subaru Hyper Suprime-Cam, XMM-Newton, and infrared source catalogs. 53 X-ray sources satisfying i band magnitude fainter than 23.5 mag and X-ray counts with EPIC-PN detector larger than 70 are selected from 9.1 deg^2, and their spectral energy distributions (SEDs) and X-ray spectra are analyzed. 44 objects with an X-ray to i-band flux ratio F_X/F_i>10 are classified as extreme X-ray-to-optical flux sources. SEDs of 48 among 53 are represented by templates of type 2 AGNs or starforming galaxies and show signature of stellar emission from host galaxies in the optical in the source rest frame. Infrared/optical SEDs indicate significant contribution of emission from dust to infrared fluxes and that the central AGN is dust obscured. Photometric redshifts determined from the SEDs are in the range of 0.6-2.5. X-ray spectra are fitted by an absorbed power law model, and the intrinsic absorption column densities are modest (best-fit log N_H = 20.5-23.5 cm^-2 in most cases). The absorption corrected X-ray luminosities are in the range of 6x10^42 - 2x10^45 erg s^-1. 20 objects are classified as type 2 quasars based on X-ray luminsosity and N_H. The optical faintness is explained by a combination of redshifts (mostly z>1.0), strong dust extinction, and in part a large ratio of dust/gas.

We present spectroscopic identification of 32 new quasars and luminous galaxies discovered at 5.7 < z < 6.8. This is the second in a series of papers presenting the results of the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The photometric candidates were selected by a Bayesian probabilistic algorithm, and then observed with spectrographs on the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous paper, we have now identified 64 HSC sources over about 430 deg2, which include 33 high-z quasars, 14 high-z luminous galaxies, 2 [O III] emitters at z ~ 0.8, and 15 Galactic brown dwarfs. The new quasars have considerably lower luminosity (M1450 ~ -25 to -22 mag) than most of the previously known high-z quasars. Several of these quasars have luminous (> 10^(43) erg/s) and narrow (< 500 km/s) Ly alpha lines, and also a possible mini broad absorption line system of N V 1240 in the composite spectrum, which clearly separate them from typical quasars. On the other hand, the high-z galaxies have extremely high luminosity (M1450 ~ -24 to -22 mag) compared to other galaxies found at similar redshift. With the discovery of these new classes of objects, we are opening up new parameter spaces in the high-z Universe. Further survey observations and follow-up studies of the identified objects, including the construction of the quasar luminosity function at z ~ 6, are ongoing.

We examine the clustering of quasars over a wide luminosity range, by utilizing 901 quasars at $\overline{z}_{\rm phot}\sim3.8$ with $-24.73<M_{\rm 1450}<-22.23$ photometrically selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) S16A Wide2 date release and 342 more luminous quasars at $3.4<z_{\rm spec}<4.6$ having $-28.0<M_{\rm 1450}<-23.95$ from the Sloan Digital Sky Survey (SDSS) that fall in the HSC survey fields. We measure the bias factors of two quasar samples by evaluating the cross-correlation functions (CCFs) between the quasar samples and 25790 bright $z\sim4$ Lyman Break Galaxies (LBGs) in $M_{\rm 1450}<-21.25$ photometrically selected from the HSC dataset. Over an angular scale of \timeform{10.0"} to \timeform{1000.0"}, the bias factors are $5.93^{+1.34}_{-1.43}$ and $2.73^{+2.44}_{-2.55}$ for the low and high luminosity quasars, respectively, indicating no luminosity dependence of quasar clustering at $z\sim4$. It is noted that the bias factor of the luminous quasars estimated by the CCF is smaller than that estimated by the auto-correlation function (ACF) over a similar redshift range, especially on scales below \timeform{40.0"}. Moreover, the bias factor of the less-luminous quasars implies the minimal mass of their host dark matter halos (DMHs) is $0.3$-$2\times10^{12}h^{-1}M_{\odot}$, corresponding to a quasar duty cycle of $0.001$-$0.06$.

We present the luminosity function of z=4 quasars based on the Hyper Suprime-Cam Subaru Strategic Program Wide layer imaging data in the g, r, i, z, and y bands covering 339.8 deg^2. From stellar objects, 1666 z~4 quasar candidates are selected by the g-dropout selection down to i=24.0 mag. Their photometric redshifts cover the redshift range between 3.6 and 4.3 with an average of 3.9. In combination with the quasar sample from the Sloan Digital Sky Survey in the same redshift range, the quasar luminosity function covering the wide luminosity range of M1450=-22 to -29 mag is constructed. It is well described by a double power-law model with a knee at M1450=-25.36+-0.13 mag and a flat faint-end slope with a power-law index of -1.30+-0.05. The knee and faint-end slope show no clear evidence of redshift evolution from those at z~2. The flat slope implies that the UV luminosity density of the quasar population is dominated by the quasars around the knee, and does not support the steeper faint-end slope at higher redshifts reported at z>5. If we convert the M1450 luminosity function to the hard X-ray 2-10keV luminosity function using the relation between UV and X-ray luminosity of quasars and its scatter, the number density of UV-selected quasars matches well with that of the X-ray-selected AGNs above the knee of the luminosity function. Below the knee, the UV-selected quasars show a deficiency compared to the hard X-ray luminosity function. The deficiency can be explained by the lack of obscured AGNs among the UV-selected quasars.

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.

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.

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.

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.

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.

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 $\mu$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 $\mu$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 $\mu$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 $\phi$ = -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 $\rho_{\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.

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.

We investigate the UV-optical (longward of Ly$\alpha$ 1216\AA) 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 $\alpha_{\nu}^{\text{dif}}\sim +1/3$ (in the form of $f_{\nu}\propto \nu^{\alpha_{\nu}}$), which is significantly bluer than that of the composite spectrum $\alpha_{\nu}^{\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$\AA. These results are compared with previous studies, and the physical mechanisms of the variability of the continuum and emission lines are discussed.

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.

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).

The Ultra Luminous X-ray (ULX) source HLX-1 in the galaxy ESO 243-49 has an observed maximum unabsorbed X-ray luminosity of 1.3e42 erg/s (0.2-10.0 keV). From the conservative assumption that this value exceeds the Eddington limit by at most a factor of 10, the minimum mass is then 500 solar masses. The X-ray luminosity varies by a factor of 40 with an apparent recurrence timescale of approximately one year. This X-ray variability is associated with spectral state transitions similar to those seen in black hole X-ray binaries. Here we discuss our recent modelling of all the X-ray data for HLX-1 and show that it supports the idea that this ULX is powered by sub- and near Eddington accretion onto an intermediate mass black hole. We also present evidence for transient radio emission which is consistent with a discrete jet ejection event as well as comment on the nature of the environment around HLX-1 in light of recent Hubble Space Telescope photometry.

The hyper luminous X-ray source HLX-1 in the galaxy ESO 243-49, currently the best intermediate mass black hole candidate, displays spectral transitions similar to those observed in Galactic black hole binaries, but with a luminosity 100-1000 times higher. We investigated the X-ray properties of this unique source fitting multi-epoch data collected by Swift, XMM-Newton & Chandra with a disk model computing spectra for a wide range of sub- and super-Eddington accretion rates assuming a non-spinning black hole and a face-on disk (i = 0 deg). Under these assumptions we find that the black hole in HLX-1 is in the intermediate mass range (~2 x 10^4 M_odot) and the accretion flow is in the sub-Eddington regime. The disk radiation efficiency is eta = 0.11 +/-0.03. We also show that the source does follow the L_X ~ T^4 relation for our mass estimate. At the outburst peaks, the source radiates near the Eddington limit. The accretion rate then stays constant around 4 x 10^(-4) M_odot yr^(-1) for several days and then decreases exponentially. Such "plateaus" in the accretion rate could be evidence that enhanced mass transfer rate is the driving outburst mechanism in HLX-1. We also report on the new outburst observed in August 2011 by the Swift-X-ray Telescope. The time of this new outburst further strengthens the ~1 year recurrence timescale.

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.

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.

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.

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.

Formation of primordial black holes (PBHs) on astrophysical mass scales is a natural consequence of inflationary cosmology if the primordial perturbation spectrum has a large and negative running of the spectral index as observationally inferred today, because double inflation is required to explain it and fluctuations on some astrophysical scales are enhanced in the field oscillation regime in between. It is argued that PBHs thus produced can serve as intermediate-mass black holes (IMBHs) which act as the observed ultraluminous X-ray sources (ULXs) by choosing appropriate values of the model parameters in their natural ranges. Our scenario can be observationally tested in near future because the mass of PBHs is uniquely determined once we specify the values of the spectral index and its running on large scales.

In order to investigate how the growth of galactic bulges progresses with the growth of central black holes (BHs), we observed molecular gas (fuel for the coming star formation) in possibly young active galaxies, narrow-line Seyfert 1 galaxies (NLS1s). We present the results of radio observations of 12CO(1--0) using the Nobeyama Millimeter Array (with 2--4 kpc spatial resolution) for two FIR-bright NLS1s, yielding the first detection of their CO emission. Corresponding molecular--gas masses M(H2) of (1-3) X 109 Msun are the 2nd and 4th largest ones among NLS1s. By estimating dynamical masses and bulge masses M(bulge) for these two NLS1s using CO channel map and CO line widths, we found M(H2) amount to 0.13--0.35 of these masses. Taking account the star formation efficiency (~ 0.1), the increase in M(bulge) in those NLS1s in the near future (~< 10^{7.5} yr) is expected not to be a huge fraction (1--5% of the preexisting stars). Bulge growth may have finished before BH growth, or bulge--BH coevolution may proceed with many, occasional discrete events, where one coevolution event produces only a small amount of mass growth of BHs and of bulges. We also discuss the ratios of star-formation rate--to--gas accretion rate onto BHs, finding that two NLS1s have very small ratios (~ 1) compared with the M(bulge)/M(BH) ratios found in active and inactive galaxies (~ 700). This huge difference suggests either the non-overlapped coevolution, long star formation duration or temporarily high accretion rate during NLS1 phase.

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.

We have analyzed the archival XMM-Newton data of the bright Ultra-Luminous X-ray Source (ULX) M82 X-1 with an 105 ksec exposure when the source was in the steady state. Thanks to the high photon statistics from the large effective area and long exposure, we were able to discriminate different X-ray continuum spectral models. Neither the standard accretion disk model (where the radial dependency of the disk effective temperature is T(r) \propto r^-3/4) nor a power-law model gives a satisfactory fit. In fact, observed curvature of the M82 X-1 spectrum was just between those of the two models. When the exponent of the radial dependence (p in T(r) \propto r^-p) of the disk temperature is allowed to be free, we obtained p =0.61^+0.03_-0.02. Such a reduction of p from the standard value 3/4 under extremely high mass accretion rates is predicted from the accretion disk theory as a consequence of the radial energy advection. Thus, the accretion disk in M82 X-1 is considered to be in the Slim disk state, where an optically thick Advection Dominant Accretion Flow (ADAF) is taking place. We have applied a theoretical slim disk spectral model to M82 X-1, and estimated the black hole mass ~ 19-32 M_odot. We propose that M82 X-1 is a relatively massive stellar black hole which has been produced through evolution of an extremely massive star, shining at a super-Eddington luminosity by several times the Eddington limit.

We present optical images of nearby 50 narrow-line Seyfert 1 galaxies (NLS1s) which cover all the NLS1s at z<0.0666 and $\delta \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$\beta$ 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.

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$\beta$ 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$\beta$ 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$\beta$ profile is sensitive primarily to Eddington ratio, but that inclination effects play a role in some cases. (abridged)

A comparative XMM-Newton archival data spectral study of the ultraluminous X-ray source (ULX) M33 X-8 has been performed by using both the standard disk model and the newly developed slim disk models. The results of this analysis tend to confirm the hypothesis that M33 X-8 is an X-ray binary with a stellar-mass black hole accreting at super-Eddington rate.

We have obtained optical intermediate resolution spectra (R = 3000) of the narrow-line quasars DMS 0059-0055 and PG 1543+489. The [O III] emission line in DMS 0059-0055 is blueshifted by 880 km/s relative to Hbeta. We also confirm that the [O III] emission line in PG 1543+489 has a relative blueshift of 1150 km/s. These two narrow-line quasars show the largest [O III] blueshifts known to date among type 1 active galactic nuclei (AGNs). The [O III] emission lines in both objects are broad (1000 - 2000 km/s) and those in DMS 0059-0055 show strong blue asymmetry. We interpret the large blueshift and the profile of the [O III] lines as the result of an outflow interacting with circumnuclear gas. Among type 1 AGNs with large blueshifted [O III], there is no correlation between the Eddington ratios and the amount of [O III] blueshifts. Combining our new data with published results, we confirm that the Eddington ratios of the such AGNs are the highest among AGNs with the same black hole masses. These facts suggest that the Eddington ratio is a necessary condition or the [O III] blueshifts weakly depend on the Eddington ratio. Our new sample suggests that there are possible necessary conditions to produce an outflow besides a high Eddington ratio: large black hole mass (> 10^7 M_solar) or high mass accretion rate (> 2 M_solar/yr) or large luminosity (lambda L_{lambda} (5100A) > 10^44.6 erg/s).

Using the BH masses deduced from the empirical relation of Kaspi et al. (2000) and assuming that the optical luminosity is provided by the accretion disc, we show that Narrow Line Seyfert Galaxies 1 (NLS1s) accrete at super-Eddington rates, while their luminosity stays of the order of the Eddington limit. We take into account the possibility of a non-viscous energy release in the gravitationally unstable region of the disc. It leads to a smaller accretion rate and to a redder continuum than a standard disc, which agrees better with the observations. The observed bolometric luminosities appear to saturate at a few times the Eddington luminosity for super-Eddington accretion rates, as predicted by slim disc models. The accretion rate stays always of the order of a few M$_{\odot}$/yr, indicating that the growing of the BH is mass supply limited . Since the masses of the BH increases by one order of magnitude in a few 10$^7$ years, it could explain why NLS1s appear to not follow the same BH - bulge relation as other galaxies. NLS1s should thus play an important role in shaping the mass function of local BHs. We discuss the possibility that the masses could be systematically underestimated due to an inclination effect, and we conclude that the accretion rates could thus be strongly overestimated, but only in a small proportion of objects.

Narrow-Line Seyfert 1 galaxies (NLS1s) and Narrow-Line quasars (NLQs) seem to amount to ~ 10-30 % of active galactic nuclei (AGNs) in the local universe. Together with their average accretion rate, we argue that the black hole (BH) growth by factor of 8-800 happens in these super-Eddington accretion phase of AGNs. Moreover, there is a possible, systematic underestimation of accretion rates (in the Eddington unit) due to an overestimation of BH mass by massive accretion discs for super-Eddington objects. If it is true, the factor of BH growth above may be larger by order(s) of magnitude. In contrast, the growth factor expected in sub-Eddington phase is only ~ 2. Therefore, the cosmic BH growth by accretion is likely dominated by super-Eddington phase, rather than sub-Eddington phase which is the majority among AGNs. This analysis is based on the fraction and the average accretion rate of NLS1s and NLQs obtained for z ~< 0.5. If those numbers are larger at higher redshift (where BHs were probably less grown), super-Eddington accretion would be even more important in the context of cosmic BH growth history.

Self-gravitating accretion discs have only been studied in a few nearby objects using maser spots at the parsec-scale. We find a new spectral window for observing the self-gravitating accretion disc in super-Eddington accreting Active Galactic Nuclei (AGNs). This window is determined by calculating the outermost radius (rsg) of a non self-gravitating disc and the corresponding emission wavelength (lsg) as a function of various disc parameters. We find that lsg reaches ~ 4000 AA for alpha=0.1, when Mdot > 70 (M_BH / 10^7 M_Sun)^-1 L_Edd/c^2 (where alpha, Mdot, M_BH and L_Edd are, respectively, the viscosity parameter, gas accretion rate onto the central black hole (BH), the BH mass and the Eddington luminosity). Moreover, lsg is as small as 1500 AA for alpha = 0.001, which is the smallest alpha case in this study. Therefore, the window for observing the self-gravitating part of an AGN accretion disc is from 2 micron to lsg. Incidentally, rsg can be less than the photon trapping radius for Mdot >~ 10^3.3 L_Edd/c^2. Namely, a self-gravitating, optically-thick, advection-dominated accretion disc is expected to appear in the extremely high accretion rate regime. Next, we demonstrate that the Mid-Infrared to X-ray spectrum of a bright, well-studied Narrow-Line Seyfert 1 galaxy, Ton S 180, is indeed well fitted by the spectrum arising from the following three components: an inner slim disc (with a corona), an outer, self-gravitating non-Keplerian disc and a dusty torus. The total mass, BH mass plus the entire disc mass, is found to be about (1.4 - 8.0) M_BH. If the surface density varies with radius r in proportion to r^-0.6, the total mass is consistent with the central mass estimated by H_beta and [O III] widths.

We present new interpretation of Ton S180 spectrum obtained by {\it Chandra} Spectrometer (Low Energy Transmission Grating). Several narrow absorption lines and a few emission disk lines have been successfully fitted to the data. We have not found any significant edges accompanying line emission. We propose the interpretation of narrow lines consistent with the paper recently written by Krolik (2002), where warm absorber is strongly inhomogeneous. Such situation is possible in so called multi-phase medium, where regions with different ionization states, densities and temperatures may coexist in thermal equilibrium under constant pressure. We illustrate this scenario with theoretical spectra of radiation transfered through a stratified cloud with constant pressure (instead of constant density) computed by code {\sc titan} in plane parallel approximation. Detected spectral features are faint and their presence do not alter the broad band continuum. We model the broad band continuum of Ton S180 assuming an irradiated accretion disk with a dissipative warm skin. The set of parameters appropriate for the data cannot be determined uniquely but models with low values of the black hole mass have too hot and radially extended warm skin to explain the formation of soft X-ray disk lines seen in the data.

Super-Eddington accretion onto black holes (BHs) may occur at ultraluminous compact X-ray sources in nearby galaxies, Galactic microquasars and narrow-line Seyfert 1 galaxies (NLS1s). Effects of electron scattering (opacity and Comptonization) and the relativistic correction (gravitational redshift and transverse Doppler effect) on the emergent spectra from super-Eddington accretion flows onto non-rotating BHs are examined for 10^{1.5} and 10^{6.5} M_sun BH masses (M_BH). With mdot [= Mdot / (L_Edd / c^2), where Mdot is the accretion rate] > 100, the spectral hardening factor via electron scattering is \lesssim 2.3 - 6.5. Due to the mdot-sensitive hardening factor, the color temperature of the innermost radiation is not proportional to L^{0.25}, differing from the simplest standard accretion disk. The model is applied to optical--soft X-ray emission from NLS1s. We pick up one NLS1, namely PG 1448+273 with an inferred M_BH of 10^{6.4} M_sun, among the highest mdot candidates. The broadband spectral distribution is successfully reproduced by the model with an extremely high mdot (= 1000) and the viscosity parameter alpha of 0.01. This implies that this object, as well as some other highest mdot systems, is really young: the inferred age, M_BH / Mdot, is about 10^6 years. We also briefly discuss the distribution of mdot for transient and highly variable NLS1s, finding that those are located at 3 \lesssim mdot \lesssim 300. Such a moderately high accretion rate is indicative of thermal instability. Furthermore, mdot for a possible type-2 counterpart of NLS1s, NGC 1068, is found to be similar to mdot for NLS1s.

We present preliminary results from a long Chandra HETG observation of the luminous Narrow-line Seyfert 1 galaxy 1H 0707-495. We find a complex X-ray spectrum comprised of a two-component continuum with superimposed emission and absorption lines. The short time scale X-ray variability is different than observed in other AGN: the soft X-rays vary markedly less than the hard X-rays. This behavior is similar to that of high-state Galactic black holes. We also investigated the long time scale variability, and discovered an apparent bimodal flux distribution. We postulate that the bimodality is the signature of the radiation pressure instability, and note that this instability may be expected in luminous NLS1s.

We extended the disk corona model (Meyer & Meyer-Hofmeister 1994; Meyer, Liu, & Meyer-Hofmeister 2000a) to the inner region of galactic nuclei by including different temperatures in ions and electrons as well as Compton cooling. We found that the mass evaporation rate and hence the fraction of accretion energy released in the corona depend strongly on the rate of incoming mass flow from outer edge of the disk, a larger rate leading to more Compton cooling, less efficient evaporation and a weaker corona. We also found a strong dependence on the viscosity, higher viscosity leading to an enhanced mass flow in the corona and therefore more evaporation of gas from the disk below. If we take accretion rates in units of the Eddington rate our results become independent on the mass of the central black hole. The model predicts weaker contributions to the hard X-rays for objects with higher accretion rate like narrow-line Seyfert 1 galaxies (NLS1s), in agreement with observations. For luminous active galactic nuclei (AGN) strong Compton cooling in the innermost corona is so efficient that a large amount of additional heating is required to maintain the corona above the thin disk.

Cosmological evolution of the QSO luminosity functions (LFs) at NIR/optical/X-ray bands for 1.3 < z < 3.5 is investigated based on the realistic QSO spectra. The accretion-disk theory predicts that although QSO luminosities only depend on mass-accretion rate, \Mdot, QSO spectra have a dependence on black-hole mass, M_{BH}, as well. The smaller M_{BH} is and/or the larger \Mdot is, the harder becomes the QSO NIR/optical/UV spectrum. We model disk spectra which can reproduce these features and calculated LFs for redshift z ~ 3 with the assumption of new-born QSOs being shining at the Eddington luminosity. The main results are: (i) the observed LFs at optical and X-rays can be simultaneously reproduced. (ii) LFs at optical and X-ray bands are not sensitive to M_{BH}, while LFs at NIR bands are; about one order of magnitude difference is expected in volume number densities at L_{I, J} ~ 10^{46} erg s^{-1} between the case that all QSOs would have the same spectral shape as that of M_{BH} = 10^{9} M_sun and the case with M_{BH} = 10^{11} M_sun. (iii) The resultant LFs at NIR are dominated by 10^{7} M_sun black-holes at L_{I, J} ~ 10^{44} erg s^{-1}, and by 10^{11} M_sun black-holes at L_{I, J} \~ 10^{46} erg s^{-1}. Future infrared observations from space(e.g.NGST) will probe cosmological evolution of black hole masses. For redshift z < 3, on the other hand, the observed optical/X-ray LFs can be fitted, if the initial QSO luminosity L_0 is below the Eddington luminosity. Interestingly, the best fitting values of l = L_0/L_{Edd} are different in B- and X-ray bands; l_B ~ 2.5 l_X. The reason for this discrepancy is briefly discussed.

We construct a disk-corona model which accounts for the optical-to-X-ray spectral properties of Seyfert nuclei and QSOs. We study emission spectrum emerging from a vertical disk-corona structure composed of two-temperature plasma by solving hydrostatic equilibrium and radiative transfer self-consistently. A fraction $f$ of viscous heating is assumed to be dissipated in a corona with a Thomson optical depth of $\tau_c$, where advective cooling is also included, and a remaining fraction, $1-f$, dissipates within a main body of the disk. Our model can nicely reproduce the soft X-ray excess with a power-law shape and the hard tail extending to $\sim$ 50 keV. The different spectral slopes ($\alpha \sim$ 1.5 below 2keV and $\sim$ 0.5 above, where $\Fnu \propto \nu^{- \alpha}$) are the results of different emission mechanisms and different sites; the former slope is due to unsaturated Comptonization from the innermost zone and the latter is due to a combination of the Comptonization, bremsstrahlung and a reflection of the coronal radiation at the disk-corona boundary from the inner to surrounding zone ($\leq$ 300 Schwarzschild radii). The emergent optical spectrum is redder ($\alpha \sim 0.3$) than that of the standard disk ($\alpha \sim -0.3$), being consistent with observations, due to the different efficiencies of spectral distortion of disk emission at different radii. Further, we find that the cut-off frequency of the hard X-ray ($\sim$ coronal electron temperature) and broad-band spectral shape are insensitive to the black-hole mass, while the peak frequency of the big blue bump is sensitive to the mass as the peak frequency $\propto \Mbh^{-1/4}$.

To explain the broad-band spectral energy distributions (SED) of Seyfert nuclei and QSOs, we study the emission spectrum emerging from a vertical disk-corona structure composed of a two-temperature plasma by solving hydrostatic equilibrium and radiative transfer self-consistently. Our model can nicely reproduce the soft X-ray excess with alpha (L_frequency proportional to the frequency to the power of -alpha) of about 1.5 and the hard tail extending to about 50 keV with alpha about 0.5. The different spectral slopes (alpha about 1.5 below 2 keV and about 0.5 above) are the results of different emission mechanisms: unsaturated Comptonization in the former and a combination of Comptonization, bremsstrahlung, and reflection of the coronal radiation at the disk-corona boundary in the latter.

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.

Rapid fluctuation with a frequency dependence of $1/f^{\alpha}$ (with $\alpha \simeq 1 - 2$) is characteristic of radiation from black-hole objects. Its origin remains poorly understood. We examine the three-dimensional magnetohydrodynamical (MHD) simulation data, finding that a magnetized accretion disk exhibits both $1/f^\alpha$ fluctuation (with $\alpha \simeq 2$) and a fractal magnetic structure (with the fractal dimension of $D \sim 1.9$). The fractal field configuration leads reconnection events with a variety of released energy and of duration, thereby producing $1/f^\alpha$ fluctuations.

Number of monitoring observations of continuum emission from Active Galactic Nuclei (AGNs) have been made in optical--X-ray bands. The results obtained so far show (i) random up and down on timescales longer than decades, (ii) no typical timescales of variability on shorter timescales and (iii) decreasing amplitudes as timescales become shorter. The second feature indicates that any successful model must produce a wide variety of shot-amplitudes and -durations over a few orders in their light curves. In this sense, we conclude that the disk instability model is favored over the starburst model, since fluctuations on days are hard to produce by the latter model. Inter-band correlations and time lags also impose great constraints on models. Thus, constructing wavelength and time dependent models remains as a future work.

Aperiodic optical variability is a common property of Active Galactic Nuclei (AGNs), though its physical origin is still open to question. To study the origin of the optical -- ultraviolet variability in AGN, we compare light curves of two models to observations of quasar 0957+561 in terms of a structure function analysis. In the starburst (SB) model, random superposition of supernovae in the nuclear starburst region produce aperiodic luminosity variations, while in the disk-instability (DI) model, variability is caused by instabilities in the accretion disk around a supermassive black hole. We calculate fluctuating light curves and structure functions, $V(\tau)$, by simple Monte-Carlo simulations on the basis of the two models. Each resultant $V(\tau)$ possesses a power-law portion, $[V(\tau)]^{1/2} \propto \tau^{\beta}$, at short time lags ($\tau$). The two models can be distinguished by the logarithmic slope, $\beta$; $\beta \sim$ 0.74--0.90 in the SB model and $\beta \sim$ 0.41--0.49 in the DI model, while the observed light curves exhibit $\beta \sim$ 0.35. Therefore, we conclude that the DI model is favored over the SB model to explain the slopes of the observational structure function, in the case of 0957+561, though this object is a radio-loud object and thus not really a fair test for the SB model. In addition, we examine the time-asymmetry of the light curves by calculating $V(\tau)$ separately for brightening and decaying phases. The two models exhibit opposite trends of time-asymmetry to some extent, although the present observation is not long enough to test this prediction.