A hot plasma is the dominant phase of the interstellar medium of early-type galaxies. Its origin can reside in stellar mass losses, residual gas from the formation epoch, and accretion from outside of the galaxies. Its evolution is linked to the dynamical structure of the host galaxy, to the supernova and AGN feedback, and to (late-epoch) star formation, in a way that has yet to be fully understood. Important clues about the origin and evolution of the hot gas come from the abundances of heavy metals, that have been studied with increasing detail with XMM-Newton and Chandra. We present recent high resolution hydrodynamical simulations of the hot gas evolution that include the above processes, and where several chemical species, originating in AGB stars and supernovae of type Ia and II, have also been considered. The high resolution, of few parsecs in the central galactic region, allows us to track the metal enrichment, transportation and dilution throughout the galaxy. The comparison of model results with observed abundances reveals a good agreement for the region enriched by the AGN wind, but also discrepancies for the diffuse hot gas; the latter indicate the need for a revision of standard assumptions, and/or the importance of neglected effects as those due to the dust, and/or residual uncertainties in deriving abundances from the X-ray spectra.

What is the expected infrared output of elliptical galaxies? Here we report the latest findings obtained in this high time resolution (~10 years) and high spatial resolution (2.5 parsec at center) study. We add a set of grain physics to the MACER code, including (a) dust grains made in passive stellar evolution; (b) dust grain growth due to collision and sticking; (c) grain destruction due to thermal sputtering; (d) dust cooling of hot gas via inelastic collisions; and (e) radiation pressure on dust grains. The code improvements enable us to analyze the effects of dust on metal depletion and AGN obscuration, and also to assess the infrared output of the modeled galaxies. We simulate a representative massive elliptical galaxy of a central stellar velocity dispersion ~ 260 km/s and modest rotation. We find that: (1) the circumnuclear disk (of a size <~ 1 kpc) is dusty in its outer region where most of the metals are in dust grains, while in the inner disk most of the dust grains are destroyed by the AGN irradiation; (2) the dusty disk is optically thick to both the starlight within the disk and the radiation from the central AGN. Thus the AGN will be obscured behind the disk, and the latter is of a covering factor ~ 0.2; (3) the dust infrared emission is mainly due to the AGN irradiation. The median infrared luminosity is ~ 2e44 erg/s, and it can reach >~ 1e46 erg/s during outbursts; (4) the duty cycles of the AGN activities, star formation, and the dust infrared luminosity roughly match observations, e.g., in most of its lifetime, the simulated galaxy is a stereotypical "quiescent" elliptical galaxy with L_{IR} ~ 1e11*L_{solar}, while the star formation rate can exceed 250 M_{solar}/yr during central outbursts.

NGC 3115 is known as the low-luminosity active galactic nucleus which hosts the nearest ($z\sim0.002$) billion solar mass supermassive black hole ($\sim1.5\times10^9~M_\odot$). Its Bondi radius $r_\mathrm{B}$ ($\sim3\farcs6$) can be readily resolved with Chandra, which offers us an excellent opportunity to investigate the accretion flow onto a supermassive black hole. In this paper, we perform two-dimensional hydrodynamical numerical simulations, tailored for NGC 3115, on the mass flow across the Bondi radius. Our best fittings for the density and temperature agree well with the observations of the hot interstellar medium in the centre of NGC 3115. We find that the flow properties are solely determined by the local galaxy properties in the galaxy centre: (1) stellar winds (including supernova ejecta) supply the mass and energy sources for the accreting gas; (2) similar to the one-dimensional calculations, a stagnation radius $r_\mathrm{st}\sim0.1~r_\mathrm{B}$ is also found in the two-dimensional simulations, which divides the mass flow into an inflow-outflow structure; (3) the radiatively inefficient accretion flow theory applies well inside the stagnation radius, where the gravity is dominated by the supermassive black hole and the gas is supported by rotation; (4) beyond the stagnation radius, the stellar gravity dominates the spherical-like fluid dynamics and causes the transition from a steep density profile outside to a flat density profile inside the Bondi radius.

