We show the best current simulations of the absorption and emission features predicted from thermal-radiative winds produced from X-ray illumination of the outer accretion disc in binary systems. We use the density and velocity structure derived from a radiation hydrodynamic code as input to a Monte-Carlo radiation transport calculation. The initial conditions are matched to those of the black hole binary system H1743-322 in its soft, disc dominated state, where wind features are seen in Chandra grating data. Our simulation fits well to the observed line profile, showing that these physical wind models can be the origin of the absorption features seen, rather than requiring a magnetically driven wind. We show how the velocity structure is the key observable discriminator between magnetic and thermal winds. Magnetic winds are faster at smaller radii, whereas thermal winds transition to a static atmosphere at smaller radii. New data from XRISM (due for launch Jan 2022) will give an unprecedented view of the physics of the wind launch and acceleration processes, but the existence of static atmospheres in small disc systems already rules out magnetic winds which assume self-similar magnetic fields from the entire disc as the origin of the absorption features seen.

Based on recent X-ray observations, ultra-fast outflows from supermassive black holes are expected to have enough energy to dramatically affect their host galaxy but their launch and acceleration mechanisms are not well understood. We perform two-dimensional radiation hydrodynamics simulations of UV line-driven disc winds in order to calculate the mass loss rates and kinetic power in these models. We develop a new iterative technique which reduces the mass accretion rate through the inner disc in response to the wind mass loss. This makes the inner disc is less UV bright, reducing the wind power compared to previous simulations which assumed a constant accretion rate with radius. The line-driven winds in our simulations are still extremely powerful, with around half the supplied mass accretion rate being ejected in the wind for black holes with mass $10^8-$$10^{10}\, M_\odot$ accreting at $L/L_{\rm Edd}=0.5-$$0.9$. Our results open up the way for estimating the growth rate of supermassive black hole and evaluating the kinetic energy ejected into the inter stellar medium (active galactic nuclei feedback) based on a physical model of line-driven disc winds.

We calculate the pulsed fraction (PF) of the super-critical column accretion flows onto magnetized neutron stars (NSs), of which the magnetic axis is misaligned with the rotation axis, based on the simulation results by Kawashima et al.(2016, PASJ, 68, 83). Here, we solve the geodesic equation for light in the Schwarzschild spacetime in order to take into account the light bending effect. The gravitational redshift and the relativistic doppler effect from gas motions of the accretion columns are also incorporated. The pulsed emission appears since the observed luminosity, which exceeds the Eddington luminosity for the stellar-mass black holes, periodically changes via precession of the column caused by the rotation of the NS. The PF tends to increase as $\theta_{\rm obs}$ approaching to $\theta_{\rm B}$, where $\theta_{\rm obs}$ and $\theta_{\rm B}$ are the observer's viewing angle and the polar angle of the magnetic axis measured from the rotation axis. The maximum PF is around 50 %. Also, we find that the PF becomes less than 5 % for $\theta_{\rm obs} \lesssim 5^\circ$ or for $\theta_{\rm B} \lesssim 5^\circ$. Our results are consistent with observations of ultraluminous X-ray pulsars (ULXPs) with few exceptions, since the ULXPs mostly exhibit the PF of $\lesssim$ 50 %. Our present study supports the hypothesis that the ULXPs are powered by the super-critical column accretion onto NSs.

UltraFast Outflows (UFO) are observed in some active galactic nuclei (AGN), with blueshifted and highly ionised Fe-K absorption features. AGN typically have a UV bright accretion flow, so UV line driving is an obvious candidate for launching these winds, but this mechanism requires material with UV opacity, in apparent conflict with the observed high ionisation state of the wind. In this paper we synthesise the X-ray energy spectra resulting from different lines of sight through a state of the art radiation hydrodynamics UV line driven disc wind simulation. We demonstrate that there are some lines of sight which only intercept highly ionised and fast outflowing material. The cooler material required for the UV line driving acceleration is out of the line of sight, close to the disc, shielded from the X-rays by a failed wind. We fit these simulated wind spectra to data from the archetypal UFO source PG 1211+143 and show that they broadly reproduce the depth and velocity of the iron absorption lines seen. This directly demonstrates that UV line driving is a viable mechanism to launch UFOs. We simulate microcalorimeter observations of this wind and show that their high energy resolution can resolve the detailed structure in the wind and recover the wind energetics when combined with models which correctly estimate the wind launch radius. New data from microcalorimeters will pave the way for physical predictions of AGN wind feedback in cosmological simulations.

The spectacular images of the M87 black hole taken by the Event Horizon Telescope (EHT) have opened a new era of black hole research. One of the next issues is to take polarization images around the central black hole (BH). Since radio emission is produced by synchrotron process, polarization properties should vividly reflect the magnetic field structure at the jet base and thus provide good information regarding the magnetic mechanism of jet formation. With this kept in mind we perform general relativistic (GR) radiative transfer calculations of polarized light based on the GR magnetohydrodynamic (MHD) simulation data of accretion flow and outflow in M87, to obtain their linear and circular polarization images in the horizon-scale. We found that the linear polarization components from the jet base and inner accretion flow should experience Faraday rotation when passing through magnetized plasmas around the BH, thus sensitively depending on the BH spin. Through the comparison with intensity image at 1.3mm by EHT and the polarized degree and the rotation measure (RM) at 1.3mm with the Submillimeter Array, the model with the spin parameter of a=0.9M_BH (with M_BH being the BH mass) is favored over other models with a=0.5M_BH or 0.99M_BH, though we need further systematic studies for confirmation. We also find in low-temperature models clear ring-like image in the circular polarization map, which arises because of Faraday conversion of the linearly polarized synchrotron emission and is thus indicative of magnetic field. This occurs only when the emission region is threaded with well-ordered magnetic fields and hence no clear images are expected in high-temperature disk models, in which disk emission is appreciable. We will be able to elucidate the field configuration through the comparison between the simulated polarization images and future polarimetry with EHT and other VLBI observations.

Blueshifted absorption lines are seen in high inclination black hole binary systems in their disc dominated states, showing these power an equatorial disc wind. While some contribution from magnetic winds remain a possibility, thermal and thermal-radiative winds are expected to be present. We show results from radiation hydrodynamic simulations which show that the additional radiation force from atomic features (bound-free and lines) are important along with electron scattering. Together, these increase the wind velocity at high inclinations, so that they quantitatively match the observations in H 1743-322, unlike purely thermal winds which are too slow. We highlight the role played by shadowing of the outer disc from the (sub grid) inner disc Compton heated layer, and show that the increase in shadow from the higher Compton temperature after the spectral transition to the hard state leads to strong suppression of the wind. Thermal-radiative winds explain all of the spectral features (and their disappearance) in this simplest wind system and magnetic winds play only a minor role. We speculate that thermal-radiative winds can explain all the spectral features seen in the more complex (larger disc size) binaries, GRO J1655-40 and GRS 1915+105, without requiring magnetic winds.

We investigate the properties of accretion flows onto a black hole (BH) with a mass of $M_{\rm BH}$ embedded in an initially uniform gas cloud with a density of $n_{\infty}$ in order to study rapid growth of BHs in the early Universe. In previous work, the conditions for super-Eddington accretion from outside the Bondi radius were studied by assuming that radiation produced at the vicinity of the central BH has a single-power-law spectrum $\nu^{-\alpha}$ at $h\nu \geq 13.6~{\rm eV}$ ($\alpha \sim 1.5$). However, radiation spectra depends on the BH mass and accretion rate. Here, we perform two-dimensional multi-frequency radiation hydrodynamical simulations taking into account more realistic radiation spectra associated with the properties of nuclear accretion disks. We find that the condition for a transition to super-Eddington accretion is alleviated for a wide range of masses ($10\lesssim M_{\rm BH}/M_{\odot} \lesssim 10^6$) because photoionization for accretion disk spectra are less efficient than those for single-power-law spectra. For disk spectra, the transition to super-Eddington is more likely to occur for lower BH masses because the radiation spectra become too hard to ionize the gas. Even when accretion flows are exposed to anisotropic radiation, the effect due to radiation spectra shrinks the ionized region and likely leads to the transition to a wholly neutral accretion phase. Finally, by generalizing our simulation results, we construct a new analytical criterion required for super-Eddington accretion; $(M_{\rm BH}/10^5~M_{\odot}) (n_{\infty}/10^4~{\rm cm^{-3}}) \gtrsim 2.4~ (\langle\epsilon\rangle /100~{\rm eV})^{-5/9}$, where $\langle\epsilon\rangle$ is the mean energy of ionizing radiation from the central BH.

