The supermassive black holes (SMBHs) of narrow-line Seyfert 1 galaxies (NLS1s) are at the lowest end of mass function of active galactic nuclei (AGNs) and preferentially reside in late-type host galaxies with pseudobulges, which are thought to be formed by internal secular evolution. On the other hand, the population of radio-loud NLS1s presents a challenge for the relativistic jet paradigm that powerful radio jets are exclusively associated with very high mass SMBHs in elliptical hosts, which are built-up through galaxy mergers. We investigated distorted radio structures associated with the nearest gamma-ray emitting, radio-loud NLS1 1H 0323+342. This provides supporting evidence for the merger hypothesis based on the past optical/near-infrared observations of its host galaxy. The anomalous radio morphology consists of two different structures, the inner curved structure of currently active jet and the outer linear structure of low-brightness relics. Such a coexistence might be indicative of the stage of an established black hole binary with precession before the black holes coalesce in the galaxy merger process. 1H 0323+342 and other radio-loud NLS1s under galaxy interactions may be extreme objects on the evolutionary path from radio-quiet NLS1s to normal Seyfert galaxies with larger SMBHs in classical bulges through mergers and merger-induced jet phases.

The nearby bright radio galaxy 3C 84 at the center of Perseus cluster is one of the ideal targets to explore the jet in active galactic nuclei (AGNs) and its parsec-scale environment. A recent research of Fujita & Nagai revealed the existence of the northern counter-jet component (N1) located at $\sim$ 2 mas north from the central core in VLBI images at 15 and 43 GHz and they are explained by the free-free absorption (FFA) due to an ionized plasma foreground. Here we report a new quasi-simultaneous observation of 3C 84 with the Korean VLBI Network (KVN) at 86 GHz and the KVN and VERA (VLBI Exploration of Radio Astrometry) Array (KaVA) at 43 GHz in 2016 February. We succeeded the first detection of N1 at 86 GHz and the data show that N1 still has an inverted spectrum between 43 and 86 GHz with its spectral index $\alpha$ ($S_{\nu} \propto \nu^{\alpha}$) of $1.19 \pm 0.43$, while the approaching lobe component has the steep spectrum with the index of $-0.54 \pm 0.30$. Based on the measured flux asymmetry between the counter and approaching lobes, we constrain the averaged number density of the FFA foreground $n_{\rm e}$ as $1.8 \times 10^{4}~{\rm cm^{-3}} \lesssim n_{\rm e} \lesssim 1.0 \times 10^{6}~{\rm cm^{-3}}$. Those results suggest that the observational properties of the FFA foreground can be explained by the dense ionized gas in the circumnuclear disk and/or assembly of clumpy clouds at the central $\sim 1$ pc region of 3C 84.

We study the obscuring structure of circumnuclear disks (CNDs) by considering supernova (SN) feedbacks from nuclear starburst and the effect of anisotropic radiative pressure from AGNs. We suppose that the mass accretion onto a central supermassive black hole (SMBH) is triggered by SN-driven turbulence within CNDs, and we explore how the structures of CNDs depend on the BH mass ($M_{\rm BH}$) and AGN luminosity ($L_{\rm AGN}$). We find that the obscuring fraction ($f_{\rm obs}$) peaks at $\sim10\%$ of the Eddington luminosity ($L_{\rm Edd}$), and its maximal value is $f_\mathrm{obs} \sim 0.6$ for less massive SMBHs (e.g., $M_{\rm BH} < 10^{8}M_{\odot}$). This is because the scale height of CNDs is determined by the SN-driven accretion for a smaller $L_{\rm AGN}$, while the dusty molecular gas in CNDs is blown away by the radiation pressure from AGNs beyond the critical luminosity. On the other hand, for massive SMBHs (e.g., $M_{\rm BH} > 10^{8}M_{\odot}$), $f_{\rm obs}$ is always smaller than $0.2$, and it is almost independent of $L_{\rm AGN}$ because the scale height of CNDs is mainly controlled by the maximal star-formation efficiency ($C_{\rm *, max}$) in CNDs. By comparison with the obscuring fractions suggested from the mid-infrared observations of nearby AGNs, the SN plus radiative feedback model with $C_{\rm *, max} = 10^{-7}\, {\rm yr}^{-1}$ well reproduces the observations for $M_\mathrm{BH} = 10^8 M_\odot$. We also find that the intense starburst or the existence of dust-free absorbers inside CNDs are necessary, to explain X-ray observations.

We present ALMA Band 6 observations of the CO(2-1), HCN(3-2), and HCO$^{+}$(3-2) lines in the nearby radio galaxy / brightest cluster galaxy (BCG) of NGC 1275 with the spatial resolution of $\sim20$ pc. In the previous observations, CO(2-1) emission was detected as radial filaments lying in the east-west direction. We resolved the inner filament and found that the filament cannot be represented by a simple infalling stream both morphologically and kinematically. The observed complex nature of the filament resembles the cold gas structure predicted by recent numerical simulations of cold chaotic accretion. A crude estimate suggests that the accretion rate of the cold gas can be higher than that of hot gas. Within the central 100 pc, we detected a rotational disk of the molecular gas whose mass is $\sim10^{8} M_{\sun}$. This is the first evidence of the presence of massive cold gas disk on this spatial scale for BCGs. The disk rotation axis is approximately consistent with the axis of the radio jet on subpc scales. This probably suggests that the cold gas disk is physically connected to the innermost accretion disk which is responsible for jet launching. We also detected absorption features in the HCN(3-2) and HCO$^{+}$(3-2) spectra against the radio continuum emission mostly radiated by $\sim1.2$-pc size jet. The absorption features are blue-shifted from the systemic velocity by $\sim$300-600~km~s$^{-1}$, which suggests the presence of outflowing gas from the active galactic nucleus (AGN). We discuss the relation of the AGN feeding with cold accretion, the origin of blue-shifted absorption, and estimate of black hole mass using the molecular gas dynamics.

The narrow-line Seyfert 1 (NLS1) galaxy PMN J0948+0022 is an archetype of gamma-ray-emitting NLS1s in active galactic nuclei (AGNs). In this study, we investigate its radio structures using archival data obtained using the Karl G.~Jansky very large array (VLA) and the very long baseline array (VLBA). The new VLA images reveal an outermost radio emission separated by 9.1 arcsec. Its resolved structure and steep spectrum suggest a terminal shock in a radio lobe energized by the jet from the PMN J0948+0022 nucleus. This large-scale radio component is found at almost the same position angle as that of the pc-scale jet, indicating a stable jet direction up to ~1 Mpc. Its apparent one-sidedness implies a moderate advancing speed ($\beta>0.27$). The kinematic age of 1 x 10^7 year of the jet activity is consistent with the expected NLS1 phase of ~10^7 year in the AGN lifetime. The VLBA image reveals the jet structure at distances ranging from r=0.82 milliarcsec to 3.5 milliarcsec, corresponding to approximately 100 pc, where superluminal motions were found. The jet width profile ($\propto r^{1.12}$) and flux-density profile ($\propto r^{-1.44}$) depending on the distance $r$ along the jet suggest that the jet kinetic energy is converted into internal energy in this region. The jet is causally connected in a nearly conical streamline, which is supported by ambient pressure at 100 pc scales in the host galaxy of PMN J0948+0022.