We add a suite of chemical abundances to the MACER (Massive AGN Controlled Ellipticals Resolved) 2D code, by solving 12 additional continuity equations for H, He, C, N, O, Ne, Mg, Si, S, Ca, Fe and Ni respectively with sources from AGB stars and supernovae of type Ia and II with metal yields based on standard stellar physics. New stars, formed in Toomre unstable circumnuclear disks (of a size $\lesssim150$ parsec), are assumed to have a top-heavy initial mass function with a power index of 1.65. The metal dilution effects due to cosmic accretion are also included. With the high resolution of few parsecs in central regions, resolved black hole accretion and AGN feedback, we can track the metal enrichment, transportation and dilution throughout the modeled massive elliptical galaxy of velocity dispersion $\sim280$ km/s. We retrieve the chemical composition of the BAL winds launched by the central AGN, synthesize the X-ray features of the hot ISM, and find that (1) the simulated metallicity in the BAL winds could be up to $\sim 8 Z_\odot$, while that of the hot ISM in the host galaxy is $\sim 2.3 Z_\odot$, matching well with SDSS observations of BLR gas; (2) the X-ray emitting hot gas is metal enriched with a typical value $\sim2.5 Z_\odot$; (3) the circumunuclear cold gas disk, where the metals are condensed, further enriched and recycled, plays a critical role in the metal enrichment; (4) the black hole accretion rate $\dot{M}_{\rm BH}$ linearly correlates with the star formation rate $\dot{M}_\star^+$ in the circumnuclear disk, i.e, $\dot{M}_\star^+ \sim 7.7 \dot{M}_{\rm BH}$, but lagged in time by roughly $10^6$ years.

Active galactic nuclei (AGN) with relativistic jets, powered by gas accretion onto their central supermassive black hole (SMBH), are unique laboratories for studying the physics of matter and elementary particles in extreme conditions that cannot be realized on Earth. For a long time since the discovery of AGN, photons were the only way to probe the underlying physical processes. The recent discovery of a very high energy neutrino, IceCube-170922A, coincident with a flaring blazar, TXS 0506+056, provides the first evidence that AGN jets are multi-messenger sources; they are capable of accelerating hadrons to very high energies, while producing non-thermal EM radiation and high-energy neutrinos. This new era of multi-messenger astronomy, which will mature in the next decade, offers us the unprecedented opportunity to combine more than one messenger to solve some long-standing puzzles of AGN jet physics: How do jets dissipate their energy to accelerate particles? What is the jet total kinetic power? Where and how do jets produce high-energy emission and neutrinos? What physical mechanisms drive the particle acceleration?

Based on our previous modeling of AGN feedback in isolated elliptical galaxies (Gan et al. 2014) using the MACER (Massive AGN Controlled Ellipticals Resolved) code, we extend and improve the model to include rotation, to limit star formation to regions of high density and low temperature, to facilitate angular momentum transfer via the Toomre instability in gaseous disks, and to improve the treatment of hot mode (low accretion rate) AGN feedback. The model galaxy now has an extended dark matter profile that matches with standard observations, but has a resolution of parsecs in the inner region and resolves the Bondi radius. We find that the results agree reasonably well with a panoply of observations: (1) both AGN activity and star formation are primarily in central cold gaseous disks, are bursty and mainly driven by the Toomre instability; (2)AGN duty cycle agrees well with the Soltan argument, i.e., the AGN spends most of its lifetime when it is in low luminosity (half of time with $L/L_{Edd}<7\times10^{-5}$), while emitting most of its energy when it is in high luminosity (half of radiant energy emitted with $L/L_{Edd}>0.06$); (3) the total star formation is $\sim$ few percents of the initial stellar mass, occurring in the bursts that would be associated with the observed E+A phenomenon. Most of the star formation occurs in the circumnuclear disk of a size <1 kpc, which is in agreement with recent observations; (4) the ISM X-ray luminosity varies within a reasonable range (median $L_{\rm X, ISM}=9.1\times10^{39}$ erg/s) in agreement with observations.