To what extent can the one-dimensional slim disk model reproduce the multi-dimensional results of global radiation-hydrodynamic simulations of super-Eddington accretion? With this question in mind, we perform systematic simulation study of accretion flow onto a non-spinning black hole for a variety of black hole masses of $(10 - 10^7) M_{\odot}$ and mass accretion rates of $(1.4 \times 10^2 - 5.6 \times 10^3)L_{\rm Edd}/c^2$ (with $L_{\rm Edd}$ and $c$ being the Eddington luminosity and speed of light). In order to adequately resolve large-scale outflow structure, we extensively expand a simulation box to cover the space of $3000 r_{\rm S}$ (with $r_{\rm S}$ being the Schwarzschild radius), larger than those in most previous studies, so that we can put relatively large angular momentum to the gas injected from the outer simulation boundary. The adopted Keplerian radius, at which the centrifugal force balances with the gravitational force, is $r_{\rm K}=300 r_{\rm S}$. The injected mass first falls and is accumulated at around this radius and then slowly accretes towards the central black hole via viscosity. We simulate such accretion processes, taking inverse and bulk Compton scattering into account. The simulated accretion flow is in a quasi-steady state inside $r_{\rm qss}\sim 200 r_{\rm S}$. Within this radius the flow properties are, as a whole, in good agreement with those described by the slim disk model except that the radial density profile of the underlying disk is much flatter, $\rho\propto r^{-0.73}$ (cf. $\rho\propto r^{-3/2}$ in the slim disk model), due probably to efficient convection. We find very weak outflow from inside $r\sim 200 r_{\rm S}$ unlike the past studies.

Using two-dimensional magnetohydrodynamic simulations, we investigated the gas dynamics around a black hole plunging into a molecular cloud. In these calculations, we assumed a parallel-magnetic-field layer in the cloud. The size of the accelerated region is far larger than the Bondi-Hoyle-Lyttleton radius, being approximately inversely proportional to the Alfv\'en Mach number for the plunging black hole. Our results successfully reproduce the "Y" shape in position velocity maps of the "Bullet" in the W44 molecular cloud. The size of the Bullet is also reproduced within an order of magnitude using a reasonable parameter set. This consistency supports the shooting model of the Bullet, according to which an isolated black hole plunged into a molecular cloud to form a compact broad-velocity-width feature.

We perform global three-dimensional (3D) radiation-hydrodynamic (RHD) simulations of out- flow from supercritical accretion flow around a 10 Msun black hole. We only solve the outflow part, starting from the axisymmetric 2D simulation data in a nearly steady state but with small perturbations in a sinusoidal form being added in the azimuthal direction. The mass accretion rate onto the black hole is ~10^2 L_E/c^2 in the underlying 2D simulation data and the outflow rate is ~10 L_E/c^2 (with LE and c being the Eddington luminosity and speed of light, respectively). We first confirm the emergence of clumpy outflow, which was discovered by the 2D RHD simulations, above the photosphere located at a few hundreds of Schwarzschild radii (r_S) from the central black hole. As prominent 3D features we find that the clumps have the shape of a torn sheet, rather than a cut string, and that they are rotating around the central black hole with a sub-Keplerian velocity at a distance of ~10^3 r_S from the center. The typical clump size is ~30 r_S or less in the radial direction, and is more elongated in the angular directions, ~hundreds of r_S at most. The sheet separation ranges from 50 to 150 r_S. We expect stochastic time variations when clumps pass across the line of the sight of a distant observer. Variation timescales are estimated to be several seconds for a black hole with mass of ten to several tens of Msun, in rough agreement with the observations of some ultra-luminous X-ray sources.

To understand why supercritical accretion is feasible onto a neutron star, we carefully examine the accretion flow dynamics by 2.5-dimensional general relativistic (GR) radiation magnetohydrodynamic (RMHD) simulations, comparing the cases of accretion onto a non-magnetized neutron star (NS) and that onto a black hole (BH). Supercritical BH accretion is relatively easy, since BH can swallow excess radiation energy, so that radiation flux can be inward in its vicinity. This mechanism can never work for NS which has a solid surface. In fact, we find that the radiation force is always outward. Instead, we found significant reduction in the mass accretion rate due to strong radiation-pressure driven outflow. The radiation flux $F_\mathrm{rad}$ is self-regulated such that the radiation force balances with the sum of gravity and centrifugal forces. Even when the radiation energy density much exceeds that expected from the Eddington luminosity $E_\mathrm{rad} \simeq F_\mathrm{rad}\tau/c> 10^2 L_\mathrm{Edd}/(4\pi r^2 c)$, the radiation flux is always kept below the certain value which makes it possible not to blow all the gas away from the disk. These effects make supercritical accretion feasible. We also find that a settling region, where accretion is significantly decelerated by radiation cushion, is formed around the NS surface. In the settling region, the radiation temperature and mass density roughly follow $T_\mathrm{rad} \propto r^{-1}$ and $\rho \propto r^{-3}$, respectively. No settling region appears around the BH so that matter can be directly swallowed by the BH with supersonic speed.

X-ray continuum spectra of super-Eddington accretion flow are studied by means of Monte Carlo radiative transfer simulations based on the radiation hydrodynamic simulation data, in which both of thermal and bulk Compton scatterings are taken into account. We compare the calculated spectra of accretion flow around black holes with masses of $M_{\rm BH} = 10, 10^2, 10^3$, and $10^4 M_\odot$ for a fixed mass injection rate (from the computational boundary at $10^3 r_{\rm s}$) of $10^3 L_{\rm Edd}/c^2$ (with $r_{\rm s}$, $L_{\rm Edd}$, and $c$ being the Schwarzschild radius, the Eddington luminosity, and the speed of light, respectively). The soft X-ray spectra exhibit mass dependence in accordance with the standard-disk relation; the maximum surface temperature is scaled as $T \propto M_{\rm BH}^{-1/4}$. The spectra in the hard X-ray bands, by contrast, look quite similar among different models, if we normalize the radiation luminosity by $M_{\rm BH}$. This reflects that the hard component is created by thermal and bulk Compton scattering of soft photons originating from an accretion flow in the over-heated and/or funnel regions, the temperatures of which have no mass dependence. The hard X-ray spectra can be reproduced by a Wien spectrum with temperature of $T\sim 3$ keV accompanied by a hard excess at photon energy above several keV. The excess spectrum can be well fitted with a power law with a photon index of $\Gamma \sim 3$. This feature is in good agreement with that of the recent NuSTAR observations of ULX (Ultra-Luminous X-ray sources).

By performing 2.5-dimensional general relativistic radiation magnetohydrodynamic simulations, we demonstrate supercritical accretion onto a non-rotating, magnetized neutron star, where the magnetic field strength of dipole fields is $10^{10}$ G on the star surface. We found the supercritical accretion flow consists of two parts; the accretion columns and the truncated accretion disk. The supercritical accretion disk, which appears far from the neutron star, is truncated at around $\sim 3R_*$ ($R_*=10^6$ cm is the neutron star radius), where the magnetic pressure via the dipole magnetic fields balances with the radiation pressure of the disks. The angular momentum of the disk around the truncation radius is effectively transported inward through magnetic torque by dipole fields, inducing the spin up of a neutron star. The evaluated spin up rate, $\sim -10^{-11}$ s s$^{-1}$, is consistent with the recent observations of the ultra luminous X-ray pulsars. Within the truncation radius, the gas falls onto neutron star along dipole fields, which results in a formation of accretion columns onto north and south hemispheres. The net accretion rate and the luminosity of the column are $\sim 66L_{\rm Edd}/c^2$ and $\lesssim 10L_{\rm Edd}$, where $L_{\rm Edd}$ is the Eddington luminosity and c is the light speed. Our simulations support a hypothesis whereby the ultra luminous X-ray pulsars are powered by the supercritical accretion onto the magnetized neutron stars.

We propose that superluminous transients that appear at central regions of active galactic nuclei (AGNs) such as CSS100217:102913+404220 (CSS100217) and PS16dtm, which reach near or super-Eddington luminosities of the central black holes, are powered by the interaction between accretion disk winds and clouds in broad-line regions (BLRs) surrounding them. If the disk luminosity temporary increases by, e.g., limit-cycle oscillations, leading to a powerful radiatively driven wind, strong shock waves propagate in the BLR. Because the dense clouds in the AGN BLRs typically have similar densities to those found in Type IIn supernovae, strong radiative shocks emerge and efficiently convert the ejecta kinetic energy to radiation. As a result, transients similar to Type IIn supernovae can be observed at AGN central regions. Since a typical black-hole disk wind velocity is ~ 0.1c where c is the speed of light, the ejecta kinetic energy is expected to be ~ 1e52 erg when ~ 1 Msun is ejected. This kinetic energy is transformed to radiation energy in a timescale for the wind to sweep up a similar mass to itself in the BLR, which is a few hundred days. Therefore, both luminosities (~ 1e44 erg/s) and timescales (~ 100 days) of the superluminous transients from AGN central regions match to those expected in our interaction model. If CSS100217 and PS16dtm are related to the AGN activities triggered by limit-cycle oscillations, they become bright again in coming years or decades.

