Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

We propose a compact binary model with an eccentric orbit to explain periodically active fast radio burst (FRB) sources, where the system consists of a neutron star (NS) with strong dipolar magnetic fields and a magnetic white dwarf (WD). In our model, the WD fills its Roche lobe at periastron, and mass transfer occurs from the WD to the NS around this point. The accreted material may be fragmented into a number of parts, which arrive at the NS at different times. The fragmented magnetized material may trigger magnetic reconnection near the NS surface. The electrons can be accelerated to an ultra-relativistic speed, and therefore the curvature radiation of the electrons can account for the burst activity. In this scenario, the duty cycle of burst activity is related to the orbital period of the binary. We show that such a model may work for duty cycles roughly from ten minutes to two days. For the recently reported 16.35-day periodicity of FRB 180916.J0158+65, our model does not naturally explain such a long duty cycle, since an extremely high eccentricity ($e>0.95$) is required.

When the matter from a companion star is accreted towards the central compact accretor, i.e. a black hole (BH) or a neutron star (NS), an accretion disc and a jet outflow will form, providing bight X-ray and radio emission, which is known as X-ray binaries (XRBs). In the low/hard state, there exist disc-jet couplings in XRBs, but it remains uncertain whether the jet power comes from the disc or the central accretor. Moreover, BHXRBs have different properties compared with NSXRBs: quiescent BHXRBs are typically two to three orders of magnitude less luminous than NSXRBs in X-ray, whereas BHXRBs are more radio loud than NSXRBs. In observations, an empirical correlation has been established between radio and X-ray luminosity, $L_{\rm R} \propto L_{\rm X}^b$, where $b\sim 0.7$ for BHXRBs and $b \sim 1.4$ for non-pulsating NSXRBs. However, there are some outliers of BHXRBs showing unusually steep correlation as NSXRBs at higher luminosities. In this work, under the assumption that the origin of jet power is related to the internal energy of the inner disc, we apply our magnetized, radiatively efficient thin disc model and the well-known radiatively inefficient accretion flow model to NSXRBs and BHXRBs. We find that the observed radio/X-ray correlations in XRBs can be well understood by the disc-jet couplings.

Active galactic nuclei (AGNs) have long been observed to "twinkle" (i.e., their brightness varies with time) on timescales from days to years in the UV/optical bands. Such AGN UV/optical variability is essential for probing the physics of supermassive black holes (SMBHs), the accretion disk, and the broad-line region. Here we show that the temperature fluctuations of an AGN accretion disk, which is magnetically coupled with the corona, can account for observed high-quality AGN optical light curves. We calculate the temperature fluctuations by considering the gas physics of the accreted matter near the SMBH. We find that the resulting simulated AGN UV/optical light curves share the same statistical properties as the observed ones as long as the dimensionless viscosity parameter $\alpha$, which is widely believed to be controlled by magnetohydrodynamic (MHD) turbulence in the accretion disk, is about $0.01$---$0.2$. Moreover, our model can simultaneously explain the larger-than-expected accretion disk sizes and the dependence of UV/optical variability upon wavelength for NGC 5548. Our model also has the potential to explain some other observational facts of AGN UV/optical variability, including the timescale-dependent bluer-when-brighter color variability and the dependence of UV/optical variability on AGN luminosity and black hole mass. Our results also demonstrate a promising way to infer the black-hole mass, the accretion rate, and the radiative efficiency, thereby facilitating understanding of the gas physics and MHD turbulence near the SMBH and its cosmic mass growth history by fitting the AGN UV/optical light curves in the era of time-domain astronomy.

All stellar mass black holes have hitherto been identified by X-rays emitted by gas that is accreting onto the black hole from a companion star. These systems are all binaries with black holes below 30 M$_{\odot}$$^{1-4}$. Theory predicts, however, that X-ray emitting systems form a minority of the total population of star-black hole binaries$^{5,6}$. When the black hole is not accreting gas, it can be found through radial velocity measurements of the motion of the companion star. Here we report radial velocity measurements of a Galactic star, LB-1, which is a B-type star, taken over two years. We find that the motion of the B-star and an accompanying H$\alpha$ emission line require the presence of a dark companion with a mass of $68^{+11}_{-13}$ M$_{\odot}$, which can only be a black hole. The long orbital period of 78.9 days shows that this is a wide binary system. The gravitational wave experiments have detected similarly massive black holes$^{7,8}$, but forming such massive ones in a high-metallicity environment would be extremely challenging to current stellar evolution theories$^{9-11}$.

We study the prospects of searching for black hole (BH) binary systems with a stellar-mass BH and a non-compact visible companion, by utilizing the spectroscopic data of Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST). We simulate the Galactic BH binary population and determine its optical visibility by considering the stellar synthetic population model and the distributions of binary orbital parameters. By convolving the visibility of BH binaries with the LAMOST detection sensitivity, we predict that $\gtrsim$ 400 candidate BH binaries can be found by the low-resolution, non-time-domain survey, and $\sim$ 50-350 candidates by the LAMOST ongoing medium-resolution, time-domain spectroscopic survey. Most of the candidates are short-period (0.2-2 days) binaries with M-, K-, G-, or F-type companions, in which $\sim$ 47% have a mass function (the lower limit of the BH mass) larger than 3 $M_{\odot}$. By complementing the LAMOST spectroscopic data with other photometric/spectroscopic surveys or follow-up observations, these candidates could be confirmed. Therefore, by exploring the LAMOST data, we can enlarge the sample of dynamically confirmed BH binaries significantly, which can improve our understanding of the mass distribution of BHs and the stellar evolution model.

The thickness of a slim disk is determined by the balance between the radiation force and the vertical component of the gravity of the black hole (BH). It was found that vertical gravity increases with the disk height, and it will decrease with the disk height if the disk thickness is above a critical value, which implies that gas at the disk surface may be driven into outflows by radiation force when the disk thickness surpasses the critical value. In this work, we derive a global solution to a slim disk with radiation-driven outflows. We find that the outflows are driven from the disk surface if the mass accretion rate $\dot{m} \gtrsim 1.78-1.91$ ($\dot{m}=\dot{M}/\dot{M}_{\rm Edd}$, and $\dot{M}_{\rm Edd}=L_{\rm Edd}/0.1c^2$) depending on BH mass, while the outflows are suppressed when the mass accretion rate is lower than this critical value. The mass accretion rate decreases with decreasing radius in the disk with outflows, and the rate of the gas swallowed by the BH is always limited to $\dot{m}_{\rm in}\sim 1.78-1.91$ even if the mass accretion rate at the outer edge of the disk is very high. This may set constraints on the massive BH growth through accretion in the early Universe. Due to the presence of outflows, there is an upper limit on the radiation flux of the disk, which leads to saturation of the continuum spectra of the disk with outflows at the high energy band. This may be tested by the observations of quasars or/and BH X-ray binaries.

Most dynamically confirmed stellar-mass black holes and the candidates were originally selected from X-ray outbursts. In the present work, we search for black hole candidates in the LAMOST survey by using the spectra along with photometry from the ASAS-SN survey, where the orbital period of the binary may be revealed by the periodic light curve, such as the ellipsoidal modulation type. Our sample consists of 9 binaries, where each source contains a giant star with large radial velocity variation ($\Delta V_{\rm R} > 70~{\rm km~s^{-1}}$) and periods known from light curves. We focus on the 9 sources with long periods ($T_{\rm ph} > 5$ days) and evaluate the mass $M_2$ of the optically invisible companion. Since the observed $\Delta V_{\rm R}$ from only a few repeating spectroscopic observations is a lower limit of the real amplitude, the real mass $M_2$ can be significantly higher than the current evaluation. It is likely an efficient method to place constraints on $M_2$ by combining $\Delta V_{\rm R}$ from LAMOST and $T_{\rm ph}$ from ASAS-SN, particularly by the ongoing LAMOST Medium Resolution Survey.

We investigated effects on flow variables of transonic advection-dominated accretion flows (ADAFs) for different outer boundary locations (BLs) with a changing energy constant ($E$) of the flow. We used the ADAF solutions and investigated a general power index rule of a radial bulk velocity $(\vr\propto r^{-p})$ with different BLs, but the power index with radius for a rotation velocity and sound speed is unchanged. Here, $p\geq0.5$ is a power index. This power rule gives two types of self-similar solutions; first, when $p=0.5$ gives a self-similar solution of a first kind and exists for infinite length, which has already been discovered for the ADAFs by Narayan \& Yi, and second, when $p>0.5$ gives a self-similar solution of a second kind and exists for finite length, which corresponds to our new solutions for the ADAFs. By using this index rule in fluid equations, we found that the Mach number ($M$) and advection factor ($\fadv$) vary with the radius when $p>0.5$. The local energies of the ADAFs and the Keplerian disk are matched very well at the BLs. So, this theoretical study is supporting a two-zone configuration theory of the accretion disk, and we also discussed other possible hybrid disk geometries. The present study can have two main implications with a variation of the $p$; first, one that can help with the understanding of outflows and non-thermal spectrum variations in black hole candidates, and second, one that can help with solving partial differential equations for any sized advective disk.