Understanding how galaxies maintain the inefficiency of star formation with physically self-consistent models is a central problem for galaxy evolution. Although numerous theoretical models have been proposed in recent decades, the debate still exists. By means of high-resolution two-dimensional hydrodynamical simulations, we study the three feedback effects (the stellar wind heating, SNe feedback, and AGN feedback) in suppressing star formation activities on the evolution of early-type galaxies with different stellar masses. AGN feedback models are updated based on \citet{Yuan2018}. The gas sources comes exclusively from the mass losses of dying low-mass stars for most of our models. We find that SNe feedback can keep star formation at a significantly low level for low mass elliptical galaxies for a cosmological evolution time. For the high mass galaxies, AGN feedback can efficiently offset the radiative cooling and thus regulate the star formation activities. Such a suppression of star formation is extremely efficient in the inner region of the galaxies. AGB heating cannot account for this suppression for low and high mass galaxies. The X-ray temperature $T_{\rm X}$ and luminosity $L_{\rm X}$ of hot plasma can be in agreement with the observed data with the inclusion of effective feedback processes. These results thus suggest that we can use $T_{\rm X}$ and $L_{\rm X}$ to probe the role of different feedback processes. The inclusion of additional gas sources can make the mass scale between SNe and AGN feedback dominating in suppressing star formation decrease to an observationally inferred value of a few $10^{10}~M_{\odot}$.

This is the second paper of our series of works of studying the effects of active galactic nuclei (AGN) feedback on the cosmological evolution of an isolated elliptical galaxy by performing two-dimensional high-resolution hydrodynamical numerical simulations. In these simulations, the inner boundary is chosen so that the Bondi radius is resolved. Physical processes like star formation, SNe Ia and II are taken into account. Compared to previous works, the main improvements is that we adopt the most updated AGN physics, which is described in detail in the first paper of this series (Yuan et al. 2018, Paper I). These improvements include the discrimination of the two accretion modes of the central AGN and the most updated descriptions of the wind and radiation in the two modes. In Paper I, we consider the case that the specific angular momentum of the gas in the galaxy is very low. In this paper, we consider the case that the specific angular momentum of the gas is high. In the galactic scale, we adopt the gravitational torques raised due to non-axisymmetric structure in the galaxy as the mechanism of the transfer of angular momentum of gas, as proposed in some recent works. Since our simulations are axisymmetric, we make use of a parameterized prescription to mimic this mechanism. Same as Paper I, we investigate the AGN light curve, typical AGN lifetime, growth of the black hole mass, AGN duty-cycle, star formation, and the X-ray surface brightness of the galaxy. Special attention is paid to the effects of specific angular momentum of the galaxy on these properties. We find that some results are qualitatively similar to those shown in Paper I, while some results such as star formation and black hole growth do show a significant difference due to the mass concentration in the galactic disk as a consequence of galactic rotation.

This paper summarizes our recent works of studying AGN feedback in an isolated elliptical galaxy by performing high-resolution hydrodynamical numerical simulations. Bondi radius is resolved and the mass accretion rate of the black hole is calculated. The most updated AGN physics, namely the discrimination of cold and hot accretion modes and the exact descriptions of the AGN radiation and wind for a given accretion rate are adopted and their interaction with the gas in the host galaxy is calculated. Physical processes such as star formation and SNe feedback are taken into account. Consistent with observation, we find the AGN spends most of the time in the low-luminosity regime. AGN feedback overall suppresses the star formation; but depending on location in the galaxy and time, it can also enhance it. The light curve of specific star formation rate is not synchronous with the AGN light curve. We find that wind usually plays a dominant role in controlling the AGN luminosity and star formation, but radiation also cannot be neglected.

We investigate the effects of AGN feedback on the cosmological evolution of an isolated elliptical galaxy by performing two-dimensional high-resolution hydrodynamical numerical simulations. The inner boundary of the simulation is chosen so that the Bondi radius is resolved. Compared to previous works, the two accretion modes, namely hot and cold, which correspond to different accretion rates and have different radiation and wind outputs, are carefully discriminated and the feedback effects by radiation and wind in each mode are taken into account. The most updated AGN physics, including the descriptions of radiation and wind from the hot accretion flows and wind from cold accretion disks, are adopted. Physical processes like star formation, Type Ia and Type II supernovae are taken into account. We study the AGN light curve, typical AGN lifetime, growth of the black hole mass, AGN duty-cycle, star formation, and the X-ray surface brightness of the galaxy. We compare our simulation results with observations and find general consistency. Comparisons with previous simulation works find significant differences, indicating the importance of AGN physics. The respective roles of radiation and wind feedbacks are examined and it is found that they are different for different problems of interest such as AGN luminosity and star formation. We find that it is hard to neglect any of them, so we suggest to use the names of "cold feedback mode" and "hot feedback mode" to replace the currently used ones.