We perform two-dimensional radiation hydrodynamical simulations of accretion flows onto a black hole (BH) with a mass of $10^3\leq M_{\rm BH}/M_{\odot} \lesssim 10^6$ in order to study rapid growth of BHs in the early Universe. For spherically symmetric flows, hyper-Eddington accretion onto the BH from outside the Bondi radius can occur unimpeded by radiation feedback only when the BH mass is higher than $\simeq 10^4~M_{\odot}(n_\infty/10^5~{\rm cm}^{-3})^{-1}(T_\infty/10^4~{\rm K})^{3/2}$, where $n_\infty$ and $T_\infty$ are the density and temperature of ambient gas. Here, we study the properties of accretion flows exposed to anisotropic radiation from a nuclear accretion disk with a luminosity higher than the Eddington value ($L_{\rm Edd}$) due to collimation toward the bipolar directions. We find that, unlike the spherically symmetric case, even less massive BHs with $M_{\rm BH} < 10^4~M_{\odot}$ can be fed by surrounding gas at high accretion rates of $\gtrsim L_{\rm Edd}/c^2$ through the equatorial plane, while ionized regions expand to the polar directions producing hot outflows with $T\sim 10^5$K. For more massive BHs with $M_{\rm BH}\gtrsim 5\times 10^5~M_{\odot}$, neutral gas through the equatorial plane totally covers the central radiating region due to the non-radial gas motions, and thus the emergent radiation in all directions is blocked. Because of efficient recombination by hydrogen, the entire flow results in neutral and warm gas with $T \simeq 8000~{\rm K}$ . The central BH is fed through the equator at the averaged rate of $\sim 5\times 10^4~L_{\rm Edd}/c^2$, which corresponds to $\sim 50~\%$ of the inflow rate from the Bondi radius. Moreover, radiation momentum absorbed by neutral hydrogen produces warm outflows toward the bipolar directions at $\sim 30~\%$ of the BH feeding rate and with a typical velocity of $\simeq 50~{\rm km~s}^{-1}$.

In order to reveal the origin of the ultra-fast outflows (UFOs) that are frequently observed in active galactic nuclei (AGNs), we perform two-dimensional radiation hydrodynamics simulations of the line-driven disk winds, which are accelerated by the radiation force due to the spectral lines. The line-driven winds are successfully launched for the range of $M_{\rm BH}=10^{6-9}M_\odot$ and $\varepsilon$=0.1-0.5, and the resulting mass outflow rate ($\dot{M_{\rm w}}$), momentum flux ($\dot{p_{\rm w}}$), and kinetic luminosity ($\dot{E_{\rm w}}$) are in the region containing 90% of the posterior probability distribution in the $\dot M_{\rm w}$-$L_{\rm bol}$ plane, $\dot p_{\rm w}$-$L_{\rm bol}$ plane, and $\dot E_{\rm w}$-$L_{\rm bol}$ plane shown in Gofford et al., where $M_{\rm BH}$ is the black hole mass, $\varepsilon$ is the Eddington ratio, and $L_{\rm bol}$ is the bolometric luminosity. The best-fitting relations in Gofford et al., $d\log\dot{M_{\rm w}}/d\log{L_{\rm bol}}\sim 0.9$, $d\log\dot{p_{\rm w}}/d\log{L_{\rm bol}}\sim 1.2$, and $d\log\dot{E_{\rm w}}/d\log{L_{\rm bol}}\sim 1.5$, are roughly consistent with our results, $d\log\dot{M_{\rm w}}/d\log{L_{\rm bol}}\sim 9/8$, $d\log\dot{p_{\rm w}}/d\log{L_{\rm bol}}\sim 10/8$, and $d\log\dot{E_{\rm w}}/d\log{L_{\rm bol}}\sim 11/8$. In addition, our model predicts that no UFO features are detected for the AGNs with $\varepsilon \lesssim 0.01$, since the winds do not appear. Also, only AGNs with $M_{\rm BH} \lesssim 10^8 M_\odot$ exhibit the UFOs when $\varepsilon \sim 0.025$. These predictions nicely agree with the X-ray observations (Gofford et al.). These results support that the line-driven disk wind is the origin of the UFOs.

Prompted by the recent discovery of pulsed emission from an ultra-luminous X-ray source, M82 X-2 ("ULX-pulsar"), we perform a two-dimensional radiation-hydrodynamic simulation of a super-critical accretion flow onto a neutron star through a narrow accretion column. We set an accretion column with a cone shape filled with tenuous gas with density of $10^{-4} {\rm g}~ {\rm cm}^{-3}$ above a neutron star and solve the two dimensional gas motion and radiative transfer within the column. The side boundaries are set such that radiation can freely escape, while gas cannot. Since the initial gas layer is not in a hydrostatic balance, the column gas falls onto the neutron-star surface, thereby a shock being generated. As a result, the accretion column is composed of two regions: an upper, nearly free-fall region and a lower settling region, as was noted by Basko \& Sunyaev (1976). The average accretion rate is very high; ${\dot M}\sim 10^{2-3} L_{\rm E}/c^2$ (with $L_{\rm E}$ being the Eddington luminosity), and so radiation energy dominates over gas internal energy entirely within the column. Despite the high accretion rate, the radiation flux in the laboratory frame is kept barely below $L_{\rm E}/(4\pi r^2)$ at a distance $r$ in the settling region so that matter can slowly accrete. This adjustment is made possible, since large amount of photons produced via dissipation of kinetic energy of matter can escape through the side boundaries. The total luminosity can greatly exceed $L_{\rm E}$ by several orders of magnitude, whereas the apparent luminosity observed from the top of the column is much less. Due to such highly anisotropic radiation fields, observed flux should exhibit periodic variations with the rotation period, provided that the rotation and magnetic axes are misaligned.

Using three-dimensional general relativistic radiation magnetohydrodynamics simulations of accretion flows around stellar mass black holes, we report that the relatively cold disk ($\gtrsim 10^{7}$K) is truncated near the black hole. Hot and less-dense regions, of which the gas temperature is $ \gtrsim 10^9$K and more than ten times higher than the radiation temperature (overheated regions), appear within the truncation radius. The overheated regions also appear above as well as below the disk, and sandwich the cold disk, leading to the effective Compton upscattering. The truncation radius is $\sim 30 r_{\rm g}$ for $\dot{M} \sim L_{\rm Edd}/c^2$, where $r_{\rm g}, \dot M, L_\mathrm{Edd}, c$ are the gravitational radius, mass accretion rate, Eddington luminosity, and light speed. Our results are consistent with observations of very high state, whereby the truncated disk is thought to be embedded in the hot rarefied regions. The truncation radius shifts inward to $\sim 10 r_{\rm g}$ with increasing mass accretion rate $\dot{M} \sim 100 L_{\rm Edd}/c^2$, which is very close to an innermost stable circular orbit. This model corresponds to the slim disk state observed in ultra luminous X-ray sources. Although the overheated regions shrink if the Compton cooling effectively reduces the gas temperature, the sandwich-structure does not disappear at the range of $\dot{M} \lesssim 100L_{\rm Edd}/c^2$. Our simulations also reveal that the gas temperature in the overheated regions depends on black hole spin, which would be due to efficient energy transport from black hole to disks through the Poynting flux, resulting gas heating.

Using two-dimensional radiation hydrodynamic simulations, we investigate origin of the ultra fast outflows (UFOs) that are often observed in luminous active galactic nuclei (AGNs). We found that the radiation force due to the spectral lines generates strong winds (line-driven disk winds) that are launched from the inner region of accretion disks (~30 Schwarzschild radii). A wide range of black hole masses ($M_{\rm BH}$) and Eddington ratios ($\varepsilon$) was investigated to study conditions for causing the line-driven winds. For $M_{\rm BH} = 10^6-10^9 M_\odot$ and $\varepsilon = 0.1-0.7$, funnel-shaped disk winds appear, in which dense matter is accelerated outward with an opening angle of 70-80 deg and with 10% of the light speed. If we observe the wind along its direction, the velocity, the column density, and the ionization state are consistent with those of the observed UFOs. As long as the obscuration by the torus does not affect the observations of X-ray bands, the UFOs could be statistically observed in about 13-28% of the luminous AGNs, which is not inconsistent with the observed ratio (~40%). We also found that the results are insensitive to the X-ray luminosity and the density of the disk surface. Thus, we can conclude that the UFOs could exist in any luminous AGNs, such as narrow-line Seyfert 1s (NLS1s) and quasars with $\varepsilon > 0.1$, in which fast line-driven winds are associated.

