The star formation in molecular clouds is inefficient. The ionizing EUV radiation ($h \nu \geq 13.6$ eV) from young clusters has been considered as a primary feedback effect to limit the star formation efficiency (SFE). We here focus on effects of the stellar FUV radiation (6 eV $\leq h \nu \leq$ 13.6 eV) during the cloud disruption stage. The FUV radiation may further reduce the SFE via photoelectric heating, and it also affects the chemical states of the gas that is not converted to stars ("cloud remnants") via photodissociation of molecules. We have developed a one-dimensional semi-analytic model which follows the evolution of both the thermal and chemical structure of a photodissociation region (PDR) during the dynamical expansion of an HII region. We investigate how the FUV feedback limits the SFE, supposing that the star formation is quenched in the PDR where the temperature is above a threshold value (e.g., 100K). Our model predicts that the FUV feedback contributes to reduce the SFEs for the massive ($M_{\rm cl} \gtrsim 10^5 M_{\odot}$) clouds with the low surface densities ($\Sigma_{\rm cl} \lesssim 100$ M$_{\odot}$pc$^{-2}$). Moreover, we show that a large part of the H$_2$ molecular gas contained in the cloud remnants should be "CO-dark" under the FUV feedback for a wide range of cloud properties. Therefore, the dispersed molecular clouds are potential factories of the CO-dark gas, which returns into the cycle of the interstellar medium.

We perform two-dimensional radiation hydrodynamical simulations of accretion flows onto black holes (BHs) at the nuclei of protogalaxies, and study the impact of mechanical and radiative feedback on rapid growth of BHs. The outflows deposit mass, momentum and energy into the surrounding medium and prevent mass accretion onto the BH, resulting in the reduction of radiative output. We find that when the BH is embedded in a dense gas core, ionizing radiation attenuated by inefficient BH feeding owing to mechanical feedback hardly affects the gas dynamics at the BH gravitational sphere of influence, from which intense inflows of neutral gas occur at rates substantially exceeding the Eddington limit without impeded by photoionization and heating. Since mechanical power of outflows driven by the rapidly accreting BH is sufficiently strong, bipolar outflows completely evacuate the surrounding gas in the polar region but mass inflows through the equatorial region maintain the BH accretion rate as high as $\sim 300-10^3~\dot{M}_{\rm Edd}$, which is reduced by one order of magnitude from those with radiative feedback alone. Furthermore, we find that the critical gas density required for rapid accretion is lower by a factor of $\sim 3$ nearly independently of BH mass, when mechanical feedback is considered. By studying the dependence on outflow model parameters (e.g., opening angle, mass-loading degree into outflows, velocity), we conclude that contrary to naive expectation, the existence of stronger outflow leads to the transition to rapid accretion phases more efficiently. Rapidly growing BHs inject mechanical power with $\sim 0.1-1\%$ of the radiative luminosity into their host galaxy scales, which is used for 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.

IW And stars are a subgroup of dwarf novae characterized by repetitive light variations of the intermediate-brightness state with oscillations, which is terminated by brightening. This group of dwarf novae is also known to exhibit a wide variety even within one system in long-term light curves including usual dwarf-nova outbursts, Z Cam-type standstills, and so on, besides the typical IW And-type variations mentioned above. Following the recent observations suggesting that some IW And stars seem to have tilted disks, we have investigated how the thermal-viscous instability works in tilted accretion disks in dwarf novae and whether it could reproduce the essential features of the light curves in IW And stars. By adopting various simplifying assumptions for tilted disks, we have performed time-dependent one-dimensional numerical simulations of a viscous disk by taking into account various mass supply patterns to the disk; that is, the gas stream from the secondary star flows not only to the outer edge of the disk but also to the inner portions of the disk. We find that tilted disks can achieve a new kind of accretion cycle, in which the inner disk almost always stays in the hot state while the outer disk repeats outbursts, thereby reproducing alternating mid-brightness interval sometimes with dips and brightening, which are quite reminiscent of the most characteristic observational light variations of IW And stars. Further, we have found that our simulations produce diverse light variations, depending on different mass supply patterns even without time variations in mass transfer rates. This could explain the wide variety in long-term light curves of IW And stars.

The black hole spacetime is described by general relativity and characterized by two quantities: the black hole mass and spin. Black hole spin measurement requires information from the vicinity of the event horizon, which is spatially resolved for the Galactic center SagittariusA* (SgrA*) and nearby radio galaxy M87 by means of very long baseline interferometry (VLBI) observations with the Event Horizon Telescope (EHT). In this paper, we simulate EHT observations for a gas cloud intermittently falling onto a black hole, and construct a method for spin measurement based on its relativistic flux variation. The light curve of the infalling gas cloud is composed of peaks formed by photons which directly reach a distant observer and by secondary ones reaching the observer after more than one rotation around the black hole. The time interval between the peaks is determined by a period of photon rotation near the photon circular orbit which uniquely depends on the spin. We perform synthetic EHT observations for SgrA* under a more realistic situation that a number of gas clouds intermittently fall towards the black hole with various initial parameters. Even for this case, the black hole spin dependence is detectable in correlated flux densities which are accurately calibrated by baselines between sites with redundant stations. The synthetic observations indicate that our methodology can be applied to EHT observations of Sgr A* since April 2017.

The cataclysmic variable MV Lyr was present in the Kepler field yielding a light curve with the duration of almost 1500 days with 60 second cadence. Such high quality data of this nova-like system with obvious fast optical variability show multicomponent power density spectra. Our goal is to study the light curve from different point of view, and perform a shot profile analysis. We search for characteristics not discovered with standard power density spectrum based methods. The shot profile method identifies individual shots in the light curve, and averages them in order to get all substructures with typical time scales. We also tested the robustness of our analysis using simple shot noise model. We obtained mean profiles with multicomponent features. The shot profile method distinguishes substructures with similar time scales which appear as a single degenerate feature in power density spectra. Furthermore, this method yields the identification of another high frequency component in the power density spectra of Kepler and XMM-Newton data not detected so far. Moreover, we found side-lobes accompanied with the central spike, making the profile very similar to another Kepler data of blazar W2R 1926+42, and Ginga data of Cyg X-1. All three objects show similar time scale ratios of the rising vs. declining part of the central spikes, while the two binaries have also similar rising profiles of the shots described by a power-law function. The similarity of both binary shot profiles suggests that the shots originate from the same origin, e.g. aperiodic mass accretion in the accretion disc. Moreover, the similarity with the blazar may imply that the ejection fluctuations in the blazar jet are connected to accretion fluctuations driving the variability in binaries. This points out to connection between jet and the accretion disc.

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.

We present the results from the Hitomi Soft Gamma-ray Detector (SGD) observation of the Crab nebula. The main part of SGD is a Compton camera, which in addition to being a spectrometer, is capable of measuring polarization of gamma-ray photons. The Crab nebula is one of the brightest X-ray / gamma-ray sources on the sky, and, the only source from which polarized X-ray photons have been detected. SGD observed the Crab nebula during the initial test observation phase of Hitomi. We performed the data analysis of the SGD observation, the SGD background estimation and the SGD Monte Carlo simulations, and, successfully detected polarized gamma-ray emission from the Crab nebula with only about 5 ks exposure time. The obtained polarization fraction of the phase-integrated Crab emission (sum of pulsar and nebula emissions) is (22.1 $\pm$ 10.6)% and, the polarization angle is 110.7$^o$ + 13.2 / $-$13.0$^o$ in the energy range of 60--160 keV (The errors correspond to the 1 sigma deviation). The confidence level of the polarization detection was 99.3%. The polarization angle measured by SGD is about one sigma deviation with the projected spin axis of the pulsar, 124.0$^o$ $\pm$0.1$^o$.

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.

Many Seyfert galaxies are known to exhibit Fe-K broad emission line features in their X-ray energy spectra. The observed lines have three distinct features; (1) the line profiles are skewed and show significant low-energy tails, (2) the Fe-K band have low variability, which produces a broad and deep dip in the root-mean-square (rms) spectra, and (3) photons in this band have time lags behind those in the adjacent energy bands with amplitudes of several $R_g/c$, where $R_g$ is the gravitational radius. The "relativistic light bending model" is proposed to explain these observed features, where a compact X-ray source (lamp post) above an extreme Kerr black hole illuminates the innermost area of the accretion disc. In this paper, we critically examine the relativistic light bending model by computing the rms spectra and the lag features using a ray-tracing technique, when a lamp post moves vertically on the black hole spin axis. As a result, we found that the observed deep rms dip requires that the iron is extremely overabundant ($\gtrsim10$ solar), whereas the observed lag amplitude is consistent with the normal iron abundance. Furthermore, disappearance of the lag in the high-flux state requires a source height as high as $\sim40\,R_g$, which contradicts the relativistically broad emission line feature. Our simulations agree with the data that the reverberation feature moves to lower frequencies with larger source height, however, if this scenario is correct, the simulations predict detection of a clear Fe-K lag at low frequencies, which is not constrained in the data. Therefore, we conclude that the relativistic light bending model may not explain the characteristic Fe-K spectral variations in Seyfert galaxies.

We present results from the Hitomi X-ray observation of a young composite-type supernova remnant (SNR) G21.5$-$0.9, whose emission is dominated by the pulsar wind nebula (PWN) contribution. The X-ray spectra in the 0.8-80 keV range obtained with the Soft X-ray Spectrometer (SXS), Soft X-ray Imager (SXI) and Hard X-ray Imager (HXI) show a significant break in the continuum as previously found with the NuSTAR observation. After taking into account all known emissions from the SNR other than the PWN itself, we find that the Hitomi spectra can be fitted with a broken power law with photon indices of $\Gamma_1=1.74\pm0.02$ and $\Gamma_2=2.14\pm0.01$ below and above the break at $7.1\pm0.3$ keV, which is significantly lower than the NuSTAR result ($\sim9.0$ keV). The spectral break cannot be reproduced by time-dependent particle injection one-zone spectral energy distribution models, which strongly indicates that a more complex emission model is needed, as suggested by recent theoretical models. We also search for narrow emission or absorption lines with the SXS, and perform a timing analysis of PSR J1833$-$1034 with the HXI and SGD. No significant pulsation is found from the pulsar. However, unexpectedly, narrow absorption line features are detected in the SXS data at 4.2345 keV and 9.296 keV with a significance of 3.65 $\sigma$. While the origin of these features is not understood, their mere detection opens up a new field of research and was only possible with the high resolution, sensitivity and ability to measure extended sources provided by an X-ray microcalorimeter.

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.

The present paper investigates the temperature structure of the X-ray emitting plasma in the core of the Perseus cluster using the 1.8--20.0 keV data obtained with the Soft X-ray Spectrometer (SXS) onboard the Hitomi Observatory. A series of four observations were carried out, with a total effective exposure time of 338 ks and covering a central region $\sim7'$ in diameter. The SXS was operated with an energy resolution of $\sim$5 eV (full width at half maximum) at 5.9 keV. Not only fine structures of K-shell lines in He-like ions but also transitions from higher principal quantum numbers are clearly resolved from Si through Fe. This enables us to perform temperature diagnostics using the line ratios of Si, S, Ar, Ca, and Fe, and to provide the first direct measurement of the excitation temperature and ionization temperature in the Perseus cluster. The observed spectrum is roughly reproduced by a single temperature thermal plasma model in collisional ionization equilibrium, but detailed line ratio diagnostics reveal slight deviations from this approximation. In particular, the data exhibit an apparent trend of increasing ionization temperature with increasing atomic mass, as well as small differences between the ionization and excitation temperatures for Fe, the only element for which both temperatures can be measured. The best-fit two-temperature models suggest a combination of 3 and 5 keV gas, which is consistent with the idea that the observed small deviations from a single temperature approximation are due to the effects of projection of the known radial temperature gradient in the cluster core along the line of sight. Comparison with the Chandra/ACIS and the XMM-Newton/RGS results on the other hand suggests that additional lower-temperature components are present in the ICM but not detectable by Hitomi SXS given its 1.8--20 keV energy band.

The Hitomi SXS spectrum of the Perseus cluster, with $\sim$5 eV resolution in the 2-9 keV band, offers an unprecedented benchmark of the atomic modeling and database for hot collisional plasmas. It reveals both successes and challenges of the current atomic codes. The latest versions of AtomDB/APEC (3.0.8), SPEX (3.03.00), and CHIANTI (8.0) all provide reasonable fits to the broad-band spectrum, and are in close agreement on best-fit temperature, emission measure, and abundances of a few elements such as Ni. For the Fe abundance, the APEC and SPEX measurements differ by 16%, which is 17 times higher than the statistical uncertainty. This is mostly attributed to the differences in adopted collisional excitation and dielectronic recombination rates of the strongest emission lines. We further investigate and compare the sensitivity of the derived physical parameters to the astrophysical source modeling and instrumental effects. The Hitomi results show that an accurate atomic code is as important as the astrophysical modeling and instrumental calibration aspects. Substantial updates of atomic databases and targeted laboratory measurements are needed to get the current codes ready for the data from the next Hitomi-level mission.