GRB 170817A, detected by Fermi-GBM 1.7\,s after the merger of a neutron star (NS) binary, provides the first direct evidence for a link between such a merger and a short-duration gamma-ray burst. The X-ray observations after GRB 170817A indicate a possible X-ray flare with a peak luminosity $L_{\rm peak} \sim 2\times 10^{39}\,{\rm erg\,s}^{-1}$ near day 156. Here we show that this X-ray flare may be understood based on a slim disc around a compact object. On the one hand, there exists the maximal accretion rate $\dot M_{\rm max}$ for the slim disc, above which an optically thick outflow is significant and radiation from the disc is obscured. Based on the energy balance analysis, we find that $\dot M_{\rm max}$ is in the range of $\sim 4\dot M_{\rm Edd}$ and $\sim 21\dot M_{\rm Edd}$ when the angular velocity of the slim disc is between $\rm (1/5)^{1/2}\Omega_K$ and $\rm \Omega_K$ (where $\dot M_{\rm Edd}$ is the Eddington accretion rate and $\Omega_K$ is the Keplerian angular velocity). With $\dot M_{\rm max}$, the slim disc can provide a luminosity $\sim L_{\rm peak}$ for a compact object of $2.5 M_{\sun}$. On the other hand, if the merger of two NSs forms a typical neutrino-dominated accretion disc whose accretion rate $\dot M$ follows a power-law decline with an index $-1.8$ , then the system must pass through the outflow regime and enter the slim disc in $\sim 11-355$ days. These results imply that a post-merger slim accretion disc could account for the observed late-time $L_{\rm peak}$.

We propose the thick-disc model of Gu et al. 2016 to interpret the transition between soft ultraluminous state (SUL) and supersoft ultraluminous (SSUL) state in NGC 247. As accretion rate increases, the inner disc will puff up and act as shield to block the innermost X-ray emission regions and absorb both soft and hard X-ray photons. The absorbed X-ray emission will be re-radiated as a much softer blackbody X-ray spectrum. Hence NGC 247 shows flux dips in the hard X-ray band and transits from the SUL state to the SSUL state. The $\sim 200$s transition timescale can be explained by the viscous timescale. According to our model, the inner disc in the super-soft state is thicker and has smaller viscous timescale than in the soft state. X-ray flux variability, which is assumed to be driven by accretion rate fluctuations, might be viscous time-scale invariant. Therefore, in the SSUL state, NGC 247 is more variable. The bolometric luminosity is saturated in the thick disc; the observed radius-temperature relation can therefore be naturally explained.

By taking into account photon absorption, we investigate the vertical structure of accretion flows with comparable radiation and gas pressures. We consider two separate energy equations for matter and radiation in the diffusion limit. In order to solve the set of radiation hydrodynamic equations in steady state and axisymmetric configuration, we employ self-similar technique in the radial direction. We need the reflection symmetry about the mid-plane to find gas density at the equator. For a typical solution, we assume that the gas pressure has 10-50\% portion of the total pressure. In this paper, since the radiation energy is involved directly, we are able to estimate how much energy of viscous heating is transported in the radial direction and advected towards the central object. Our results show that although the mass accretion rate does not approach the Eddington limit, the energy advection is rather high. Moreover, in a disc with greater accretion rate and less portion of gas pressure at the total pressure, more energy is advected to its center. In addition, as we expect the accretion flow becomes thicker with greater values of gas pressure. Based on Solberg-H\~A\v{z}iland conditions, we notice that the flow is convectively stable in all parts of such a disc.

iPTF14hls is an unusually bright, long-lived II-P supernova (SN), whose light curve has at least five peaks. We propose that the outflows from the black hole hyperaccretion systems in the center of the collapsars should continuously inject into the envelope. For a jet-driven core-collapsar model, the outflow feedback results in prolonging the accretion timescale and fluctuating accretion rates in our analytic solutions. Thus, the long period of luminous, varying SN iPTF14hls might originate from the choked jets, which are regulated by the feedback of the strong disk outflows in a massive core-collapsar. One can expect that jet-driven iPTF14hls may last no more than approximately 3,000 days, and the luminosity may quickly decrease in the later stages. Moreover, the double-peak light curves in some SNe might be explained by the outflow feedback mechanism.

We propose a method to search for stellar-mass black hole (BH) candidates with giant companions from spectroscopic observations. Based on the stellar spectra of LAMOST Data Release 6, we obtain a sample of seven giants in binaries with large radial velocity variation $\Delta V_R > 80~{\rm km~s^{-1}}$. With the effective temperature, surface gravity, and metallicity provided by LAMOST, and the parallax given by {\it Gaia}, we can estimate the mass and radius of the giant, and therefore evaluate the possible mass of the optically invisible star in the binary. We show that the sources in our sample are potential BH candidates, and are worthy of dynamical measurement by further spectroscopic observations. Our method may be particularly valid for the selection of BH candidates in binaries with unknown orbital periods.