We perform two-dimensional simulations to study how the wind strength changes with accretion rate. We take into account bremsstrahlung, synchrotron radiation and the Comptonization. We find that when the accretion rate is low, radiative cooling is not important, the accretion flow is hot accretion flow. For the hot accretion flow, wind is very strong. The mass flux of wind can be $\sim 50\%$ of the mass inflow rate. When the accretion rate increases to a value at which radiative cooling rate is roughly equal to or slightly larger than viscous heating rate, cold clumps can form around the equatorial plane. In this case, the gas pressure gradient force is small and wind is very weak. Our results may be useful for the sub-grid model of active galactic nuclear feedback study.

The distinctive morphology of head-tail radio galaxies reveals strong interactions between the radio jets and their intra-cluster environment, the general consensus on the morphology origin of head-tail sources is that radio jets are bent by violent intra-cluster weather. We demonstrate in this paper that such strong interactions provide a great opportunity to study the jet properties and also the dynamics of intra-cluster medium (ICM). By three-dimensional magnetohydrodynamical simulations, we analyse the detailed bending process of a magnetically dominated jet, based on the magnetic tower jet model. We use stratified atmospheres modulated by wind/shock to mimic the violent intra-cluster weather. Core sloshing is found to be inevitable during the wind-cluster core interaction, which induces significant shear motion and could finally drive ICM turbulence around the jet, making it difficult for jet to survive. We perform detailed comparison between the behaviour of pure hydrodynamical jets and magnetic tower jet, and find that the jet-lobe morphology could not survive against the violent disruption in all of our pure hydrodynamical jet models. On the other hand, the head-tail morphology is well reproduced by using a magnetic tower jet model bent by wind, in which hydrodynamical instabilities are naturally suppressed and the jet could always keep its integrity under the protection of its internal magnetic fields. Finally, we also check the possibility for jet bending by shock only. We find that shock could not bend jet significantly, so could not be expected to explain the observed long tails in head-tail radio galaxies.

We investigate state transitions in black hole X-ray binaries through different parameters by using two-dimensional axisymmetric hydrodynamical simulation method. For radiative cooling in hot accretion flow, we take into account the bremsstrahlung, synchrotron and synchrotron-self Comptonization self-consistently in the dynamics. Our main result is that the state transitions occur when the accretion rate reaches a critical value $\dot M \sim 3\alpha\ \dot M_{\rm Edd}$, above which cold and dense clumpy/filamentary structures are formed, embedded within the hot gas. We argued this mode likely corresponds to the proposed two-phase accretion model, which may be responsible for the intermediate state of black hole X-ray binaries. When the accretion rate becomes sufficiently high, the clumpy/filamentary structures gradually merge and settle down onto the mid-plane. Eventually the accretion geometry transforms to a disc-corona configuration. In summary our results are consistent with the truncated accretion scenario for the state transition.

Numerical simulations of black hole hot accretion flows have shown the existence of strong wind. Those works focus only on the region close to black hole thus it is unknown whether or where the wind production stops at large radii. To address this question, Bu et al. (2016) have performed hydrodynamic (HD) simulations by taking into account the gravitational potential of both the black hole and the nuclear star clusters. The latter is assumed to be $\propto \sigma^2 \ln(r)$, with $\sigma$ being the velocity dispersion of stars and $r$ be the distance from the center of the galaxy. It was found that when the gravity is dominated by nuclear stars, i.e., outside of radius $R_A\equiv GM_{\rm BH}/\sigma^2$, winds can no longer be produced. That work, however, neglects the magnetic field, which is believed to play a crucial dynamical role in the accretion and thus must be taken into account. In this paper, we revisit this problem by performing magneto-hydrodynamical (MHD) simulations. We confirm the result of Bu et al. (2016), namely wind can't be produced at the region of $R>R_A$. Our result, combined with the results of Yuan et al. (2015), indicates that the formula describing the mass flux of wind $\dot{M}_{\rm wind}=\dot{M}_{\rm BH}(r/20r_s)$ can only be applied to the region where the black hole potential is dominant. Here $\dot{M}_{\rm BH}$ is the mass accretion rate at the black hole horizon and the value of $R_A$ is similar to the Bondi radius.