We present Hitomi observations of N132D, a young, X-ray bright, O-rich core-collapse supernova remnant in the Large Magellanic Cloud (LMC). Despite a very short observation of only 3.7 ks, the Soft X-ray Spectrometer (SXS) easily detects the line complexes of highly ionized S K and Fe K with 16-17 counts in each. The Fe feature is measured for the first time at high spectral resolution. Based on the plausible assumption that the Fe K emission is dominated by He-like ions, we find that the material responsible for this Fe emission is highly redshifted at ~800 km/s compared to the local LMC interstellar medium (ISM), with a 90% credible interval of 50-1500 km/s if a weakly informative prior is placed on possible line broadening. This indicates (1) that the Fe emission arises from the supernova ejecta, and (2) that these ejecta are highly asymmetric, since no blue-shifted component is found. The S K velocity is consistent with the local LMC ISM, and is likely from swept-up ISM material. These results are consistent with spatial mapping that shows the He-like Fe concentrated in the interior of the remnant and the S tracing the outer shell. The results also show that even with a very small number of counts, direct velocity measurements from Doppler-shifted lines detected in extended objects like supernova remnants are now possible. Thanks to the very low SXS background of ~1 event per spectral resolution element per 100 ks, such results are obtainable during short pointed or slew observations with similar instruments. This highlights the power of high-spectral-resolution imaging observations, and demonstrates the new window that has been opened with Hitomi and will be greatly widened with future missions such as the X-ray Astronomy Recovery Mission (XARM) and Athena.

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.

We report a Hitomi observation of IGR J16318-4848, a high-mass X-ray binary system with an extremely strong absorption of N_H~10^{24} cm^{-2}. Previous X-ray studies revealed that its spectrum is dominated by strong fluorescence lines of Fe as well as continuum emission. For physical and geometrical insight into the nature of the reprocessing material, we utilize the high spectroscopic resolving power of the X-ray microcalorimeter (the soft X-ray spectrometer; SXS) and the wide-band sensitivity by the soft and hard X-ray imager (SXI and HXI) aboard Hitomi. Even though photon counts are limited due to unintended off-axis pointing, the SXS spectrum resolves Fe K{\alpha_1} and K{\alpha_2} lines and puts strong constraints on the line centroid and width. The line width corresponds to the velocity of 160^{+300}_{-70} km s^{-1}. This represents the most accurate, and smallest, width measurement of this line made so far from any X-ray binary, much less than the Doppler broadening and shift expected from speeds which are characteristic of similar systems. Combined with the K-shell edge energy measured by the SXI and HXI spectra, the ionization state of Fe is estimated to be in the range of Fe I--IV. Considering the estimated ionization parameter and the distance between the X-ray source and the absorber, the density and thickness of the materials are estimated. The extraordinarily strong absorption and the absence of a Compton shoulder component is confirmed. These characteristics suggest reprocessing materials which are distributed in a narrow solid angle or scattering primarily with warm free electrons or neutral hydrogen.

The origin of the narrow Fe-K{\alpha} fluorescence line at 6.4 keV from active galactic nuclei has long been under debate; some of the possible sites are the outer accretion disk, the broad line region, a molecular torus, or interstellar/intracluster media. In February-March 2016, we performed the first X-ray microcalorimeter spectroscopy with the Soft X-ray Spectrometer (SXS) onboard the Hitomi satellite of the Fanaroff-Riley type I radio galaxy NGC 1275 at the center of the Perseus cluster of galaxies. With the high energy resolution of ~5 eV at 6 keV achieved by Hitomi/SXS, we detected the Fe-K{\alpha} line with ~5.4 {\sigma} significance. The velocity width is constrained to be 500-1600 km s$^{-1}$ (FWHM for Gaussian models) at 90% confidence. The SXS also constrains the continuum level from the NGC 1275 nucleus up to ~20 keV, giving an equivalent width ~20 eV of the 6.4 keV line. Because the velocity width is narrower than that of broad H{\alpha} line of ~2750 km s$^{-1}$, we can exclude a large contribution to the line flux from the accretion disk and the broad line region. Furthermore, we performed pixel map analyses on the Hitomi/SXS data and image analyses on the Chandra archival data, and revealed that the Fe-K{\alpha} line comes from a region within ~1.6 kpc from the NGC 1275 core, where an active galactic nucleus emission dominates, rather than that from intracluster media. Therefore, we suggest that the source of the Fe-K{\alpha} line from NGC 1275 is likely a low-covering fraction molecular torus or a rotating molecular disk which probably extends from a pc to hundreds pc scale in the active galactic nucleus system.