We apply our two-dimensional (2D), radially self-similar steady-state accretion flow model to the analysis of hydrodynamic simulation results of supercritical accretion flows. Self-similarity is checked and the input parameters for the model calculation, such as advective factor and heat capacity ratio, are obtained from time-averaged simulation data. Solutions of the model are then calculated and compared with the simulation results. We find that in the converged region of the simulation, excluding the part too close to the black hole, the radial distribution of azimuthal velocity $v_\phi$, density $\rho$ and pressure $p$ basically follows the self-similar assumptions, i.e. they are roughly proportional to $r^{-0.5}$, $r^{-n}$, and $r^{-(n+1)}$, respectively, where $n\sim0.85$ for the mass injection rate of $1000L_\mathrm{E}/c^2$, and $n\sim0.74$ for $3000L_\mathrm{E}/c^2$. The distribution of $v_r$ and $v_\theta$ agrees less with self-similarity, possibly due to convective motions in the $r\theta$ plane. The distribution of velocity, density and pressure in $\theta$ direction obtained by the steady model agrees well with the simulation results {within the calculation boundary of the steady model}. Outward mass flux in the simulations is overall directed toward polar angle of 0.8382 rad ($\sim 48.0^\circ$) for $1000L_\mathrm{E}/c^2$, and 0.7852 rad ($\sim 43.4^\circ$) for $3000L_\mathrm{E}/c^2$, and $\sim$94\% of the mass inflow are driven away as outflow, while outward momentum and energy fluxes are focused around the polar axis. Part of these fluxes lie in the region that are not calculated by the steady model, and special attention should be paid when the model is applied.

By performing two-dimensional radiation hydrodynamics simulations with large computational domain of 5000 Schwarzschild radius, we revealed that wide-angle outflow is launched via the radiation force from the super-critical accretion flows around black holes. The angular size of the outflow, of which the radial velocity (v_r) is over the escape velocity (v_esc), increases with an increase of the distance from the black hole. As a result, the mass is blown away with speed of v_r > v_esc in all direction except for the very vicinity of the equatorial plane, theta=0-85^circ, where theta is the polar angle. The mass ejected from the outer boundary per unit time by the outflow is larger than the mass accretion rate onto the black hole, ~150L_Edd/c^2, where L_Edd and c are the Eddington luminosity and the speed of light. Kinetic power of such wide-angle high-velocity outflow is comparable to the photon luminosity and is a few times larger than the Eddington luminosity. This corresponds to ~10^39-10^40 erg/s for the stellar mass black holes. Our model consistent with the observations of shock excited bubbles observed in some ultra-luminous X-ray sources (ULXs), supporting a hypothesis that ULXs are powered by the super-critical accretion onto stellar mass black holes.

By performing 2.5-dimensional special relativistic radiation magnetohydrodynamics simulations, we study the super-critical accretion disks and the outflows launched via the radiation force. We find that the outflow is accelerated by the radiation flux force, but the radiation drag force prevents the outflow velocity from increasing. The outflow velocity saturates around 30-40% of the light speed around the rotation axis, since then the flux force balances with the drag force. Our simulations show that the outflow velocity is kept nearly constant in the regime of \dot{M}_{BH} ~ 100-1000 L_{Edd}/c^2, where \dot{M}_{BH} is the mass accretion rate, L_{Edd} is the Eddington luminosity, and c is the light speed. Such a faster outflow is surrounded by a slower outflow of ~ 0.1c. This velocity is also determined by force balance between the radiation flux force and the radiation drag. The radiation drag works to collimate the slower outflow in cooperation with the Lorentz force, although the faster outflow is mainly collimated by the Lorentz force. The kinetic energy is carried by the slower outflow rather than by the faster outflow. The total kinetic luminosity of the outflow as well as the photon luminosity is ~ L_{Edd}, almost independent of the mass accretion rate.

We present our discovery of dramatic variability in SDSS J1100+4421 by the high-cadence transient survey Kiso Supernova Survey (KISS). The source brightened in the optical by at least a factor of three within about half a day. Spectroscopic observations suggest that this object is likely a narrow-line Seyfert 1 galaxy (NLS1) at z=0.840, however with unusually strong narrow emission lines. The estimated black hole mass of ~ 10^7 Msun implies bolometric nuclear luminosity close to the Eddington limit. SDSS J1100+4421 is also extremely radio-loud, with a radio loudness parameter of R ~ 4 x 10^2 - 3 x 10^3, which implies the presence of relativistic jets. Rapid and large-amplitude optical variability of the target, reminiscent of that found in a few radio- and gamma-ray loud NLS1s, is therefore produced most likely in a blazar-like core. The 1.4 GHz radio image of the source shows an extended structure with a linear size of about 100 kpc. If SDSS J1100+4421 is a genuine NLS1, as suggested here, this radio structure would then be the largest ever discovered in this type of active galaxies.

We study the dynamics of super-Eddington accretion flows by performing two-dimensional radiation-hydrodynamic simulations. Compared with previous works, in this paper we include the $T_{\theta\phi}$ component of the viscous stress and consider various values of the viscous parameter $\alpha$. We find that when $T_{\theta\phi}$ is included, the rotational speed of the high-latitude flow decreases, while the density increases and decreases at the high and low latitudes, respectively. We calculate the radial profiles of inflow and outflow rates. We find that the inflow rate decreases inward, following a power law form of $\dot{M}_{\rm in}\propto r^s$. The value of $s$ depends on the magnitude of $\alpha$ and is within the range of $\sim 0.4-1.0$. Correspondingly, the radial profile of density becomes flatter compared with the case of a constant $\dot{M}(r)$. We find that the density profile can be described by $\rho(r)\propto r^{-p}$, and the value of $p$ is almost same for a wide range of $\alpha$ ranging from $\alpha=0.1$ to $0.005$. The inward decrease of inflow accretion rate is very similar to hot accretion flows, which is attributed to the mass loss in outflows. To study the origin of outflow, we analyze the convective stability of slim disk. We find that depending on the value of $\alpha$, the flow is marginally stable (when $\alpha$ is small) or unstable (when $\alpha$ is large). This is different from the case of hydrodynamical hot accretion flow where radiation is dynamically unimportant and the flow is always convectively unstable. We speculate that the reason for the difference is because radiation can stabilize convection. The origin of outflow is thus likely because of the joint function of convection and radiation, but further investigation is required.

In order to understand the physical processes underlying clump formation in outflow from supercritical accretion flow, we performed two-dimensional radiation hydrodynamic (RHD) simulations. We focus our discussion on the nature of RHD instability in marginally optically thick, plane-parallel, super-Eddington atmosphere. Initially we set two-layered atmosphere with a density contrast of 100 exposed to strong, upward continuum-radiation force; the lower layer is denser than the upper one, condition for an RHD instability. We assume non-zero but negligible gravitational force, compared with the radiation force. We find that short wavelength perturbations first grow, followed by growth of longer wavelength patterns, which lead to the formation of clumpy structure. The typical size of clumps (clouds) corresponds to about one optical depth. An anti-correlation between the radiation pressure and the gas pressure is confirmed: this anti-correlation provides a damping mechanism of longer wavelength perturbations than the typical clump size. Matter and radiation energy densities are correlated. These features are exactly what we found in the radiation-magnetohydrodynamic (radiation-MHD) simulations of supercritical outflow.

We develop a numerical scheme for solving a fully special relativistic resistive radiation magnetohydrodynamics. Our code guarantees conservations of total mass, momentum and energy. Radiation energy density and radiation flux are consistently updated using the M-1 closure method, which can resolve an anisotropic radiation fields in contrast to the Eddington approximation as well as the flux-limited diffusion approximation. For the resistive part, we adopt a simple form of the Ohm's law. The advection terms are explicitly solved with an approximate Riemann solver, mainly HLL scheme, and HLLC and HLLD schemes for some tests. The source terms, which describe the gas-radiation interaction and the magnetic energy dissipation, are implicitly integrated, relaxing the Courant-Friedrichs-Lewy condition even in optically thick regime or a large magnetic Reynolds number regime. Although we need to invert $4\times 4$ (for gas-radiation interaction) and $3\times 3$ (for magnetic energy dissipation) matrices at each grid point for implicit integration, they are obtained analytically without preventing massive parallel computing. We show that our code gives reasonable outcomes in numerical tests for ideal magnetohydrodynamics, propagating radiation, and radiation hydrodynamics. We also applied our resistive code to the relativistic Petschek type magnetic reconnection, revealing the reduction of the reconnection rate via the radiation drag.