Previous works show strong winds exist in hot accretion flows around black holes. Those works focus only on the region close to the black hole thus it is unknown whether or where the wind production stops at large radii. In this paper, we investigate this problem by hydrodynamical simulations. We take into account gravities of both the black hole and the nuclear star clusters. For the latter, we assume that the velocity dispersion of stars is a constant and its gravitational potential $\propto \sigma^2 \ln (r)$, where $\sigma$ is the velocity dispersion of stars and $r$ is the distance from the center of the galaxy. We focus on the region where the gravitational potential is dominated by the star cluster. We find, same as the accretion flow at small radii, the mass inflow rate decreases inward and the flow is convectively unstable. However, trajectory analysis shows that there is very few wind launched from the flow. Our result, combined with the results of Yuan et al. (2015), indicates that the mass flux of wind launched from hot accretion flow $\dot{M}_{\rm wind}=\dot{M}_{\rm BH}(r/20r_s)$, with $r\la R_A\equiv GM_{\rm BH}/\sigma^2$. Here $\dot{M}_{\rm BH}$ is accretion rate at black hole horizon. $R_A$ is similar to Bondi radius. We argue that the inward decrease of inflow rate is not because of mass loss via wind, but because of convective motion. The disappearance of wind outside of $R_A$ must be because of the change of the gravitational potential, but the exact reason remains to be probed.

In a previous work, we have shown that the formation of the Fermi bubbles can be due to the interaction between winds launched from the hot accretion flow in Sgr A* and the interstellar medium (ISM). In that work, we focus only on the morphology. In this paper we continue our study by calculating the gamma-ray radiation. Some cosmic ray protons (CRp) and electrons must be contained in the winds, which are likely formed by physical processes such as magnetic reconnection. We have performed MHD simulations to study the spatial distribution of CRp, considering the advection and diffusion of CRp in the presence of magnetic field. We find that a permeated zone is formed just outside of the contact discontinuity between winds and ISM, where the collisions between CRp and thermal nuclei mainly occur. The decay of neutral pions generated in the collisions, combined with the inverse Compton scattering of background soft photons by the secondary leptons generated in the collisions and primary CR electrons can well explain the observed gamma-ray spectral energy distribution. Other features such as the uniform surface brightness along the latitude and the boundary width of the bubbles are also explained. The advantage of this accretion wind model is that the adopted wind properties come from the detailed small scale MHD numerical simulation of accretion flows and the value of mass accretion rate has independent observational evidences. The success of the model suggests that we may seriously consider the possibility that cavities and bubbles observed in other contexts such as galaxy clusters may be formed by winds rather than jets.

Previous MHD simulations have shown that wind must exist in black hole hot accretion flows. In this paper, we continue our study by investigating the detailed properties of wind, such as mass flux and poloidal speed, and the mechanism of wind production. For this aim, we make use of a three dimensional GRMHD simulation of hot accretion flows around a Schwarzschild black hole. The simulation is designed so that the magnetic flux is not accumulated significantly around the black hole. To distinguish real wind from turbulent outflows, we track the trajectories of the virtual Largrangian particles from simulation data. We find two types of real outflows, i.e., a quasi-relativistic jet close to the axis and a sub-relativistic wind subtending a much larger solid angle. Most of the wind originates from the surface layer of the accretion flow. The poloidal wind speed almost remains constant once they are produced, but the flux-weighted wind speed roughly follows $v_{\rm p, wind}(r)\approx 0.25 v_k(r)$. The mass flux of jet is much lower but the speed is much higher, $v_{\rm p,jet}\sim (0.3-0.4) c$. Consequently, both the energy and momentum fluxes of the wind are much larger than those of the jet. We find that the wind is produced and accelerated primarily by the combination of centrifugal force and magnetic pressure gradient, while the jet is mainly accelerated by magnetic pressure gradient. Finally, we find that the wind production efficiency $\epsilon_{\rm wind}\equiv\dot{E}_{\rm wind}/\dot{M}_{\rm BH}c^2\sim 1/1000$, in good agreement with the value required from large-scale galaxy simulations with AGN feedback.