Extending the earlier measurements reported in Hitomi collaboration (2016, Nature, 535, 117), we examine the atmospheric gas motions within the central 100~kpc of the Perseus cluster using observations obtained with the Hitomi satellite. After correcting for the point spread function of the telescope and using optically thin emission lines, we find that the line-of-sight velocity dispersion of the hot gas is remarkably low and mostly uniform. The velocity dispersion reaches maxima of approximately 200~km~s$^{-1}$ toward the central active galactic nucleus (AGN) and toward the AGN inflated north-western `ghost' bubble. Elsewhere within the observed region, the velocity dispersion appears constant around 100~km~s$^{-1}$. We also detect a velocity gradient with a 100~km~s$^{-1}$ amplitude across the cluster core, consistent with large-scale sloshing of the core gas. If the observed gas motions are isotropic, the kinetic pressure support is less than 10\% of the thermal pressure support in the cluster core. The well-resolved optically thin emission lines have Gaussian shapes, indicating that the turbulent driving scale is likely below 100~kpc, which is consistent with the size of the AGN jet inflated bubbles. We also report the first measurement of the ion temperature in the intracluster medium, which we find to be consistent with the electron temperature. In addition, we present a new measurement of the redshift to the brightest cluster galaxy NGC~1275.

We investigate dynamical and radiative properties of arc-shaped gas clouds falling onto a stellar mass black hole based on the three-dimensional general relativistic radiation-magnetohydrodynamics (3D-GRRMHD) simulation data. Assuming that the gas clouds radiate due mainly to the free-free emission and/or optically thick, inverse Compton scattering, we calculate how the emissivity distributions develop with time. We find that (1) gas clouds, each of which has a ring-like or arc shape, are intermittently formed, and that (2) they slowly fall to the black hole, keeping nearly the Keplerian orbital velocity. These features support the dynamical properties of the gas clouds assumed in the spin measurement method proposed by Moriyama & Mineshige 2015, but the radius of the inner edge of the accretion disk is larger than that of the marginally stable orbit (ISCO). Next, we examine how each gas cloud is observed by a distant observer by calculating the photon trajectories in the black hole space-time. The luminosity of the accretion flow exhibits significant time variations on different timescales, reflecting the time evolution of the gas density distributions. The relatively slow variations on the time durations of 0.08-0.10 sec is due to the formation and fall of gas clouds, while quasi-periodic flux peaks with short time intervals (0.01 sec) are due to the quasi-periodic enhancement of light from the non-axisymmetric arc-shaped clouds through the beaming effect. This may account for the high-frequency quasi-periodic oscillations (HF QPOs) observed in black hole binaries. The observational implications and future issues are briefly discussed.

Thanks to its high spectral resolution (~5 eV at 6 keV), the Soft X-ray Spectrometer (SXS) on board Hitomi enables us to measure the detailed structure of spatially resolved emission lines from highly ionized ions in galaxy clusters for the first time. In this series of papers, using the SXS we have measured the velocities of gas motions, metallicities and the multi-temperature structure of the gas in the core of the Perseus cluster. Here, we show that when inferring physical properties from line emissivities in systems like Perseus, the resonant scattering (RS) effect should be taken into account. In the Hitomi waveband, RS mostly affects the FeXXV He$\alpha$ line ($w$) - the strongest line in the spectrum. The flux measured by Hitomi in this line is suppressed by a factor ~1.3 in the inner ~30 kpc, compared to predictions for an optically thin plasma; the suppression decreases with the distance from the center. The $w$ line also appears slightly broader than other lines from the same ion. The observed distortions of the $w$ line flux, shape and distance dependence are all consistent with the expected effect of the resonant scattering in the Perseus core. By measuring the ratio of fluxes in optically thick ($w$) and thin (FeXXV forbidden, He$\beta$, Ly$\alpha$) lines, and comparing these ratios with predictions from Monte Carlo radiative transfer simulations, the velocities of gas motions have been obtained. The results are consistent with the direct measurements of gas velocities from line broadening described elsewhere in this series, although the systematic and statistical uncertainties remain significant. Further improvements in the predictions of line emissivities in plasma models, and deeper observations with future X-ray missions will enable RS measurements to provide powerful constraints on the amplitude and anisotropy of clusters gas motions.

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

To search for giant X-ray pulses correlated with the giant radio pulses (GRPs) from the Crab pulsar, we performed a simultaneous observation of the Crab pulsar with the X-ray satellite Hitomi in the 2 -- 300 keV band and the Kashima NICT radio observatory in the 1.4 -- 1.7 GHz band with a net exposure of about 2 ks on 25 March 2016, just before the loss of the Hitomi mission.The timing performance of the Hitomi instruments was confirmed to meet the timing requirement and about 1,000 and 100 GRPs were simultaneously observed at the main and inter-pulse phases, respectively, and we found no apparent correlation between the giant radio pulses and the X-ray emission in either the main or inter-pulse phases.All variations are within the 2 sigma fluctuations of the X-ray fluxes at the pulse peaks, and the 3 sigma upper limits of variations of main- or inter- pulse GRPs are 22\% or 80\% of the peak flux in a 0.20 phase width, respectively, in the 2 -- 300 keV band.The values become 25\% or 110\% for main or inter-pulse GRPs, respectively, when the phase width is restricted into the 0.03 phase.Among the upper limits from the Hitomi satellite, those in the 4.5-10 keV and the 70-300 keV are obtained for the first time, and those in other bands are consistent with previous reports.Numerically, the upper limits of main- and inter-pulse GRPs in the 0.20 phase width are about (2.4 and 9.3) $\times 10^{-11}$ erg cm$^{-2}$, respectively. No significant variability in pulse profiles implies that the GRPs originated from a local place within the magnetosphere and the number of photon-emitting particles temporally increases.However, the results do not statistically rule out variations correlated with the GRPs, because the possible X-ray enhancement may appear due to a $>0.02$\% brightening of the pulse-peak flux under such conditions.

The Crab nebula originated from a core-collapse supernova (SN) explosion observed in 1054 A.D. When viewed as a supernova remnant (SNR), it has an anomalously low observed ejecta mass and kinetic energy for an Fe-core collapse SN. Intensive searches were made for a massive shell that solves this discrepancy, but none has been detected. An alternative idea is that the SN1054 is an electron-capture (EC) explosion with a lower explosion energy by an order of magnitude than Fe-core collapse SNe. In the X-rays, imaging searches were performed for the plasma emission from the shell in the Crab outskirts to set a stringent upper limit to the X-ray emitting mass. However, the extreme brightness of the source hampers access to its vicinity. We thus employed spectroscopic technique using the X-ray micro-calorimeter onboard the Hitomi satellite. By exploiting its superb energy resolution, we set an upper limit for emission or absorption features from yet undetected thermal plasma in the 2-12 keV range. We also re-evaluated the existing Chandra and XMM-Newton data. By assembling these results, a new upper limit was obtained for the X-ray plasma mass of <~ 1Mo for a wide range of assumed shell radius, size, and plasma temperature both in and out of the collisional equilibrium. To compare with the observation, we further performed hydrodynamic simulations of the Crab SNR for two SN models (Fe-core versus EC) under two SN environments (uniform ISM versus progenitor wind). We found that the observed mass limit can be compatible with both SN models if the SN environment has a low density of <~ 0.03 cm-3 (Fe core) or <~ 0.1 cm-3 (EC) for the uniform density, or a progenitor wind density somewhat less than that provided by a mass loss rate of 10-5 Mo yr-1 at 20 km s-1 for the wind environment.

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}$.

High-resolution X-ray spectroscopy with Hitomi was expected to resolve the origin of the faint unidentified E=3.5 keV emission line reported in several low-resolution studies of various massive systems, such as galaxies and clusters, including the Perseus cluster. We have analyzed the Hitomi first-light observation of the Perseus cluster. The emission line expected for Perseus based on the XMM-Newton signal from the large cluster sample under the dark matter decay scenario is too faint to be detectable in the Hitomi data. However, the previously reported 3.5 keV flux from Perseus was anomalously high compared to the sample-based prediction. We find no unidentified line at the reported high flux level. Taking into account the XMM measurement uncertainties for this region, the inconsistency with Hitomi is at a 99% significance for a broad dark-matter line and at 99.7% for a narrow line from the gas. We do not find anomalously high fluxes of the nearby faint K line or the Ar satellite line that were proposed as explanations for the earlier 3.5 keV detections. We do find a hint of a broad excess near the energies of high-n transitions of Sxvi (E=3.44 keV rest-frame) -- a possible signature of charge exchange in the molecular nebula and another proposed explanation for the unidentified line. While its energy is consistent with XMM pn detections, it is unlikely to explain the MOS signal. A confirmation of this interesting feature has to wait for a more sensitive observation with a future calorimeter experiment.

We present spectral and timing analyses of simultaneous X-ray and UV observations of the VY Scl system MV Lyr taken by XMM-Newton, containing the longest continuous X-ray+UV light curve and highest signal-to-noise X-ray (EPIC) spectrum to date. The RGS spectrum displays emission lines plus continuum, confirming model approaches to be based on thermal plasma models. We test the sandwiched model based on fast variability that predicts a geometrically thick corona that surrounds an inner geometrically thin disc. The EPIC spectra are consistent with either a cooling flow model or a 2-T collisional plasma plus Fe emission lines in which the hotter component may be partially absorbed which would then originate in a central corona or a partially obscured boundary layer, respectively. The cooling flow model yields a lower mass accretion rate than expected during the bright state, suggesting an evaporated plasma with a low density, thus consistent with a corona. Timing analysis confirms the presence of a dominant break frequency around log(f/Hz) = -3 in the X-ray Power Density Spectrum (PDS) as in the optical PDS. The complex soft/hard X-ray light curve behaviour is consistent with a region close to the white dwarf where the hot component is generated. The soft component can be connected to an extended region. We find another break frequency around log(f/Hz) = -3.4 that is also detected by Kepler. We compared flares at different wavelengths and found that the peaks are simultaneous but the rise to maximum is delayed in X-rays with respect to UV.

We analyze the ${\it Kepler}$ monitoring light curve of a blazar W2R 1926$+$42 to examine features of microvariability by means of the "shot analysis" technique. We select 195 intra-day, flare-like variations (shots) for the continuous light curve of Quarter 14 with a duration of 100 d. In the application of the shot analysis, an averaged profile of variations is assumed to converge with a universal profile which reflects a physical mechanism generating the microvariability in a blazar jet, although light-variation profiles of selected shots show a variety. A mean profile, which is obtained by aligning the peaks of the 195 shots, is composed of a spiky-shape shot component at $\pm$0.1 d (with respect to the time of the peak), and two slow varying components ranging from $-$0.50 d to $-$0.15 d and from 0.10 d to 0.45 d of the peak time. The former spiky feature is well represented by an exponential rise of 0.043$\pm$0.001 d and an exponential decay of 0.061$\pm$0.002 d. These timescales are consistent with that corresponding to a break frequency of a power spectrum density calculated from the obtained light curve. After verification with the Monte-Carlo method, the exponential shape, but not the observed asymmetry, of the shot component can be explained by noise variation. The asymmetry is difficult to explain through a geometrical effect (i.e. changes of the geometry of the emitting region), but is more likely to be caused by the production and dissipation of high-energy accelerated particles in the jet. Additionally, durations of the detected shots show a systematic variation with a dispersion caused by a statistical randomness. A comparison with the variability of Cygnus X-1 is also briefly discussed.

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.

Clusters of galaxies are the most massive gravitationally-bound objects in the Universe and are still forming. They are thus important probes of cosmological parameters and a host of astrophysical processes. Knowledge of the dynamics of the pervasive hot gas, which dominates in mass over stars in a cluster, is a crucial missing ingredient. It can enable new insights into mechanical energy injection by the central supermassive black hole and the use of hydrostatic equilibrium for the determination of cluster masses. X-rays from the core of the Perseus cluster are emitted by the 50 million K diffuse hot plasma filling its gravitational potential well. The Active Galactic Nucleus of the central galaxy NGC1275 is pumping jetted energy into the surrounding intracluster medium, creating buoyant bubbles filled with relativistic plasma. These likely induce motions in the intracluster medium and heat the inner gas preventing runaway radiative cooling; a process known as Active Galactic Nucleus Feedback. Here we report on Hitomi X-ray observations of the Perseus cluster core, which reveal a remarkably quiescent atmosphere where the gas has a line-of-sight velocity dispersion of 164+/-10 km/s in a region 30-60 kpc from the central nucleus. A gradient in the line-of-sight velocity of 150+/-70 km/s is found across the 60 kpc image of the cluster core. Turbulent pressure support in the gas is 4% or less of the thermodynamic pressure, with large scale shear at most doubling that estimate. We infer that total cluster masses determined from hydrostatic equilibrium in the central regions need little correction for turbulent pressure.

We propose a new method for black hole spin measurement. In this method, we consider a gas blob or ring falling onto a black hole from the marginally stable orbit, keeping its initial orbital angular momentum. We calculate the gas motion and photon trajectories in the Kerr space-time and, assuming that the gas blob or ring emits monochromatic radiation, carefully examine how it is observed by a distant observer. The light curve of the orbiting gas blob is composed of many peaks because of periodic enhancement of the flux due to the gravitational lensing and beaming effects. Further, the intensity of each peak first gradually increases with time due to the focusing effect around the photon circular orbit and then rapidly decreases due to the gravitational redshift, as the gas blob approaches the event horizon. The light curve of the gas ring is equivalent to a superposition of those of the blobs with various initial orbital phases, and so it is continuous and with no peaks. The flux first gradually increases and then rapidly decays, as in the blob model. The flux variation timescale depends on the black hole spin and is independent from the inclination angle, while time averaged frequency shift have dependences of both effects. We can thus, in principle, determine spin and inclination angle from observations. The observational implications and future issues are briefly discussed.

Hyperaccretion flows with mass accretion rates far above the Eddington rate have an N-shaped equilibrium curve on the $\Sigma$-$\dot{M}$ plane (with $\Sigma$ and $\dot{M}$ being surface density and mass accretion rate, respectively). The accretion flow on the lower $\Sigma$ branch of the N-shape is optically thick, advection-dominated accretion flow (ADAF) while that on the upper one is neutrino-dominated accretion flow (NDAF). The middle branch has a negative slope on the $\Sigma$-$\dot{M}$ plane, meaning that the flow on this branch is secularly unstable. To investigate how the instability affects the flow structure and what observable signatures are produced, we study the time evolution of the unstable hyperaccretion flow in response to variable mass injection rates by solving the height-averaged equations for viscous accretion flows. When a transition occurs from the lower branch to the upper branch (or from the upper branch to the lower branch), the surface density rapidly increases (decreases) around that transition region, which induces locally enhanced mass flow (referred to as non-steady mass flow) into (out of) that region. We confirm that the non-steady flow can create a kind of disturbance and that it propagates over the whole disk. However, the non-steady mass flow is not strong enough to induce coherent transition over the whole disk, unless the mass injection rate varies with time. When the injection rate continuously changes, the neutrino luminosity varies intermittently, thus producing step-function-like light curves, as the radiation efficiency discontinuously changes every time the local transition occurs. The effects of changing the N-shape and possible observational consequences on the gamma-ray burst emission are briefly discussed in relation to gamma-ray bursts.

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.

Blazars show rapid and violent variabilities, which timescale are often less than a day. We studied intraday variations by applying a shot analysis technique to Kepler monitoring of blazar W2R 1926+42 in Quarter 14. We obtained a mean profile calculated from 195 rapid variations. The mean profile shows three components; one is a sharp structure distributed within $\pm$0.1 day of the peak, and two slow-varying components. This spiky-peak component reflects features of rapid variations directly. The profile of peak component shows an exponential rise and decay of which timescales are different, 0.0416 and 0.0588 day respectively. This component is too sharp to represent a standard function which is often used to express blazar variations. This asymmetric profile at the peak is difficult to be explained by a simple variation of the Doppler factor by changing a geometry of the emitting region. This result indicates that intraday variations arise from a production of high-energy accelerated particles in the jet.

The stochastic variability (flickering) of the nova-like system (subclass of cataclysmic variable) MV Lyr yields a complicated power density spectrum with four break frequencies. Scaringi et al. (2012) analysed high-cadence Kepler data of MV Lyr, taken almost continuously over 600 days, giving the unique opportunity to study multicomponent Power Density Spectra (PDS) over a wide frequency range. We modelled this variability with our statistical model based on disc angular momentum transport via discrete turbulent bodies with an exponential distribution of the dimension scale. Two different models were used, a full disc (developed from the white dwarf to the outer radius of ~10^10 cm) and a radially thin disc (a ring at a distance of ~10^10 cm from the white dwarf) that imitates an outer disc rim. We succeed in explaining the two lowest observed break frequencies assuming typical values for a disc radius of 0.5 and 0.9 times the primary Roche lobe and an alpha parameter of 0.1-0.4. The highest observed break frequency was also modelled, but with a rather small accretion disc with radius of 0.3 times the primary Roche lobe and a high alpha value of 0.9 consistent with previous findings by Scaringi (2014). Furthermore, the simulated light curves exhibit the typical linear rms-flux proportionality linear relation and the typical log-normal flux distribution. As the turbulent process is generating fluctuations in mass accretion that propagate through the disc, this confirms the general knowledge that the typical rms-flux relation is mainly generated by these fluctuations. In general a higher rms is generated by a larger amount of superposed flares which is compatible with a higher mass accretion rate expressed by a larger flux.

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions developed by the Institute of Space and Astronautical Science (ISAS), with a planned launch in 2015. The ASTRO-H mission is equipped with a suite of sensitive instruments with the highest energy resolution ever achieved at E > 3 keV and a wide energy range spanning four decades in energy from soft X-rays to gamma-rays. The simultaneous broad band pass, coupled with the high spectral resolution of Delta E < 7 eV of the micro-calorimeter, will enable a wide variety of important science themes to be pursued. ASTRO-H is expected to provide breakthrough results in scientific areas as diverse as the large-scale structure of the Universe and its evolution, the behavior of matter in the gravitational strong field regime, the physical conditions in sites of cosmic-ray acceleration, and the distribution of dark matter in galaxy clusters at different redshifts.

There is still 10-20% uncertainty on the neutron star (NS) mass-radius relation. These uncertainties could be reduced by an order of magnitude through an unambiguous measure of M/R from the surface redshift of a narrow line, greatly constraining the Equation of State for ultra-dense material. It is possible that the SXS on ASTRO-H can detect this from an accreting neutron star with low surface velocity in the line of sight i.e. either low inclination or low spin. Currently there is only one known low inclination LMXB, Ser X-1, and one known slow spin LMXB, J17480-2446 in Terzan 5. Ser X-1 is a persistent source which is always in the soft state (banana branch), where the accreting material should form a equatorial belt around the neutron star. A pole-on view should then allow the NS surface to be seen directly. A 100 ks observation should allow us to measure M/R if there are any heavy elements in the photosphere at the poles. Conversely, J17480-2446 in Terzan 5 is a transient accretion powered millisecond pulsar, where the accreting material is collimated onto the magnetic pole in the hard (island) state (L_x < 0.1 L_Edd). The hotspot where the shock illuminates the NS surface is clearly seen in this state. A 100 ks ToO observation of this (or any other similarly slow spin system) in this state, may again allow the surface redshift to be directly measured. (abstract continues)

Thanks to extensive observations with X-ray missions and facilities working in other wavelengths, as well as rapidly--advancing numerical simulations of accretion flows, our knowledge of astrophysical black holes has been remarkably enriched. Rapid progress has opened new areas of enquiry, including measurements of black hole spin, the properties and driving mechanisms of jets and disk winds, the impact of feedback into local environments, the origin of periodic and aperiodic X-ray variations, and the nature of super-Eddington accretion flows, among others. The goal of this White Paper is to illustrate how ASTRO-H can make dramatic progress in the study of astrophysical black holes, particularly the study of black hole X-ray binaries.

Shortly after the seminal paper {\sl "Self-Organized Criticality: An explanation of 1/f noise"} by Bak, Tang, and Wiesenfeld (1987), the idea has been applied to solar physics, in {\sl "Avalanches and the Distribution of Solar Flares"} by Lu and Hamilton (1991). In the following years, an inspiring cross-fertilization from complexity theory to solar and astrophysics took place, where the SOC concept was initially applied to solar flares, stellar flares, and magnetospheric substorms, and later extended to the radiation belt, the heliosphere, lunar craters, the asteroid belt, the Saturn ring, pulsar glitches, soft X-ray repeaters, blazars, black-hole objects, cosmic rays, and boson clouds. The application of SOC concepts has been performed by numerical cellular automaton simulations, by analytical calculations of statistical (powerlaw-like) distributions based on physical scaling laws, and by observational tests of theoretically predicted size distributions and waiting time distributions. Attempts have been undertaken to import physical models into the numerical SOC toy models, such as the discretization of magneto-hydrodynamics (MHD) processes. The novel applications stimulated also vigorous debates about the discrimination between SOC models, SOC-like, and non-SOC processes, such as phase transitions, turbulence, random-walk diffusion, percolation, branching processes, network theory, chaos theory, fractality, multi-scale, and other complexity phenomena. We review SOC studies from the last 25 years and highlight new trends, open questions, and future challenges, as discussed during two recent ISSI workshops on this theme.

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 analysed an X-ray light curve of the dwarf nova RU Peg taken by XMM-Newton with a duration of 46300 s. The power density spectrum has a complicated shape with two red noise and two white noise components, indicating the presence of two turbulent regions. We developed a statistical "toy model" to study light curves with variability produced by an unstable turbulent accretion flow from the inner disc. Our results are consistent with a disc truncation radius maximally 0.8 x 10^9 cm. We found that any fluctuation of the viscous mass accretion at the inner disc are visible as UV and X-ray variations with the same break frequency in the power density spectrum. This process is generating low frequency variability. A second break suggests the presence of a faster X-ray variability component which must be generated by another process likely localised between the inner disc and the white dwarf.

A hyperaccretion flow around a stellar mass black hole is thought to be the most plausible engine that powers gamma-ray bursts (GRBs). The flow efficiently cools via neutrino emission at >~ 0.003-0.01 M_sun s^{-1} (corresponding to a luminosity of ~10^50 erg s^{-1}), while neither neutrino nor photon emission is efficient below this rate, so the flow should be advection dominated. We carefully solve how a transition occurs from the advection-dominated to the neutrino-dominated branches, and find that the slope of the thermal equilibrium curve is negative in the surface density - accretion rate (Sigma-Mdot) plane, a condition for viscous instability, at radii smaller than ~12R_g (with R_g being the gravitational radius). We also confirm that the flow is thermally stable. The consequence of this instability is the formation of a clumpy structure in the flow. This is because the larger (respectively smaller) surface density is, the smaller (respectively larger) the mass accretion rate from the region in question becomes, leading to growth of the density contrast. The timescale for clump formation is estimated to be shorter than 0.1s. The observational implication is discussed in the context of GRBs. We suggest that this might explain the origin of the large variability observed in the prompt emission of GRBs.

We study the observational signatures of the relativistic slim disc of 10 M_sun black hole, in a wide range of mass accretion rate, mdot, dimensionless spin parameter, a_ast, and viewing angle, i. In general, the innermost temperature, T_in increases with the increase of i for a fixed value of mdot and a_ast, due to the Doppler effect. However, for i > 50 and mdot > mdot_turn, T_in starts to decrease with the increase of mdot. This is a result of self-obscuration -- the radiation from the innermost hot part of the disc is blocked by the surrounding cooler part. The value of mdot_turn and the corresponding luminosities depend on a_ast and i. Such obscuration effects cause an interesting behavior on the disc luminosity (L_disc) -- T_in plane for high inclinations. In addition to the standard-disc branch which appears below mdot_turn and which obeys L_disc propto T_in^4 -relation, another branch above mdot_turn, which is nearly horizontal, may be observed at luminosities close to the Eddington luminosity. We show that these features are likely observed in a Galactic X-ray source, GRS 1915+105. We support a high spin parameter (a_ast > 0.9) for GRS 1915+105 since otherwise the high value of T_in and small size of the emitting region (r_in < 1r_S) cannot be explained.

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.

Rapid spectral changes in the hard X-ray on a time scale down to ~0.1 s are studied by applying "shot analysis" technique to the Suzaku observations of the black hole binary Cygnus X-1, performed on 2008 April 18 during the low/hard state. We successfully obtained the shot profiles covering 10--200 keV with the Suzaku HXD-PIN and HXD-GSO detector. It is notable that the 100-200 keV shot profile is acquired for the first time owing to the HXD-GSO detector. The intensity changes in a time-symmetric way, though the hardness does in a time-asymmetric way. When the shot-phase-resolved spectra are quantified with the Compton model, the Compton y-parameter and the electron temperature are found to decrease gradually through the rising phase of the shot, while the optical depth appears to increase. All the parameters return to their time-averaged values immediately within 0.1 s past the shot peak. We have not only confirmed this feature previously found in energies below ~60 keV, but also found that the spectral change is more prominent in energies above ~100 keV, implying the existence of some instant mechanism for direct entropy production. We discuss possible interpretations on the rapid spectral changes in the hard X-ray band.

We present observations of a transient He-like Fe K alpha absorption line in Suzaku observations of the black hole binary Cygnus X-1 on 2011 October 5 near superior conjunction during the high/soft state, which enable us to map the full evolution from the start and the end of the episodic accretion phenomena or dips for the first time. We model the X-ray spectra during the event and trace their evolution. The absorption line is rather weak in the first half of the observation, but instantly deepens for ~10 ks, and weakens thereafter. The overall change in equivalent width is a factor of ~3, peaking at an orbital phase of ~0.08. This is evidence that the companion stellar wind feeding the black hole is clumpy. By analyzing the line with a Voigt profile, it is found to be consistent with a slightly redshifted Fe XXV transition, or possibly a mixture of several species less ionized than Fe XXV. The data may be explained by a clump located at a distance of ~10^(10-12) cm with a density of ~10^((-13)-(-11)) g cm^-3, which accretes onto and/or transits the line-of-sight to the black hole, causing an instant decrease in the observed degree of the ionization and/or an increase in density of the accreting matter. Continued monitoring for individual events with future X-ray calorimeter missions such as ASTRO-H and AXSIO will allow us to map out the accretion environment in detail and how it changes between the various accretion states.

In order to elucidate the emission properties of ultraluminous X-ray sources (ULXs) during their power-law (PL) state, we examined long-term X-ray spectral data of IC342 X-1 during its PL state by using our own Suzaku data and the archival data by XMM-Newton, Chandra, and Swift observations. The PL state of this source seems to be classified into two sub-states in terms of the X-ray luminosities in 0.5-10 keV: the low luminosity PL state with 4-6*10^{39} erg/s and the high luminosity one with 1.1-1.4*10^{40} erg/s. During the Suzaku observations which were made in 2010 August and 2011 March, X-1 stayed in the low luminosity PL state. The observed X-ray luminosity (4.9-5.6*10^{39} erg/s) and the spectral shape (photon index = 1.67-1.83) slightly changed between the two observations. Using the Suzaku PIN detector, we for the first time confirmed a PL tail extending up to at least 20 keV with no signatures of a high-energy turnover in both of the Suzaku observations. In contrast, a turnover at about 6 keV was observed during the high luminosity PL state in 2004 and 2005 with XMM-Newton. Importantly, photon indices are similar between the two PL states and so is the Compton y-parameters of y ~ 1, which indicates a similar energy balance (between the corona and the accretion disk) holding in the two PL states despite different electron temperatures. From spectral similarities with recent studies about other ULXs and the Galactic black hole binary GRS1915+105, IC342 X-1 is also likely to be in a state with a supercritical accretion rate, although more sensitive higher energy observations would be necessary to conclude.

We study the basic properties of accretion flows onto binary supermassive black holes, including the cases in which a circumbinary disk is misaligned with the binary orbital plane, by means of three-dimensional Smoothed Particle Hydrodynamics simulations. We find that a circular binary system with a misaligned circumbinary disk normally produces a double peaked mass-accretion-rate variation per binary orbit. This is because each black hole passes across the circumbinary disk plane and captures gas twice in one orbital period. Even in misaligned systems, however, a single peaked mass-accretion-rate variation per binary orbit is produced, if the orbital eccentricity is moderately large (e\lesssim0.3). The number of peaks in mass accretion rates can be understood simply in terms of the orbital phase dependence of the distance between each binary black hole and its closest inner edge of the circumbinary disk. In the cases of eccentric binary black holes having different masses, the less massive black hole can get closer to the circumbinary disk than the massive one, thus tidally splitting gas from its inner edge, but the created gas flows are comparably captured by both black holes with a short time delay. As a consequence, the combined light curve shows periodic occurrence of double-peaked flares with a short interval. This may account for the observed light variations of OJ287.

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions initiated by the Institute of Space and Astronautical Science (ISAS). ASTRO-H will investigate the physics of the high-energy universe via a suite of four instruments, covering a very wide energy range, from 0.3 keV to 600 keV. These instruments include a high-resolution, high-throughput spectrometer sensitive over 0.3-2 keV with high spectral resolution of Delta E < 7 eV, enabled by a micro-calorimeter array located in the focal plane of thin-foil X-ray optics; hard X-ray imaging spectrometers covering 5-80 keV, located in the focal plane of multilayer-coated, focusing hard X-ray mirrors; a wide-field imaging spectrometer sensitive over 0.4-12 keV, with an X-ray CCD camera in the focal plane of a soft X-ray telescope; and a non-focusing Compton-camera type soft gamma-ray detector, sensitive in the 40-600 keV band. The simultaneous broad bandpass, coupled with high spectral resolution, will enable the pursuit of a wide variety of important science themes.

We investigate the accretion process in high-luminosity AGNs (HLAGNs) in the scenario of the disk evaporation model. Based on this model, the thin disk can extend down to the innermost stable circular orbit (ISCO) at accretion rates higher than $0.02\dot{M}_{\rm Edd}$; while the corona is weak since part of the coronal gas is cooled by strong inverse Compton scattering of the disk photons. This implies that the corona cannot produce as strong X-ray radiation as observed in HLAGNs with large Eddington ratio. In addition to the viscous heating, other heating to the corona is necessary to interpret HLAGN. In this paper, we assume that a part of accretion energy released in the disk is transported into the corona, heating up the electrons and thereby radiated away. We for the first time, compute the corona structure with additional heating, taking fully into account the mass supply to the corona and find that the corona could indeed survive at higher accretion rates and its radiation power increases. The spectra composed of bremsstrahlung and Compton radiation are also calculated. Our calculations show that the Compton dominated spectrum becomes harder with the increase of energy fraction ($f$) liberating in the corona, and the photon index for hard X-ray($2-10 \rm keV$) is $2.2 < \Gamma < 2.7 $. We discuss possible heating mechanisms for the corona. Combining the energy fraction transported to the corona with the accretion rate by magnetic heating, we find that the hard X-ray spectrum becomes steeper at larger accretion rate and the bolometric correction factor ($L_{\rm bol}/L_{\rm 2-10keV}$) increases with increasing accretion rate for $f<8/35$, which is roughly consistent with the observational results.

When binary black holes are embedded in a gaseous environment, a rotating disk surrounding them, the so-called circumbinary disk, will be formed. The binary exerts a gravitational torque on the circumbinary disk and thereby the orbital angular momentum is transferred to it, while the angular momentum of the circumbinary disk is transferred to the binary through the mass accretion. The binary undergoes an orbital decay due to both the gravitational wave emission and the binary-disk interaction. This causes the phase evolution of the gravitational wave signal. The precise measurement of the gravitational wave phase thus may provide information regarding the circumbinary disk. In this paper, we assess the detectability of the signature of the binary-disk interaction using the future space-borne gravitational wave detectors such as DECIGO and BBO by the standard matched filtering analysis. We find that the effect of the circumbinary disk around binary black holes in the mass range $6M_sun\le{M}\lesssim3\times10^3M_sun$ is detectable at a statistically significant level in five year observation, provided that gas accretes onto the binary at a rate greater than $\dot{M}\sim1.4\times10^{17} [gs^{-1}] j^{-1}(M/10M_sun)^{33/23}$ with 10% mass-to-energy conversion efficiency, where j represents the efficiency of the angular momentum transfer from the binary to the circumbinary disk. We show that $O(0.1)$ coalescence events are expected to occur in sufficiently dense molecular clouds in five year observation. We also point out that the circumbinary disk is detectable, even if its mass at around the inner edge is by over 10 orders of magnitude less than the binary mass.

We present a study of the flickering activity in two nova like systems KR Aur and UU Aqr. We applied a statistical model of flickering simulations in accretion discs based on turbulent angular momentum transport between two adjacent rings with an exponential distribution of the turbulence dimension scale. The model is based on a steady state disc model which is satisfied in the case of hot ionized discs of nova like cataclysmic variables. Our model successfully fits the observed power density spectrum of KR Aur with the disc parameter alpha = 0.10 - 0.40 and an inner disc truncation radius in the range R_in = 0.88 - 1.67 x 10^9 cm. The exact values depend on the mass transfer rate in the sense that alpha decreases and R_in increases with mass transfer rate. In any case, the inner disc radius found for KR Aur is considerably smaller than in quiescent dwarf novae, as predicted by the disc instability model. On the other hand, our simulations fail to reproduce the power density spectrum of UU Aqr. A tantalizing explanation involves the possible presence of spiral waves, which are expected in UU Aqr, because of its low mass ratio, but not in KR Aur. In general our model predicts the observed concentration of flickering in the central disc. We explain this by the radial dependence of the angular momentum gradient.

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.

In order to examine if the next generation space VLBI,such as VSOP-2 (VLBI Space Observatory Programme-2), will make it possible to obtain direct images of the accretion flow around the M87 black hole, we calculate the expected observed images by the relativistic ray-tracing simulations under the considerations of possible observational errors. We consider various cases of electron temperature profiles, as well as a variety of the distance, mass, and spin of the M87 black hole. We find it feasible to detect an asymmetric intensity profile around the black hole caused by rapid disk rotation, as long as the electron temperature does not steeply rises towards the black hole, as was predicted by the accretion disk theory and the three dimensional magnetohydrodynamic simulations. Further, we can detect a deficit in the observed intensity around the black hole when the apparent size of the gravitational radius is larger than $\gtrsim 1.5 \mu$ arcseconds. In the cases that the inner edge of the disk is located at the radius of the innermost stable circular orbit (ISCO), moreover, even the black hole spin will be measured. We also estimate the required signal-to-noise ratio $\mathcal{R}_{\rm SN}$ for achieving the scientific goals mentioned above, finding that it should be at least 10 at 22 GHz. To conclude, direct mapping observations by the next generation space VLBI will provide us a unique opportunity to provide the best evidence for the presence of a black hole and to test the accretion disk theory.

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions initiated by the Institute of Space and Astronautical Science (ISAS). ASTRO-H will investigate the physics of the high-energy universe by performing high-resolution, high-throughput spectroscopy with moderate angular resolution. ASTRO-H covers very wide energy range from 0.3 keV to 600 keV. ASTRO-H allows a combination of wide band X-ray spectroscopy (5-80 keV) provided by multilayer coating, focusing hard X-ray mirrors and hard X-ray imaging detectors, and high energy-resolution soft X-ray spectroscopy (0.3-12 keV) provided by thin-foil X-ray optics and a micro-calorimeter array. The mission will also carry an X-ray CCD camera as a focal plane detector for a soft X-ray telescope (0.4-12 keV) and a non-focusing soft gamma-ray detector (40-600 keV) . The micro-calorimeter system is developed by an international collaboration led by ISAS/JAXA and NASA. The simultaneous broad bandpass, coupled with high spectral resolution of Delta E ~7 eV provided by the micro-calorimeter will enable a wide variety of important science themes to be pursued.

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.

To understand the origin of hard X-ray emissions from black hole X-ray binaries during their low/hard states, we calculate the X-ray spectra of black-hole accretion flow for the following three configurations of hot and cool media: (a) an inner hot flow and a cool outer disk (inner hot flow model), (b) a cool disk sandwiched by disk coronae (disk corona model), and (c) the combination of those two (hybrid model). The basic features we require for successful models are (i) significant hard X-ray emission whose luminosity exceeds that of soft X-rays, (ii) high hard X-ray luminosities in the range of (0.4 - 30) times 10^{37} erg s^{-1}, and (iii) the existence of two power-law components in the hard X-ray band with the photon indices of Gamma_s ~ 2 > Gamma_h, where Gamma_s and Gamma_h are the photon indices of the softer (<10 keV) and the harder (>10 keV) power-law components, respectively. Contribution by non-thermal electrons nor time-dependent evolution are not considered. We find that Models (a) and (b) can be ruled out, since the spectra are always dominated by the soft component, and since only one power-law component, at most, can be reproduced. Only Model (c) can account for sufficiently strong hard X-ray emissions, as well as the existence of the two power-law components, for a large ratio of the accretion rate in the corona to that in the thin disk. The outer disk corona (where the Compton y-parameter is smaller, y < 1) produces the softer power-law component with photon index of Gamma_s ~ 2, whereas the inner hot flow (where y gtrsim 1) generates the harder component with Gamma_h < 2. This model can also account for the observed relationship between the photon index and the reflection fraction.

To understand the nature of supercritical accretion, we systematically analyze the {\it RXTE}/PCA data of GRS 1915+105 in its quasi-steady states, by choosing data with small variability during 1999 -- 2000. We apply a multicolor disk plus a thermal Comptonization model and take into consideration accurate interstellar absorption, a reflection component, and absorption features from the disk wind self-consistently. There is a strong correlation between the inner disk temperature and the fraction of the disk component. Most of the Comptonization-dominated spectra show $T_{\rm in} \sim 1$ keV with a high electron temperature of $>10$ keV, which may correspond to the very high state in canonical black hole X-ray binaries (BHBs). By contrast, the disk-dominated spectra have $T_{\rm in} \sim 2$ keV with a low temperature ($<10$ keV) and optically thick Comptonization, and show two separate branches in the ($L$ -- $T_{\rm in}$) diagram. The lower branch clearly follows the $L \propto T_{\rm in}^4$-track. Furthermore, applying the extended disk blackbody model, we find that 9 out of 12 datasets with disk luminosity above $0.3 L_{\rm E}$ prefer a flatter temperature gradient than that in the standard disk ($p < 0.7$). We interpret that, in the lower branch, the disk extends down to the innermost stable circular orbit, and is most probably in the slim disk state. A rapidly spinning black hole can explain both the lack of the $L \propto T_{\rm in}^2$-track and a high value of spectral hardening factor ($\sim 4$) that would be required for a non-rotating black hole. The spectra in the upper branch are consistent with the picture of a truncated disk with low temperature Comptonization. This state is uniquely observed from GRS 1915+105 among BHBs, which may be present at near-Eddington luminosity.

We use XMM-Newton and Swift data to study spectral variability in the ultraluminous X-ray source (ULX), Holmberg IX X-1. The source luminosity varies by a factor 3-4, giving rise to corresponding spectral changes which are significant, but subtle, and not well tracked by a simple hardness ratio. Instead, we co-add the Swift data in intensity bins and do full spectral fitting with disc plus thermal Comptonisation models. All the data are well-fitted by a low temperature, optically thick Comptonising corona, and the variability can be roughly characterised by decreasing temperature and increasing optical depth as the source becomes brighter, as expected if the corona is becoming progressively mass loaded by material blown off the super-Eddington inner disc. This variability behaviour is seen in other ULX which have similar spectra, but is opposite to the trend seen in ULX with much softer spectra. This supports the idea that there are two distinct physical regimes in ULXs, where the spectra go from being dominated by a disc-corona to being dominated by a wind.

We study mass flow rate through a disc resulting from a varying mass supply rate. Variable mass supply rate occurs, e.g., during disc state transitions, and in interacting eccentric binaries. It is, however, damped by the viscosity of the disc. Here, we calculate this damping in detail. We derive an analytical description of the propagation of the flow rate using the solution of Lynden-Bell & Pringle, in which the disc is assumed to extend to infinity. In particular, we derive the accretion-rate Green's function, and its Fourier transform, which gives the fractional damping at a given variability frequency. We then compare this model to that of a finite disc with the mass supply at its outer edge. We find significant differences with respect to the infinite disc solution, which we find to overestimate the viscous damping. In particular, the asymptotic form of the Green's function is power-law for the infinite disc and exponential for the finite one. We then find a simple fitting form for the latter, and also calculate its Fourier transform. In general, the damping becomes very strong when the viscous time at the outer edge of the disc becomes longer than the modulation time scale. We apply our results to a number of astrophysical systems. We find the effect is much stronger in low-mass X-ray binaries, where the disc size is comparable to that of the Roche lobe, than in high-mass binaries, where the wind-fed disc can have a much smaller size.

Radiatively inefficient accretion flows (RIAFs) are common feature of low-luminosity accretion flows, including quiescent states of X-ray binaries and low-lunimosity active galactic nuclei. Thermally driven winds are expected from such hot accretion flows. By assuming that the flow has self-similarity structure in the radial direction, we solve the vertical structure of the wind and accretion flows simultaneously and evaluate the mass loss rates by wind. We find that the ratio of the outflow rate to the accretion rate is approximately unity for a viscosity parameter, alpha lesssim 0.1, despite some uncertainties in the angular momentum and temperature distributions. That is, the accretion rate in the RIAFs is roughly proportional to the radius. Moreover, we elucidate the effect of cooling by wind on the underneath accretion flow, finding that this effect could be important for calculating energy spectrum of the RIAF. Observational implications are briefly discussed in the context of Sgr A*.

The source of viscosity in astrophysical accretion flows is still a hotly debated issue. We investigate the contribution of convective turbulence to the total viscosity in a self-consistent approach, where the strength of convection is determined from the vertical disc structure itself. Additional sources of viscosity are parametrized by a beta-viscosity prescription, which also allows an investigation of self-gravitating effects. In the context of accretion discs around stellar mass and intermediate mass black holes, we conclude that convection alone cannot account for the total viscosity in the disc, but significantly adds to it. For accretion rates up to 10% of the Eddington rate, we find that differential rotation provides a sufficiently large underlying viscosity. For higher accretion rates, further support is needed in the inner disc region, which can be provided by an MRI-induced viscosity. We briefly discuss the interplay of MRI, convection and differential rotation. We conduct a detailed parameter study of the effects of central masses and accretion rates on the disc models and find that the threshold value of the supporting viscosity is determined mostly by the Eddington ratio with only little influence from the central black hole mass.

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 propose a methodology to derive a black-hole mass for super-critical accretion flow. Here, we use the extended disk blackbody (extended DBB) model, a fitting model in which the effective temperature profile obeys the relation $T_{\rm eff} \propto r^{-p}$, with $r$ being the disk radius and $p$ being treated as a fitting parameter. We first numerically calculate the theoretical flow structure and its spectra for a given black-hole mass, $M$, and accretion rate, $\dot{M}$. Through fitting to the theoretical spectra by the extended DBB model, we can estimate the black-hole mass, $M_{\rm x}$, assuming that the innermost disk radius is $r_{\rm in}=3r_{\rm g} (\propto M_{\rm x})$, where $r_{\rm g}$ is the Schwarzschild radius. We find, however, that the estimated mass deviates from that adopted in the spectral calculations, $M$, even for low-$\dot{M}$ cases. We also find that the deviations can be eliminated by introducing a new correction for the innermost radius. Using this correction, we calculate mass correction factors, $M/M_{\rm x}$, in the super-critical regimes for some sets of $M$ and $\dot M$, finding that a mass correction factor ranges between $M/M_{\rm x} \sim$ 1.2 -- 1.6. The higher is $\dot{M}$, the larger does the mass correction factor tend to be. Since the correction is relatively small, we can safely conclude that the black holes in ULXs which Vierdayanti et al. (2006, PASJ, 58, 915) analyzed are stellar-mass black holes with the mass being $< 100M_{\odot}$.

We have investigated the effect of initial mass segregation on the runaway merging of stars. The evolution of multi-mass, dense star clusters was followed by means of direct N-body simulations of up to 131.072 stars. All clusters started from King models with dimensionless central potentials of 3.0 <= W_0 <= 9.0. Initial mass segregation was realized by varying the minimum mass of a certain fraction of stars whose either (1) distances were closest to the cluster center or (2) total energies were lowest. The second case is more favorable to promote the runaway merging of stars by creating a high-mass core of massive, low-energy stars. Initial mass segregation could decrease the central relaxation time and thus help the formation of a high-mass core. However, we found that initial mass segregation does not help the runaway stellar merger to happen if the overall mass density profile is kept constant. This is due to the fact that the collision rate of stars is not increased due to initial mass segregation. Our simulations show that initial mass segregation is not sufficient to allow runaway merging of stars to occur in clusters with central densities typical for star clusters in the Milky Way.

Hierarchical structure formation inevitably leads to the formation of supermassive binary black holes (BBHs) with a sub-parsec separation in galactic nuclei. However, to date there has been no unambiguous detection of such systems. In an effort to search for potential observational signatures of supermassive BBHs, we performed high-resolution smoothed particle hydrodynamics (SPH) simulations of two black holes in a binary of moderate eccentricity surrounded by a circumbinary disk. Building on our previous work, which has shown that gas can periodically transfer from the circumbinary disk to the black holes when the binary is on an eccentric orbit, the current set of simulations focuses on the formation of the individual accretion disks, their evolution and mutual interaction, and the predicted radiative signature. The variation in mass transfer with orbital phase from the circumbinary disk induces periodic variations in the light curve of the two accretion disks at ultraviolet wavelengths, but not in the optical or near-infrared. Searches for this signal offer a promising method to detect supermassive BBHs.

We calculate the radiation spectrum and its time variability of the black hole accretion disk-corona system based on the three-dimensional magnetohydrodynamic simulation. In explaining the spectral properties of active galactic nuclei (AGNs), it is often assumed that they consist of a geometrically thin, optically thick disk and hot, optically thin corona surrounding the thin disk. As for a model of the corona, we adopt the simulation data of three-dimensional, non-radiative MHD accretion flows calculated by Kato and coworkers, while for a thin disk we assume a standard type disk. We perform Monte Carlo radiative transfer simulations in the corona, taking into account the Compton scattering of soft photons from the thin disk by hot thermal electrons and coronal irradiation heating of the thin disk, which emits blackbody radiation. By adjusting the density parameter of the MHD coronal flow, we can produce the emergent spectra which are consistent with those of typical Seyfert galaxies. Moreover, we find rapid time variability in X-ray emission spectra, originating from the density fluctuation produced by the magnetorotational instability in the MHD corona. The features of reflection component including iron fluorescent line emission are also briefly discussed.

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.

Hypercritical accretion flows onto stellar mass black holes (BHs) are commonly considered as a promising model of central engines of gamma-ray bursts (GRBs). In this model a certain fraction of gravitational binding energy of accreting matter is deposited to the energy of relativistic jets via neutrino annihilation and/or magnetic fields. However, some recent studies have indicated that the energy deposition rate by neutrino annihilation is somewhat smaller than that needed to power a GRB. To overcome this difficulty, Ramirez-Ruiz & Socrates (2005) proposed that high energy neutrinos from hot corona above the accretion disk might enhance the efficiency of energy deposition. We elucidate the disk corona model in the context of hypercritical accretion flows. From the energy balance in the disk and the corona, we can calculate the disk and coronal temperature, Td and Tc, and neutrino spectra, taking into account the neutrino cooling processes by neutrino-electron scatterings and neutrino pair productions. The calculated neutrino spectra consist of two peaks; one by the neutrino emission from the disk and the other by that from the corona. We find that the disk corona can enhance the efficiency of energy release but only by a factor of 1.5 or so, unless the height of the corona is very small, H<<r. This is because the neutrino emission is very sensitive to the temperature of emitting region, and then the ratio Tc/Td cannot be so large.

Black holes and neutron stars present extreme forms of matter that cannot be created as such in a laboratory on Earth. Instead, we have to observe and analyze the experiments that are ongoing in the Universe. The most telling observations of black holes and neutron stars come from dense stellar systems, where stars are crowded close enough to each other to undergo frequent interactions. It is the interplay between black holes, neutron stars and other objects in a dense environment that allows us to use observations to draw firm conclusions about the properties of these extreme forms of matter, through comparisons with simulations. The art of modeling dense stellar systems through computer simulations forms the main topic of this review.

(Simplified) We present a theoretical framework which establishes how the core radius of a star cluster varies with the mass of an assumed central black hole. Our result is that the ratio of core to half-mass radius varies as the 3/4 power of the ratio of the black hole to cluster mass, when the system is well relaxed. The theory compares favourably with a number of simulations of this problem, which extend to black hole masses of order 10% of the cluster mass. Though strictly limited as yet to clusters with stars of equal mass, our conclusion strengthens the view that clusters with large core radii are the most promising candidates in which to find a massive black hole.

We investigate the structure and stability of hypercritical accretion flows around stellar-mass black holes, taking into account neutrino cooling, lepton conservation, and firstly a realistic equation of state in order to properly treat the dissociation of nuclei. We obtain the radial distributions of physical properties, such as density, temperature and electron fraction, for various mass accretion rates $0.1\sim 10M_{\odot}{\rm s}^{-1}$. We find that, depending on mass accretion rates, different physics affect considerably the structure of the disk; most important physics is (1) the photodissociation of nuclei around $r\sim 100r_g$ for relatively low mass accretion rates ($\dot{M}\sim 0.01-0.1M_{\odot} {\rm s}^{-1}$), (2) efficient neutrino cooling around $r\sim 10-100r_g$ for moderately high mass accretion rate ($\dot{M}\sim 0.2-1.0M_{\odot}{\rm s}^{-1}$), and (3) neutrino trapping ($r\sim 3-10r_g$) for very high mass accretion rate ($\dot{M}\gtrsim 2.0M_{\odot}{\rm s}^{-1}$). We also investigate the stability of hypercritical accretion flows by drawing the thermal equilibrium curves, and find that efficient neutrino cooling makes the accretion flows rather stable against both thermal and viscous modes.

We investigate the accretion flows onto the supermassive binary black holes (SMBBHs) from the circumbinary disk with the equal mass, eccentric binary on the subparsec scale, using Smoothed Particle Hydrodynamics (SPH) code. We find that the material can be supplied from circumbinary disk, which leads to the formation of two accretion disks around the SMBBHs. The mass accretion rates significantly modulate with the binary orbital motion. These could provide the observable diagnosis of the existence of the supermassive binary black holes (e.g. OJ287) on the subparsec scale in merged galactic nuclei.

We study the accretion flows from the circumbinary disks onto the supermassive binary black holes in a subparsec scale of the galactic center, using a smoothed particles hydrodynamics (SPH) code. Simulation models are presented in four cases of a circular binary with equal and unequal masses, and of an eccentric binary with equal and unequal masses. We find that the circumblack-hole disks are formed around each black holes regardless of simulation parameters. There are two-step mechanisms to cause an accretion flow from the circumbinary disk onto supermassive binary black holes: First, the tidally induced elongation of the circumbinary disk triggers mass inflow towards two closest points on the circumbinary disk from the black holes. Then, the gas is increasingly accumulated on these two points owing to the gravitational attraction of black holes. Second, when the gas can pass across the maximum loci of the effective binary potential, it starts to overflow via their two points and freely infalls to each black hole. In circular binaries, the gas continues to be supplied from the circumbinary disk (i.e. the gap between the circumbinary disk and the binary black hole is always closed.) In eccentric binaries, the mass supply undergoes the periodic on/off transitions during one orbital period because of the variation of periodic potential. The gap starts to close after the apastron and to open again after the next periastron passage. Due to this gap closing/opening cycles, the mass-capture rates are eventually strongly phase dependent. This could provide observable diagnosis for the presence of supermassive binary black holes in merged galactic nuclei.

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

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.

We present the first study of the dynamical evolution of an isolated star cluster that combines a significant population of primordial binaries with the presence of a central black hole. We use equal-mass direct N-body simulations, with N ranging from 4096 to 16384 and a primordial binary ratio of 0-10%; the black hole mass is about one percent of the total mass of the cluster. The evolution of the binary population is strongly influenced by the presence of the black hole, which gives the cluster a large core with a central density cusp. Starting from a variety of initial conditions (Plummer and King models), we first encounter a phase, that last approximately 10 half-mass relaxation times, in which binaries are disrupted faster compared to analogous simulations without a black hole. Subsequently, however, binary disruption slows down significantly, due to the large core size. The dynamical interplay between the primordial binaries and the black hole thus introduces new features with respect to the scenarios investigated so far, where the influence of the black hole and of the binaries have been considered separately. A large core to half mass radius ratio appears to be a promising indirect evidence for the presence of a intermediate-mass black hole in old globular clusters.

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 dynamics of an injected outflow propagating in a progenitor in the context of the collapsar model for gamma-ray bursts (GRBs) through two dimensional axisymmetric relativistic hydrodynamic simulations. Initially, we locally inject an outflow near the center of a progenitor. We calculate 25 models, in total, by fixing its total input energy to be 10^{51} ergs s^{-1} and radius of the injected outflow to be $7\times 10^7$ cm while varying its bulk Lorentz factor, $\Gamma_{0} = 1.05\sim 5$, and its specific internal energy, $\epsilon_0/c^2 = 0.1\sim 30$. The injected outflow propagates in the progenitor and drives a large-scale outflow or jet. We find a smooth but dramatic transition from a collimated jet to an expanding outflow among calculated models. The maximum Lorentz factor is, on the other hand, sensitive to both of $\Gamma_0$ and $\epsilon_0$; roughly $\Gamma_{\rm max} \sim \Gamma_0 (1+\epsilon_0/c^2)$. Our finding will explain a smooth transition between the GRBs, X-ray rich GRBs (XRRs) and X-ray Flashes (XRFs) by the same model but with different $\epsilon_0$ values.

We present the first study of the dynamical evolution of a star cluster that combines a significant population of primordial binaries with the presence of a central black hole. We use direct N-body simulations, with a black hole mass of about one percent of the total mass of the cluster. The evolution of the binary population is strongly influenced by the presence of the black hole, which gives the cluster a large core with a central density cusp. Starting from a Plummer profile, we first encounter a phase, that last approximately 10 half-mass relaxation times, in which binaries are disrupted faster compared to analogous simulations without a black hole. Subsequently, however, binary disruption slows down significantly, due to the large core size. The dynamical interplay between the primordial binaries and the black hole thus introduces new features with respect to the scenarios investigated so far, where the influence of the black hole and of the binaries have been considered separately. Specifically, the pattern of binary destruction by an intermediate-mass black hole may leave a fingerprint that could be detected observationally.

We investigate the iron K$\alpha$ fluorescent line produced by hard X-ray photons from magnetic reconnection-heated corona. The hot corona with temperature being about $10^9$K can irradiate the underlying disk with a continuum X-ray spectrum produced via thermal Comptonization. Then the iron atoms in the disk photoelectrically absorb X-ray photons and radiate K$\alpha$ line photons. Therefore, the activity of corona is responsible to the iron line emission from the underlying disk. In previous studies, oversimplified X-ray photon sources are often assumed above the disk in order to compute the iron line profile or power-law line emissivity profiles are assumed with an index being a free parameter. We adopt the more realistic corona model constructed by Liu et al. in which the corona is heated by magnetic energy released through the reconnection of magnetic flux loops and which has no free parameter. Then the accretion energy is dominantly dissipated in the corona, in which X-ray photons are efficiently produced and irradiate the underlying disk. We find the local emmisivity of iron line on the disk is approximated as $F_{{\rm K}\alpha}(r)\propto r^{-5}$. The iron line profiles derived from this model give excellent fits to the observational data of MCG-6-30-15 with the profiles derived theoretically for $i\sim 30^{\circ}$ for energy band 4-7keV. Possible origins of line variability are briefly discussed.

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 compile a sample consisting of 56 radio-quiet active galactic nuclei so as to investigate statistical properties of hot corona of accretion disks from {\em ASCA} observations. The black-hole masses in the sample are estimated via several popular methods and the bolometric luminosities from the multi-wavelength continuum. This allows us to estimate the Eddington ratio (${\cal E}\equiv L_{\rm Bol}/L_{\rm Edd}$) so that the undergoing physical processes can be tested via hard X-ray data. We find a strong correlation between ${\cal F}_{\rm X}\equiv L_{\rm 2-10keV}/L_{\rm Bol}$ and ${\cal E}$ as ${\cal F}_{\rm X}\propto {\cal E}^{-0.64}$ with a multivariate regression. This indicates that the release of gravitational energy in the hot corona is controlled by the Eddington ratio. On the other hand, the correlation between the hard X-ray spectral index ($\Gamma$) and ${\cal E}$ depends critically on the types of objects: $\Gamma$ is nearly constant ($\Gamma \propto {\cal E}^0$) in broad-line Seyfert 1's (BLS1s), whereas $\Gamma \propto \log {\cal E}^{0.18}$ in narrow-line Seyfert 1's (NLS1s), although not very significant. We can set constraints on the forms of magnetic stress tensor on the condition that ${\cal F}_{\rm X}$ is proportional to the fraction $f$ of gravitational energy dissipated in the hot corona and that $f$ is proportional to magnetic energy density in the disk. We find that the shear stress tensor $t_{r\phi}\propto P_{\rm gas}$ is favored by the correlation in the present sample, where $P_{\rm gas}$ is the gas pressure.

We present three-dimensional (3-D) magnetohydrodynamical (MHD) simulations of radiatively inefficient accretion flow around black holes. General relativistic effects are simulated by using the pseudo-Newtonian potential. We start calculations with a rotating torus threaded by localized poloidal magnetic fields with plasma beta, a ratio of the gas pressure to the magnetic pressure, $\beta =10$ and 100. When the bulk of torus material reaches the innermost region close to a central black hole, a magnetically driven jet emerges. This magnetic jet is derived by vertically inflating toroidal fields (`magnetic tower') and has a two-component structure: low-$\beta$ ($\lsim 1$) plasmas threaded with poloidal (vertical) fields are surrounded by that with toroidal fields. The collimation width of the jet depends on external pressure, pressure of ambient medium; the weaker the external pressure is, the wider and the longer-lasting becomes the jet. Unless the external pressure is negligible, the bipolar jet phase ceases after several dynamical timescales at the original torus position and a subsequent quasi-steady state starts. The black hole is surrounded by quasi-spherical zone with highly inhomogeneous structure in which toroidal fields are dominant except near the rotation axis. Mass accretion takes place mainly along the equatorial plane. Comparisons with other MHD simulation results and observational implications are discussed.