The tidal disruption event by a supermassive black hole in Swift J1644+57 can trigger limit-cycle oscillations between a supercritically accreting X-ray bright state and a subcritically accreting X-ray dim state. Time evolution of the debris gas around a black hole with mass $M=10^{6} {\MO}$ is studied by performing axisymmetric, two-dimensional radiation hydrodynamic simulations. We assumed the $\alpha$-prescription of viscosity, in which the viscous stress is proportional to the total pressure. The mass supply rate from the outer boundary is assumed to be ${\dot M}_{\rm supply}=100L_{\rm Edd}/c^2$, where $L_{\rm Edd}$ is the Eddington luminosity, and $c$ is the light speed. Since the mass accretion rate decreases inward by outflows driven by radiation pressure, the state transition from a supercritically accreting slim disk state to a subcritically accreting Shakura-Sunyaev disk starts from the inner disk and propagates outward in a timescale of a day. The sudden drop of the X-ray flux observed in Swift J1644+57 in August 2012 can be explained by this transition. As long as ${\dot M}_{\rm supply}$ exceeds the threshold for the existence of a radiation pressure dominant disk, accumulation of the accreting gas in the subcritically accreting region triggers the transition from a gas pressure dominant Shakura-Sunyaev disk to a slim disk. This transition takes place at $t {\sim}~50/({\alpha}/0.1)$ days after the X-ray darkening. We expect that if $\alpha > 0.01$, X-ray emission with luminosity $\gtrsim 10^{44}$ ${\rm erg}{\cdot}{\rm s}^{-1}$ and jet ejection will revive in Swift J1644+57 in 2013--2014.

Significant fraction of matter in supercritical (or super-Eddington) accretion flow is blown away by radiation force, thus forming outflows, however, the properties of such radiation-driven outflows have been poorly understood. We have performed global two-dimensional radiaion-magnetohydrodynamic simulations of supercritical accretion flow onto a black hole with 10 or 10^8 solar masses in a large simulation box of 514 r_S x 514 r_S (with r_S being the Schwarzschild radius). We confirm that uncollimated outflows with velocities of 10 percents of the speed of light emerge from the innermost part of the accretion flow over wide angles of 10 - 50 degree from the disk rotation axis. Importantly, the outflows exhibit clumpy structure above heights of ~ 250 r_S. The typical size of the clumps is ~ 10 r_S, which corresponds to one optical depth, and their shapes are slightly elongated along the outflow direction. Since clumps start to form in the layer above which (upward) radiation force overcomes (downward) gravity force, Rayleigh-Taylor instability seems to be of primary cause. In addition, a radiation hydrodynamic instability, which arises when radiation funnels through radiation-pressure supported atmosphere, may also help forming clumps of one optical depth. Magnetic photon bubble instability seems not to be essential, since similar clumpy outflow structure is obtained in non-magnetic radiation-hydrodynamic simulations. Since the spatial covering factor of the clumps is estimated to be ~ 0.3 and since they are marginally optically thick, they will explain at least some of rapid light variations of active galactic nuclei. We further discuss a possibility of producing broad-line clouds by the clumpy outflow.

We study the origin of absorption features on the blue side of the C IV broad emission line of the large-separation lensed quasar SDSS J1029+2623 at z_em ~ 2.197. The quasar images, produced by a foreground cluster of galaxies, have a maximum separation angle of ~ 22".5. The large angular separation suggests that the sight-lines to the quasar central source can go through different regions of outflowing winds from the accretion disk of the quasar, providing a unique opportunity to study the structure of outflows from the accretion disk, a key ingredient for the evolution of quasars as well as for galaxy formation and evolution. Based on medium- and high-resolution spectroscopy of the two brightest images conducted at the Subaru telescope, we find that each image has different intrinsic levels of absorptions, which can be attributed either to variability of absorption features over the time delay between the lensed images, ~ 774 days, or to the fine structure of quasar outflows probed by the multiple sight-lines toward the quasar. While both these scenarios are consistent with the current data, we argue that they can be distinguished with additional spectroscopic monitoring observations.

We propose an explicit-implicit scheme for numerically solving Special Relativistic Radiation Hydrodynamic (RRHD) equations, which ensures a conservation of total energy and momentum (matter and radiation). In our scheme, 0th and 1st moment equations of the radiation transfer equation are numerically solved without employing a flux-limited diffusion (FLD) approximation. For an hyperbolic term, of which the time scale is the light crossing time when the flow velocity is comparable to the speed of light, is explicitly solved using an approximate Riemann solver. Source terms describing an exchange of energy and momentum between the matter and the radiation via the gas-radiation interaction are implicitly integrated using an iteration method. The implicit scheme allows us to relax the Courant-Friedrichs-Lewy condition in optically thick media, where heating/cooling and scattering timescales could be much shorter than the dynamical timescale. We show that our numerical code can pass test problems of one- and two-dimensional radiation energy transport, and one-dimensional radiation hydrodynamics. Our newly developed scheme could be useful for a number of relativistic astrophysical problems. We also discuss how to extend our explicit-implicit scheme to the relativistic radiation magnetohydrodynamics.

The disk wind, which is powered by the radiation force due to spectral lines (line force), is studied for broad absorption line (BAL) quasars. We investigate the structure of the disk wind based on the non-hydrodynamic method and compare with wind properties inferred from X-ray observations of BAL quasars. In this paper, we apply the stellar wind theory to the initial condition (the mass outflow rate at the base of the wind). We found the funnel-shaped winds with a half opening angle of 50^{circ} for the case of epsilon=0.3-0.9 and M_{BH}=10^{7-8.5}M_odot, where epsilon is the Eddington ratio and M_{BH} is the black hole mass. Thus, the absorption features are observed for an observer of which a viewing angle is around 50^{circ}. A probability of BAL quasars is 7-11%, which is roughly consistent the abundance ratio of BAL quasars, 10-15%. Here, the probability is estimated by the solid angle, that the absorbing features would be detected, divided by 4pi. In contrast, if the Eddington ratio is smaller than 0.01 or if the black hole is very massive, M_{BH} < 10^9M_{odot}, the disk wind is not launched due to the less effective line force. Then, the quasars are identified as non-BAL quasars independently of the observer's viewing angle.

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 present the detailed global structure of black hole accretion flows and outflows through newly performed two-dimensional radiation-magnetohydrodynamic simulations. By starting from a torus threaded with weak toroidal magnetic fields and by controlling the central density of the initial torus, rho_0, we can reproduce three distinct modes of accretion flow. In model A with the highest central density, an optically and geometrically thick supercritical accretion disk is created. The radiation force greatly exceeds the gravity above the disk surface, thereby driving a strong outflow (or jet). Because of the mild beaming, the apparent (isotropic) photon luminosity is ~22L_E (where L_E is the Eddington luminosity) in the face-on view. Even higher apparent luminosity is feasible if we increase the flow density. In model B with a moderate density, radiative cooling of the accretion flow is so efficient that a standard-type, cold, and geometrically thin disk is formed at radii greater than ~7R_S (where R_S is the Schwarzschild radius), while the flow is radiatively inefficient otherwise. The magnetic-pressure-driven disk wind appears in this model. In model C the density is too low for the flow to be radiatively efficient. The flow thus becomes radiatively inefficient accretion flow, which is geometrically thick and optically thin. The magnetic-pressure force, in cooperation with the gas-pressure force, drives outflows from the disk surface, and the flow releases its energy via jets rather than via radiation. Observational implications are briefly discussed.

Relativistic jets from compact objects are ubiquitous phenomena in the Unvierse, but their driving mechanism has been an enigmatic issue over many decades. Two basic models have been extensively discussed: magnetohydrodynamic (MHD) jets and radiation-hydrodynamic (RHD) jets. Currently, the former is more widely accepted, since magnetic field is expected to provide both the acceleration and collimation mechanisms, whereas radiation field cannot collimate outflow. Here, we propose a new type of jets, radiation-magnetohydrodynamic (RMHD) jets, based on our global RMHD simulation of luminous accretion flow onto a black hole shining above the Eddington luminosity. The RMHD jet can be accelerated up to the relativistic speed by the radiation-pressure force and is collimated by the Lorentz force of a magnetic tower, inflated magnetic structure made by toroidal magnetic field lines accumulated around the black hole, though radiation energy greatly dominates over magnetic energy. This magnetic tower is collimated by a geometrically thick accretion flow supported by radiation-pressure force. This type of jet may explain relativistic jets from Galactic microquasars, appearing at high luminosities.