Based on two-dimensional high resolution hydrodynamic numerical simulation, we study the mechanical and radiative feedback effects from the central AGN on the cosmological evolution of an isolated elliptical galaxy. Physical processes such as star formation and supernovae are considered. The inner boundary of the simulation domain is carefully chosen so that the fiducial Bondi radius is resolved and the accretion rate of the black hole is determined self-consistently. In analogy to previous works, we assume that the specific angular momentum of the galaxy is low. It is well-known that when the accretion rates are high and low, the central AGNs will be in cold and hot accretion modes, which correspond to the radiative and kinetic feedback modes, respectively. The emitted spectrum from the hot accretion flows is harder than that from the cold accretion flows, which could result in a higher Compton temperature accompanied by a more efficient radiative heating, according to previous theoretical works. Such a difference of the Compton temperature between the two feedback modes, the focus of this study, has been neglected in previous works. Significant differences in the kinetic feedback mode are found as a result of the stronger Compton heating and accretion becomes more chaotic. More importantly, if we constrain models to correctly predict black hole growth and AGN duty cycle after cosmological evolution, we find that the favored model parameters are constrained: mechanical feedback efficiency diminishes with decreasing luminosity (the maximum efficiency being $\simeq 10^{-3.5}$) and X-ray Compton temperature increases with decreasing luminosity, although models with fixed mechanical efficiency and Compton temperature can be found that are satisfactory as well. We conclude that radiative feedback in the kinetic mode is much more important than previously thought.

We perform two dimensional hydrodynamic numerical simulations to study the positive active galactic nucleus feedback which triggers, rather than suppresses, star formation. Recently, it was shown by Nayakshin et al. and Ishibashi et al. that star formation occurs when the cold interstellar medium is squeezed by the impact of mass outflow or radiation pressure, respectively. Mass outflow is ubiquitous in this astrophysical context, and radiation pressure is also important if the AGN is luminous. For the first time in this subject, we incorporate both mass outflow feedback and radiative feedback into our model. Consequently, the ISM is shocked into shells by the AGN feedback, and these shells soon fragment into clumps and filaments because of Rayleigh-Taylor and thermal instabilities. We have two major findings: (1) the star formation rate can indeed be very large in the clumps and filaments. However, the resultant star formation rate density is too large compared with previous works, which is mainly because we ignore the fact that most of the stars that are formed would be disrupted when they move away from the galactic center. (2) Although radiation pressure feedback has a limited effect, when mass outflow feedback is also included, they reinforce each other. Specifically, in the gas-poor case, mass outflow is always the dominant contributor to feedback.

We perform time-dependent, 2DHD numerical simulations to study the dynamics of a slowly rotating accretion flow from sub-pc to pc scales under the irradiation from the central AGN. Compared to previous work, we improve the calculation of the radiative force due to X-rays. More importantly, in addition to radiative pressure and radiative heating/cooling directly from the central AGN, in the momentum equation we also include the force due to the scattered and reprocessed photons. We find that the accretion flow properties change significantly due to this "re-radiation" effect. The inflow rate at the inner boundary is reduced, while the outflow rate at the outer boundary is enhanced by about one order of magnitude. This effect is more significant when the density at the outer boundary is higher. The properties of outflows such as velocity, momentum and energy fluxes, and the ratio of outflow rate and the accretion rate, are calculated. We find that the efficiency of transferring the radiation power into the kinetic power of outflow is typically $10^{-3}$, far below the value of $\sim 0.05$ which is assumed in some cosmological simulations. The effect of the temperature of the gas at the outer boundary ($T_0$) is investigated. When $T_0$ is high, the emitted luminosity of the accretion flow oscillates. This is because in this case the gas around the Bondi radius can be more easily heated to be above the virial temperature due to its high internal energy. Another question we hope to address is the so-called "sub-Eddington" puzzle. Observationally, the luminosity of almost all AGNs are sub-Eddington, while theoretically the luminosity of an accretion flow can easily be super-Eddington. We find that even when the re-radiation effect is included and outflow does become much stronger, the luminosity, while reduced, can still be super-Eddington.