Impressed by the widespread misunderstanding of the issue, we return to the old question of the location of the inner edge of accretion disk around black hole. We recall the fundamental results obtained in the 1970's and 1980's by Warsaw and Kyoto research groups that proved, in particular, that the inner edge does not coincide with the location of the innermost stable Keplerian circular orbit. We give some novel illustrations of this particular point and of some other fundamental results obtained by Warsaw and Kyoto groups. To investigate the flow dynamics of the inner edge of accretion disk, we carefully solve the structure of the transonic flow and plot the effective potential profile based on the angular-momentum distribution calculated numerically. We show that the flow does not have a potential minimum for accretion rates, {\dot M} > 10 L_E/c^2 (with L_E being the Eddington luminosity and $c$ being the speed of light). This property is realized even in relatively small viscosity parameters (i.e., \alpha ~ 0.01), because of the effect of pressure gradient. In conclusion, the argument based on the last circular orbit of a test particle cannot give a correct inner boundary of the super-critical flow and the inner edge should be determined in connection with radiation efficiency. The same argument can apply to optically thin ADAF. The interpretation of the observed QPO frequencies should be re-considered, since the assumption of Kepler rotation velocity can grossly over- or underestimate the disk rotation velocity, depending on the magnitude of viscosity.

To understand the bursting behavior of the microquasar GRS 1915+105, we calculate time evolution of a luminous, optically thick accretion disk around a stellar mass black hole undergoing limit-cycle oscillations between the high- and low- luminosity states. We, especially, carefully solve the behavior of the innermost part of the disk, since it produces significant number of photons during the burst, and fit the theoretical spectra with the multi-color disk model. The fitting parameters are $\Tin$ (the maximum disk temperature) and $\Rin$ (the innermost radius of the disk). We find an abrupt, transient increase in $\Tin$ and a temporary decrease in $\Rin$ during a burst, which are actually observed in GRS 1915+105. The precise behavior is subject to the viscosity prescription. We prescribe the radial-azimuthal component of viscosity stress tensor to be $T_{r \phi}=-\alpha \Pi (p_{\rm gas}/p)^{\mu}$, with $\Pi$ being the height integrated pressure, $\alpha$ and $\mu$ being the parameter, and $p$ and $p_{\rm gas}$ being the total pressure and gas pressure on the equatorial plane, respectively. Model with $\mu=0.1$ can produce the overall time changes of $\Tin$ and $\Rin$, but cannot give an excellent fit to the observed amplitudes. Model with $\mu=0.2$, on the other hand, gives the right amplitudes, but the changes of $\Tin$ and $\Rin$ are smaller. Although precise matching is left as future work, we may conclude that the basic properties of the bursts of GRS 1915+105 can be explained by our ``limit-cycle oscillation'' model. It is then required that the spectral hardening factor at high luminosities should be about 3 at around the Eddington luminosity instead of less than 2 as is usually assumed.

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.