We perform a three-dimensional nested-grid radiation magneto-hydrodynamics (RMHD) simulation with self-gravity to study the early phase of the low-mass star formation process from a rotating molecular cloud core to a first adiabatic core just before the second collapse begins. Radiation transfer is handled with the flux-limited diffusion approximation, operator-splitting and implicit time-integrator. In the RMHD simulation, the outer region of the first core attains a higher entropy and the size of first core is larger than that in the magnetohydrodynamics simulations with the barotropic approximation. Bipolar molecular outflow consisting of two components is driven by magnetic Lorentz force via different mechanisms, and shock heating by the outflow is observed. Using the RMHD simulation we can predict and interpret the observed properties of star-forming clouds, first cores and outflows with millimeter/submillimeter radio interferometers, especially the Atacama Large Millimeter/submillimeter Array (ALMA).

By solving radiative transfer equations, we examine three-dimensional radiative properties of a magnetohydrodynamic accretion flow model confronting with the observed spectrum of Sgr A*, in the vicinity of supermassive black hole at the Galactic centre. As a result, we find that the core of radio emission is larger than the size of the event horizon shadow and its peak location is shifted from the gravitational centre. We also find that the self-absorbed synchrotron emissions by the superposition of thermal electrons within a few tens of the Schwartzschild radius can account for low-frequency spectra below the critical frequency $\nu_{c}\approx 10^{12}$ Hz. Above the critical frequency, the synchrotron self-Compton emission by thermal electrons can account for variable emissions in recent near-infrared observations. In contrast to the previous study by Ohsuga et al. (2005), we found that the X-ray spectra by Bremsstrahlung emission of thermal electrons for the different mass accretion rates can be consistent with both the flaring state and the quiescent state of Sgr A* observed by {\it Chandra}.

Photon trapping and outflow are two key physics associated with the supercritical accretion flow. We investigate the conflict between these two processes based on two-dimensional radiation-hydrodynamic (RHD) simulation data and construct a simplified (radially) one-dimensional model. Mass loss due to outflow, which is not considered in the slim-disk model, will reduce surface density of the flow, and if very significant, it will totally suppress photon trapping effects. If the photon trapping is very significant, conversely, outflow will be suppressed because radiation pressure force will be reduced. To see what actually occurs, we examine the RHD simulation data and evaluate the accretion rate and outflow rate as functions of radius. We find that the former monotonically decreases, while the latter increases, as the radius decreases. However, the former is kept constant at small radii, inside several Schwarzschild radii, since the outflow is suppressed by the photon trapping effects. To understand the conflict between the photon trapping and outflow in a simpler way, we model the radial distribution of the accretion rate from the simulation data and build up a new (radially) one-dimensional model, which is similar to the slim-disk model but incorporates the mass loss effects due to the outflow. We find that the surface density (and, hence, the optical depth) is much reduced even inside the trapping radius, compared with the case without outflow, whereas the effective temperature distribution hardly changes. That is, the emergent spectra do not sensitively depend on the amount of mass outflow. We conclude that the slim-disk approach is valid for interpreting observations, even if the outflow is taken into account.

Supercritical accretion flows inevitably produce radiation-pressure driven outflows, which will Compton up-scatter soft photons from the underlying accretion flow, thereby making hard emission. We perform two dimensional radiation hydrodynamic simulations of supercritical accretion flows and outflows, incorporating such Compton scattering effects, and demonstrate that there appears a new hard spectral state at higher photon luminosities than that of the slim-disk state. In this state, as the photon luminosity increases, the photon index decreases and the fraction of the hard emission increases. The Compton $y$-parameter is of the order of unity (and thus the photon index will be $\sim 2$) when the apparent photon luminosity is ${\sim}30L_{\rm E}$ (with $L_{\rm E}$ being the Eddington luminosity) for nearly face-on sources. This explains the observed spectral hardening of the ULX NGC1313 X-2 in its brightening phase and thus supports the model of supercritical accretion onto stellar mass black holes in this ULX.

Black-hole accretion systems are known to possess several distinct modes (or spectral states), such as low/hard state, high/soft state, and so on. Since the dynamics of the corresponding flows is distinct, theoretical models were separately discussed for each state. We here propose a unified model based on our new, global, two-dimensional radiation-magnetohydrodynamic simulations. By controlling a density normalization we could for the first time reproduce three distinct modes of accretion flow and outflow with one numerical code. When the density is large (model A), a geometrically thick, very luminous disk forms, in which photon trapping takes place. When the density is moderate (model B), the accreting gas can effectively cool by emitting radiation, thus generating a thin disk, i.e., the soft-state disk. When the density is too low for radiative cooling to be important (model C), a disk becomes hot, thick, and faint; i.e., the hard-state disk. The magnetic energy is amplified within the disk up to about twice, 30%, and 20% of the gas energy in models A, B, and C, respectively. Notably, the disk outflows with helical magnetic fields, which are driven either by radiation pressure force or magnetic pressure force, are ubiquitous in any accretion modes. Finally, our simulations are consistent with the phenomenological alpha-viscosity prescription, that is, the disk viscosity is proportional to the pressure.

We investigate the stability of black hole accretion disks in a primordial environment (POP III disks for short), by solving the vertical structure of optically thick disks, including convective energy transport, and by employing a one-zone model for optically thin isothermal disks. Because of the absence of metals in POP III disks, we find significant differences in stability associated with ionization between POP III disks and the disks of solar metallicity. An unstable branch in S-shaped equilibrium curves on the Mdot-Sigma (mass accretion rate - surface density) plane extends to a larger surface density compared with the case of disks of solar metallicity. The resulting equilibrium loci indicate that quasi-periodic oscillations in luminosity can also be driven in POP III disks, and their maximal luminosity is typically by an order of magnitude larger than that of the disks of solar metallicity. Such a strong outburst of POP III disks can be observed by future huge telescopes, in case that the mass is supplied onto the disks at the Bondi accretion rates in typical virialized small dark halos.

Although the occurrence of steady supercritical disk accretion onto a black hole has been speculated about since the 1970s, it has not been accurately verified so far. For the first time, we previously demonstrated it through two-dimensional, long-term radiation-hydrodynamic simulations. To clarify why this accretion is possible, we quantitatively investigate the dynamics of a simulated supercritical accretion flow with a mass accretion rate of ~10^2 L_E/c^2 (with L_E and c being, respectively, the Eddington luminosity and the speed of light). We confirm two important mechanisms underlying supercritical disk accretion flow, as previously claimed, one of which is the radiation anisotropy arising from the anisotropic density distribution of very optically thick material. We qualitatively show that despite a very large radiation energy density, E_0>10^2L_E/(4 pi r^2 c) (with r being the distance from the black hole), the radiative flux F_0 cE_0/tau could be small due to a large optical depth, typically tau 10^3, in the disk. Another mechanism is photon trapping, quantified by vE_0, where v is the flow velocity. With a large |v| and E_0, this term significantly reduces the radiative flux and even makes it negative (inward) at r<70r_S, where r_S is the Schwarzschild radius. Due to the combination of these effects, the radiative force in the direction along the disk plane is largely attenuated so that the gravitational force barely exceeds the sum of the radiative force and the centrifugal force. As a result, matter can slowly fall onto the central black hole mainly along the disk plane with velocity much less than the free-fall velocity, even though the disk luminosity exceeds the Eddington luminosity. Along the disk rotation axis, in contrast, the strong radiative force drives strong gas outflows.

We performed two-dimensional radiation hydrodynamic simulations of supercritical accretion flows around neutron stars (NSs). In contrast with the accretion flows onto black holes (BHs), we find that the shell-shaped high-density regions form around the NSs, since the radiation force is enhanced in the innermost regions. The enhanced radiation force drives strong outflows above and below the disk. The mass-accretion rate onto the NS exceeds the critical rate, $L_{\rm E}/c^2$, with $L_{\rm E}$ being the Eddington luminosity. However it is about $20-30%$ of that onto the BH, under the condition that we employ the same mass-input rate, $\dot{M}_{\rm input}$, which is mass injected from the outer disk boundary per unit time. The mass-outflow rate is a few-times larger in flows around NSs than in flows around BHs. The supercritical NS accretion flows mainly release the accretion energy as the kinetic energy of the outflows, though the disk luminosity is predominant over the kinetic energy output rate in the BH accretion flows. The resulting velocity and mass-outflow rate of the outflows are $0.2-0.3c$ and $150-700L_{\rm E}/c^2$, respectively, for the mass-input rate of $3\times 10^2\lsim \dot{M}_{\rm input}/(L_{\rm E}/c^2)\lsim 3\times 10^3$. This implies that the SS433 jets can be roughly explained by the supercritical accretion onto a NS. However, the collimation angle of the outflows in our simulations ($\sim 20^\circ$) is larger than that of the SS433 jets (a few degrees).