We apply epicyclic resonances to the magnetic connection (MC) of a black hole (BH) with a relativistic accretion disc, interpreting the high frequency quasi-periodic oscillations (HFQPOs) with 3:2 pairs observed in three BH X-ray binaries. It turns out that the 3:2 HFQPO pairs are associated with the steep power-law states, and the severe damping can be overcome by transferring energy and angular momentum from a spinning BH to the inner disc in the MC process.

We discuss transportation and redistribution of energy and angular momentum in the magnetic connection (MC) process and Blandford-Payne (BP) process. The MC results in readjusting the interior viscous torque, and its effects are operative not only in but also beyond the MC region. The BP process is invoked to transfer the "excessive" angular momentum from an accretion disc. In addition, we derive a criterion for the interior viscous torque to resolve the puzzle of the overall equilibrium of angular momentum in disc accretion. It turns out that the BP efficiency of extracting angular momentum and the intensity of the outflow are required to be greater than some critical values.

Recent observations and theoretical work on gamma-ray bursts (GRBs) favor the central engine model of a Kerr black hole (BH) surrounded by a magnetized neutrino-dominated accretion flow (NDAF). The magnetic coupling between the BH and disk through a large-scale closed magnetic field exerts a torque on the disk, and transports the rotational energy from the BH to the disk. We investigate the properties of the NDAF with this magnetic torque. For a rapid spinning BH, the magnetic torque transfers enormous rotational energy from BH into the inner disk. There are two consequences: (i) the luminosity of neutrino annihilation is greatly augmented; (ii) the disk becomes thermally and viscously unstable in the inner region, and behaves S-Shape of the surface density versus accretion rate. It turns out that magnetically torqued NDAF can be invoked to interpret the variability of gamma-ray luminosity. In addition, we discuss the possibility of restarting the central engine to produce the X-ray flares with required energy.

We propose a simplified model of outflow/jet driven by the Blandford-Payne (BP) process from advection dominated accretion flows (ADAF) and derive the expressions of the BP power and disk luminosity based on the conservation laws of mass, angular momentum and energy. We fit the 2--10 keV luminosity and kinetic power of 15 active galactic nucleus (AGNs) of sub-Eddington luminosity. It is found that there exists an anti-correlation between the accretion rate and the advection parameter, which could be used to explain the correlation between Eddington-scaled kinetic power and bolometric luminosity of the 15 samples. In addition, the Ledlow-Owen relation for FR I/II dichotomy is re-expressed in a parameter space consisting of logarithm of dimensionless accretion rate versus that of the BH mass. It turns out that the FR I/II dichotomy is determined mainly by the dimensionless accretion rate, being insensitive to the BH mass. And the dividing accretion rate is less than the critical accretion rate for ADAFs, suggesting that FR I sources are all in the ADAF state.

We discuss one of the possible origins of large-scale magnetic fields based on a continuous distribution of toroidal electric current flowing in the inner region of the disc around a Kerr black hole (BH) in the framework of general relativity. It turns out that four types of configuration of the magnetic connection (MC) are generated, i.e., MC of the BH with the remote astrophysical load (MCHL), MC of the BH with the disc (MCHD), MC of the plunging region with the disc (MCPD) and MC of the inner and outer disc regions (MCDD). It turns out that the Blandford-Znajek (BZ) process can be regarded as one type of MC, i.e., MCHL. In addition, we propose a scenario for fitting the quasi-periodic oscillations in BH binaries based on MCDD associated with the magnetic reconnection.