Recently, Chartas et al. (2001) detected a rapid X-ray flare in the gravitationally lensed, multiple image quasar RX J0911.4+0551. Dramatic events, such as rapid X-ray flares, are useful in providing high precision measurements of the time delays between multiple images. In this paper, we argue that there is a new possibility in measurements of time delays between multiple images of gravitationally lensed quasars; constrain the locations of putative flares that give rise to the intrinsic rapid variabilities of quasars. The realization, however, of these goals cannot be presently achieved due to the limited accuracy of the current measurements. We predict that timing flares with accuracies of the order of a few seconds will be needed to probe the location of the flares. Our proposing method will work with better instruments in near future, such as XEUS.

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.

The detectability of gravitational wave (GW) radiation from accretion disks is discussed based on various astrophysical contexts. In order to emit GW radiation, the disk shape should lose axial symmetry. We point out that a significant deformation is plausible in non-radiative hot accretion disks because of enhanced magnetic activity, whereas it is unlikely for standard-type cool disks. We have analyzed the 3D magnetohydrodynamical (MHD) simulation data of magnetized accretion flow, finding non-axisymmetric density patterns. The corresponding ellipticity is $\epsilon \sim 0.01$. The expected time variations of GW radiation are overall chaotic, but there is a hint of quasi-periodicity. GW radiation has no interesting consequence, however, in the case of close binaries, because of very tiny disk masses. GW radiation is not significant, either, for AGN because of very slow rotation velocities. The most promising case can be found in gamma-ray bursts or supernovae, in which a massive torus (or disk) with a solar mass or so may be formed around a stellar-mass compact object as the result of a merger of compact objects, or by the fallback of exploded material towards the center in a supernova. Although much more intense GW radiation is expected before the formation of the torus, the detection of GW radiation in the subsequent accretion phase is of great importance, since it will provide a good probe to investigating their central engines.

