The authors of this white paper met on 16-17 January 2020 at the New Jersey Institute of Technology, Newark, NJ, for a 2-day workshop that brought together a group of heliophysicists, data providers, expert modelers, and computer/data scientists. Their objective was to discuss critical developments and prospects of the application of machine and/or deep learning techniques for data analysis, modeling and forecasting in Heliophysics, and to shape a strategy for further developments in the field. The workshop combined a set of plenary sessions featuring invited introductory talks interleaved with a set of open discussion sessions. The outcome of the discussion is encapsulated in this white paper that also features a top-level list of recommendations agreed by participants.

Absorption spectroscopy is widely used in sensing and astronomy to understand molecular compositions on microscopic to cosmological scales. However, typical dispersive spectroscopic techniques require multichannel detection, fundamentally limiting the ability to detect extremely weak signals when compared to direct photometric methods. We report the realization of direct spectral molecular detection using a silicon nanophotonic waveguide resonator, obviating dispersive spectral acquisition. We use a thermally tunable silicon ring resonator with a transmission spectrum matched and cross-correlated to the quasi-periodic vibronic absorption lines of hydrogen cyanide. We show that the correlation peak amplitude is proportional to the number of overlapping ring resonances and gas lines, and that molecular specificity is obtained from the phase of the correlation signal in a single detection channel. Our results demonstrate on-chip correlation spectroscopy that is less restricted by the signal-to-noise penalty of other spectroscopic approaches, enabling the detection of faint spectral signatures.

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{\gamma,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{\gamma,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB.

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.

Sgr A* represents a unique laboratory for the detailed study of accretion processes around a low-luminosity supermassive black hole (SMBH). Recent X-ray observations have allowed for spatially resolved modeling of the emission from the accretion flow around the SMBH, placing tight constraints on the flux and spectral shape of the accretion from the inner region with $r<10^3~r_g$, where $r_g$ is the gravitational radius of the black hole. We present here the first modeling of the multi-band spectral energy distribution (SED) of this inner region to better constrain the physical condition of the innermost accretion flow. Our modeling uses the Markov chain Monte Carlo (MCMC) method to fit the SED, accounting for the limitations on the accretion rate at the outer radius of $10^3~r_g$ from the earlier works and the domination of the accretion flow within $30~r_g$ to the sub-mm bump. It is found that the fitting results of the outflow index could be very different. If only the most luminous part of the SED, the sub-mm bump, is considered, the outflow index is about 0, while if low-frequency radio data and X-ray data are also included, the outflow index could be 0.37 or even higher. The great difference of the fitting results indicates that the outflow index should be variable along radius, with a strong outflow in the outer region and a weak outflow in the innermost region. Such weak outflow agrees with numerical simulations and is possible to explain the multi-band SED even better.

The detection of highly ionized metal absorption lines in the X-ray spectra of the Galactic X-ray binaries (XRBs) implies the distribution of hot gas along the sightline toward the background sources. However, the origin of this hot gas is still unclear: it can arise in the hot interstellar medium (ISM), or is intrinsic to the XRBs. In this paper, we present an XMM-Newton survey of the OVII absorption lines in the spectra of Galactic XRBs. A total of 33 XRBs were selected, with 29 low mass XRBs and 4 high mass XRBs. At more than 3$\sigma$ threshold, OVII absorption line was detected in 16 targets, among which 4 were newly discovered in this work. The average line equivalent width (EW) is centered around $\sim$ 20 m\AA. Additionally we do not find strong correlations between the OVII EWs and the Galactic neutral absorption N$_{\rm H}$, the Galactic coordinates or the distance of background targets. Such non-correlation may suggest the contamination of the circumstellar material, or the lack of constrains on the line Doppler-b parameter. We also find regardless the direction of the XRBs, the OVII absorption lines were always detected when the flux of the background XRBs reach a certain level, suggesting a uniform distribution of this hot gas. We estimate a ratio of 0.004 --- 0.4 between hot and neutral phase of the interstellar medium. This is the second paper in the series following Fang et al. 2015, in which we focused on the local OVII absorption lines detected in the background AGN spectra. Detailed modeling of the hot ISM distribution will be investigated in a future paper.

The clumpy disc, or inhomogeneous accretion flow, has been proposed to explain the properties of accreting black hole systems. However, the observational evidences remain to be explored. In this work, we calculate the profiles of Fe K{\alpha} lines emitted from the inhomogeneous accretion flow through the ray-tracing technique, in order to find possible observable signals of the clumps. Compared with the skewed doublepeaked profile of the continuous standard accretion disc, the lines show a multi-peak structure when the emissivity index is not very steep. The peaks and wings are affected by the position and size of the cold clumps. When the clump is small and is located in the innermost region, due to the significant gravitational redshift, the blue wing can overlap with the red wing of the outer cold disc/clump, forming a fake peak or greatly enhancing the red peak. Given high enough resolution, it is easier to constrain the clumps around the supermassive black holes than the clumps in stellar mass black holes due to the thermal Doppler effect.

The model of inhomogeneous accretion flow, in which cold clumps are surrounded by hot gas or corona, has been proposed to explain the spectral features of black hole X-ray binaries (BHXBs). In this work, we try to find possible observational features in the continuum that can indicate the existence of clumps. The spectra of inhomogeneous accretion flow are calculated via the Monte Carlo method. Since the corresponding accretion flow is unsteady and complex, the accretion flow is described by a set of free parameters, the ranges of which can include the real cases. The influences of the parameters are investigated. It is found that the thermal component of the spectra deviates from the multi-color black body spectra in the middle power-law part. On the one hand, a warp appears due to the gap region between the clumps and the outer cold disk, and on the other hand, the slope of the line connecting the thermal peaks deviates from 4/3. The warp feature, as well as the correlation between the thermal peak at higher frequency and the spectral index, are possible to indicate the existence of clumps, and are worthy of further investigation with more self-consistent models.

The active galactic nucleus (AGN) NGC 7213 shows a complex correlation between the monochromatic radio luminosity $L_R$ and the 2--10 keV X-ray luminosity $L_X$, i.e. the correlation is unusually weak with $p\sim 0$ (in the form $L_R\propto L_X^p$) when $L_X$ is below a critical luminosity, and steep with $p>1$ when $L_X$ is above that luminosity. Such a hybrid correlation in individual AGNs is unexpected as it deviates from the Fundamental Plane of AGN activity. Interestingly, a similar correlation pattern is observed in the black hole X-ray binary H1743--322, where it has been modelled by switching between different modes of accretion. We propose that the flat $L_R$--$L_X$ correlation of NGC 7213 is due to the presence of a luminous hot accretion flow, an accretion model whose radiative efficiency is sensitive to the accretion rate. Given the low luminosity of the source, $L_X\sim 10^{-4}$ of the Eddington luminosity, the viscosity parameter is determined to be small, $\alpha\approx 0.01$. We also modelled the broad-band spectrum from radio to $\gamma$-rays, the time lag between the radio and X-ray light curves, and the implied size and the Lorentz factor of the radio jet. We predict that NGC 7213 will enter into a two-phase accretion regime when $L_X > 1.5 \times 10^{42}\, {\rm erg\,s^{-1}}$. When this happens, we predict a softening of the X-ray spectrum with the increasing flux and a steep radio/X-ray correlation.

Neutrino-dominated accretion flows (NDAFs) around rotating stellar-mass black holes (BHs) have been theorized as the central engine of relativistic jets launched in massive star core collapse events or compact star mergers. In this work, we calculate the electron neutrino/anti-neutrino spectra of NDAFs by fully taking into account the general relativistic effects, and investigate the effects of viewing angle, BH spin, and mass accretion rate on the results. We show that even though a typical NDAF has a neutrino luminosity lower than that of a typical supernova (SN), it can reach $10^{50}-10^{51}~{\rm erg~s^{-1}}$ peaking at $\sim 10$ MeV, making them potentially detectable with the upcoming sensitive MeV neutrino detectors if they are close enough to Earth. Based on the observed GRB event rate in the local universe and requiring that at least 3 neutrinos are detected to claim a detection, we estimate a detection rate up to $\sim$ (0.10-0.25) per century for GRB-related NDAFs by the Hyper-Kamiokande (Hyper-K) detector if one neglects neutrino oscillation. If one assumes that all Type Ib/c SNe have an engine-driven NDAF, the Hyper-K detection rate would be $\sim$ (1-3) per century. By considering neutrino oscillations, the detection rate may decrease by a factor of 2-3. Detecting one such event would establish the observational evidence of NDAFs in the universe.

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.

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.

Highly ionized, z=0 metal absorption lines detected in the X-ray spectra of background active galactic nuclei (AGNs) provide an effective method to probe the hot ($T\sim10^6$ K) gas and its metal content in and around the Milky Way. We present an all-sky survey of the $K_{\alpha}$ transition of the local O VII absorption lines obtained by Voigt-profile fitting archival XMM-Newton observations. A total of 43 AGNs were selected, among which 12 are BL Lac-type AGNs, and the rest are Seyfert 1 galaxies. At above the $3\sigma$ level the local O VII absorption lines were detected in 21 AGNs, among which 7 were newly discovered in this work. The sky covering fraction, defined as the ratio between the number of detections and the sample size, increases from at about 40% for all targets to 100% for the brightest targets, suggesting a uniform distribution of the O VII absorbers. We correlate the line equivalent width with the Galactic coordinates and do not find any strong correlations between these quantities. Some AGNs have warm absorbers that may complicate the analysis of the local X-ray absorber since the recession velocity can be compensated by the outflow velocity, especially for the nearby targets. We discuss the potential impact of the warm absorbers on our analysis. A comprehensive theoretical modelling of the X-ray absorbers will be presented in a later paper.

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.

The dynamical and radiative properties of the quiescent state (X-ray luminosity $<10^{34}\ {\rm erg s^{-1} cm^{-2}}$) of black hole X-ray transients (BHXTs) remains unclear, mainly because of low-luminosity and poor data quantity. We demonstrate that, the simultaneous multi-wavelength (including radio, optical, ultraviolet and X-ray bands) spectrum of V404 Cyg in its bright quiescent state can be well described by the radiation from the companion star and more importantly, the compact jet. Neither the outer thin disc nor the inner hot accretion flow is important in the total spectrum. Together with several additional recent findings, i.e. the power-law X-ray spectrum and the constant X-ray spectral shape (or constant photon index) in contrast to the dramatic change in the X-ray luminosity, we argue the quiescent state spectrum of BHXTs is actually jet-dominated. Observational features consistent with this jet model are also discussed as supporting evidences.

The correlation between the kinetic power $P_{\rm jet}$ and intrinsic bolometric luminosity $L_{\rm jet}$ of jets may reveal the underlying jet physics in various black hole systems. Based on the recent work by Nemmen et al. (2012), we re-investigate this correlation with additional sources of black-hole X-ray binaries (BXBs) in hard/quiescent states and low-luminosity active galactic nuclei (LLAGNs). The new sample includes 29 sets of data from 7 BXBs and 20 LLAGNs, with $P_{\rm jet}$ and $L_{\rm jet}$ being derived from spectral modeling of the quasi-simultaneous multi-band spectra under the accretion-jet scenario. Compared to previous works, the range of luminosity is now enlarged to more than $20$ decades, i.e. from $10^{31}{\rm erg/s}$ to $10^{52}{\rm erg/s}$, which allows for better constraining of the correlation. One notable result is that the jets in BXBs and LLAGNs almost follow the same $P_{\rm jet} - L_{\rm jet}$ correlation that was obtained from blazars and gamma-ray bursts (GRBs). The slope indices we derived are $1.03\pm0.01$ for the whole sample, $0.85\pm0.06$ for BXB subsample, $0.71\pm0.11$ for LLAGN subsample, and $1.01\pm0.05$ for the LLAGN-blazar subsample, respectively. The correlation index around unit implies the independence of jet efficiency on the luminosity or kinetic power. Our results may further support the hypothesis that similar physical processes exist in the jets of various black hole systems.

We propose a disc-corona model in which a geometrically thin, optically thick disc surrounds a Kerr black hole, and magnetic fields exert a time-steady torque on the inner edge of the accretion disc. The analytical expression of the total gravitational power is derived from the thin-disc dynamics equations by using this new boundary condition. It is shown that the magnetic torque can considerably enhance the amount of energy released in the disc-corona system. Furthermore, the global solutions of this disc-corona system are obtained numerically. We find that the fraction of the power dissipated into the corona in the total for such disc-corona system increases with the increasing dimensionless black hole spin parameter $a_\ast $, but is insensitive on the $\Delta \varepsilon $ which is the additional radiative efficiency parameter relevant to magnetic torque, for $\Delta\varepsilon > 1$. In addition, the emerged spectra from this disc-corona system are simulated by using Monte-Carlo method, and the effect of the different parameters on the output spectra is discussed.

The correlation between radio and X-ray luminosity in the hard state of black hole X-ray binaries is important for unveiling the relation between the accretion flow and the jets. In this paper, we have modeled the quasi-simultaneous multi-band observations of a recently discovered transient X-ray source, IGR J17177-3656. It is found that the source is probably an outlier following the steep radio/X-ray correlation rather than an outlier in the transition region as suggested by Paizis et al. (2011). It is also found that the multi-band spectrum can be successfully modeled by the luminous hot accretion flow (LHAF) but less likely by the advection dominated accretion flow (ADAF). Our results support the point that LHAF can explain the steep radio/X-ray correlation.

We present calculations of expected continuum emissions from Sedov-Taylor phase Type Ia supernova remnants (SNRs), using the energy spectra of cosmic ray (CR) electrons and protons from nonlinear diffusive shock acceleration (DSA) simulations. A new, general-purpose radiative process code, Cosmicp, was employed to calculate the radiation expected from CR electrons and protons and their secondary products. These radio, X-ray and gamma-ray emissions are generally consistent with current observations of Type Ia SNRs. The emissions from electrons in these models dominate the radio through X-ray bands. Decays of \pi^0 s from p-p collisions mostly dominate the gamma-ray range, although for a hot, low density ISM case (n_{ISM}=0.003 cm^{-3}), the pion decay contribution is reduced sufficiently to reveal the inverse Compton contribution to TeV gamma-rays. In addition, we present simple scalings for the contributing emission processes to allow a crude exploration of model parameter space, enabling these results to be used more broadly. We also discuss the radial surface brightness profiles expected for these model SNRs in the X-ray and gamma-ray bands.

When we model black hole accretion sources such as active galactic nuclei and black hole X-ray binaries as advection-dominated accretion flows (ADAFs), it is neccesary to use the global solution to the equations rather than the simpler self-similar solution, since the latter is inaccurate in the region near the black hole where most of the radiation is emitted. However, technically, it is a difficult task to calculate the global solution because of the transonic nature of the flow, which makes it a two-point boundary value problem. In this paper we propose a simplified approach for calculating the global ADAF solution. We replace the radial momentum equation by a simple algebraic relation between the angular velocity of the gas and the Keplerian angular velocity, while keeping all other equations unchanged. It is then easy to solve the differential energy equations to obtain an approximate global solution. By adjusting the free parameters, we find that for almost any accretion rate and for $\alpha=0.1-0.3$ we can get good simplified global solutions. The predicted spectra from the approximate solutions are very close to the spectra obtained from the true global solutions.

A large-scale closed magnetic field can transfer angular momentum and energy between a black hole (BH) and its surrounding accretion flow. We investigate the effects of this magnetic coupling (MC) process on the dynamics of a hot accretion flow (e.g., an advection dominated accretion flow, hereafter ADAF). The energy and angular momentum fluxes transported by the magnetic field are derived by an equivalent circuit approach. For a rapidly rotating BH, it is found that the radial velocity and the electron temperature of the accretion flow decrease, whereas the ion temperature and the surface density increase. The significance of the MC effects depends on the value of the viscous parameter \alpha. The effects are obvious for \alpha=0.3 but nearly ignorable for \alpha=0.1. For a BH with specific angular momentum, a_*=0.9, and \alpha=0.3, we find that for reasonable parameters the radiative efficiency of a hot accretion flow can be increased by about 30%.

A model of magnetic coupling (MC) of a rotating black hole (BH) with advection- dominated accretion flow (MCADAF) is proposed. It turns out that MCADAF providers a natural explanation for the transition radius between ADAF and SSD, and could be used to interpret the highest luminosity of GX 339-4 in hard-state. A very steep emissivity index can be produced in the innermost part of the MCADAF,which is consistent with the recent XMM-Newton observations of the nearby bright Seyfert 1 galaxy MCG-6-30-15 and with two X-ray binaries (XRBs): XTE J1655-500 and GX 339-4. In addition, we estimate the BH spins in Seyfert 1 galaxy MCG-6-30-15 and in the two XRBs based on this model.

Recently, much attention has been paid to the magnetic coupling (MC) process, which is supported by very high emissivity indexes observed in Seyfert 1 galaxy MCG-6-30-15 and GBHC XTE J1650-500. But the rotational energy transferred from a black hole is simply assumed to be radiated away from the surrounding accretion disk in black-body spectrum, which is obviously not consistent with the observed hard power-law X-ray spectra. We intend to introduce corona into the MC model to make it more compatible with the observations. We describe the model and the procedure of a simplified Monte Carlo simulation, compare the output spectra in the cases with and without the MC effects, and discuss the influences of three parameters involved in the MC process on the output spectra. It is shown that the MC process augments radiation fluxes in the UV or X-ray band. The emergent spectrum is affected by the BH spin and magnetic field strength at the BH horizon, while it is almost unaffected by the radial profile of the magnetic field at the disk. Introducing corona into the MC model will improve the fitting of the output spectra from AGNs and GBHCs.

We propose a model to explain a puzzling 3:2 frequency ratio of high frequency quasi-periodic oscillations (HFQPOs) in black hole (BH) X-ray binaries, GRO J1655-40, GRS 1915+105 and XTE J1550-564. In our model a non-axisymmetric magnetic coupling (MC) of a rotating black hole (BH) with its surrounding accretion disc coexists with the Blandford-Znajek (BZ) process. The upper frequency is fitted by a rotating hotspot near the inner edge of the disc, which is produced by the energy transferred from the BH to the disc, and the lower frequency is fitted by another rotating hotspot somewhere away from the inner edge of the disc, which arises from the screw instability of the magnetic field on the disc. It turns out that the 3:2 frequency ratio of HFQPOs in these X-ray binaries could be well fitted to the observational data with a much narrower range of the BH spin. In addition, the spectral properties of HFQPOs are discussed. The correlation of HFQPOs with jets from microquasars is contained naturally in our model.

A toy model for gamma-ray burst-supernovae (GRB-SN) is discussed by considering the coexistence of baryon poor outflows from black holes (BHs) and a powerful spin-connection to surrounding disk, giving rise to consistent calorimetry as described by van Putten (2001a) in a variant of the Blandford-Znajek process (BZ, 1977). In this model the half-opening angle of the magnetic flux tube on the horizon is determined by the mapping relation between the angular coordinate on the BH horizon and the radial coordinate on the surrounding accretion disk. The GRB is powered by the baryon poor outflows in the BZ process, and the associated SN is powered by very small fraction of the spin energy transferred from the BH to the disk in the magnetic coupling (MC) process. The timescale of the GRB is fitted by the duration of the open magnetic flux on the horizon. It turns out that the data of several GRB-SNe are well fitted with our model.

A toy model for magnetic extraction of energy from black hole (BH) accretion disk is discussed by considering the restriction of the screw instability to the magnetic field configuration. Three mechanisms of extracting energy magnetically are involved. (1) The Blandford-Znajek (BZ) process is related to the open magnetic field lines connecting the BH with the astrophysical load; (2) the magnetic coupling (MC) process is related to the closed magnetic field lines connecting the BH with its surrounding disk; and (3) a new scenario (henceforth the DL process) for extracting rotational energy from the disk is related to the open field lines connecting the disk with the astrophysical load. The expressions for the electromagnetic powers and torques are derived by using the equivalent circuits corresponding to the above energy mechanisms. It turns out that the DL power is comparable with the BZ and MC powers as the BH spin approaches unity. The radiation from a quasi-steady thin disk is discussed in detail by applying the conservation laws of mass, energy and angular momentum to the regions corresponding to the MC and DL processes. In addition, the poloidal currents and the current densities in BH magnetosphere are calculated by using the equivalent circuits.

Screw instability of the magnetic field connecting a rotating black hole (BH) with its surrounding disk is discussed based on the model of the coexistence of the Blandford-Znajek (BZ) process and the magnetic coupling (MC) process (CEBZMC). A criterion for the screw instability with the state of CEBZMC is derived based on the calculations of the poloidal and toroidal components of the magnetic field on the disk. It is shown by the criterion that the screw instability will occur, if the BH spin and the power-law index for the variation of the magnetic field on the disk are greater than some critical values. It turns out that the instability occurs outside some critical radii on the disk. It is argued that the state of CEBZMC always accompanies the screw instability. In addtition, we show that the screw instability contributes only a small fraction of magnetic extraction of energy from a rotating BH.

Effects of magnetic coupling (MC) of a rotating black hole (BH) with its surrounding accretion disk are discussed in detail in the following aspects: (i) The mapping relation between the angular coordinate on the BH horizon and the radial coordinate on the disk is modified based on a more reasonable configuration of magnetic field, and a condition for coexistence of the Blandford-Znajek (BZ) and the MC process is derived. (ii) The transfer direction of energy and angular momentum in MC process is described equivalently by the co-rotation radius and by the flow of electromagnetic angular momentum and redshifted energy, where the latter is based on an assumption that the theory of BH magnetosphere is applicable to both the BZ and MC processes. (iii) The profile of the current on the BH horizon and that of the current density flowing from the magnetosphere onto the horizon are given in terms of the angular coordinate of the horizon. It is shown that the current on the BH horizon varies with the latitude of the horizon and is not continuous at the angular boundary between the open and closed magnetic field lines. (iv) The MC effects on disk radiation are discussed, and a very steep emissivity is produced by MC process, which is consistent with the recent XMM-Newton observation of the nearby bright Seyfert 1 galaxy MCG-6-30-15 by a variety of parameters of the BH-disk system.

An analytic model of rotating hot spot is proposed to explain kilohertz quasi-periodic oscillations (kHz QPO) in X-ray binaries, which is based on the magnetic coupling (MC) of a rotating black hole (BH) with its surrounding accretion disc. The hot spot in the inner region of the disc is produced by energy transferred from a spinning BH with non-axisymmetric magnetic field. The frequency and energy band of the hot spot turns out to be related to six parameters, by which the strength and the azimuthal angular region of the bulging magnetic field on the BH horizon, the mass and spin of the BH, and the power-law index of the magnetic field varying with the radial coordinate on the disc are described. In addition, the correlation of the fluctuation of the bulging magnetic field with the widths of QPO frequency is discussed.

The transfer of energy and angular momentum in the magnetic coupling (MC) of a rotating black hole (BH) with its surrounding accretion disc is discussed based on a mapping relation derived by considering the conservation of magnetic flux with two basic assumptions: (i)the magnetic field on the horizon is constant, (ii) the magnetic field on the disc surface varies as a power law with the radial coordinate of the disc. The following results are obtained: (i) the transfer direction of energy and angular momentum between the BH and the disc depends on the position of a co-rotation radius relative to the MC region on the disc, which is eventually determined by the BH spin; (ii) the evolution characteristics of a rotating BH in the MC process without disc accretion are depicted in a parameter space, and a series of values of the BH spin are given to indicate the evolution characteristics; (iii) the efficiency of converting accreted mass into radiation energy of a BH-disc system is discussed by considering coexistence of disc accretion and the MC process; (iv) the MC effects on disc radiation and emissivity index are discussed and it is concluded that they are consistent with the recent XMM-Newton observation of the nearby bright Seyfert 1 galaxy MCG-6-30-15 with reference to a variety of parameters of the BH-disc system.

An analytic model of the evolution of a rotating black hole (BH) is proposed by considering the coexistence of disk accretion with the Blandford-Znajek process. The evolutionary characteristics of the BH are described in terms of three parameters: the BH spin $a_*$, the ratio $k$ of the angular velocity of the magnetic field lines to that of the BH horizon and the parameter $\lambda$ indicating the position of the inner edge of the disk. It is shown that the ratio $k$ being a little greater than 0.5 affects the evolutionary characteristics of the BH significantly, and the BH spin increases rather than decreases in its evolutionary process provided that the initial value of the BH spin is located in an appropriate value range determined by ratio $k$. Our calculations show that the system of a BH accretion disk with $k=0.6$ might provide a much higher output energy in a shorter timescale for gamma-ray bursts than the same system with $k=0.5$.