Radio jets in active galaxies have been expected to interact with circumnuclear environments in their early phase evolutions. By performing the multi-epoch monitoring observation with the KVN and VERA Array (KaVA) at 43~GHz, we investigate the kinematics of the notable newborn bright component C3 located at the tip of the recurrent jet of 3C~84. During 2015 August-September, we discover the flip of C3 and the amount of the flip is about 0.4~milli-arcsecond in angular scale, which corresponds to 0.14 parsec in physical scale. After the flip of C3, it wobbled at the same location for a few months and then it restarted to propagate towards the southern direction. The flux density of C3 coherently showed the monotonic increase during the observation period. The flip is in good agreement with hydrodynamical simulations of jets in clumpy ambient medium. We estimate the number density of the putative clump based on the momentum balance between the jet thrust and the ram pressure from the clump and it is about $10^{3-5}~{\rm cm^{-3}}$. We briefly discuss possible origins of the clump.

We present the result of our spectroscopic follow-up observation for faint quasar candidates at z~5 in a part of the Canada-France-Hawaii Telescope Legacy Survey wide field. We select nine photometric candidates and identify three z~5 faint quasars, one z~4 faint quasar, and a late-type star. Since two faint quasar spectra show Civ emission line without suffering from a heavy atmospheric absorption, we estimate the black hole mass (M$_{BH}$) and Eddington ratio (L/L$_{Edd}$) of them. The inferred log M$_{BH}$ are 9.04+/-0.14 and 8.53+/-0.20, respectively. In addition, the inferred log (L/L$_{Edd}$) are -1.00+/-0.15 and -0.42+/-0.22, respectively. If we adopt that L/L$_{Edd}$= constant or $\propto$ (1+z)^2, the seed black hole masses (M$_{seed}$) of our z~5 faint quasars are expected to be >10^5 M$_\odot$ in most cases. We also compare the observational results with a mass accretion model where angular momentum is lost due to supernova explosions (Kawakatu & Wada 2008). Accordingly, M$_{BH}$ of the z~5 faint quasars in our sample can be explained even if M$_{seed}$ is ~10^3M$_\odot$. Since z~6 luminous qusars and our z~5 faint quasars are not on the same evolutionary track, z~6 luminous quasars and our z~5 quasars are not the same populations but different populations, due to the difference of a period of the mass supply from host galaxies. Furthermore, we confirm that one can explain M$_{BH}$ of z~6 luminous quasars and our z~5 faint quasars even if their seed black holes of them are formed at z~7.

We used the Atacama Large Millimeter/submillimeter Array (ALMA) to map the CO(3-2) and the underlying continuum emissions around the type 1 low-luminosity active galactic nucleus (LLAGN; bolometric luminosity $\lesssim 10^{42}$ erg~s$^{-1}$) of NGC 1097 at $\sim 10$ pc resolution. These observations revealed a detailed cold gas distribution within a $\sim 100$ pc of this LLAGN. In contrast to the luminous Seyfert galaxy NGC 1068, where a $\sim 7$ pc cold molecular torus was recently revealed, a distinctively dense and compact torus is missing in our CO(3-2) integrated intensity map of NGC 1097. Based on the CO(3-2) flux, the gas mass of the torus of NGC 1097 would be a factor of $\gtrsim 2-3$ less than that found for NGC 1068 by using the same CO-to-H$_2$ conversion factor, which implies less active nuclear star formation and/or inflows in NGC 1097. Our dynamical modeling of the CO(3-2) velocity field implies that the cold molecular gas is concentrated in a thin layer as compared to the hot gas traced by the 2.12 $\mu$m H$_2$ emission in and around the torus. Furthermore, we suggest that NGC 1097 hosts a geometrically thinner torus than NGC 1068. Although the physical origin of the torus thickness remains unclear, our observations support a theoretical prediction that geometrically thick tori with high opacity will become deficient as AGNs evolve from luminous Seyferts to LLAGNs.

We explore physical properties of the shocked external medium (i.e., a shell) in 3C 84 associated with the recurrent radio lobe born around 1960. In the previous work of Ito et al., we investigated a dynamical and radiative evolution of such a shell after the central engine stops the jet launching and we found that a fossil shell emission overwhelms that of the rapidly fading radio lobe. We apply this model to 3C 84 and find the followings: (i) The fossil shell made of shocked diffuse ambient matter with the number density of 0.3 cm$^{-3}$ radiates bright Inverse-Compton (IC) emission with the seed photons of the radio emission from the central compact region and the IC emission is above the sensitivity threshold of the Cherenkov Telescope Array (CTA). (ii) When the fossil shell is produced in a geometrically thick ionized plasma with the number density of $10^{3}$ cm$^{-3}$ and the field strength in the shell may reach about 17 mG in the presence of magnetic fields amplification and the radio emission becomes comparable to the sensitivity of deep imaging VLBI observations. A possible production of ultra high energy cosmic-rays (UHECRs) in the dense shocked plasma is also argued.

We present a positive correlation between the mass of dense molecular gas ($M_{\rm dense}$) of $\sim 100$ pc scale circumnuclear disks (CNDs) and the black hole mass accretion rate ($\dot{M}_{\rm BH}$) in total 10 Seyfert galaxies, based on data compiled from the literature and an archive (median aperture $\theta_{\rm med}$ = 220 pc). A typical $M_{\rm dense}$ of CNDs is 10$^{7-8}$ $M_\odot$, estimated from the luminosity of the dense gas tracer, the HCN($1-0$) emission line. Because dense molecular gas is the site of star formation, this correlation is virtually equivalent to the one between nuclear star formation rate and $\dot{M}_{\rm BH}$ revealed previously. Moreover, the $M_{\rm dense}-\dot{M}_{\rm BH}$ correlation was tighter for CND-scale gas than for the gas on kpc or larger scales. This indicates that CNDs likely play an important role in fueling black holes, whereas $>$kpc scale gas does not. To demonstrate a possible approach for studying the CND-scale accretion process with the Atacama Large Millimeter/submillimeter Array (ALMA), we used a mass accretion model where angular momentum loss due to supernova explosions is vital. Based on the model prediction, we suggest that only the partial fraction of the mass accreted from the CND ($\dot{M}_{\rm acc}$) is consumed as $\dot{M}_{\rm BH}$. However, $\dot{M}_{\rm acc}$ agrees well with the total nuclear mass flow rate (i.e., $\dot{M}_{\rm BH}$ + outflow rate). Although these results are still tentative with large uncertainties, they support the view that star formation in CNDs can drive mass accretion onto supermassive black holes in Seyfert galaxies.

We investigate the ability of jets in active galactic nuclei (AGNs) to break out of the ambient gas with sufficiently large advance velocities. Using observationally estimated jet power, we analyze 28 bright elliptical galaxies in nearby galaxy clusters. Because the gas density profiles in the innermost regions of galaxies have not been resolved so far, we consider two extreme cases for temperature and density profiles. We also follow two types of evolution for the jet cocoons: being driven by the pressure inside the cocoon (Fanaroff-Riley [FR] I type), and being driven by the jet momentum (FR II type). Our main result is that regardless of the assumed form of density profiles, jets with observed powers of <~ 10^{44} erg s^{-1} are not powerful enough to evolve as FR II sources. Instead, they evolve as FR I sources and appear to be decelerated below the buoyant velocities of the cocoons when jets were propagating through the central dense regions of the host galaxies. This explains the reason why FR I sources are more frequent than FR II sources in clusters. Furthermore, we predict the sizes of X-ray cavities from the observed jet powers and compare them with the observed ones --- they are consistent within a factor of two if the FR I type evolution is realized. Finally, we find that the jets with a power >~ 10^{44} erg s^{-1} are less affected by the ambient medium, and some of them, but not all, could serve as precursors of the FR II sources.

