A significant amount 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} M_{\odot}$ in a large simulation box of 514 $r_{\rm S} \times 514 r_{\rm S}$ (with $r_{\rm S}$ being the Schwarzschild radius). We confirm that uncollimated outflows with velocities of 10 percent of the speed of light emerge from the innermost part of the accretion flow at a wide angle of 10$^{\circ}$ –50$^{\circ}$ from the disk rotation axis. Importantly, the outflows exhibit clumpy structures above heights of $\sim 250 r_{\rm S}$. The typical size of the clumps is $\sim 10 r_{\rm 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 the (upward) radiation is superior in force to the (downward) gravity, the Rayleigh–Taylor instability seems to be a primary cause. In addition, a radiation-hydrodynamic instability, which arises when radiation funnels through a radiation-pressure-supported atmosphere, may also help to form clumps of one optical depth. A magnetic photon bubble instability does not seem to be essential, since a similar clumpy outflow structure is obtained in nonmagnetic radiation-hydrodynamic simulations. Since the spatial covering factor of the clumps is estimated to be $\sim$ 0.3, and since they are marginally optically thick, they will explain at least some of the rapid light variations of active galactic nuclei. We further discuss a possibility of producing broad-line region clouds by the clumpy outflow.

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.

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.

Photon trapping and outflow are two key physics associated with the supercritical accretion flow. We investigate the conflict between these two processes based on twodimensional 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. The observational implications of our findings are briefly discussed in relation to ultra-luminous X-ray sources.

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 timeaveraged 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 distributions of azimuthal velocity f v , density ρ and pressure p basically follow the self-similar assumptions, i.e., they are roughly proportional to 2 , and ∼94% of the mass inflow is driven away as outflow, while outward momentum and energy fluxes are focused around the polar axis. Parts of these fluxes lie in the region that is not calculated by the steady model, and special attention should be paid when the model is applied.