We consider photon emission from a supercritical accretion disk in which photons in the inner regions are trapped and advected in towards the center. Such disks are believed to be present in failed supernovae and gamma-ray bursts, and possibly also in ultra-luminous X-ray sources. We show that the luminosity from a supercritical accretion disk is proportional to the logarithm of the mass accretion rate when the vertical profile of the matter density is exponential. We obtain analytical results also for other density profiles, but the dependence of the luminosity on the accretion rate is no longer logarithmic.

The supercritical disk accretion flow with radiatively driven outflows is studied based on two-dimensional radiation-hydrodynamic simulations for a wide range of the mass input rate, $\dot{M}_{\rm input}$, which is the mass supplied from the outer region to the disk per unit time. The $\alpha$-prescription is adopted for the viscosity. We employ $\alpha=0.5$, as well as $\alpha=0.1$ for $\dot{M}_{\rm input}\ge 3\times 10^2L_{\rm E}/c^2$ and only $\alpha=0.5$ for $\dot{M}_{\rm input}\le 10^2L_{\rm E}/c^2$, where $L_{\rm E}$ is the Eddington luminosity and $c$ is the speed of light. The quasi-steady disk and radiately driven outflows form in the case in which the mass input rate highly exceeds the critical rate, $\dot{M}_{\rm input}>3\times 10^2 L_{\rm E}/c^2$. Then, the disk luminosity as well as the kinetic energy output rate by the outflow exceeds the Eddington luminosity. The moderately supercritical disk, $\dot{M}_{\rm input}\sim 10-10^2 L_{\rm E}/c^2$, exhibits limit-cycle oscillations. The disk luminosity goes up and down across the Eddington luminosity, and the radiatively driven outflows intermittently appear. The time averaged mass, momentum, and kinetic energy output rates by the outflow as well as the disk luminosity increase with increase of the mass input rate, $\propto \dot{M}_{\rm input}^{0.7-1.0}$ for $\alpha=0.5$ and $\propto \dot{M}_{\rm input}^{0.4-0.6}$ for $\alpha=0.1$. Our numerical simulations show that the radiatively driven outflow model for the correlation between black hole mass and bulge velocity dispersion proposed by \citeauthor{SR98} and \citeauthor{King03} is successful if $\dot{M}_{\rm input}c^2/L_{\rm E} \sim$ a few 10 ($\alpha=0.5$) or $\gsim$ a few ($\alpha=0.1$).

Recently, the issue of the role of the Eddington limit in accretion discs became a matter of debate. While the classical (spherical) Eddington limit is certainly an over-simplification, it is not really clear how to treat it in a flattened structure like an accretion disc. We calculate the critical accretion rates and resulting disc luminosities for various disc models corresponding to the classical Eddington limit by equating the attractive and repulsive forces locally. We also discuss the observational appearance of such highly accreting systems by analyzing their spectral energy distributions. Our calculations indicate that the allowed mass accretion rates differ considerably from what one expects by applying the Eddington limit in its classical form, while the luminosities only weakly exceed their classical equivalent. Depending on the orientation of the disc relative to the observer, mild relativistic beaming turns out to have an important influence on the disc spectra. Thus, possible super-Eddington accretion, combined with mild relativistic beaming, supports the idea that ultraluminous X-ray sources host stellar mass black holes and accounts partially for the observed high temperatures of these objects.

Spectral properties of super-Eddington accretion flows are investigated by means of a parallel line-of-sight calculation. The subjacent model, taken from two-dimensional radiation hydrodynamic simulations by Ohsuga et al. (2005), consists of a disc accretion region and an extended atmosphere with high velocity outflows. The non-gray radiative transfer equation is solved, including relativistic effects, by applying the FLD approximation. The calculated spectrum is composed of a thermal, blackbody-like emission from the disc which depends sensitively on the inclination angle, and of high energy X-ray and gamma-ray emission from the atmosphere. We find mild beaming effects in the thermal radiation for small inclination angles. If we compare the face-on case with the edge-on case, the average photon energy is larger by a factor of ~1.7 due mainly to Doppler boosting, while the photon number density is larger by a factor of ~3.7 due mainly to anisotropic matter distribution around the central black hole. This gives an explanation for the observed X-ray temperatures of ULXs which are too high to be explained in the framework of intermediate-mass black holes. While the main features of the thermal spectral component are consistent with more detailed calculations of slim accretion discs, the atmosphere induces major changes in the high-energy part, which cannot be reproduced by existing models. In order to interpret observational data properly, simple approaches like the Eddington-Barbier approximation cannot be applied.

We investigate the time evolution of luminous accretion disks around black holes, conducting the two-dimensional radiation-hydrodynamic simulations. We adopt the alpha prescription for the viscosity. The radial-azimuthal component of viscous stress tensor is assumed to be proportional to the total pressure in the optically thick region, while the gas pressure in the optically thin regime. The viscosity parameter, alpha, is taken to be 0.1. We find the limit-cycle variation in luminosity between high and low states. When we set the mass input rate from the outer disk boundary to be 100 L_E/c^2, the luminosity suddenly rises from 0.3L_E to 2L_E, where L_E is the Eddington luminosity. It decays after retaining high value for about 40 s. Our numerical results can explain the variation amplitude and duration of the recurrent outbursts observed in microquasar, GRS 1915+105. We show that the multi-dimensional effects play an important role in the high-luminosity state. In this state, the outflow is driven by the strong radiation force, and some part of radiation energy dissipated inside the disk is swallowed by the black hole due to the photon-trapping effects. This trapped luminosity is comparable to the disk luminosity. We also calculate two more cases: one with a much larger accretion rate than the critical value for the instability and the other with the viscous stress tensor being proportional to the gas pressure only even when the radiation pressure is dominant. We find no quasi-periodic light variations in these cases. This confirms that the limit-cycle behavior found in the simulations is caused by the disk instability.

We investigate the dependence of the viewing angle in supercritical accretion flows and discuss the observational implications of galactic black-hole candidates and ultraluminous X-ray sources. When the mass accretion rate exceeds the critical rate, then the shape of the disk is geometrically thick due to the enhanced radiation pressure. The model spectra of supercritical accretion flows strongly depend on the inclination angle. Because the outer disk blocks the emission from the disk inner region for high inclination angle. We also find that the spectral properties of low-inclination angle and low accretion-rate disks are very similar to those of high-inclination and high accretion rate disks. That is, if an object has a high inclination and high accretion rate, such a system suffers from self-occultation and the spectrum will be extremely soft. Therefore, we cannot discriminate these differences from spectrum shapes only. Conversely, if we use the self-occultation properties, we could constrain the inclination angle of the system. We suggest that some observed high temperature ultraluminous X-ray sources have near face-on geometry, i < 40, and Galactic black hole candidate, XTE J1550-564, possesses relatively high-inclination angles, i > 60.

The quasi-steady structure of super-critical accretion flows around a black hole is studied based on the two-dimensional radiation-hydrodynamical (2D-RHD) simulations. The super-critical flow is composed of two parts: the disk region and the outflow regions above and below the disk. Within the disk region the circular motion as well as the patchy density structure are observed, which is caused by Kelvin-Helmholtz instability and probably by convection. The mass-accretion rate decreases inward, roughly in proportion to the radius, and the remaining part of the disk material leaves the disk to form outflow because of strong radiation pressure force. We confirm that photon trapping plays an important role within the disk. Thus, matter can fall onto the black hole at a rate exceeding the Eddington rate. The emission is highly anisotropic and moderately collimated so that the apparent luminosity can exceed the Eddington luminosity by a factor of a few in the face-on view. The mass-accretion rate onto the black hole increases with increase of the absorption opacity (metalicity) of the accreting matter. This implies that the black hole tends to grow up faster in the metal rich regions as in starburst galaxies or star-forming regions.