We propose a model of magnetic connection (MC) of a black hole with its surrounding accretion disc based on large-scale magnetic field. The MC gives rise to transport of energy and angular momentum between the black hole and the disc, and the closed field lines pipe the hot matter evaporated from the disc, and shape it in the corona above the disc to form a magnetically induced disc-corona system, in which the corona has the same configuration as the large-scale magnetic field. We numerically solve the dynamic equations in the context of the Kerr metric, in which the large-scale magnetic field is determined by dynamo process and equipartition between magnetic pressure and gas pressure. Thus we can obtain a global solution rather than assuming the distribution of large-scale magnetic field beforehand. The main MC effects lie in three aspects. (1) The rotational energy of a fast-spinning black hole can be extracted, enhancing the dissipation in the accretion disc, (2) the closed field lines provide a natural channel for corona matter escaping from disc and finally falling into black hole, and (3) the scope of the corona can be bounded by the conservation of magnetic flux. We simulate the high-energy spectra of this system by using Monte-Carlo method, and find that the relative hardness of the spectra decreases as accretion rate or black hole spin $a_*$ increases. We fit the typical X-ray spectra of three black-hole binaries (GRO J1655-40, XTE 1118+480 and GX 339-4) in the low/hard or very high state.

Association of the high-frequency quasi-periodic oscillation (HFQPO) pairs with the broad Fe K line in XTE J1550-564 and GRO J1655-40 is discussed based on the magnetic coupling (MC) of a rotating black hole (BH) with its surrounding disc. The 3:2 HFQPO pairs are interpreted by virtue of the inner and outer hotspots arising from non-axisymmetric magnetic field, where the inner hotspot is produced by a torque exerted at the inner edge of the disc, and the outer hotspot is created by the screw instability of the large-scale magnetic field. The very steep emissivity index is created predominantly by the torque exerted at the inner edge of the disc. It turns out that the 3:2 HFQPO pairs observed in the two sources can be fitted by tuning several model parameters, such as the BH spin, and the main features of this model lie in three aspects. (1) The condition for only one HFQPO is discussed based on the two mechanisms for producing the 3:2 HFQPO pairs, (2) an explanation is given for a systematic shift away from disc dominated flux with the increasing power-law flux as the HFQPO pairs shift from the higher to lower frequencies, which is consistent with the analysis given by Remillard et al. (2002), and (3) the BH spin in XTE J1550-564 and GRO J1655-40 can be estimated by combining the 3:2 HFQPO pairs with the very steep emissivity index required for fitting the broad Fe K emission line.

Aims. A simplified model of jet power from active galactic nuclei is proposed in which the relationship between jet power and disk luminosity is discussed by combining disk accretion with two mechanisms of extracting energy magnetically from a black hole accretion disk, i.e., the Blandford-Payne (BP) and the Blandford-Znajek (BZ) processes. Methods. By including the BP process into the conservation laws of mass, angular momentum and energy, we derive the expressions of the BP power and disk luminosity, and the jet power is regarded as the sum of the BZ and BP powers. Results. We find that the disk radiation flux and luminosity decrease because a fraction of the accretion energy is channelled into the outflow/jet in the BP process. It is found that the dominant cooling mode of the accretion disk is determined mainly by how the poloidal magnetic field decreases with the cylindrical radius of the jet. By using the parameter space we found, which consists of the black hole spin and the self-similar index of the configuration of the poloidal magnetic field frozen in the disk, we were able to compare the relative importance of the following quantities related to the jet production: (1) the BP power versus the disk luminosity, (2) the BP power versus the BZ power, and (3) the jet power versus the disk luminosity. In addition, we fit the jet power and broad-line region luminosity of 11 flat-spectrum radio quasars (FSRQs) and 17 steep-spectrum radio quasars (SSRQs) based on our model.

The effects of magnetic coupling (MC) process on the inner edge of the disc are discussed in detail. It is shown that the inner edge can deviate from the innermost stable circular orbit (ISCO) due to the magnetic transfer of energy and angular momentum between a Kerr black hole (BH) and its surrounding accretion disc. It turns out that the inner edge could move inward and outward for the BH spin $a_{*}$ being greater and less than 0.3594, respectively. The MC effects on disc radiation are discussed based on the displaced inner edge. A very steep emissivity can be provided by the MC process, which is consistent with the observation of MCG-6-30-15. In addition, the BH spins of GRO J1655-40 and GRS 1915+105 are detected by X-ray continuum fitting based on this model.