(Abridged) We propose a new method to estimate the Bondi (hot gas) accretion rates, \dot M_B, onto the supermassive black holes (SMBHs) at the centres of elliptical galaxies. It can be applied even if the Bondi radius is not well-resolved in X-ray observations. This method is based on two simple assumptions: (1) hot gas outside the Bondi radius is in nearly a hydrostatic equilibrium, and (2) the gas temperature near the galaxy centre is close to the virial temperature. We apply this method to 28 bright elliptical galaxies in nearby galaxy clusters (27 of them are the brightest cluster galaxies; BCGs). We find a correlation between the Bondi accretion rates and the power of jets associated with the SMBHs over four orders of magnitude in \dot M_B. For most galaxies, the accretion rates are large enough to account for the jet powers. Our results indicate that the feedback from the active galactic nuclei (AGN) correlates with the properties of the hot gas surrounding the SMBHs. We also find that more massive SMBHs in BCGs tend to have larger specific growth rates. This may explain the hyper masses (~ 10^10 M_sun) of some of the SMBHs. Comparison between the accretion rates and the X-ray luminosities of the AGN suggests that the AGN in the BCGs are extremely radiatively inefficient compared with X-ray binaries in the Milky Way, even when their Eddington accretion ratio, \dot M_B/\dot M_Edd, exceeds 0.01. Lastly, we find a tight correlation between the Bondi accretion rates and the X-ray luminosities of cool cores. Their relation is linear and the power generated by the Bondi accretion is large enough to compensate the radiative cooling of the cool cores. Although the `classical' Bondi accretion model is a greatly oversimplified one, the correlations we find here demonstrate that the accretion onto the SMBHs reflects broadly the properties of the Bondi accretion.

We examine the plasma composition of relativistic jets in four FRII radio galaxies by analyzing the total cocoon pressure in terms of partial pressures of thermal and non-thermal electrons/positrons and protons. The total cocoon pressure is determined by cocoon dynamics via comparison of theoretical model with the observed cocoon shape. By inserting the observed number density of non-thermal electrons/positrons and the upper limit of thermal electron/positron number density into the equation of state, the number density of protons is constrained. We apply this method to four FRII radio galaxies (Cygnus A, 3C 219, 3C 223 and 3C 284), for which the total cocoon pressures have been already evaluated. We find that the positron-free plasma comprising of protons and electrons is ruled out, when we consider plausible particle distribution functions. In other words, the mixture of positrons is required for all four FRII radio galaxies; the number density ratio of electrons/positrons to protons is larger than two. Thus, we find that the plasma composition is independent of the jet power and the size of cocoons. We also investigate the additional contribution of thermal electrons/positrons and protons on the cocoon dynamics. When thermal electrons/positrons are absent, the cocoon is supported by the electron/ proton plasma pressure, while both electron/positron pressure supported and electron/proton plasma pressure supported cocoons are allowed if the number density of thermal electrons/positrons is about 10 times larger than that of non-thermal ones.

We analyse the environment of the supermassive black hole (SMBH) in the centre of a massive elliptical galaxy NGC 1275 in the Perseus cluster, hosting the radio source 3C 84. We focus on the young radio lobe observed inside the estimated Bondi accretion radius. We discuss the momentum balance between the jet associated with the lobe and the surrounding gas. The results are compared with the proper motion of the radio lobe obtained with the very long baseline interferometry. We find that under assumption of a high-density environment >~ 100 cm^-3), the jet power must be comparable to the Eddington luminosity --- this is clearly inconsistent with the current moderate activity of 3C 84, which indicates instead that the jet is expanding in a very low density region (<~1 cm^-3), along the rotation axis of the accretion flow. The power required for the jet to expand in the low-density environment is comparable to the past average jet power estimated from the X-ray observations. We estimate the classical Bondi accretion rate, assuming that (1) gas accretion is spherically symmetric, (2) accretion is associated with the jet environment, and (3) the medium surrounding the jet is representative of the properties of the dominant accreting gas. We find that Bondi accretion is inconsistent with the estimated jet power. This means that either accretion of the cold gas in the NGC 1275 is more efficient than that of the hot gas, or the jets are powered by the SMBH spin.

We investigate shell emission associated with dying radio loud AGNs. First, based on our recent work by Ito et al. (2015), we describe the dynamical and spectral evolutions of shells after stopping the jet energy injection. We find that the shell emission overwhelms that of the radio lobes soon after stopping the jet energy injection because fresh electrons are continuously supplied into the shell via the forward shock while the radio lobes rapidly fade out without jet energy injection. We find that such fossil shells can be a new class of target sources for SKA telescope. Next, we apply the model to the nearby radio source 3C84. Then, we find that the fossil shell emission in 3C84 is less luminous in radio band while it is bright in TeV gamma-ray band and it can be detectable by CTA.

We examine the fate of a dead radio source in which jet injection from the central engine has stopped at an early stage of its evolution ($t = t_j \lesssim 10^5$ yr). To this aim, we theoretically evaluate the evolution of the emission from both the lobe and the shell, which are composed of shocked jet matter and a shocked ambient medium, respectively. Based on a simple dynamical model of expanding lobe and shell, we clarify how the broadband spectrum of each component evolves before and after the cessation of the jet activity. It is shown that the spectrum is strongly dominated by the lobe emission while the jet is active ($t \leq t_j$). On the other hand, once the jet activity has ceased ($t > t_j$), the lobe emission fades out rapidly, since fresh electrons are no longer supplied from the jet. Meanwhile, shell emission only shows a gradual decrease, since accelerated electrons are continuously supplied from the bow shock that is propagating into the ambient medium. As a result, overall emission from the shell overwhelms that from the lobe at wide range of frequencies from radio up to gamma-ray soon after the jet activity has ceased. Our result predicts a new class of dead radio sources that are dominated by shell emission. We suggest that the emission from the shell can be probed in particular at a radio wavelengths with the Square Kilometer Array (SKA) phase 1.

We explore non-thermal emission from a shocked interstellar medium, which is identified as an expanding shell, driven by a relativistic jet in active galactic nuclei (AGNs). In this work, we particularly focus on parsec-scale size mini-shells surrounding mini radio lobes. From radio to X-ray band, the mini radio lobe emission dominates the faint emission from the mini-shell. On the other hand, we find that inverse-Compton (IC) emission from the shell can overwhelm the associated lobe emission at very high energy (VHE; E>100 GeV) gamma-ray range, because energy densities of synchrotron photons from the lobe and/or soft photons from the AGN nucleus are large and IC scattering effectively works. The predicted IC emission from nearby mini-shells can be detected with the Cherencov Telescope Array (CTA) and they are potentially a new class of VHE gamma-ray emitters.

We report the finding of kiloparsec (kpc)-scale radio structures in three radio-loud narrow-line Seyfert 1 (NLS1) galaxies from the Faint Images of the Radio Sky at Twenty-centimeters (FIRST) of the Very Large Array (VLA), which increases the number of known radio-loud NLS1s with kpc-scale structures to six, including two gamma-ray emitting NLS1s (PMN J0948+0022 and 1H 0323+342) detected by the Fermi Gamma-ray Space Telescope. The detection rate of extended radio emissions in NLS1s is lower than that in broad-line active galactic nuclei (AGNs) with a statistical significance. We found both core-dominated (blazar-like) and lobe-dominated (radio-galaxy-like) radio structures in these six NLS1s, which can be understood in the framework of the unified scheme of radio-loud AGNs that considers radio galaxies as non-beamed parent populations of blazars. Five of the six NLS1s have (i) extended radio luminosities suggesting jet kinetic powers of >~10^44 erg/s, which is sufficient to make jets escape from hosts' dense environments, (ii) black holes of >~10^7 solar mass, which can generate the necessary jet powers from near-Eddington mass accretion, and (iii) two-sided radio structures at kpc scales, requiring expansion rates of ~0.01c--0.3c and kinematic ages of >~10^7 years. On the other hand, most typical NLS1s would be driven by black holes of <~10^7 solar mass in a limited lifetime of ~10^7 years. Hence the kpc-scale radio structures may originate in a small window of opportunity during the final stage of the NLS1 phase just before growing into broad-line AGNs.

Young radio galaxies (YRGs) provide an ideal laboratory to explore the connection between accretion disk and radio jet thanks to their recent jet formation. We investigate the relationship between the emission-line properties, the black hole accretion rate, and the radio properties using a sample of 34 low-redshift (z < 0.4) YRGs. We classify YRGs as high-excitation galaxies (HEGs) and low-excitation galaxies (LEGs) based on the flux ratio of high-ionization to low-ionization emission lines. Using the H{\alpha} luminosities as a proxy of accretion rate, we find that HEGs in YRGs have \sim1 dex higher Eddington ratios than LEGs in YRGs, suggesting that HEGs have higher mass accretion rate or higher radiative efficiency than LEGs. In agreement with previous studies, we find that the luminosities of emission lines, in particular H{\alpha}, are correlated with radio core luminosity, suggesting that accretion and young radio activities are fundamentally connected.

We examine plasma composition of jets in active galactic nuclei through the comparison of the total pressure ($P$) with partial pressures of electrons and protons in a cocoon. The total pressure is estimated from the analysis of an expanding cocoon dynamics. We determine the average kinetic energy per particle for several representative cases of particle energy distribution such as one- and two-temperature thermal plasmas and non-thermal electrons by evaluating the dissipation of total kinetic energy of the jet into the internal energy of cocoon plasma. The number density of the total electrons/positrons ($n_{\pm}$) in the cocoon is constrained by using the particle supply from hot spots and the absence of thermal bremsstrahlung emission from radio lobes. By inserting $P$, $n_{\pm}$ and the particle energy of each population into the equation of state, the number density ($n_{p}$) and pressure ($P_{p}$) of protons in the cocoon can be constrained. Applying this method to Cygnus A, we find that (i) electron/positron ($e^{\pm}$) pairs always dominate in terms of number density, but that (ii) either an "$e^{\pm}$-supported cocoon (i.e., $P_{\pm} >P_{p}$)" or "proton-supported one (i.e, $P_{\pm} <P_{p}$)" is possible.

We propose a new method to explore the candidate super-Eddington active galactic nuclei (AGNs). We examine the properties of infrared (IR) emission from the inner edge of the dusty torus in AGNs, which are powered by super- or sub-Eddington accretion flows around black holes, by considering the dependence of the polar angle on the radiation flux of accretion flows (Watarai et al. 2005). We find that for super-Eddington AGNs, of which the mass accretion rate is more than 10^2 times larger than the Eddington rate, the ratio of the AGN IR luminosity and the disc bolometric luminosity is less than 10^-2, unless the half opening angle of the torus (theta_torus) is small (theta_torus <65 degree). This is due to the self-occultation effect, whereby the self-absorption at the outer region of the super-Eddington flow dilutes the illumination of the torus. Such a small luminosity ratio is not observed in sub-Eddington AGNs, whose mass accretion rate is comparable to or no more than 10 times larger than the Eddington mass accretion rate, except for extremely thin tori (theta_torus >85 degree). We also consider the properties of the near-IR (NIR) emission radiated from hot dust >1000 K. We find that super-Eddington AGNs have a ratio of the NIR luminosity to the bolometric luminosity, L_NIR,AGN/L_bol,disc, at least one order of magnitude smaller than for sub-Eddington AGNs for a wide range of half opening angle (theta_torus > 65 degree), for various types of dusty torus model. Thus, a relatively low L_NIR,AGN/L_bol,disc is a property that allows identification of candidate super-Eddington AGNs. Lastly, we discuss the possibility that NIR-faint quasars at redshift z=6 discovered by a recent deep SDSS survey may be young quasars whose black holes grow via super-Eddington accretion.

We explore the evolution of the emissions by accelerated electrons in shocked shells driven by jets in active galactic nuclei (AGNs). Focusing on powerful sources which host luminous quasars, we evaluated the broadband emission spectra by properly taking into account adiabatic and radiative cooling effects on the electron distribution. The synchrotron radiation and inverse Compton (IC) scattering of various photons that are mainly produced in the accretion disc and dusty torus are considered as radiation processes. We show that the resultant radiation is dominated by the IC emission for compact sources (< 10kpc), whereas the synchrotron radiation is more important for larger sources. We also compare the shell emissions with those expected from the lobe under the assumption that a fractions of the energy deposited in the shell and lobe carried by the non-thermal electrons are $\epsilon_e \sim 0.01$ and $\epsilon_{e, lobe} \sim 1$, respectively. Then, we find that the shell emissions are brighter than the lobe ones at infra-red and optical bands when the source size is > 10kpc, and the IC emissions from the shell at > 10 GeV can be observed with the absence of contamination from the lobe irrespective of the source size. In particular, it is predicted that, for most powerful nearby sources ($L_j \sim 10^{47} ergs s^{-1}$), TeV gamma-rays produced via the IC emissions can be detected by the modern Cherenkov telescopes such as MAGIC, HESS and VERITAS.

We present the results of infrared L-band (3-4 micron) slit spectroscopy of 30 PG QSOs at z < 0.17, the representative sample of local high-luminosity, optically selected AGNs. The 3.3 micron polycyclic aromatic hydrocarbon (PAH) emission feature is used to probe nuclear (< a few kpc) starburst activity and to investigate the connections between AGNs and nuclear starbursts in PG QSOs. The 3.3 micron PAH emission is detected in the individual spectra of 5/30 of the observed PG QSOs. We construct a composite spectrum of PAH-undetected PG QSOs and discern the presence of the 3.3 micron PAH emission therein. We estimate the nuclear-starburst and AGN luminosities from the observed 3.3 micron PAH emission and 3.35 micron continuum luminosities, respectively, and find that the nuclear-starburst-to-AGN luminosity ratios in PG QSOs are similar to those of previously studied AGN populations with lower luminosities, suggesting that AGN-nuclear starburst connections are valid over the wide luminosity range of AGNs in the local universe. The observed nuclear-starburst-to-AGN luminosity ratios in PG QSOs with available supermassive black hole masses are comparable to a theoretical prediction based on the assumption that the growth of a supermassive black hole is controlled by starburst-induced turbulence.

We explore the emissions by accelerated electrons in shocked shells driven by jets in active galactic nuclei (AGNs). Focusing on powerful sources which host luminous quasars, the synchrotron radiation and inverse Compton (IC) scattering of various photons that are mainly produced in the core are considered as radiation processes. We show that the radiative output is dominated by the IC emission for compact sources (< 30kpc), whereas the synchrotron radiation is more important for larger sources. It is predicted that, for powerful sources ($L_j \sim 10^{47} ergs s^{-1}$), GeV-TeV gamma-rays produced via the IC emissions can be detected by the Fermi satellite and modern Cherenkov telescopes such as MAGIC, HESS and VERITAS if the source is compact.

Based on the physical model of a supermassive black hole (SMBH) growth via gas accretion in a circumnuclear disk (CND) proposed by Kawakatu & Wada (2008), we describe the formation of high-$z$ ($z > 6$) quasars (QSOs) whose BH masses are M_{BH}> 10^{9} M_{\odot}. We derive the necessary conditions to form QSOs at z > 6 by only gas accretion: (i) A large mass supply with M_{sup} > 10^{10}M_{\odot} from host galaxies to CNDs, because the final BH mass is only 1-10% of the total supplied mass from QSO hosts. (ii) High star formation efficiency for a rapid BH growth. We also find that if the BH growth is limited by the Eddington accretion, the final BH mass is greatly suppressed. Thus, the super-Eddington growth is required for the QSO formation. The evolution of the QSO luminosity depends on the redshift z_{i} at which accretion onto a seed BH is initiated. In other words, the brighter QSOs at z >6 favor the late growth of SMBHs (i.e., z_{i}=10) rather than early growth (i.e., z_{i}=30). Moreover, we predict the observable properties and the evolution of QSOs at z >6. In a QSO phase, there should exist a stellar rich massive CND, whose gas mass is about 10% of the dynamical mass inside 0.1-1 kpc}. On the other hand, in a phase where the BH grows (i.e., a proto-QSO phase), the proto-QSO has a gas rich massive CNDs whose gas mass is comparable to the dynamical mass (abridged).

We investigated the relationship between supermassive black hole (SMBH) mass and host starburst luminosity in Seyfert galaxies and Palomar-Green QSOs, focusing on the host galaxy morphology. Host starburst luminosity was derived from the 11.3 micron polycyclic aromatic hydrocarbon luminosity. We found that the SMBH masses of elliptical-dominated host galaxies are more massive than those of disk-dominated host galaxies statistically. We also found that the SMBH masses of disk-dominated host galaxies seem to be suppressed even under increasing starburst luminosity. These findings imply that final SMBH mass is strongly regulated by host galaxy morphology. This can be understood by considering the radiation drag model as the SMBH growth mechanism, taking into account the radiation efficiency of the host galaxy.

We argue that the origin of "FRI/FRI{-.1em}I dichotomy" -- the division between Fanaroff-Riley class I (FRI) with subsonic lobes and class I{-.1em}I (FRI{-.1em}I) radio sources with supersonic lobes is sharp in the radio-optical luminosity plane (Owen-White diagram) -- can be explained by the deceleration of advancing radio lobes. The deceleration is caused by the growth of the effective cross-sectional area of radio lobes. We derive the condition in which an initially supersonic lobe turns into a subsonic lobe, combining the ram-pressure equilibrium between the hot spots and the ambient medium with the relation between "the hot spot radius" and "the linear size of radio sources" obtained from the radio observations. We find that the dividing line between the supersonic lobes and subsonic ones is determined by the ratio of the jet power $L_{\rm j}$ to the number density of the ambient matter at the core radius of the host galaxy $\bar{n}_{\rm a}$. It is also found that there exists the maximal ratio of $(L_{\rm j}/\bar{n}_{\rm a})$ and its value resides in $(L_{\rm j}/\bar{n}_{\rm a})_{\rm max}\approx 10^{44-47} {\rm erg} {\rm s}^{-1} {\rm cm}^{3}$, taking account of considerable uncertainties. This suggests that the maximal value $(L_{\rm j}/\bar{n}_{\rm a})_{\rm max}$ separates between FRIs and FRI{-.1em}Is.

We present a new prediction of GeV $\gamma$-ray emission from radio lobes of young AGN jets. In the previous work of Kino et al. (2007), MeV $\gamma$-ray bremsstrahlung emission was predicted from young cocoons/radio-lobes in the regime of no coolings. In this study, we include cooling effects of bremsstrahlung emission and adiabatic loss. With the initial conditions determined by observed young radio lobes, we solve a set of equations describing the expanding lobe evolution. Then we find that the lobes initially have electron temperature of $\sim$GeV, and they cool down to $\sim$MeV by the adiabatic loss. Correspondingly, the lobes initially yield bright bremsstrahlung luminosity in $\sim$GeV range and they fade out. We estimate these $\gamma$-ray emissions and show that nearby young radio lobes could be detected with Fermi Gamma-ray Space Telescope.

Gamma-ray emission from cocoons of young radio galaxies is predicted. Considering the process of adiabatic injection of the shock dissipation energy and mass of the relativistic jet into the cocoon, we find that the thermal electron temperature of the cocoon is typically predicted to be of the order of $\sim$ MeV, and is determined only by the bulk Lorentz factor of the jet. Together with the time-dependent dynamics of the cocoon expansion, we find that young cocoons can yield thermal Bremsstrahlung emissions at energies $\sim$MeV. Hotter cocoons (i.e., GeV) for younger sources are also discussed.

In order to reveal the long-term evolution of relativistic jets in active galactic nuclei (AGNs), we examine the dynamical evolution of variously-sized radio galaxies [i.e., compact symmetric objects (CSOs), medium-size symmetric objects (MSOs), Fanaroff-Riley type II radio galaxies (FRIIs)]. By comparing the observed relation between the hot spot size and the linear size of radio source with a coevolution model of hot spot and cocoon, we find that the advance speed of hot spots and lobes inevitably show the deceleration phase (CSO-MSO phase) and the acceleration phase (MSO-FRII phase). The deceleration is caused by the growth of the cross-sectional area of the cocoon head. Moreover, by comparing the hot spot speed with the sound speed of the ambient medium, we predict that only CSOs whose initial advance speed is higher than 0.3-0.5c can evolve into FRIIs.

We investigate the optical emission-line flux ratios of narrow-line regions, in order to determine whether the formation of AGN jets requires specific accretion conditions. We find that bright compact radio galaxies, which are powerful radio galaxies in the early stage of the jet activity, exhibit systematically larger flux ratios of [O{\sc i}]$\lambda$6300/[O{\sc iii}]$\lambda$5007 and smaller flux ratios of [O{\sc iii}]$\lambda$5007/[O{\sc iii}]$\lambda$4363 than radio-quiet (RQ) Seyfert 2 galaxies. Comparing the observed line ratios with photoionization models, it is found that the difference in the flux ratio of low- to high-ionization lines (e.g., [O{\sc i}]$\lambda$6300/[O{\sc iii}]$\lambda$5007) can be well understood by the difference in the spectral energy distribution (SED) of ionizing sources. Powerful young radio-loud (YRL) AGNs favor SED without a strong big blue bump, i.e., a radiatively inefficient accretion flow (RIAF), while RQ AGNs are consistent with the models adopting SED with a strong big blue bump, i.e., a geometrically thin, optically thick disk. These findings imply that the formation of powerful AGN jets requires the accretion disk with harder ionizing SED (i.e., a RIAF). We discuss the obscuring structure of YRL AGNs as a plausible origin of the difference in flux ratios of [O{\sc iii}]$ \lambda$5007/[O{\sc iii}]$\lambda$4363.

We investigate the total kinetic powers (L_j) and ages (t_age) of powerful jets in FR II radio galaxies by comparison of the dynamical model of expanding cocoons with observations. We select four FR II radio sources (Cygnus A, 3C 223, 3C 284, and 3C 219), for which the mass-density profiles of intracluster medium (ICM) are known in the literature. It is found that large fractions > 0.02 - 0.7 of the Eddington luminosity (L_Edd) are carried away as a kinetic power of jet. The upper limit of estimated 2 L_j / L_Edd are larger than unity (< 10) for some sources, suggesting a possibility of super-Eddington mass accretions. As a consequence of the large powers, we also find that the total energy stored in the cocoon (E_c) exceeds the energy derived from the minimum energy condition for the energy of radiating non-thermal electrons and magnetic fields (E_min): 4< E_c /E_min <310. This implies that most of the energy in cocoon is carried by invisible components such as thermal leptons (electron and positron) and/or protons.

We examine the evolution of variously-sized radio galaxies [i.e., compact symmetric objects (CSOs), medium-size symmetric objects (MSOs), Fanaroff-Riley type II radio galaxies (FRIIs)], by comparing the relation between the hot spot size and the projected linear size with a coevolution model of hot spots and a cocoon. We take account of the deceleration effect by the cocoon head growth. We find that the advance speed of hot spots and lobes inevitably show the deceleration phase (CSO-MSO phase) and the acceleration phase (MSO-FRII phase). This is ascribed to the change of the power-law index of ambient density profile in the MSO phase ($\sim $1 kpc). It is also found that the cocoon shape becomes nearly spherical or disrupted for MSOs, while an elongated morphology is predicted for CSOs and FRIIs. This seems to be consistent with the higher fraction of distorted morphology of MSOs than that of CSOs and FRI. Finally, we predict that only CSOs whose initial advance speed is higher than about 0.1c can evolve into FRIIs, comparing the hot spot speed with the sound speed of the ambient medium.

Motivated by the recent observational and theoretical evidence that long Gamma-Ray Bursts (GRBs) are likely associated with low metallicity, rapidly rotating massive stars, we examine the cosmological star formation rate (SFR) below a critical metallicity Z_crit Z_sun/10 - Z_sun/5, to estimate the event rate of high-redshift long GRB progenitors. To this purpose, we exploit a galaxy formation scenario already successfully tested on a wealth of observational data on (proto)spheroids, Lyman break galaxies, Lyman alpha emitters, submm galaxies, quasars, and local early-type galaxies. We find that the predicted rate of long GRBs amounts to about 300 events/yr/sr, of which about 30 per cent occur at z>~6. Correspondingly, the GRB number counts well agree with the bright SWIFT data, without the need for an intrinsic luminosity evolution. Moreover, the above framework enables us to predict properties of the GRB host galaxies. Most GRBs are associated with low mass galaxy halos M_H<~10^11 M_sun, and effectively trace the formation of small galaxies in such halos. The hosts are young, with age smaller than 5*10^7 yr, gas rich, but poorly extincted (A_V<~0.1) because of their chemical immaturity; this also implies high specific SFR and quite extreme alpha-enhancement. Only the minority of hosts residing in large halos with M_H>~10^12 M_sun have larger extinction (A_V~0.7-1), SFRs exceeding 100 M_sun/yr and can be detected at submm wavelengths. Most of the hosts have UV magnitudes in the range -20 <~M_1350<~ -16, and Lyman alpha luminosity in the range 2*10^40 <~L_Lya<~2*10^42 erg/s. GRB hosts are thus tracing the faint end of the luminosity function of Lyman break galaxies and Lyman alpha emitters.

We propose a new evolutionary model of a supermassive black hole (SMBH) and a circumnuclear disk (CND), taking into account the mass-supply from a host galaxy and the physical states of CND. In the model, two distinct accretion modes depending on gravitational stability of the CND play a key role on accreting gas to a SMBH. (i) If the CMD is gravitationally unstable, energy feedback from supernovae (SNe) supports a geometrically thick, turbulent gas disk. The accretion in this mode is dominated by turbulent viscosity, and it is significantly larger than that in the mode (ii), i.e., the CMD is supported by gas pressure. Once the gas supply from the host is stopped, the high accretion phase ($\sim 0.01- 0.1 M_{\odot} {\rm yr}^{-1}$) changes to the low one (mode (ii), $\sim 10^{-4} M_{\odot} {\rm yr}^{-1}$), but there is a delay with $\sim 10^{8}$ yr. Through this evolution, the gas-rich CND turns into the gas poor stellar disk. We found that not all the gas supplied from the host galaxy accrete onto the SMBH even in the high accretion phase (mode (i)), because the part of gas is used to form stars. As a result, the final SMBH mass ($M_{\rm BH,final}$) is not proportional to the total gas mass supplied from the host galaxy ($M_{\rm sup}$); $M_{\rm BH,final}/M_{\rm sup}$ decreases with $M_{\rm sup}$.This would indicate that it is difficult to form a SMBH with $\sim 10^{9} M_{\odot}$ observed at high-$z$ QSOs. The evolution of the SMBH and CND would be related to the evolutionary tracks of different type of AGNs.

We investigated the correlation between nuclear/circumnuclear starbursts around the active galactic nuclei (AGNs) and the AGN activities for 43 Seyfert galaxies in the CfA and 12 micron samples. We found that circumnuclear starburst luminosity as well as nuclear starburst luminosity are positively correlated with AGN luminosity. Moreover, nuclear starburst luminosity is more strongly correlated with the AGN luminosity normalized with AGN Eddington luminosity than is circumnuclear starburst luminosity. This implies that starbursts nearer the AGN could have a greater effect on AGN mass accretion. We also discuss these results from the viewpoint of the radiation effects from starbursts and sequential starbursts.

We investigate the total kinetic powers (L_{j}) and ages (t_{age}) of powerful jets of four FR II radio sources (Cygnus A, 3C 223, 3C 284, and 3C 219) by the detail comparison of the dynamical model of expanding cocoons with observed ones. It is found that these sources have quite large kinetic powers with the ratio of L_{j} to the Eddington luminosity (L_{Edd}) resides in $0.02 <L_{j}/L_{Edd} <10$. Reflecting the large kinetic powers, we also find that the total energy stored in the cocoon (E_{c}) exceed the energy derived from the minimum energy condition (E_{min}): $2< E_{c}/E_{min} <160$. This implies that a large amount of kinetic power is carried by invisible components such as thermal leptons (electron and positron) and/or protons.

Massive tori with $\approx 10^{8-9}M_{\odot}$ are predicted to extend on $\sim $100 pc scale around the centre of elliptical galaxy progenitors by a model of a supermassive black hole (SMBH) growth coeval to the spheroidal population of the host galaxy. Direct detection of such massive tori would cast light on a key physical condition that allows the rapid growth of SMBHs and the appearance of QSOs at high redshift. For this reason, we examine the detectability of such structures at substantial redshift with the Atacama Large Millimeter Array (ALMA). We propose that submillimeter galaxies (SMGs) are the best targets to test our predictions. In order to assess the observational feasibility, we estimate the expected number counts of SMGs with massive tori and check the detectability with the ALMA instrument, the unique facility which can resolve the central region of high redshift objects. Our work shows that ALMA will be able to resolve and detect high-$J$ ($J >$ 4) CO emissions from $\sim$100 pc scale extended massive tori up to $z\approx2$. Observations of lensed SMGs will yield excellent spatial resolution, allowing even to resolve their massive tori at higher redshift. We discuss further the detectability of the HCN molecule, as a better tracer of the high density gas expected in such tori. The final goal of these kind of observations is to pinpoint possible physical mechanisms that storage in the very central galactic regions very large amount gas on timescale of several 10$^{8}$ yr.

We discovered a significant anti-correlation between the mass of a supermassive black hole (SMBH), $M_{\rm BH}$, and the luminosity ratio of infrared to active galactic nuclei (AGN) Eddington luminosity, $L_{\rm IR}/L_{\rm Edd}$, over four orders of magnitude for ultraluminous infrared galaxies with type I Seyfert nuclei (type I ULIRGs) and nearby QSOs. This anti-correlation ($M_{\rm BH}$ vs. $L_{\rm IR}/L_{\rm Edd}$) can be interpreted as the anti-correlation between the mass of a SMBH and the rate of mass accretion onto a SMBH normalized by the AGN Eddington rate, $\dot{M}_{\rm BH}/\dot{M}_{\rm Edd}$. In other words, the mass accretion rate $\dot{M}_{\rm BH}$ is not proportional to that of the central BH mass. Thus, this anti-correlation indicates that BH growth is determined by the external mass supply process, and not the AGN Eddington-limited mechanism. Moreover, we found an interesting tendency for type I ULIRGs to favor a super-Eddington accretion flow, whereas QSOs tended to show a sub-Eddington flow. On the basis of our findings, we suggest that a central SMBH grows by changing its mass accretion rate from super-Eddington to sub-Eddington. According to a coevolution scenario of ULIRGs and QSOs based on the radiation drag process, it has been predicted that a self-gravitating massive torus, whose mass is larger than a central BH, exists in the early phase of BH growth (type I ULIRG phase) but not in the final phase of BH growth (QSO phase). At the same time, if one considers the mass accretion rate onto a central SMBH via a turbulent viscosity, the anti-correlation ($M_{\rm BH}$ vs. $L_{\rm IR}/L_{\rm Edd}$) is well explained by the positive correlation between the mass accretion rate $\dot{M}_{\rm BH}$ and the mass ratio of a massive torus to a SMBH.

Strong $\gamma$-ray emission from cocoons of young radio galaxies is newly predicted. Considering the process of adiabatic injection of the shock dissipation energy and mass of the relativistic jet in active nuclei (AGNs) into the cocoon, while assuming thermalizing electron plasma interactions, we find that the thermal electron temperature of the cocoon is typically predicted in $\sim$MeV, which is determined only by the bulk Lorentz factor of the relativistic jet. Together with the time-dependent dynamics of the cocoon expansion, we find that young cocoons can yield thermal bremsstrahlung emissions at energies $\sim$ MeV.

We describe the dynamical evolution of hot spots velocity, pressure and mass density in radio loud active galactic nuclei (AGNs), taking proper account of (1) the conservations of the mass, momentum, and kinetic energy flux of the unshocked jet, (2) the deceleration process of the jet by shocks, and (3) the cocoon expansion without assuming the constant aspect ratio of the cocoon. By the detailed comparison with two dimensional relativistic hydrodynamic simulations, we show that our model well reproduces the whole evolution of relativistic jets. Our model can explain also the observational trends of the velocity, the pressure, the size, and mass density of hot spots in compact symmetric objects (CSOs) and FR II radio galaxies.

We examine the constraints imposed on the total kinetic power ($L_{\rm j}$) and the age ($t_{\rm age}$) of relativistic jets in FR II radio sources in a new way. We solve the dynamical expansion of its cocoon embedded in the intra-cluster medium (ICM) and obtain the analytic solution of its physical quantities in terms of $L_{\rm j}$ and $t_{\rm age}$. The estimate of $L_{\rm j}$ and $t_{\rm age}$ is done by the comparison of the model and the observed shape of the cocoon. The analysis is focused on Cygnus A and we find that (i) the source age is $3 {\rm Myr} <t_{\rm age}<30 {\rm Myr}$ and (ii) the total kinetic power of the jet is estimated as $0.2\times 10^{46}{\rm erg s^{-1}} <L_{\rm j}<1\times 10^{48} {\rm erg s^{-1}}$ which is larger than 1% of the Eddington luminosity of Cygnus A.

We examine whether the ultraluminous infrared galaxies that contain a type I Seyfert nucleus (a type I ULIRG) are in the transition stage from ULIRGs to quasi-stellar objects (QSOs). To inspect this issue, we compare the black hole (BH) mass, the bulge luminosity and the far infrared luminosity among type I ULIRGs, QSOs and elliptical galaxies. As a result, we find the following results; (1) The type I ULIRGs have systematically smaller BH masses in spite of the comparable bulge luminosity relative to QSOs and elliptical galaxies. (2) The far-infrared luminosity of most type I ULIRGs is larger than the Eddington luminosity. We show that above results do not change significantly for 3 type I ULIRGs that we can estimate the visual extinction from the column density. Also, for all 8 type I ULIRGs, we investigate the effect of uncertainties of BH mass measurments and our sample bias, so that it turns out that our results do not alter even if we consider above two effects. In addition, Anabuki (2004) revealed that their X-ray properties are similar to those of the narrow line Seyfert 1 galaxies. These would indicate that active galactic nuclei (AGNs) with a high mass accretion rate exist in the type I ULIRGs. Based on all of these findings, we conclude that it would be a natural interpretation that type I ULIRGs are the early phase of BH growth, namely the missing link between ULIRGs and QSOs. Moreover, by comparing our results with a theoretical model of a coevolution scenario of a QSO BH and a galactic bulge, we show clearly that this explanation would be valid.

The matter content of jets in active galactic nuclei is examined in a new way. We model the dynamical expansion of its cocoon embedded in the intra-cluster medium (ICM). By comparing the observed shape of the cocoon with that expected from the theoretical model, we estimate the total pressure ($P_{\rm c}$) and electron temperature ($T_{e}$) of the cocoon. The number density of the total electrons ($n_{e^{-}}$) is constrained by using the non-thermal spectrum of the hot spot and the analysis of the momentum balance between the jet thrust and the rum pressure of ICM. Together with the obtained $P_{\rm c}$, $T_{e}$ and $n_{e^{-}}$, we constrain the matter content in the jets. We find that, in the case of Cygnus A, the ratio of number density of protons to that of electrons is of order of $10^{-3}$ . This implies the existence of a large number of positron in the jet.

he formation and growth of supermassive black holes (SMBHs) physically linked with bulges are considered. We focus on the radiation hydrodynamic process for the growth of SMBH in the optically thick starburst phase, where radiation from bulge stars drives the mass accretion on to a galactic center through radiation drag effect. In the present scenario, the AGN luminosity-dominant phase (QSO phase) is preceded by the host luminosity-dominat phase, which is called ``proto-QSO phase''. In this phase, there exists the massive dusty disks within younger bulges. Also, the proto-QSO phase is anticipated by an optically-thick ultraluminous infrared galaxy (ULIRG) phase. Furthermore, such radiation hydrodynamic model has been also applied to disk galaxies. It turns out that the mass of a SMBH primarily correlates with a bulge component even in a disk galaxy. Thus, by analogy to proto-QSOs, the BH growing phase in disk galaxies may have massive dusty disks within younger bulges.

Combining a theoretical model of mass accretion onto a galactic center with a high-resolution $N$-body/SPH simulation, we investigate the formation of an intermediate massive black hole (IMBH) during the hierarchical formation of a small spiral galaxy (with a total mass of $10^{10}M_{\odot}$) in the high-$z$ universe. We found that the rate of average mass accretion to the nucleus due to the radiation drag exerted by newly formed stars in the forming galaxy is $\approx 10^{-5}M_{\odot}$yr$^{-1}$. As a result of this accretion, an IMBH with $\approx 10^{4}M_{\odot}$ can be formed in the center of the spiral galaxy at $z\sim 4$. We found that a central BH coevolves with the dark matter halo from $z\sim 15$ to $z\sim 2$. The mass ratio of the BH to the dark matter halo is nearly constant $\approx (1-3) \times 10^{-6}$ from $z\sim 10$ to $z\sim 2$. This is because that change in the dark matter potential enhances star formation in the central part of the galaxy, and as a result the BH evolves due to mass accretion via the radiation drag. Therefore, our model naturally predicts a correlation between massive BHs and dark matter halos. Moreover, it is found that the final BH-to-bulge mass ratio ($\approx 5\times 10^{-5}$) in a small spiral galaxy at high-$z$ is much smaller than that in the large galaxies ($\approx 10^{-3}$). Our results also suggest that the scatter in the observed scaling relations between the bulge mass and black hole mass are caused by a time lag between BH growth and growth of bulge. We also predict that the X-ray luminosity of AGN is positively correlated with the CO luminosity in the central region. By comparing our results with the properties of Lyman break galaxies (LBGs), it is predicted that some LBGs have massive BHs of $\approx 10^{6}-10^{7}M_{\odot}$.

The hydrodynamic formation of massive black holes (BHs) in globular clusters is considered. In particular, we examine the possibility of the BH formation induced by the radiation drag that is exerted on the interstellar matter by stellar radiation in globular clusters. The radiation drag extracts angular momentum from interstellar gas and thus allows the gas to accrete on to the cluster center.By incorporating the realistic chemical evolution of globular clusters, we scrutinize the efficiency of the radiation drag to assess the total accreted mass. As a result, we find that if a globular cluster is more massive than $\approx 6\times 10^{6}M_{\odot}$ in the present-day stellar component, a massive BH with $>260M_{\odot}$ can form within it. But, the BH-to-cluster mass ratio s considerably smaller than the BH-to-bulge mass ratio ($\approx 10^{-3}$) found in galactic bulges. The results are not sensitive to the assumed stellar initial mass function and star formation rate in the cluster, as long as the resultant color-magnitude relation, metallicity, and mass-to-luminosity ratio satisfy those observed in globular clusters. Hence, the putative linear relation between BH mass and bulge mass($M_{\rm BH}-M_{\rm bulge}$ relation) cannot be extrapolated to globular cluster systems. In the present regime, we discuss the BH formation in M15, G1, $\omega$ Cen, the M33 nucleus, and the compact X-ray sources in M82. Finally, we argue observational indications of the formation process of massive BHs in globular clusters. We find that the final phase of BH growth due to the radiation drag can be observed as ultraluminous X-ray sources (ULXs) with $\sim 10^{41}{\rm erg} {\rm s}^{-1}$.

A potential mechanism is proposed to account for the fact that supermassive black holes (SMBHs) in disk galaxies appear to be smaller than those in elliptical galaxies in the same luminosity range. We consider the formation of SMBHs by radiation drag (Poynting-Robertson effect), which extracts angular momentum from interstellar matter and thereby drives the mass accretion onto a galactic center. In particular, we quantitatively scrutinize the efficiency of radiation drag in a galaxy composed of bulge and disk, and elucidate the relation between the final mass of SMBH and the bulge-to-disk ratio of the galaxy. As a result, it is found that the BH-to-galaxy mass ratio, $M_{\rm BH}/M_{\rm galaxy}$, decreases with a smaller bulge-to-disk ratio, due to the effects of geometrical dilution and opacity, and is reduced maximally by two orders of magnitude, resulting in $M_{\rm BH}/M_{\rm galaxy}\approx 10^{-5}$. In contrast, if only the bulge components in galaxies are focused, the BH-to-bulge mass ratio becomes $M_{\rm BH}/M_{\rm bulge}\approx 10^{-3}$, which is similar to that found in elliptical galaxies. Thus, it turns out that the mass of SMBH primarily correlates with a bulge, not with a disk, consistently with observational data.

Based on the radiation hydrodynamical model for the black hole(BH) growth,incorporated with the chemical evolution of the early-type host galaxy, we construct the coevolution model of a QSO BH and the host galaxy. As a result, it is found that after a galactic wind epoch,the luminosity is shifted from the host-dominant phase to the AGN-dominant phase (QSO phase) in the timescale of a few $10^{8}$ years.The former phase corresponds to the early stage of growing BH, and can be regarded as a ``proto-QSO'' phase. It has observable characteristic properties (detail inthis paper).By comparing these predictions with recent observations, radio galaxies are a possible candidate for proto-QSOs.Also, it is anticipated that the proto-QSO phase is preceded by an optically thick phase, which may correspond to ULIRGs.In this phase, $M_{\rm BH}/M_{\rm bulge}$ is predicted to be much less than $10^{-3}$ and grow with metallicity.Moreover, as precursors of ULIRGs, optically-thin star-forming galaxies are predicted. These may be in the assembly phase of Lyman break galaxies (LBGs) or Ly$\alpha$ emitters.

We quantitatively scrutinize the effects of the radiation drag arising from the radiation fields in a galactic bulge in order to examine the possibility that the radiation drag could be an effective mechanism to extract angular momentum in a spheroidal system like a bulge and allow plenty of gas to accrete onto the galactic center. For this purpose, we numerically solve the relativistic radiation hydrodynamical equation coupled with the accurate radiative transfer and quantitatively assess the radiation drag efficiency. As a result, we find that in an optically thick regime the radiation drag efficiency is sensitively dependent on the density distributions of interstellar medium (ISM). The efficiency drops according to $\tau_{\rm T}^{-2}$ in an optically thick {\it uniform} ISM, where $\tau_{\rm T}$ is the total optical depth of the dusty ISM, whereas the efficiency remains almost constant at a high level if the ISM is {\it clumpy}. Hence, if the bulge formation begins with a star formation event in a clumpy ISM, the radiation drag will effectively work to remove the angular momentum and the accreted gas may form a supermassive black hole. As a natural consequence, this mechanism reproduces a putative linear relation between the mass of a supermassive black hole and the mass of a galactic bulge, although further detailed modeling for stellar evolution is required for the more precise prediction.