In spite of a large number of global three-dimensional (3D) magneto-hydrodynamical (MHD) simulations of accretion flows and jets being made recently, their astrophysical relevance for realistic situations is not well known. In order to examine to what extent the simulated MHD flows can account for the observed spectral energy distribution (SED) of Sagittarius A* (Sgr A*), for the first time we calculate the emergent spectra from 3D MHD flows in a wide range of wavelengths (from radio to X-ray) by solving the 3D radiative transfer equations. We use the simulation data by Kato, Mineshige, and Shibata (2004) and perform Monte Carlo radiative transfer simulations, in which synchrotron emission/absorption, free-free emission/absorption, and Compton/inverse Compton scattering are taken into account. We assume two temperature plasmas and calculate electron temperatures by solving the electron energy equation. Only thermal electrons are considered. It is found that the 3D MHD flow generally over-produces X-rays by means of bremsstrahlung radiation from the regions at large radii. A flatter density profile, r^{-a} with a<1, than that of the advection-dominated accretion flow (ADAF), r^{-3/2}, is the main reason for this. If we restrict the size of the emission region to be as small as 10r_s, where r_s is the Schwarzschild radius, the MHD model can reproduce the basic features of the observed SED of Sgr A* during its flaring state. Yet, the spectrum in the quiescent state remains to be understood. We also calculate the time-dependent spectral changes, finding that the fluxes fluctuate in a wide range of the frequency and the flux at each wavelength does not always vary coherently.

We investigate the spectral energy distribution (SED) of the supercritical disk accretion flows around black holes by solving the multi-frequency zeroth moment equation of the radiation transfer equation under the flux-limited diffusion approximation as well as the energy equation of gas. Special attention is paid to the photon-trapping effects, the effects that photons are trapped within accretion flow and are swallowed by a black hole, with little being radiated away. It is found that when flow luminosity is below about twice the Eddington luminosity (L<2L_E) photon trapping is ineffective and the SED shifts to the higher-energy side as L increases. The peak frequency at which the SED reaches its peak becomes at maximum three times higher than that given by the standard-disk model, which might resolve the so-called 'too hot accretion disk problem.' When L>2L_E, in contrast the peak frequency of the SED tends to decrease with increase of mass-accretion rate. This is due to enhanced photon-trapping. Since high-energy photons are generated near the equatorial plane, they can be more effectively trapped in flow than low-energy ones and, hence, the high energy part of radiation is suppressed. Finally, the observed time variation of IC 342 S1, which is an ultraluminous X-ray source, in the X-ray HR diagram can be explained by the modulation of the mass accretion rate. In comparison with the observational data, we estimate the mass of black hole in IC 342 S1 is around 100 solar mass.

We investigate a corona coupled with underlying disk through magnetic field and radiation field, and present emergent spectra calculated from Monte Carlo simulations. The spectra are determined for given black-hole mass and accretion rate. We find two types of solutions corresponding for hard spectrum and soft spectrum. In the hard-spectrum solution, the accretion energy is dominantly dissipated in the corona, supporting a strong corona above a cool disk; The hard X-ray spectral indices are around 1.1, same for different accretion rates. In the soft-spectrum solution, the accretion energy is mainly dissipated in the disk. The coronal temperature and density are quite low. Consequently, the spectra are dominated by the disk radiation peaking at UV and soft X-rays. For low-luminosity systems, there exists only the solution of hard spectra; While for high-luminosity systems, there exist both solutions of hard and soft spectra. For middle-luminosity systems, besides the hard spectra, moderately soft spectra composed of an inner soft-spectrum solution and an outer hard-spectrum solution may occur, the softness of which increases with increasing luminosity. The hard spectra are close to the observed spectra in Seyfert galaxies and radio-quiet QSOs. The composite spectra may account for the diversity of broad band spectra observed in narrow-line Seyfert 1 galaxies.

We investigate the photon-trapping effects in the super-critical black hole accretion flows by solving radiation transfer as well as the energy equations of radiation and gas. It is found that the slim-disk model generally overestimates the luminosity of the disk at around the Eddington luminosity (L_E) and is not accurate in describing the effective temperature profile, since it neglects time delay between energy generation at deeper inside the disk and energy release at the surface. Especially, the photon-trapping effects are appreciable even below L ~ L_E, while they appear above ~ 3L_E according to the slim disk. Through the photon-trapping effects, the luminosity is reduced and the effective temperature profile becomes flatter than r^{-3/4} as in the standard disk. In the case that the viscous heating is effective only around the equatorial plane, the luminosity is kept around the Eddington luminosity even at very large mass accretion rate, Mdot>>L_E/c^2. The effective temperature profile is almost flat, and the maximum temperature decreases in accordance with rise in the mass accretion rate. Thus, the most luminous radius shifts to the outer region when Mdot/(L_E/c^2) >> 10^2. In the case that the energy is dissipated equally at any heights, the resultant luminosity is somewhat larger than in the former case, but the energy-conversion efficiency still decreases with increase of the mass accretion rate, as well. The most luminous radius stays around the inner edge of the disk in the latter case. Hence, the effective temperature profile is sensitive to the vertical distribution of energy production rates, so is the spectral shape. Future observations of high L/L_E objects will be able to test our model.

We explore the formation of dusty gas walls induced by a circumnuclear starburst around an AGN. We concentrate our attention on the role of the radiation force by a starburst as well as an AGN, where the effects of optical depth of dusty gas are taken into consideration. In two-dimensional axisymmetric space, we analyze the configuration and the stability of geometrically thin walls which are in balance between radiation pressure and gravity. As a result, it is shown that the radiation force by the circumnuclear starburst works to stabilize optically thick walls surrounding the nucleus. In the case of a brighter starburst with a fainter AGN, there form double walls, an inner one of which is located between the nucleus and the circumnuclear starburst, and an outer one of which enshrouds both the starburst regions and the nucleus. The total extinction (Av) of both walls turns out to be ~10 for a brighter starburst. As a consequence, double walls could heavily obscure the nucleus to make the AGN type 2. The outer wall may provide an explanation for the recent indications for large-scale obscuring materials in Seyfert 2's. Also, it is predicted that the AGN type is shifts from type 2 to type 1 in several times 10^7 yr according to the stellar evolution in the starburst. In contrast, if the AGN itself is much brighter than the starburst as a quasar is, then neither wall forms regardless of the starburst activity and the nucleus is likely to be identified as type 1. To conclude, the radiatively-supported gas walls could be responsible for the putative correlation between AGN type and the starbursts, whereby Seyfert 2 galaxies are more frequently associated with circumnuclear starbursts than type 1, whereas quasars are mostly observed as type 1 regardless of star-forming activity in the host galaxies.

We consider the radiation-hydrodynamical formation of a dusty wall in the circumnuclear regions of the Circinus Seyfert 2 galaxy. We focus on the radiative flux force due to the circumnuclear starburst, the nuclear starburst and the active galactic nucleus (AGN), and analyze the equilibrium configuration and stability of the dusty gas in the circumnuclear regions. The mass distributions used to model the gravitational fields are determined by the observed rotation velocities of the stars and of the gas. Moreover, by using the simple stellar evolution in the circumnuclear starburst, the bolometric luminosity of the starburst is estimated as a function of time. As a result, it is shown that the radiation force by the circumnuclear starburst does play an important role to build up a radiatively-supported dusty wall of $\la$ 100 pc, which obscures the circumnuclear regions with $A_V$ of a few magnitudes. It is also found that the age of a circumnuclear starburst when the circumnuclear regions are enshrouded by the dusty wall is constrained to be $\la 10^8$ yr, which is consistent with the age observationally inferred for the starburst, $4 10^7$ yr $-1.5 10^8$ yr.

By considering the radiative force by a circumnuclear starburst as well as an AGN, we analyze the equilibrium configuration and the stability of dusty gas in the circumnuclear regions. It is found that the radiative force by an intensive starburst can support a stable gaseous wall with a scale-height of several hundred parsecs. Moreover, by taking the simple stellar evolution in the starburst into account, we find that the covering factor of the wall decreases on a time-scale of several $10^7$ yr. The large-scale wall, if formed, works to obscure the nucleus due to the dust opacity. Hence, it is anticipated that the index of AGN type tends to shift from higher to lower in several $10^7$ yr according as the circumnuclear starburst becomes dimmer. On the other hand, if the AGN itself is brighter than the circumnuclear starburst (e.g. quasar case), no stable large-scale wall forms. In that case, the AGN is highly probably identified as type 1. The present mechanism may provide a physical explanation for the putative correlation between AGN type and host properties that Sy2's are more frequently associated with circumnuclear starbursts than Sy1's, whereas quasars are mostly observed as type 1 regardless of star-forming activity in the host galaxies.