Jet launching in radio loud (RL) quasars is one of the fundamental problems in astrophysics. Exploring the differences in the inner accretion disk properties between RL and radio quiet (RQ) quasars might yield helpful clues to this puzzle. We previously discovered that the shorter term UV/optical variations of quasars are bluer than the longer term ones, i.e., the so-called timescale-dependent color variation. This is consistent with the scheme that the faster variations come from the inner and hotter disk regions, thus providing a useful tool to map the accretion disk which is otherwise unresolvable. In this work we compare the UV/optical variations of RL quasars in SDSS Stripe 82 to those of several RQ samples, including those matched in redshift-luminosity-black hole mass and/or color-magnitude. We find that while both RL and RQ populations appear bluer when they brighten, RL quasars potentially show a weaker/flatter dependence on timescale in their color variation. We further find that while both RL and RQ populations on average show similar variation amplitudes at long timescales, fast variations of RL sources appear weaker/smaller (at timescales of ~ 25 -- 300 days in the observer's frame), and the difference is more prominent in the g-band than in the r-band. Inhomogeneous disk simulations can qualitatively reproduce these observed differences if the inner accretion disk of RL quasars fluctuates less based on simple toy models. Though the implications are likely model dependent, the discovery points to an interesting diagram that magnetic fields in RL quasars may be prospectively stronger and play a key role in both jet launching and the stabilization of the inner accretion disk.

Recent multi-band variability studies have revealed that active galactic nucleus (AGN) accretion disc sizes are generally larger than the predictions of the classical thin disc by a factor of $2\sim 3$. This hints at some missing key ingredient in the classical thin disc theory: here, we propose an accretion disc wind. For a given bolometric luminosity, in the outer part of an accretion disc, the effective temperature in the wind case is higher than that in the no-wind one; meanwhile, the radial temperature profile of the wind case is shallower than the no-wind one. In presence of winds, for a given band, blackbody emission from large radii can contribute more to the observed luminosity than the no-wind case. Therefore, the disc sizes of the wind case can be larger than those of the no-wind case. We demonstrate that a model with the accretion rate scaling as $\dot{M}_0 (R/R_{\mathrm{S}})^{\beta}$ (i.e., the accretion rate declines with decreasing radius due to winds) can match both the inter-band time lags and the spectral energy distribution of NGC 5548. Our model can also explain the inter-band time lags of other sources. Therefore, our model can help decipher current and future continuum reverberation mapping observations.

Observations have shown that UV/optical variation amplitude of quasars depend on several physi- cal parameters including luminosity, Eddington ratio, and likely also black hole mass. Identifying new factors which correlate with the variation is essential to probe the underlying physical processes. Combining ~ten years long quasar light curves from SDSS stripe 82 and X-ray data from Stripe 82X, we build a sample of X-ray detected quasars to investigate the relation between UV/optical variation amplitude ($\sigma_{rms}$) and X-ray loudness. We find that quasars with more intense X-ray radiation (com- pared to bolometric luminosity) are more variable in UV/optical. Such correlation remains highly significant after excluding the effect of other parameters including luminosity, black hole mass, Ed- dington ratio, redshift, rest-frame wavelength (i.e., through partial correlation analyses). We further find the intrinsic link between X-ray loudness and UV/optical variation is gradually more prominent on longer timescales (up to 10 years in the observed frame), but tends to disappear at timescales < 100 days. This suggests a slow and long-term underlying physical process. The X-ray reprocessing paradigm, in which UV/optical variation is produced by a variable central X-ray emission illuminating the accretion disk, is thus disfavored. The discovery points to an interesting scheme that both the X-ray corona heating and UV/optical variation is quasars are closely associated with magnetic disc turbulence, and the innermost disc turbulence (where corona heating occurs) correlates with the slow turbulence at larger radii (where UV/optical emission is produced).

GRS 1915+105 has been active for more than 26 years since it was discovered in 1992. There are hundreds of RXTE pointed observations on this source, and the quasi-regular flares with a slow rise and a sharp decrease (i.e. the "heartbeat" state) were recorded in more than 200 observations. The connections among the disk/corona, jet, and the disk wind at the heartbeat state have been extensively studied. In this work, we firstly perform a statistical analysis of the light curves and the X-ray spectra to investigate this peculiar state. We calculate the parameters for heartbeat cycles, including the recurrence time, the maximum and the minimum count rate, the flare amplitude, and the cumulative radiation for each cycle. The recurrence time has a bimodal distribution ranging from $\sim 20$ to $\sim 200$ s. The minimum count rate increases with increasing recurrence time; while the maximum count rate remains nearly constant around 2 Crab. Fitting the averaged spectrum for each observation, we find the strong correlations among the recurrence time, the apparent inner radius of the accretion disk (or the color correction factor), and the (nonthermal) X-ray luminosity. We suggest that the true inner edge of the accretion disk might always extend to the marginally stable orbit, while the change in corona size should result in the observed correlations.

The occurrence rates of bright X-ray flares in z<1 gamma-ray bursts (GRBs) with or without observed supernovae (SNe) association were compared. Our Sample I: the z<1 long GRBs (LGRBs) with SNe association (SN-GRBs) and with early Swift/X-Ray Telescope (XRT) observations, consists of 18 GRBs, among which only two GRBs have bright X-ray flares. Our Sample II: for comparison, all the z<1 LGRBs without observed SNe association and with early Swift/XRT observations, consists of 45 GRBs, among which 16 GRBs present bright X-ray flares. Thus, the study indicates a lower occurrence rate of bright X-ray flares in Sample I (11.1%) than in Sample II (35.6%). In addition, if dim X-ray fluctuations are included as flares, then 16.7% of Sample I and 55.6% of Sample II are found to have flares, again showing the discrepancy between these two samples. We examined the physical origin of these bright X-ray flares and found that most of them are probably related to the central engine reactivity. To understand the discrepancy, we propose that such a lower occurrence rate of flares in the SN-GRB sample may hint at an energy partition among the GRB, SNe, and X-ray flares under a saturated energy budget of massive star explosion.

We solved analytically viscous two-dimensional (2D) fluid equations for accretion and outflows in spherical polar coordinates ($r, \theta, \phi$) and obtained explicitly flow variables in $r-$ and $\theta -$directions around black holes (BHs). We investigated global transonic advection-dominated accretion flow (ADAF) solutions in $r-$direction on an equatorial plane with using Paczy\'nski-Wiita potential. We used radial flow variables of ADAFs with symmetric conditions on the equatorial plane, as initial values for integration in $\theta-$direction. In the study of 2D disk structure, we used two-azimuthal components of viscous stress tensors namely, $\tau_{\rm r\phi}$ and $\tau_{\rm\theta\phi}$. Interestingly, we found that the whole advective disk is not participating in outflow generation and the outflows form close to the BHs. Normally, outflow strength increased with increasing viscosity parameter ($\alpha_1$), mass-loss parameter ($s$) and decreasing gas pressure ratio ($\beta$). Outflow region increased with increasing $s$, $\alpha_1$ for $\tau_{\rm r\phi}$ and decreasing $\alpha_2$ for $\tau_{\rm\theta\phi}$. The $\tau_{\rm\theta\phi}$ is effective in angular momentum transportation at high latitude and outflows collimation along an axis of symmetry since it changes polar velocity ($v_{\rm\theta}$) of the flow. The outflow emission is also affected by the ADAF size and decreased with decreasing it. Transonic surfaces formed for both inflows ($v_{\rm r}<0$, very close to BH) and outflows ($v_{\rm r}>0$). We also explored no outflows, outflows and failed outflows regions, which mainly depend on the viscosity parameters.

Short gamma-ray bursts may originate from the merger of double neutron stars (NS) or that of a black hole (BH) and an NS. We propose that the bright X-ray flare related to the central engine reactivity may hint a BH-NS merger, since such a merger can provide more fall-back materials and therefore a more massive accretion disk than the NS-NS merger. Based on the observed 49 short bursts with Swift/X-ray Telescope follow-up observations, we find that three bursts have bright X-ray flares, among which three flares from two bursts are probably related to the central engine reactivity. We argue that these two bursts may originate from the BH-NS merger rather than the NS-NS merger. Our suggested link between the central engine-powered bright X-ray flare and the BH-NS merger event can be checked by the future gravitational wave detections from advanced LIGO and Virgo.

We study the first electromagnetic pulse after the gravitational wave chirp signal, focusing on the profile and duration. It is found that the light curve, especially the steep decay (SD) phase, can be very different by adopting different viewing angle $\theta_{\rm view}$ on the jet shell. For an on-axis jet with a power-law radiation spectrum, the observed flux in the SD is proportional to $t_{\rm{obs}}^{-2-\beta}$ with $\beta$ being the spectral index and $t_{\rm{obs}}$ being the observer time. Here, $t_{\rm{obs}}=0$ is set at the observed time of the jet ejected from the central engine. The SD may become steep by increasing $\theta_{\rm view}$. We also study the bolometric luminosity $L$ from a jet shell with a non-power-law radiation spectrum. For an on-axis jet, $L{\propto}t_{\rm{obs}}^{-3}$ is found in the SD. However, the SD is steeper than $L{\propto}t_{\rm{obs}}^{-3}$ for the radiation from an off-axis jet. The higher value of $\theta_{\rm view}$ is, the steeper of SD would be. Then, we suggest that the SD phase can be used to discriminate an off-axis jet from an on-axis jet. The reason for above behaviors is discussed. In addition, we find that the duration of first electromagnetic pulse is close to its peak time, especially for $\theta_{\rm{view}}\sim20^\circ$. This result is consistent with that found in GW~170817/GRB~170817A. Thus, the jet corresponding to the prompt emission of GRB~170817A should be immediately ejected after the merger. Our results also reveal that the duration of the first electromagnetic pulse can provide the information of the time to search gravitational waves.

The tight inter-band correlation and the lag-wavelength relation among UV/optical continua of active galactic nuclei have been firmly established. They are usually understood within the widespread reprocessing scenario, however, the implied inter-band lags are generally too small. Furthermore, it is challenged by new evidences, such as the X-ray reprocessing yields too much high frequency UV/optical variations as well as it fails to reproduce the observed timescale-dependent color variations among {\it Swift} lightcurves of NGC 5548. In a different manner, we demonstrate that an upgraded inhomogeneous accretion disk model, whose local {\it independent} temperature fluctuations are subject to a speculated {\it common} large-scale temperature fluctuation, can intrinsically generate the tight inter-band correlation and lag across UV/optical, and be in nice agreement with several observational properties of NGC 5548, including the timescale-dependent color variation. The emergent lag is a result of the {\it differential regression capability} of local temperature fluctuations when responding to the large-scale fluctuation. An average speed of propagations as large as $\gtrsim 15\%$ of the speed of light may be required by this common fluctuation. Several potential physical mechanisms for such propagations are discussed. Our interesting phenomenological scenario may shed new light on comprehending the UV/optical continuum variations of active galactic nuclei.

Two plateaus and one following bump in the X-ray light curve of GRB 170714A have been detected by the \textit{Swift}/X-Ray Telescope, which could be very meaningful for the central engine of gamma-ray bursts (GRBs), implying that the origin of this burst might be different from that of other ultra-long GRBs. We propose that merging two neutron stars into a hyper-massive quark star (QS) and then collapsing into a black hole (BH), with a delay time around $10^4$~s, could be responsible for those X-ray components. The hyper-massive QS is initially in a fluid state, being turbulent and differentially rotating, but would be solidified and release its latent heat injected into the GRB fireball (lasting about $10^3$~s during the liquid-solid phase transition). Magnetic field as high as $\sim 10^{15}$~G could be created by dynamo action of the newborn liquid QS, and a magnetar-like central engine (after solidification) supplies significant energy for the second plateau. More energy could be released during a fall-back accretion after the post-merger QS collapses to a BH, and the X-ray bump forms. This post-merger QS model might be tested by future observations, with either advanced gravitational wave detectors (e.g., advanced LIGO and VIRGO) or X-ray/optical telescopes.

Outflows launched from a remnant disc of compact binary merger may have essential contribution to the kilonova emission. Numerical calculations are conducted in this work to study the structure of accretion flows and outflows. By the incorporation of limited-energy advection in the hyper-accretion discs, outflows occur naturally from accretion flows due to imbalance between the viscous heating and the sum of the advective and radiative cooling. Following this spirit, we revisit the properties of the merger outflow ejecta. Our results show that around $10^{-3} \sim 10^{-1} M_\odot$ of the disc mass can be launched as powerful outflows. The amount of unbound mass varies with the disc mass and the viscosity. The outflow-contributed peak luminosity is around $10^{40} \sim 10^{41}$ erg/s. Such a scenario can account for the observed kilonovae associated with short gamma-ray bursts, including the recent event AT2017gfo (GW170817).

Gamma-ray bursts (GRBs) are luminous and violent phenomena in the universe. Traditionally, long GRBs are expected to be produced by the collapse of massive stars and associated with supernovae. However, some low-redshift long GRBs have no detection of supernova association, such as GRBs 060505, 060614 and 111005A. It is hard to classify these events convincingly according to usual classifications, and the lack of the supernova implies a non-massive star origin. We propose a new path to produce long GRBs without supernova association, the unstable and extremely violent accretion in a contact binary system consisting of a stellar-mass black hole and a white dwarf, which fills an important gap in compact binary evolution.

Central engine of gamma-ray bursts (GRBs) may be intermittent and launch several episodes of ejecta separated by a long quiescent interval. In this scenario, an external shock is formed due to the propagation of the first launched ejecta into the circum-burst medium and the later launched ejecta may interact with the external shock at later period. Owing to the internal dissipation, the later launched ejecta may be observed at a later time ($t_{\rm{jet}}$). In this paper, we study the relation of $t_{\rm{b}}$ and $t_{\rm{jet}}$, where $t_{\rm{b}}$ is the collision time of the later launched ejecta with the formed external shock. It is found that the relation of $t_{\rm{b}}$ and $t_{\rm{jet}}$ depends on the bulk Lorentz factor ($\Gamma_{\rm{jet}}$) of the later launched ejecta and the density ($\rho$) of the circum-burst medium. If the value of $\Gamma_{\rm{jet}}$ or $\rho$ is low, the $t_{\rm{b}}$ would be significantly larger than $t_{\rm{jet}}$. However, the $t_{\rm{b}}\sim t_{\rm{jet}}$ can be found if the value of $\Gamma_{\rm{jet}}$ or $\rho$ is significantly large. Our results can explain the large lag of the optical emission relative to the $\gamma$-ray/X-ray emission in GRBs, e.g., GRB~111209A. For GRBs with a precursor, our results suggest that the energy injection into the external shock and thus more than one external-reverse shock may appear in the main prompt emission phase. According to our model, we estimate the Lorentz factor of the second launched ejecta in GRB~160625B.

Long-duration gamma-ray bursts (LGRBs) and ultra-LGRBs (ULGRBs) originate from collapsars, in the center of which a newborn rotating stellar-mass black hole (BH) surrounded by a massive accretion disk may form. In the scenario of BH hyperaccretion inflow-outflow model and Blandford-Znajek (BZ) mechanism to trigger gamma-ray bursts (GRBs), the real accretion rate to power a BZ jet is far lower than the mass supply rate from the progenitor star. The characteristics of the progenitor stars can be constrained by GRB luminosity observations, and the results exceed usual expectations. LGRBs lasting from several seconds to tens of seconds in the rest frame may originate from solar-metallicity ($Z \sim 1~ Z_\odot$, where $Z$ and $Z_\odot$ are the metallicities of progenitor stars and the Sun), massive ($M \gtrsim 34 ~M_\odot$, where $M$ and $M_\odot$ are the masses of progenitor stars and the Sun) stars or some zero-metallicity ($Z \sim 0$) stars. A fraction of low-metallicity ($Z \lesssim 10^{-2}~Z_\odot$) stars, including Population III stars, can produce ULGRBs such as GRB 111209A. The fraction of LGRBs lasting less than tens of seconds in the rest frame is more than 40$\%$, which cannot conform to the fraction of the demanded type of progenitor star. It possibly implies that the activity timescale of central engine may be much longer than the observed timescale of prompt emission phase, as indicated by X-ray late-time activities. Alternatively, LGRBs and ULGRBs may be powered by a millisecond magnetar central engine.

Neutrino-dominated accretion flows (NDAFs) around rotating stellar-mass black holes (BHs) are plausible candidates for the central engines of gamma-ray bursts (GRBs). NDAFs are hyperaccretion disks with accretion rates in the range of around 0.001-10 $M_\odot ~\rm s^{-1}$, which have high density and temperature and therefore are extremely optically thick and geometrically slim or even thick. We review the theoretical progresses in studying the properties of NDAFs as well as their applications to the GRB phenomenology. The topics include: the steady radial and vertical structure of NDAFs and the implications for calculating neutrino luminosity and annihilation luminosity, jet power due to neutrino-antineutrino annihilation and Blandford-Znajek mechanism and their dependences on parameters such as BH mass, spin, and accretion rate, time evolution of NDAFs, effect of magnetic fields, applications of NDAF theories to the GRB phenomenology such as lightcurve variability, extended emission, X-ray flares, kilonovae, etc., as well as probing NDAFs using multi-messenger signals such as MeV neutrinos and gravitational waves.

We derive a simple analytical formula to describe the evolution of spectral index $\beta$ in the steep decay phase shaped by the curvature effect with assumption that the spectral parameters and Lorentz factor of jet shell is the same for different latitude. Here, the value of $\beta$ is estimated in 0.3$-$10keV energy band. For a spherical thin shell with a cutoff power law (CPL) intrinsic radiation spectrum, the spectral evolution can be read as a linear function of observer time. For the situation with Band function intrinsic radiation spectrum, the spectral evolution may be complex. If the observed break energy of radiation spectrum is larger than 10keV, the spectral evolution is the same as that shaped by jet shells with a CPL spectrum. If the observed break energy is less than 0.3keV, the value of $\beta$ would be a constant. Others, the spectral evolution can be approximated as a logarithmal function of the observer time in generally.

The curvature effect may be responsible for the steep decay phase observed in gamma-ray bursts. For testing the curvature effect with observations, the zero time point $t_0$ adopted to plot observer time and flux on a logarithmic scale should be appropriately selected. In practice, however, the true $t_0$ cannot be directly constrained from the data. Then, we move $t_0$ to a certain time in the steep decay phase, which can be easily identified. In this situation, we derive an analytical formula to describe the flux evolution of the steep decay phase. The analytical formula is read as $F_\nu\propto (1+\tilde t_{\rm obs}/{\tilde t_c})^{-\alpha}$ with $\alpha(\tilde{t}_{\rm obs})=2+{\int_{0}^{\log (1+\tilde{t}_{\rm obs}/{\tilde{t}_c})} {\beta(\tau)d[\log(1+\tau/{\tilde{t}_c})]}}/{\log (1 + {\tilde t}_{\rm obs}/{{\tilde t}_c})}$, where $F_\nu$ is the flux observed at frequency $\nu$, $\tilde t_{\rm obs}$ is the observer time by setting zero time point $t_0$ at a certain time in the steep decay phase, $\beta$ is the spectral index estimated around $\nu$, and ${\tilde t}_c$ is the decay timescale of the phase with $\tilde{t}_{\rm obs}{\geqslant}0$. We test the analytical formula with the data from numerical calculations. It is found that the analytical formula presents a well estimation about the evolution of flux shaped by the curvature effect. Our analytical formula can be used to confront the curvature effect with observations and estimate the decay timescale of the steep decay phase.

In the classic picture of standard thin accretion disks, the viscous heating is balanced by the radiative cooling through the diffusion process, and the radiation-pressure-dominated inner disk suffers convective instability. However, recent simulations have shown that the vertical advection process owing to the magnetic buoyancy can make significant contribution to the energy transport. In addition, no convective instability has been found by comparing the simulation results with the local convective stability criterion. In this work, following the spirit of simulations, we revisit the vertical structure of radiation-pressure-dominated thin disks by including the vertical advection process. Our study indicates a link between the additional energy transport and the convectively stable property. Thus, the vertical advection not only has significant contribution to the energy transport, but also plays an important role to make the disk convectively stable. Our analyses may be helpful to understand the discrepancy between the classic theory and simulations on standard thin disks.

Recent simulations on super-Eddington accretion flows have shown that, apart from the diffusion process, the vertical advection based on magnetic buoyancy can be a more efficient process to release the trapped photons in the optically thick disk. As a consequence, the radiative luminosity from the accretion disk can be far beyond the Eddington value. Following this spirit, we revisit the structure and radiation of hyper-accretion disks with mass accretion rates in the range $10^{-3}\sim 10~M_{\sun}~{\rm s}^{-1}$. Our results show that, due to the strong cooling through the vertical advection, the disk temperature becomes lower than that in the classic model without the vertical advection process, and therefore the neutrino luminosity from the disk is lower. On the other hand, the gamma-ray photons released through the vertical advection can be extremely super-Eddington. We argue that the large amount of escaped gamma-ray photons may have more significant contribution to the primordial fireball than the neutrino annihilation, and may hint a link between gamma-ray bursts and kilonovae in the black hole hyper-accretion scenario.

The hyperluminous X-ray source (HLX-1, the peak X-ray luminosity $\sim 10^{42}\rm erg\ s^{-1}$) near the spiral galaxy ESO 243-49 is possibly the best candidate for intermediate mass black hole (IMBH), which underwent recurrent outbursts with a period of $\sim 400$ days. The physical reason for this quasi-periodic variability is still unclear. We explore the possibility of radiation-pressure instability in accretion disk by modeling the light curve of HLX-1, and find that it can roughly reproduce the duration, period and amplitude of the recurrent outbursts HLX-1 with an IMBH of ~10^5Msun. Our result provides a possible mechanism to explain the recurrent outbursts in HLX-1. We further find a universal correlation between the outburst duration and the bolometric luminosity for the BH sources with a very broad mass range (e.g., X-ray binaries, XRBs, HLX-1 and active galactic nuclei, AGNs), which is roughly consistent with the prediction of radiation-pressure instability of the accretion disk. These results imply that "heartbeat" oscillations triggered by radiation-pressure instability may appears in different-scale BH systems.

This work focuses on a sample of seven extremely late-time X-ray flares with peak time $t_{\rm p} > 10^4 {\rm s}$, among which two flares can be confirmed as the late-time activity of central engine. The main purpose is to investigate the mechanism of such late-time flares based on the internal origin assumption. In the hyper-accreting black hole (BH) scenario, we study the possibility of two well-known mechanisms as the central engine to power such X-ray flares, i.e., the neutrino-antineutrino annihilation and the Blandford-Znajek (BZ) process. Our results show that the annihilation luminosity is far below the observational data. Thus, the annihilation mechanism cannot account for such late-time flares. For the BZ process, if the role of outflows is taken into consideration, the inflow mass rate near the horizon will be quite low such that the magnetic field will probably be too weak to power the observed X-ray flares. We therefore argue that, for the late-time flares with internal origin, the central engine is unlikely to be associated with BHs. On the contrary, a fast rotating neutron star with strong bipolar magnetic fields may be responsible for such flares.

HLX-1, currently the best intermediate-mass black hole candidate, has undergone seven violent outbursts, each with a peak X-ray luminosity of $L_{\mathrm{peak},\mathrm{X}}\sim 10^{42}\ \rm{erg\ s^{-1}}$. Interestingly, the properties of the HLX-1 outbursts evolve with time. In this work, we aim to constrain the physical parameters of the central engine of the HLX-1 outbursts in the framework of the black hole accretion. We find that the physical properties of the HLX-1 outbursts are consistent with being driven by the radiation pressure instability. This scenario can explain the evolution of the recurrent timescales of the HLX-1 outbursts as a function of the durations.

The emission from black hole binaries (BHBs) and active galactic nuclei (AGNs) displays significant aperiodic variabilities. The most promising explanation for these variabilities is the propagating fluctuations in the accretion flow. It is natural to expect that the mechanism driving variabilities in BHBs and AGNs may operate in a black hole hyper-accretion disk, which is believed to power gamma-ray bursts (GRBs). We study the variabilities of jet power in GRBs based on the model of propagating fluctuations. It is found that the variabilities of jet power and the temporal profile of erratic spikes in this scenario are similar to those in observed light curves of prompt gamma-ray emission of GRBs. Our results show that the mechanism driving X-ray variabilities in BHBs and AGNs may operate in the central engine to drive the variabilities of GRBs.

The UV/optical variability of active galactic nuclei and quasars is useful for understanding the physics of the accretion disk and is gradually attributed to the stochastic fluctuations over the accretion disk. Quasars generally appear bluer when they brighten in the UV/optical, the nature of which remains controversial. Recently \citeauthor{Sun2014} discovered that the color variation of quasars is timescale dependent, in the way that faster variations are even bluer than longer term ones. While this discovery can directly rule out models that simply attribute the color variation to contamination from the host galaxies, or to changes in the global accretion rates, it favors the stochastic disk fluctuation model as fluctuations in the innermost hotter disk could dominate the short-term variations. In this work, we show that a revised inhomogeneous disk model, where the characteristic timescales of thermal fluctuations in the disk are radius-dependent (i.e., $\tau \sim r$; based on the one originally proposed by \citeauthor{DexterAgol2011}), can well reproduce a timescale dependent color variation pattern, similar to the observed one and unaffected by the un-even sampling and photometric error. This demonstrates that one may statistically use variation emission at different timescales to spatially resolve the accretion disk in quasars, thus opens a new window to probe and test the accretion disk physics in the era of time domain astronomy. Caveats of the current model, which ought to be addressed in future simulations, are discussed.

We propose a compact binary model for the fast radio burst (FRB) repeaters, where the system consists of a magnetic white dwarf (WD) and a neutron star (NS) with strong bipolar magnetic fields. When the WD fills its Roche lobe, mass transfer will occur from the WD to the NS through the inner Lagrange point. The accreted magnetized materials may trigger magnetic reconnection when they approach the NS surface, and therefore the electrons can be accelerated to an ultra-relativistic speed. In this scenario, the curvature radiation of the electrons moving along the NS magnetic field lines can account for the characteristic frequency and the timescale of an FRB. Owing to the conservation of angular momentum, the WD may be kicked away after a burst, and the next burst may appear when the system becomes semi-detached again through the gravitational radiation. By comparing our analyses with the observations, we show that such an intermittent Roche lobe overflow mechanism can be responsible for the observed repeating behavior of FRB 121102.

Based on the relativistic mean-field effective interactions theory, and Lai dong model \citep{b37, b38, b39}, we discuss the influences of superstrong magnetic fields (SMFs) on electron Fermi energy, nuclear blinding energy, and single-particle level structure in magnetars surface. By using the method of Shell-Model Monte Carlo (SMMC), and the Random Phase Approximation (RPA) theory, we detailed analyze the neutrino energy loss rates(NELRs) by electron capture (EC) for iron group nuclei in SMFs.

Long-duration gamma-ray bursts (LGRBs) are generally considered to originate from the massive collapsars. It is believed that the central engine of gamma-ray bursts (GRBs) is a neutrino-dominated accretion flow (NDAF) around a rotating stellar-mass black hole (BH). The neutrino annihilation above the NDAF is a feasible mechanism to power GRB. In this work, we analyse the distributions of the isotropic gamma-ray radiated energy and jet kinetic energy of 48 LGRBs. According to the NDAF and fireball models, we estimate the mean accreted masses of LGRBs in our sample to investigate whether the NDAFs can power LGRBs with the reasonable BH parameters and conversion efficiency of neutrino annihilation. The results indicate that most of the values of the accreted masses are less than $5~M_\odot$ for the extreme Kerr BHs and high conversion efficiency. It suggests that the NDAFs may be suitable for most of LGRBs except for some extremely high energy sources.

We propose a geometrically thick, super-Eddington accretion disk model, where an optically thick wind is not necessary, to understand ultraluminous supersoft sources (ULSs). For high mass accretion rates $\dot M \ga 30\dot M_{\rm Edd}$ and not small inclination angles $\theta \ga 25^{\circ}$, where $\dot M_{\rm Edd}$ is the Eddington accretion rate, the hard photons from the hot inner region may be shaded by the geometrically thick inner disk, and therefore only the soft photons from the outer thin disk and the outer photosphere of the thick disk can reach the observer. Our model can naturally explain the approximate relation between the typical thermal radius and the thermal temperature, $R_{\rm bb} \propto T_{\rm bb}^{-2}$. Moreover, the thick disk model can unify ULSs and normal ultraluminous X-ray sources, where the different observational characteristics are probably related to the inclination angle and the mass accretion rate. By comparing our model with the optically thick outflow model, we find that less mass accretion rate is required in our model.

A neutrino-dominated accretion disk around a stellar-mass black hole (BH) can power a gamma-ray burst (GRB) via annihilation of neutrinos launched from the disk. For the BH hyperaccretion system, high accretion rate should trigger the violent evolution of the BH's characteristics, which further leads to the evolution of the neutrino annihilation luminosity. In this paper, we consider the evolution of the accretion system to analyze the mean time-dependent neutrino annihilation luminosity with the different mean accretion rates and initial BH parameters. By time-integrating the luminosity, the total neutrino annihilation energy with the reasonable initial disk mass can satisfy the most of short-duration GRBs and about half of long-duration GRBs. Moreover, the extreme Kerr BH should exist in the cental engines of some high-luminosity GRBs. GRBs with higher energy have to request the alternative magnetohydrodynamics processes in the centers, such as the Blandford-Znajek jet from the accretion system or the millisecond magnetar.

The formation of relativistic jets by an accreting compact object is one of the fundamental mysteries of astrophysics. While the theory is poorly understood, observations of relativistic jets from systems known as microquasars have led to a well-established phenomenology. Relativistic jets are not expected from sources with soft or supersoft X-ray spectra, although two such systems are known to produce relatively low-velocity bipolar outflows. Here we report optical spectra of an ultraluminous supersoft X-ray source (ULS) in the nearby galaxy M81 (M81 ULS-1) showing blueshifted broad H\alpha\ emission lines, characteristic of baryonic jets with relativistic speeds. The time variable jets have projected velocities ~17 per cent of the speed of light, and seem similar to those in the prototype microquasar SS 433. Such relativistic jets are not expected to be launched from white dwarfs, but an origin from a black hole or neutron star in M81 ULS-1 is hard to reconcile with its constant soft X-rays. The completely unexpected presence of relativistic jets in a ULS challenges the canonical theories for jet formation, but may possibly be explained by a long speculated super-critically accreting black hole with optically thick outflows

The variability of accretion rate, which is believed to induce the aperiodic variability of X-ray emission from disk, may affect the energy injection into the jet. In this spirit, a correlation between disk emission and jet emission can be formed even if the mean luminosity of disk emission keeps constant. In this work, these correlations are found in the situation that the luminosity of disk emission is variable and kept with a constant mean value. The obtained correlations may be shallower than that of the fundamental plane of black hole activity. In addition, the slope of correlation may increase with increasing observed frequency of jet emission. For the luminosities spacing with three days, the slope of correlation decreases with increasing black hole mass. The deviation of our found correlations from that of the fundamental plane is related to the suppression of variability in the jet emission in comparison with that in the disk emission. This mechanism may work in some sources in which shallower correlations have been reported. Moreover, it implies that luminosities used to estimate the relation of fundamental plane should cover an appropriate timescale, in which the variability of jet emission is not significantly suppressed.

By proposing a pure leptonic radiation model, we study the potential gamma-ray emissions from jets of the low-mass X-ray binaries. In this model, the relativistic electrons that are accelerated in the jets are responsible for radiative outputs. Nevertheless, dynamics of jets are dominated by the magnetic and proton-matter kinetic energies. The model involves all kinds of related radiative processes and considers the evolution of relativistic electrons along the jet by numerically solving the kinetic equation. Numerical results show that the spectral energy distributions can extend up to TeV bands, in which synchrotron radiation and synchrotron self-Compton scattering are dominant components. As an example, we apply the model to the low-mass X-ray binary GX 339-4. The results can not only reproduce the currently available observations from GX 339-4, but also predict detectable radiation at GeV and TeV bands by Fermi and CTA telescopes. The future observations with Fermi and CTA can be used to test our model, which could be employed to distinguish the origin of X-ray emissions.

A neutrino-dominated accretion flow (NDAF) around a rotating stellar-mass black hole (BH) is one of the plausible candidates for the central engine of gamma-ray bursts (GRBs). Two mechanisms, i.e., Blandford-Znajek (BZ) mechanism and neutrino annihilation process, are generally considered to power GRBs. Using the analytic solutions from Xue et al. (2013) and ignoring the effects of the magnetic field configuration, we estimate the BZ and neutrino annihilation luminosities as the functions of the disk masses and BH spin parameters to contrast the observational jet luminosities of GRBs. The results show that, although the neutrino annihilation processes could account for most of GRBs, the BZ mechanism is more effective, especially for long-duration GRBs. Actually, if the energy of afterglows and flares of GRBs is included, the distinction between these two mechanisms is more significant. Furthermore, massive disk mass and high BH spin are beneficial to power high luminosities of GRBs. Finally, we discuss possible physical mechanisms to enhance the disk mass or the neutrino emission rate of NDAFs and relevant difference between these two mechanisms.

Stellar-mass black holes (BHs) surrounded by neutrino-dominated accretion flows (NDAFs) are the plausible candidates to power gamma-ray bursts (GRBs) via neutrinos emission and their annihilation. The progenitors of short-duration GRBs (SGRBs) are generally considered to be compact binaries mergers. According to the simulation results, the disk mass of the NDAF has been limited after merger events. We can estimate such disk mass by using the current SGRB observational data and fireball model. The results show that the disk mass of a certain SGRB mainly depends on its output energy, jet opening angle, and central BH characteristics. Even for the extreme BH parameters, some SGRBs require massive disks, which approach or exceed the limits in simulations. We suggest that there may exist alternative magnetohydrodynamic processes or some mechanisms increasing the neutrino emission to produce SGRBs with the reasonable BH parameters and disk mass.

We present the effects of the vertical convection on the structure and luminosity of the neutrino-dominated accretion flow (NDAF) around a stellar-mass black hole in spherical coordinates. We found that the convective energy transfer can suppress the radial advection in the NDAF, and that the density, temperature and opening angle are slightly changed. As a result, the neutrino luminosity and annihilation luminosity are increased, which is conducive to achieve the energy requirement of gamma-ray bursts.

There is a maximum for the gravity of a black hole in the vertical direction in the accretion disc. Outflows may probably be driven from the disc if the radiation flux of the disc is greater than a critical value corresponding to the maximal vertical gravity. We find that outflows are driven by the radiation force from the disc if the accretion rate is greater than the Eddington rate. The radiation of the disc is therefore limited by such outflows. The disc luminosity, L=L_Edd\propto ln mdot, at large-mdot cases. The Eddington ratio of the disc is ~3 for mdot~100, which is significantly lower than that of a conventional slim disc without outflows. This implies that the emission from some ultra-luminous X-ray sources with highly super Eddington luminosity should be Doppler beamed, or intermediate mass black holes are in these sources instead of stellar mass black holes. The spectra of the discs with outflows are saturated in the high frequency end provided mdot>2. We suggest that the saturated emission can be observed to estimate the masses of the black holes accreting at high rates, such as the narrow-line Seyfert galaxies, with the model calculations. The rate of the mass accreted by the black hole is always around the Eddington rate even if the mass accretion rate at the outer radius is very high, because most of the gas is removed into the outflows. This implies that the luminous quasars at high redshifts z>6 should have grown up through persistent accretion at a rate close to the Eddington rate.

We investigate the thermal stability of optically thin, two-temperature, radiative cooling-dominated accretion disks. Our linear analysis shows that the disk is thermally unstable without magnetic fields, which agrees with previous stability analysis on the Shapiro-Lightman-Eardley disk. By taking into account the effects of magnetic fields, however, we find that the disk can be or partly be thermally stable. Our results may be helpful to understand the outflows in optically thin flows. Moreover, such radiative cooling-dominated disks may provide a new explanation of the different behaviors between black hole and neutron star X-ray binaries on the radio/X-ray correlation.

Based on no-outflow assumption, we investigate steady state, axisymmetric, optically thin accretion flows in spherical coordinates. By comparing the vertically integrated advective cooling rate with the viscous heating rate, we find that the former is generally less than 30% of the latter, which indicates that the advective cooling itself cannot balance the viscous heating. As a consequence, for radiatively inefficient flows with low accretion rates such as $\dot M \la 10^{-3} \dot M_{Edd}$, where $\dot M_{Edd}$ is the Eddington accretion rate, the viscous heating rate will be larger than the sum of the advective cooling rate and the radiative cooling one. Thus, no thermal equilibrium can be established under the no-outflow assumption. We therefore argue that in such case outflows ought to occur and take away more than 70% of the thermal energy generated by viscous dissipation. Similarly, for optically thick flows with extremely large accretion rates such as $\dot M \ga 10 \dot M_{Edd}$, outflows should also occur due to the limited advection and the low efficiency of radiative cooling. Our results may help to understand the mechanism of outflows found in observations and numerical simulations.

We present the effects of the self-gravity on the vertical structure and neutrino luminosity of the neutrino-dominated accretion disks in cylindrical coordinates. It is found that significant changes of the structure appear in the outer region of the disk, especially for high accretion rates (e.g., $\ga 1 M_\odot~\rm s^{-1}$), and thus causes the slight increase of the neutrino luminosity. Furthermore, the gravitational instability of the disk is reviewed by the vertical distribution of the Toomre parameter \citep{Toomre1964}, which may account for the late-time flares in gamma-ray bursts and the extended emission in short-duration gamma-ray bursts.

The particular giant X-ray bump of GRB 121027A triggered by \emph{Swift} is quite different from the typical X-ray flares in gamma-ray bursts. There exhibit four parts of the observed structural variabilities in the rise and decay phase of the bump. Considering the quality of four parts of the data, we can only analyze the data from about 5300 s to about 6100 s in the bump using the stepwise filter correlation method (Gao et al. 2012), and find that the $86^{+5.9}_{-9.4}~\rm s$ periodic oscillation may exist, which is confirmed by the Lomb-Scargle method (Scargle 1982). Furthermore, a jet precession model (Liu et al. 2010) is proposed to account for such a variability.

GRB 130925A, composed of three gamma-ray emission episodes and a series of orderly flares, has been detected by $Swift$, Fermi, Konus-$Wind$, and $INTEGRAL$. If the third weakest gamma-ray episode can be considered as a giant flare, we find that after the second gamma-ray episode observed by $INTEGRAL$ located at about 2000 s, there exists a positive relation between the time intervals of the adjacent flares and the time since the episode. We suggest that the second gamma-ray episode and its flares originate from the resumption of the accretion process due to the fragments from the collapsar falling back, such a relation may be related to a hyperaccretion disk around a precessed black hole (BH). We propose that the origin and time evolution of the flares, and the approximately symmetrical temporal structure and spectral evolution of the single flare, can be explained well by the jet precession model \citep{Liu2010}. In addition, the mass and spin of the BH can be constrained, which indicates a stellar-mass, fast rotating BH located in the center of GRB 130925A.

The existence of intermediate-mass black holes (IMBHs) in globular clusters (GCs) remains a crucial problem. Searching IMBHs in GCs reveals a discrepancy between radio observations and dynamical modelings: the upper mass limits constrained by radio observations are systematically lower than that of dynamical modelings. One possibility for such a discrepancy is that, as we suggest in this work, there exist outflows in accretion flows. Our results indicate that, for most sources, current radio observations cannot rule out the possibility that IMBHs may exist in GCs. In addition, we adopt an \dot{M}-L_{R} relation to revisit this issue, which confirms the results obtained by the Fundamental Plane relation.

Rapid temporal variability has been widely observed in the light curves of gamma-ray bursts (GRBs). One possible mechanism for such variability is related to the relativistic eddies in the jet. In this paper, we include the contribution of the inter-eddy medium together with the eddies to the gamma-ray emission. We show that the gamma-ray emission can either lead or lag behind the observed synchrotron emission, where the latter originates in the inter-eddy medium and provides most of seed photons for producing gamma-ray emission through the inverse-Compton scattering. As a consequence, we argue that the lead/lag found in non-stationary short-lived light curves may not reveal the intrinsic lead/lag of different emission components. In addition, our results may explain the lead of gamma-ray emission with respect to optical emission observed in GRB 080319B.

The neutrino-cooled accretion disk, which was proposed to work as the central engine of gamma-ray bursts, encounters difficulty in interpreting the X-ray flares after the prompt gamma-ray emission. In this paper, the magnetic coupling between the inner disk and the central black hole is taken into consideration. For mass accretion rates around $0.001 \sim 0.1 M_{\sun} s^{-1}$, our results show that the luminosity of neutrino annihilation can be significantly enhanced due to the coupling effects. As a consequence, after the gamma-ray emission, a remnant disk with mass $M_{disk} \la 0.5 M_{\sun}$ may power most of the observed X-ray flares with the rest frame duration less than 100 seconds. In addition, a comparison between the magnetic coupling process and the Blandford-Znajek mechanism is shown on the extraction of black hole rotational energy.

Neutrino-dominated accretion flows (NDAFs) around rotating stellar-mass black holes are plausible candidates for the central engines of gamma-ray bursts (GRBs). We investigate one-dimensional global solutions of NDAFs, taking account of general relativity in Kerr metric, neutrino physics and nucleosynthesis more precisely than previous works. We calculate sixteen solutions with different characterized accretion rates and black hole spins to exhibit the radial distributions of various physical properties in NDAFs. \iffalse We find that the gas pressure and the neutrino cooling always become dominated in the inner region for large accretion rate, $\gtrsim 0.1 M_\odot$ $\rm s^{-1}$. \fi We confirm that the electron degeneracy has important effects in NDAFs and we find that the electron fraction is about 0.46 in the outer region for all the sixteen solutions. From the perspective of the mass fraction, free nucleons, $\rm ^4He$, and $\rm ^{56}Fe$ dominate in the inner, middle, and outer region, respectively. The influence of neutrino trapping on the annihilation is of importance for the superhigh accretion ($\dot M=10M_{\odot} \rm s^{-1}$) and most of the sixteen solutions have an adequate annihilation luminosity for GRBs.

We investigate the vertical structure and elements distribution of neutrino-dominated accretion flows around black holes in spherical coordinates with the reasonable nuclear statistical equilibrium. According our calculations, heavy nuclei tend to be produced in a thin region near the disk surface, whose mass fractions are primarily determined by the accretion rate and the vertical distribution of temperature and density. In this thin region, we find that $^{56}\rm Ni$ is dominant for the flow with low accretion rate (e.g., $0.05 $$M_{\odot}$ $\rm s^{-1}$) but $^{56}\rm Fe$ is dominant for the high counterpart (e.g., $1 M_{\odot}$ $\rm s^{-1}$). The dominant $^{56}\rm Ni$ in the special region may provide a clue to understand the bumps in the optical light curve of core-collapse supernovae.

We study the fluctuations of standard thin accretion disks by linear analysis of the time-dependent energy equation together with the vertical hydrostatic equilibrium and the equation of state. We show that some of the simulation results in Hirose et al. (2009b), such as the time delay, the relationship of power spectra, and the correlation between magnetic energy and radiation energy, can be well understood by our analytic results.

Soft extended emission (EE) following initial hard spikes up to 100 seconds was observed with {\em Swift}/BAT for about half of short-type gamma-ray bursts (SGRBs). This challenges the conversional central engine models of SGRBs, i.e., compact star merger models. In the framework of the black hole-neutron star merger models, we study the roles of the radial angular momentum transfer in the disk and the magnetic barrier around the black hole for the activity of SGRB central engines. We show that the radial angular momentum transfer may significantly prolong the lifetime of the accretion process and multiple episodes may be switched by the magnetic barrier. Our numerical calculations based on the models of the neutrino-dominated accretion flows suggest that the disk mass is critical for producing the observed EE. In case of the mass being $\sim 0.8M_{\odot}$, our model can reproduce the observed timescale and luminosity of both the main and EE episodes in a reasonable parameter set. The predicted luminosity of the EE component is lower than the observed EE with about one order of magnitude and the timescale is shorter than 20 seconds if the disk mass being $\sim 0.2M_{\odot}$. {\em Swift}/BAT-like instruments may be not sensitive enough to detect the EE component in this case. We argue that the EE component would be a probe for merger process and disk formation for compact star mergers.

By taking into account the local energy balance per unit volume between the viscous heating and the advective cooling plus the radiative cooling, we investigate the vertical structure of radiation pressure-supported accretion disks in spherical coordinates. Our solutions show that the photosphere of the disk is close to the polar axis and therefore the disk seems to be extremely thick. However, the profile of density implies that most of the accreted matter exists in a moderate range around the equatorial plane. We show that the well-known polytropic relation between the pressure and the density is unsuitable for describing the vertical structure of radiation pressure-supported disks. More importantly, we find that the energy advection is significant even for slightly sub-Eddington accretion disks. We argue that the non-negligible advection may help to understand why the standard thin disk model is likely to be inaccurate above \sim 0.3 Eddington luminosity, which was found by some works on the black hole spin measurement. Furthermore, the solutions satisfy the Solberg-Hoiland conditions, which indicates the disk to be convectively stable. In addition, we discuss the possible link between our disk model and ultraluminous X-ray sources.

The physical nature of gamma-ray bursts (GRBs) are believed to involve an ultra-relativistic jet. The observed complex structure of light curves motivate the idea of jet precession. In this work, we study the gravitational waves of jet precession based on neutrino-dominated accretion disks around black holes, which may account for the central engine of GRBs. In our model, the jet and the inner part of the disk may precess along with the black hole, which is driven by the outer part of the disk. Gravitational waves are therefore expected to be significant from this black hole-inner disk precession system. By comparing our numerical results with the sensitivity of some detectors, we find that it is possible for DECIGO and BBO to detect such gravitational waves, particularly for GRBs in the Local Group.

In accretion systems, outflows may have significant influence on the luminosity fluctuations. In this paper, following the Lyubarskii's general scheme, we revisit the power spectral density of luminosity fluctuations by taking into account the role of outflows. Our analysis is based on the assumption that the coupling between the local outflow and inflow is weak on the accretion rate fluctuations. We find that, for the inflow mass accretion rate $\dot M \propto r^{s}$, the power spectrum of flicker noise component will present a power-law distribution $p(f) \propto f^{-(1+4s/3)}$ for advection-dominated flows. We also obtain descriptions of $p(f)$ for both standard thin discs and neutrino-cooled discs, which show that the power-law index of a neutrino-cooled disc is generally larger than that of a photon-cooled disc. Furthermore, the obtained relationship between $p(f)$ and $s$ indicates the possibility of evaluating the strength of outflows by the power spectrum in X-ray binaries and gamma-ray bursts. In addition, we discuss the possible influence of the outflow-inflow coupling on our results.

We revisit the vertical structure of neutrino dominated accretion flows (NDAFs) in spherical coordinates with a new boundary condition based on the mechanical equilibrium. The solutions show that NDAF is significantly thick. The Bernoulli parameter and neutrino trapping are determined by the mass accretion rate and the viscosity parameter. According to the distribution of the Bernoulli parameter, the possible outflow may appear in the outer region of the disk. The neutrino trapping can essentially affect the neutrino radiation luminosity. The vertical structure of NDAF is like a "sandwich", and the multilayer accretion may account for the flares in gamma-ray bursts.

We find that the ambiguity term of approximate Killing vector field is responsible for the entropy production term. Without the ambiguity term, pure Einstein theory and $f(R)$ satisfy the relation of thermodynamic equilibrium. Considering such an ambiguity term of approximate Killing vector field, we can get the entropy production term and the entropy in $f(R)$ with a form defined by Jacobson. In pure Einstein theory, the shear term is the only geometric contribution of entropy production term, while in $f(R)$ it can also contribute. We believe our approach and conclusion can be generalized to other gravity theory.

The stress evolution process is taken into account in the linear stability analysis of standard thin accretion discs. We find that the growth rate of thermally unstable modes can decrease significantly owing to the stress delay, which may help to understand the quasi-periodic variability of GRS 1915+105. We also discuss possible application of stress evolution to the stability of Shapiro-Lightman-Eardley disc.

The standard thin disk model predicts that when the accretion rate is over a small fraction of the Eddington rate, which corresponds to $L \ga 0.06 L_{Edd}$, the inner region of the disk is radiation-pressure-dominated and thermally unstable. However, observations of the high/soft state of black hole X-ray binaries with luminosity well within this regime ($0.01L_{Edd} \la L \la 0.5L_{Edd}$) indicate that the disk has very little variability, i.e., quite stable. Recent radiation magnetohydrodynamic simulations of a vertically stratified shearing box have confirmed the absence of the thermal instability. In this paper, we revisit the thermal stability by linear analysis, taking into account the role of magnetic field in the accretion flow. By assuming that the field responses negatively to a positive temperature perturbation, we find that the threshold of accretion rate above which the disk becomes thermally unstable increases significantly compared with the case of not considering the role of magnetic field. This accounts for the stability of the observed sources with high luminosities. Our model also presents a possible explanation as to why only GRS 1915+105 seems to show thermally unstable behavior. This peculiar source holds the highest accretion rate (or luminosity) among the known high state sources, which is well above the accretion rate threshold of the instability.

We propose a new holographic program of gravity in which we introduce a surface stress tensor. Our proposal differs from Verlinde's in several aspects. First, we use an open or a closed screen, a temperature is not necessary but a surface energy density and pressure are introduced. The surface stress tensor is proportional to the extrinsic curvature. The energy we use is Brown-York energy and the equipartiton theorem is violated by a non-vanishing surface pressure. We discuss holographic thermodynamics of a gas of weak gravity and find a chemical potential, and show that Verlinde's program does not lead to a reasonable thermodynamics. The holographic entropy is similar to the Bekenstein entropy bound.

A model of jet precession driven by a neutrino-cooled disc around a spinning black hole is present in order to explain the temporal structure and spectral evolution of gamma-ray bursts (GRBs). The differential rotation of the outer part of a neutrino dominated accretion disc may result in precession of the inner part of the disc and the central black hole, hence drives a precessed jet via neutrino annihilation around the inner part of the disc. Both analytic and numeric results for our model are present. Our calculations show that a black hole-accretion disk system with black hole mass $M \simeq 3.66 M_\odot$, accretion rate $\dot{M} \simeq 0.54 M_\odot \rm s^{-1}$, spin parameter $a=0.9$ and viscosity parameter $\alpha=0.01$ may drive a precessed jet with period P=1 s and luminosity $L=10^{51}$ erg s$^{-1}$, corresponding to the scenario for long GRBs. A precessed jet with $P=0.1$s and $L=10^{50}$ erg s$^{-1}$ may be powered by a system with $M \simeq 5.59 M_\odot$, $\dot{M} \simeq 0.74 M_\odot \rm s^{-1}$, $a=0.1$, and $\alpha=0.01$, possibly being responsible for the short GRBs. Both the temporal and spectral evolution in GRB pulse may explained with our model. GRB central engines likely power a precessed jet driven by a neutrino-cooled disc. The global GRB lightcurves thus could be modulated by the jet precession during the accretion timescale of the GRB central engine. Both the temporal and spectral evolution in GRB pulse may be due to an viewing effect due to the jet precession.

We revisit the vertical structure of neutrino-dominated accretion flows in spherical coordinates. We stress that the flow should be geometrically thick when advection becomes dominant. In our calculation, the luminosity of neutrino annihilation is enhanced by one or two orders of magnitude. The empty funnel along the rotation axis can naturally explain the neutrino annihilable ejection.

We revisit the vertical structure of black hole accretion disks in spherical coordinates. By comparing the advective cooling with the viscous heating, we show that advection-dominated disks are geometrically thick, i.e., with the half-opening angle > 2\pi/5, rather than slim as supposed previously in the literature.

We provide a detailed analysis of 12 XMM observations of the nearest persistent extragalactic ultraluminous X-ray source (ULX), M33 X-8. No significant spectral evolution is detected between the observations, therefore we combine the individual observations to increase the signal-to-noise ratio for spectral fitting. The combined spectra are best fitted by the self-consistent p-free disk plus power-law component model with p = 0.571_{-0.030}^{+0.032}, kT_{in} = 1.38_{-0.08}^{+0.09} keV, and the flux ratio of the p-free disk component to the power-law component being 0.63:0.37 in the 0.3 -- 10 keV band. The fitting indicates that the black hole in M33 X-8 is of \sim 10 M_{\odot} and accretes at a super-Eddington rate (\sim 1.5 L_{Edd}), and the phase of the accretion disk is close to a slim disk (p = 0.5). We report, for the first time, that an extra power-law component is required in addition to the p-free disk model for ULXs. In super-Eddington cases, the power-law component may possibly result from the optically thin inner region f the disk or a comptonized corona similar to that of a standard thin disk.

We take the vertical distribution of the radial and azimuthal velocity into account in spherical coordinates, and find that the analytic relation c_{s0}/(v_K \Theta) = [(\gamma -1)/(2\gamma)]^{1/2} is valid for both geometrically thin and thick accretion flows, where c_{s0} is the sound speed on the equatorial plane, v_K is the Keplerian velocity, \Theta is the half-opening angle of the flow, and \gamma is the adiabatic index.

We show that there exists a maximal possible accretion rate, beyond which global slim disk solutions cannot be constructed because in the vertical direction the gravitational force would be unable to balance the pressure force to gather the accreted matter. The principle for this restriction is the same as that for the Eddington luminosity and the corresponding critical accretion rate, which were derived for spherical accretion by considering the same force balance in the radial direction. If the assumption of hydrostatic equilibrium is waived and vertical motion is included, this restriction may become even more serious as the value of the maximal possible accretion rate becomes smaller. Previous understanding in the literature that global slim disk solutions could stand for any large accretion rates is due to the overestimation of the vertical gravitational force by using an approximate potential. For accretion flows with large accretion rates at large radii, outflows seem unavoidable in order for the accretion flow to reduce the accretion rate and follow a global solution till the central black hole.

We present a unified description of all the three known classes of optically thick accretion disks around black holes, namely Shakura-Sunyaev disks, slim disks, and neutrino-dominated accretion flows (NDAFs). It is found that NDAFs have both a maximal and a minimal possible mass accretion rate at their each radius. This may be suggestive of an interpretation for the origin of X-ray flares observed in gamma-ray bursts.

We show that when the gravitational force is correctly calculated in dealing with the vertical hydrostatic equilibrium of black hole accretion disks, the relationship that is valid for geometrically thin disks, i.e., $c_s/\Omega_K H =$ constant, where $c_s$ is the sound speed, $\Omega_K$ is the Keplerian angular velocity, and $H$ is the half-thickness of the disk, does not hold for slim disks. More importantly, by adopting the correct vertical gravitational force in studies of thermal equilibrium solutions, we find that there exists a maximally possible accretion rate for each radius in the outer region of optically thick accretion flows, so that only the inner region of these flows can possibly take the form of slim disks, and strong outflows from the outer region are required to reduce the accretion rate in order for slim disks to be realized.

Neutrino-cooled hyperaccretion disks around stellar mass black holes are plausible candidates for the central engine of gamma-ray bursts. We calculate the one-dimensional structure and the annihilation luminosity of such disks. The neutrino optical depth is of crucial importance in determining the neutrino cooling rate and is in turn dependent on the electron fraction, the free nucleon fraction, and the electron degeneracy, with given density and temperature of the disk matter. We construct a bridging formula for the electron fraction that works for various neutrino optical depths, and give exact definitions for the free proton fraction and free neutron fraction. We show that the electron degeneracy has important effects in the sense that it enlarges the absorption optical depth for neutrinos, and it along with the neutronization processes favored by high temperature cause the electron fraction to drop to be below 0.1 in the inner region of the disk. The resulting neutrino annihilation luminosity is considerably reduced comparing with that obtained in previous works where the electron degeneracy was not considered and the electron fraction was simply taken to be 0.5, but it is still likely to be adequate for gamma-ray bursts, and it is ejected mainly from the inner region of the disk and has an anisotropic distribution.

We first refine the fixed concept in the literature that the usage of the Newtonian potential in studies of black hole accretion is invalid and the general relativistic effect must be considered. As our main results, we then show that the energy released by neutrino annihilation in neutrino-dominated accretion flows is sufficient for gamma-ray bursts when the contribution from the optically thick region of the flow is included, and that in the optically thick region advection does not necessarily dominate over neutrino cooling because the advection factor is relevant to the geometrical depth rather than the optical depth of the flow.

We investigate the linear stability of a shocked accretion flow onto a black hole in the adiabatic limit. Our linear analyses and numerical calculations show that, despite the post-shock deceleration, the shock is generally unstable to non-axisymmetric perturbations. The simulation results of Molteni, T\'oth & Kuznetsov can be well explained by our linear eigenmodes. The mechanism of this instability is confirmed to be based on the cycle of acoustic waves between the corotation radius and the shock. We obtain an analytical formula to calculate the oscillation period from the physical parameters of the flow. We argue that the quasi-periodic oscillation should be a common phenomenon in accretion flows with angular momentum.

We obtain global solutions of radiatively inefficiently accretion flows around black holes. Whether and where convection develops in a flow are self-consistently determined with the mixing-length theory. The solutions can be divided into three types according to the strength of normal viscosity. Type I solution corresponds to large viscosity parameter $\alpha \ga 0.1$, which is purely advection-dominated and with no convection, and has been extensively studied in the literature. Type II solution is for moderate $\alpha \sim 0.01$, which has a three-zone structure. The inner zone is advection-dominated, the middle zone is convection-dominated and ranges from a few tens to a few thousands of gravitational radii, and the outer zone is convectively stable and matches outward a Keplerian disc. The net energy flux throughout the flow is inward as in type I solution. Type III solution which is for small $\alpha \la 0.001$ consists of two zones as Abramowicz et al. suggested previously: an inner advection-dominated zone and an outer convection-dominated zone, separated at a radius of a few tens of gravitational radii. This type of solution has an outward net energy flux. In both type II and III solutions the radial density profile is between the 1/2 law of self-similar convection-dominated accretion flow model and the 3/2 law of self-similar advection-dominated accretion flow model, and the efficiency of energy release is found to be extremely low. Our results are in good agreement with those of recent numerical simulations.

We use the standard Runge-Kutta method to solve the set of basic equations describing black hole accretion flows composed of two-temperature plasma. We do not invoke any extra energy transport mechanism such as thermal conduction and do not specify any ad hoc outer boundary condition for the advection-dominated accretion flow (ADAF) solution. We find that in the case of high viscosity and non-zero radiative cooling, the ADAF solution can have an asymptotic approach to the Shakura-Sunyaev disk (SSD) solution, and the SSD-ADAF transition radius is close to the central black hole. Our results further prove the mechanism of thermal instability-triggered SSD-ADAF transition suggested previously by Takeuchi & Mineshige and Gu & Lu.

The linear stability of a shocked isothermal accretion flow onto a black hole is investigated in the inviscid limit. The outer shock solution, which was previously found to be stable with respect to axisymmetric perturbations, is, however, generally unstable to non-axisymmetric ones. Eigenmodes and growth rates are obtained by numerical integration of the linearized equations. The mechanism of this instability is based on the cycle of acoustic waves between their corotation radius and the shock. It is a form of the Papaloizou-Pringle instability, modified by advection and the presence of the shock. As such it can be generalized to non isothermal shocked accretion flows. Blobs and vortices are generated by the shock as a by-product of the instability.

We show that, unlike the results presented previously in the literature, the transition from an outer Shakura-Sunyaev disk (SSD) to an advection-dominated accretion flow (ADAF) is possible for large values of the viscosity parameter $\alpha >0.5$. The transition is triggered by thermal instability of a radiation-pressure-supported SSD. The transition radius is close to the central black hole. We confirm our qualitative prediction by actually constructing global bimodal SSD-ADAF solutions.

We numerically solve the set of dynamical equations describing advection-dominated accretion flows (ADAF) around black holes, using a method similar to that of Chakrabarti (1996a). We choose the sonic radius of the flow $R_s$ and the integration constant in angular momentum equation $j$ as free parameters, and integrate the equations from the sonic point inwards to see if the solution can extend supersonically to the black hole horizon, and outwards to see if and where an acceptable outer boundary of the flow can be found. We recover the ADAF-thin disk solution constructed in Narayan, Kato, & Honma (1997, NKH97), an representative paper of the previous works on global ADAF solutions, although in that paper an apparently very different procedure was adopted. We obtain a complete picture in the form of $R_s-j$ parameter space which sums up the situation of ADAF solution at a glance. For comparison we also present the distribution of global solutions for inviscid flows in the $R_s-j$ space, which supports the view that there should be some similarities between the dynamical behavior of ADAF and that of adiabatic flows, and that there should be a continuous change from the properties of viscous flows to those of inviscid ones.