Prompted by recent discoveries of log-normal distributions in gamma-ray/X-ray temporal variabilities of gamma-ray bursts (GRBs) and blazars, we re-examine the X-ray variability of Cygnus X-1 in the hard/low state using Ginga data in 1990. It was previously reported that the distributions of the time intervals between X-ray shots (flares) deviated from the Poisson distributions at short time intervals: the occurrence of shots tended to be suppressed for several seconds before and/or after a shot event. Detailed analyses show that this deviation is larger for shots with larger peaks, and that the time-interval distributions for large shots approach log-normal distributions with a peak at the interval of 7--8 s. Furthermore, we also show that the peak-intensity distribution for the shots is consistent with the log-normal distribution, though no typical peak intensity can be seen. This might indicate the presence of a physical link connecting the physics of black hole accretion flow and that of jet/GRB formation.

It is often considered that a massive torus with solar mass or so surrounding a stellar-mass black hole may be a central engine of a gamma-ray burst. We study the properties of such massive accretion tori (or disks) based on the $\alpha$ viscosity model. For surface density exceeding about $10^{20}$ g cm$^{-2}$, which realizes when about a solar-mass material is contained within a disk with a size of $\sim 5 \times 10^6$ cm, we find that (1) luminosity of photons is practically zero due to significant photon trapping, (2) neutrino cooling dominates over advective cooling, (3) pressure of degenerate electrons dominates over pressure of gas and photons, and (4) magnetic field strength exceeds the critical value of about $4 \times 10^{13}$ G, even if we take 0.1 % of the equi-partition value. The possible observable quantum electrodynamical (QED) effects arising from super-critical fields are discussed. Most interestingly, photon splitting may occur, producing significant number of photons of energy below $\sim 511$ keV, thereby possibly suppressing e$^\pm$ pair creation.

We propose a new method, emission-line eclipse mapping, to map the velocity fields of an accretion disk in position space. Quiescent dwarf novae usually exhibit double-peaked emission line profiles because of disk rotation. Since a part of the disk having a different line-of-sight velocity is successively obscured by a companion in eclipsing systems, they show time-varying line profiles. We calculate time changes of the emission line profiles, assuming Keplerian rotation fields (v_\phi\propto r^(-1/2) with r being the distance from the disk center) and the emissivity distribution of j \propto r^(-3/2). We, then, apply the usual eclipse mapping technique to the light curves at each of 12-24 wavelengths across the line center to map the region with the same line-of-sight velocity. The reconstructed images typically exhibit a `two-eye' pattern for high line-of-sight velocities and we can recover the relation, v_\phi \propto d^(-1/2), on the assumption of the axisymmetric disk, where d is the separation between the two `eyes'. We will be able to probe the Keplerian rotation law, the most fundamental assumption adopted in many disk models, by high-speed spectroscopic observations with 8-m class telescopes.

We construct a simple model for a magnetic reconnection heated corona above a thin accretion disk in AGNs and Galactic black hole candidates(GBHCs). The basic assumptions are that (1) the magnetic reconnection heat is cooled down overwhelmingly by Compton scattering in the corona and that (2) thermal conduction is dominantly cooled by evaporation of the chromospheric plasma in the disk-corona interface before Compton cooling sets in. With these two basic equations as well as equi-partition of magnetic energy with gas energy in the disk, we can consistently determine the fraction of accretion energy dissipated in the corona without free parameters, and thus determine the temperature and all other quantities in both the corona and disk for given black hole mass and accretion rate. Then, we calculate the luminosity contributed from the disk and corona and the coronal flux weighted Compton $y$ parameter. It is found that, at a low luminosity (less than 0.1 Eddington luminosity) the spectrum is hard with energy spectral index of around 1, while at a high luminosity (larger than 0.1 Eddington luminosity) the spectrum can be either soft or hard. We also find that the situation is almost same for supermassive and stellar-mass black holes. These features are consistent with observations of AGNs and GBHCs.

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

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 examine the effects of the different viscosity prescriptions and the magnitude of the viscosity parameter, $\alpha$, on the structure of the slim disk, and discuss the observational implications on accretion-flow into a stellar-mass black hole. For the range of $\alpha = 10^{-2} \sim 10^{0}$ we calculate the disk spectra and from spectral fitting we derive $\Tin$, maximum temperature of the disk, $\Rin$, the size of the region emitting blackbody radiation with $\Tin$, and $p\equiv -{\rm dln} T_{\rm eff}/{\rm dln}~r$, the slope of the effective temperature distribution. It was founded that the estimated $\Tin$ slightly increases as $\alpha$ increases. This is because the larger the magnitude of viscosity is, the larger becomes the accretion velocity and, hence, the more enhanced becomes advective energy transport, which means less efficient radiative cooling and thus higher temperatures. Furthermore we check different viscosity prescriptions with the form of the viscous stress tensor of $t_{r \varphi} = -\alpha \beta^{\mu}p_{\rm total}$, where $\beta$ is the ratio of gas pressure to total pressure, and $\mu$ is a parameter ($0 \le \mu \le 1$). For $\mu=0$ we have previously found that as luminosity approaches the Eddington, $L_{\rm E$, $\Rin$ decreases below 3$\rg$ (with $r_{\rm g}$ being Schwarzschild radius) and the effective temperature becomes flatter, $T_{\rm eff} \propto r^{-1/2}$. Such a slim-disk nature does not appear when $\mu$ is large, $\mu \sim 0.5$, even at $L_{\rm E}$. Hence, the temperature of the innermost region of the disk sensitively depends on the $\mu$ value. We can rule out the case with large $\mu~(\sim 0.5)$, since it will not be able to produce a drop in $\Rin$ with an increase in luminosity as was observed in an ultraluminous X-ray source, IC~342, source 1.

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

We calculated the expected microlens light curves for several accretion-flow models and compared the results to find what aspects of distinct flow structures can be extracted from the microlens light curves. As for the flow models, we considered the standard disk, optically thin advection-dominated flow, and the disk-corona model by Kawaguchi et al., which can nicely reproduce the observed quasar spectra. The calculated microlens light curves of the disk-corona model exhibit rapid soft X-ray and slow hard X-ray variations. This is because the former is generated by Compton up-scattering of soft (optical -- UV) photons from the innermost part of the disk body, while the latter is emitted via bremsstrahlung within the corona at relatively large radii. Furthermore, an inhomogeneous emissivity distribution over the disk produces humps in the microlens light curves because of enhanced microlens amplification of the bright areas. It is hoped that black-hole accretion theory coupled with future multi-wavelength microlens observations will be able to reveal such spatial and temporal behaviors of flows in the vicinity of black holes.

The Ultra Luminous X-ray Sources (ULXs) are unique in exhibiting moderately bright X-ray luminosities, $L_{\rm x} \sim 10^{38-40} {\rm erg~s^{-1}}$, and relatively high blackbody temperatures, $\Tin \sim 1.0-2.0 {\rm keV}$. From the constraint that $L_{\rm x}$ cannot exceed the Eddington luminosity, $L_{\rm E}$, we require relatively high black-hole masses, $M\sim 10-100 M_\odot$, however, for such large masses the standard disk theory predicts lower blackbody temperatures, $\Tin < 1.0$ keV. To understand a cause of this puzzling fact, we carefully calculate the accretion flow structure shining at $\sim L_{\rm E}$, fully taking into account the advective energy transport in the optically thick regime and the transonic nature of the flow. Our calculations show that at high accretion rate ($\dot M \ga 30~L_{\rm E}/c^2$) an apparently compact region with a size of $\Rin \simeq (1-3)\rg$ (with $\rg$ being Schwarzschild radius) is shining with a blackbody temperature of $\Tin \simeq 1.8 (M/10M_\odot)^{-1/4}$ keV even for the case of a non-rotating black hole. Further, $\Rin$ decreases as $\dot M$ increases, on the contrary to the canonical belief that the inner edge of the disk is fixed at the radius of the marginally stable last circular orbit. Accordingly, the loci of a constant black-hole mass on the "H-R diagram" (representing the relation between $L_{\rm x}$ and $\Tin$ both on the logarithmic scales) are not straight but bent towards the lower $M$ direction in the frame of the standard-disk relation.

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

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

We argue that both the extreme soft X-ray excess and the large-amplitude variability of Narrow-Line Seyfert 1 galaxies (NLS1s) can be explained in the framework of the slim disk model. When the disk luminosity approaches the Eddington luminosity, the disk becomes a slim disk, exhibiting a multi-color blackbody spectrum with a maximum temperature, T(bb), of about 0.2 (M/1e5 solar masses)e(-1/4) keV, and size of the X-ray emitting region, r(bb), of about R(S) (the Schwarzschild radius). Furthermore, magnetic energy can be amplified up to a level exceeding radiation energy emitted from the disk, causing substantial variability in X-rays by consecutive magnetic flares.

Very recently the Chandra First Light Observation discovered a point-like source in Cassiopeia A (Cas A) supernova remnant. This detection was subsequently confirmed by the analyses of the archival data from both ROSAT and Einstein observations. Here we compare the results from these observations with the scenarios involving both black holes (BH) and neutron stars (NS). If this point source is a BH we offer as a promising model a disk-corona type model with low accretion rate where a soft photon source at ~ 0.1 keV is Comptonized by higher energy electrons in the corona. If it is a NS the dominant radiation observed by Chandra most likely originates from smaller, hotter regions of the stellar surface, but we argue that it is still worthwhile to compare the cooler component from the rest of the surface with cooling theories. We emphasize that the detection of this point source itself should potentially provide enormous impacts on the theories of supernova explosion, progenitor scenario, compact remnant formation, accretion to compact objects, and neutron star thermal evolution.

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

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

We consider the propagation of hydromagnetic waves generated by a compact turbulent source in low ionized plasmas, applying the Lighthill theory. We assume the plasma to be isothermal, and adopt a uniform, stationary medium thread by ordered magnetic fields as an initial condition. Then, the distinct properties of the hydromagnetic waves originating from a source oscillating with a fixed frequency are studied in the linear regime. As is well known, in low ionized plasmas, the generated waves dissipate due to ion-neutral damping. In this paper, the dependence of the dissipation rate on the frequency of the oscillating source is investigated. The larger the frequency becomes, the more substantial is the wave dissipation. Implications of our results on the energy source in molecular clouds are also discussed. Interestingly, since the outflow lobes associated with young stellar objects act as compact turbulent sources, hydromagnetic waves are generated by them. From our order-estimations, about 70% of the energy of the outflow itself propagates as waves or turbulences, while the remaining 30% dissipates and heats the neutrals via ion-neutral damping. Then, we confirm that the outflows are significant energy sources in molecular clouds in the context of the Lighthill theory.

Some effects of irradiation on time-dependent accretion-disk models for black hole X-ray novae are presented. Two types of irradiation are considered: direct irradiation from the inner hot disk and indirect irradiation as might be reflected by a corona or chromosphere above the disk. The shadowing effect of the time-dependent evolution of the disk height and consequent blocking of the outer disk by the inner and middle portions of the disk from the direct irradiation is included. The direct irradiation of the disk by inner layers where the soft X-ray flux is generated is found to have only a small effect on the outer disk because of shadowing. Mild indirect irradiation that flattens, but otherwise does not affect the light curve substantially, still has interesting non-linear effects on the structure of the disk as heating and cooling waves propagate. The irradiated disks do not always make simple transitions between the hot and cold states, but can linger at intermediate temperatures or even return temporarily to the hot state, depending on the irradiation and the activity in adjacent portions of the disk.

Since the radiation from different portions in the central region of a quasar can be successively amplified during a microlensing event, microlensing light curves provide us with fruitful information regarding the emissivity distribution of an accretion disk located at the quasar center. We present a basic methodology how to map the emissivity distribution of the disk as a function of the radial distance from the center, $Q(r)$, from `observed' microlens light curves during a caustic crossing event. Our proposed method is based on the standard inversion technique, the so-called regularization method, and the Abel's transformation of the one-dimensional luminosity profile integrated along the line parallel to the caustics. The technique will be used to map the disk structure in Q2237+0305, for which the HST and AXAF observations are scheduled. A reconstruction of the image on length scales of several to ten AUs is quite feasible for this source, as long as measuring errors are within 0.02 mag and the observation time intervals are a week or less.

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

The optical-ultraviolet continuum from active galactic nuclei (AGN) seem to originate from optically thick and/or thin disks, and occasionally from associated circumnuclear starburst regions. These different possible origins can, in principle, be discriminated by observations of gravitational microlensing event. We performed numerical simulations of microlensing of an AGN disk by a single lensing star passing in front of the AGN. Calculated spectral variations and light curves show distinct behavior, depending on the nature of emitting region; the time variation of a few months with strong wavelength dependence is expected in the case of an optically thick disk ({\it standard disk}), while an optically thin disk ({\it advection-dominated disk}) will produce shorter, nearly wavelength-independent variation. In the case of a starburst region much slower variations (over a year) will be superposed on the shorter variations.

We examine the radiatively driven mass accretion onto galactic nuclei due to the intensive radiation from circumnuclear starbursts. The radiation from a starburst not only contracts an inner gas disk by the radiation flux force, but also extracts angular momenta due to the relativistic radiation drag, thereby inducing an avalanche of the surface layer of disk. The present radiatively driven mass accretion may provide a physical mechanism which enables mass accretion from 100 pc scales down to pc, and it can eventually link to an advection-dominated viscous accretion onto a massive black hole. The radiation-hydrodynamical and self-gravitational instabilities of the disk are briefly discussed. In particular, the radiative acceleration possibly builds up a dusty wall, which provides another version for the formation of an obscuring torus.

The innermost regions of quasars can be resolved by a gravitational-lens {\lq}telescope{\rq} on scales down to a few AU. For the purpose, X-ray observations are most preferable, because X-rays originating from the innermost regions, can be selectively amplified by microlensing due to the so-called `caustic crossing'. If detected, X-ray variations will constrain the size of the X-ray emitting region down to a few AU. The maximum attainable resolution depends mainly on the monitoring intervals of lens events, which should be much shorter than the crossing time. On the basis of this idea, we performe numerical simulations of microlensing of an optically-thick, standard-type disk as well as an optically-thin, advection-dominated accretion flow (ADAF). Calculated spectral variations and light curves show distinct behaviors, depending on the photon energy. X-ray radiation which is produced in optically thin region, exhibits intensity variation over a few tens of days. In contrast, optical-UV fluxes, which are likely to come from optically thick region, exhibit more gradual light changes, which is consistent with the microlensing events so far observed in Q2237+0305. Currently, Q2237+0305 is being monitored in the optical range at Apache Point Observatory. Simultaneous multi-wavelength observations by X-ray sattelites (e.g., ASCA, AXAF, XMM) as well as HST at the moment of a microlens event enable us to reveal an AU scale structure of the central accretion disk around black hole.

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

The global structure of optically thin advection dominated accretion flows which are composed of two-temperature plasma around black holes is calculated. We adopt the full set of basic equations including the advective energy transport in the energy equation for the electrons. The spectra emitted by the optically thin accretion flows are also investigated. The radiation mechanisms which are taken into accout are bremsstrahlung, synchrotron emission, and Comptonization. The calculation of the spectra and that of the structure of the accretion flows are made to be completely consistent by calculating the radiative cooling rate at each radius. As a result of the advection domination for the ions, the heat transport from the ions to the electrons becomes practically zero and the radiative cooling balances with the advective heating in the energy equation of the electrons. Following up on the successful work of Narayan et al. (1995), we applied our model to the spectrum of Sgr A*. We find that the spectrum of Sgr A* is explained by the optically thin advection dominated accretion flow around a black hole of the mass M_bh=10^6 M_sun. The parameter dependence of the spectrum and the structure of the accretion flows is also discussed.

We examine the possibility that ram pressure exerted by the galactic wind from the Galaxy could have stripped gas from the Local Group dwarf galaxies, thereby affecting their star formation histories. Whether gas stripping occurs or not depends on the relative magnitudes of two counteracting forces acting on gas in a dwarf galaxy: ram pressure force by the wind and the gravitational binding force by the dwarf galaxy itself. We suggest that the galactic wind could have stripped gas in a dwarf galaxy located within the distance of $R_{c}\simeq 120(r_{s}/1 {kpc})^{3/2} ({\cal E}_{b}/10^{50} {erg})^{-1/2}$ kpc (where $r_{s}$ is the surface radius and ${\cal E}_{b}$ is the total binding energy of the dwarf galaxy, respectively) from the Galaxy within a timescale of Gyr, thereby preventing star formation there. Our result based on this Galactic wind model explains the recent observation that dwarfs located close to the Galaxy experienced star formation only in the early phase of their lifetimes, whereas distant dwarfs are still undergoing star formation. The present star formation in the Large Magellanic Cloud can also be explained through our Galactic wind model.

Evolutions of X-ray clusters of galaxies are studied by N-body (shell model) + mesh code (TVD) simulations on the assumption of spherical symmetry. We consider a density perturbation of $10^{15} M_{\odot}$ composed of dark matter and gas in cold dark matter dominated universe with the cosmological density parameter, $\Omega_0 = 1$ or $0.2$. A shock front appears during its initial collapse, moving outward as ambient gas accretes towards cluster center. The shock front separates the inner X-ray emitting, hot region, where gas is almost in hydrostatic equilibrium but with small radial infall ($\sim 100$km s$^{-1}$) being left, from the outer cool region. Gas inside the shock is strongly compressed and heated by shock so that X-ray luminosity rapidly rises in the early stage. In the late stage, however, the X-ray luminosity rises only gradually due partly to the expansion of the inner high temperature region and partly to the increase of X-ray emissivity of gas as the results of continuous adiabatic compression inside the shock. The shock front structure, which was not clearly resolved in the previous SPH simulations, is clearly captured by the present simulations. In addition, we find a sound wave propagating outward, thereby producing spatial modulations with amplitudes of $\sim$ 10 % in the radial temperature and density profiles and time variations in the strength of the shock. Such modulations, if observed, could be used as a probe to investigate the structure of clusters.

A self-similar solution for time evolution of isothermal, self-gravitating viscous disks is found under the condition that $\alpha' \equiv \alpha (H/r)$ is constant in space (where $\alpha$ is the viscosity parameter and $H/r$ is the ratio of a half-thickness to radius of the disk). This solution describes a homologous collapse of a disk via self-gravity and viscosity. The disk structure and evolution is distinct in the inner and outer parts. There is a constant mass inflow in the outer portions so that the disk has flat rotation velocity, constant accretion velocity, and surface density decreasing outward as $\Sigma \propto r^{-1}$. In the inner portions, in contrast, mass is accumulated near the center owing to the boundary condition of no radial velocity at the origin, thereby a strong central concentration being produced; surface density varies as $\Sigma \propto r^{-5/3}$. Moreover, the transition radius separating the inner and outer portions increases linearly with time. The consequence of such a high condensation is briefly discussed in the context of formation of a quasar black hole.

The response of advection-dominated accretion disks to local disturbances is examined by one-dimensional numerical simulations. It is generally believed that advection-dominated disks are thermally stable. We, however, find that any disurbance added onto accretion flow at large radii does not decay so rapidly that it can move inward with roughly the free-fall velocity. Although disturbances continue to be present, the global disk structure will not be modified largely. This can account for persistent hard X-ray emission with substantial variations observed in active galactic nuclei and stellar black hole candidates during the hard state. Moreover, when the disturbance reaches the innermost parts, an acoustic wave emerges, propagating outward as a shock wave. The resultant light variation is roughly (time) symmetric and is quite reminiscent of the observed X-ray shots of Cygnus X-1.

The most plausible mechanism for triggering the outburst of black hole candidate X-ray transients is the ionization thermal instability. The disk instability models can give the observed mass flow in quiescence, but not the X-ray spectrum. Self-irradiation of the disk in outburst may not lead to X-ray reprocessing as the dominant source of optical light, but may play a role in the "reflare." The hard power-law spectrum and radio bursts may be non-thermal processes driven by the flow of pair-rich plasma from the disk at early times and due to the formation of a pair-rich plasma corona at late times. The repeated outbursts suggest some sort of clock, but it is unlikely that it has anything to do with a simple X-ray heating of the companion star.

Recent X-ray observations have established that collisions between subclusters of galaxies are rather common phenomena. Prompted by such observations, we have performed N-body simulations of two equal-mass subclusters of galaxies, which are going to merge. We first have confirmed that only a part of kinetic energy associated with the relative motion of two subclusters is converted to internal energy of each subcluster at the first encounter. We performed simulations for a variety of dark-matter (DM) distributions, and find that the time scale for washing out the double peak structures after the first encounter strongly depends on the distribution of the density and velocity distributions of the DM. It takes longer when the DM is spatially extended. According to our calculation it takes more than $4 imes 10^9$yr after the first encounter until the density contour shows only a single peak. To explain the high fraction of the clusters with substructures among nearby clusters, Richstone,Loeb Turner(1992; hereafter RLT) required recent($$ $10^9 h^{-1}$yr) cluster formation. However, our results show that the timescale for subcluster merging is still uncertain and possibly much longer than the time scale assumed by RLT. Caution should be exercised when concluding that the density of the universe is high by using RLT's method.