Source author record

De-Fu Bu

De-Fu Bu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2022arXiv

Effect of radiation drag on the line-force-driven winds

Ultra-fast outflows (UFOs) with mildly relativistic velocities are measured using the X-ray spectra of radio-quiet and -loud active galactic nuclei (AGNs). In general, UFOs are believed to be generated from the accretion disk around a black hole (BH). A line-force driving model is suggested to be the mechanism to drive UFOs from the accretion disk. In this paper, we use the non-hydrodynamic approach to examine the influences of radiation-drag effects on the line-force-driven winds generated from the accretion disk. We find that the radiation-drag effects can significantly weaken the line-force-driven winds. Compared with the case without the radiation-drag effects, when the radiation-drag effects are considered, the maximum speed of winds is reduced by $\sim$60\%--70\%, the mass outflow rate is reduced by $\sim$50\%--80\%, and the kinetic power is reduced by about an order of magnitude. The radiation-drag effects narrow the area where the winds are generated.

preprint2022arXiv

Large-scale Dynamics of Winds Driven by Line Force from a Thin Accretion Disk

Winds play a significant role in active galactic nuclei feedback process. Previous simulations studying winds only focus on a small dynamical range. Therefore, it is unknown how far the winds can go and what the properties of the winds will be if they can move to large radii. We perform simulations to study the large scale dynamics of winds driven by line force. We find that the properties of the winds depend on both black hole mass ($M_{BH}$) and accretion disk luminosity. When the accretion disk luminosity is $0.6L_{edd}$ ($L_{edd}$ being Eddington luminosity), independent of $M_{BH}$, the winds have kinetic energy flux exceeding $1\% L_{edd}$ and can escape from the black hole potential. For the case with the accretion disk luminosity equaling 0.3$L_{edd}$, the strength of the winds decreases with the decrease of $M_{BH}$. If $M_{BH}$ decreases from $10^9$ to $10^6$ solar mass ($M_\odot$), the winds kinetic energy flux decreases from $\sim 0.01 L_{edd}$ to $ \sim 10^{-6} L_{edd}$. In case of $M_{BH}\geq 10^7 M_\odot$, winds can escape from black hole potential. In the case of $M_{BH}=10^6 M_\odot$, the winds can not escape. We find that for the ultra-fast winds observed in hard X-ray bands (\citealt{Gofford et al. 2015}), the observed dependence of the mass flux and the kinetic energy flux on accretion disk luminosity can be well produced by line force driven winds model. We also find that the properties of the ultra-fast winds observed in soft X-ray bands can be explained by the line force driven winds model.

preprint2022arXiv

Magnetically driven accretion disc winds: the role of gas thermodynamics and comparison to ultra-fast outflows

Winds are commonly observed in luminous active galactic nuclei (AGNs). A plausible model of those winds is magnetohydrodynamic (MHD) disc winds. In the case of disc winds from a thin accretion disc, isothermal or adiabatic assumption is usually adopted in such MHD models. In this work we perform two-dimensional MHD simulations implementing different thermal treatments (isothermal, adiabatic and radiative) to study their effects on winds from a thin accretion disc. We find that both the isothermal model and the adiabatic model overestimate the temperature, underestimate the power of disc winds, and cannot predict the local structure of the winds, compared to the results obtained by solving the energy equation with radiative cooling and heating. Based on the model with radiative cooling and heating, the ionization parameter, the column density and the velocity of the disc winds have been compared to the observed ultrafast outflows (UFOs). We find that in our simulations the UFOs can only be produced inside hundreds of Schwarzschild radius. At much larger radii, no UFOs are found. Thus, the pure MHD winds cannot interpret all the observed UFOs.

preprint2022arXiv

Self-similar Solution of Hot Accretion Flow with Thermal Conduction and Anisotropic Pressure

We explore the effects of anisotropic thermal conduction, anisotropic pressure, and magnetic field strength on the hot accretion flows around black holes by solving the axisymmetric, steady-state magnetohydrodynamic equations. The anisotropic pressure is known as a mechanism for transporting angular momentum in weakly collisional plasmas in hot accretion flows with extremely low mass accretion rates. However, anisotropic pressure does not extensively impact the transport of the angular momentum, it leads to shrinkage of the wind region. Our results show that the strength of the magnetic field can help the Poynting energy flux overcomes the kinetic energy flux. This result may be applicable to understand the hot accretion flow in the Galactic Center Sgr A* and M87 galaxy.

preprint2022arXiv

The origin of the optical/ultraviolet emission of optical/ultraviolet tidal disruption events

One of the most prominent problems of optical/ultraviolet (UV) tidal disruption events (TDEs) is the origin of their optical/UV emission. It has been proposed that the soft X-rays produced by the stellar debris accretion disk can be reprocessed into optical/UV photons by a surrounding optically thick envelope or outflow. However, there is still no detailed models for this mechanism. In this paper, by performing hydrodynamic simulations with radiative transfer, we calculate the optical/UV emission of the circularized stellar debris accretion flow/outflow system. We find that the optical/UV photons can be generated by reprocessing the emission of the accretion flow in the optically thick outflows. The model can well interpret the observed emission properties of optical/UV TDEs, including the emission radius, the radiation temperature and the luminosity, as well as the evolution of these quantities with time, providing a strong theoretical basis for understanding the origin of optical/UV TDEs.

preprint2021arXiv

Anisotropic multimessenger signals from black hole neutrino-dominated accretion flows with outflows in binary compact object mergers

A black hole (BH) hyperaccretion system might be born after the merger of a BH and a neutron star (NS) or a binary NS (BNS). In the case of a high mass accretion rate, the hyperaccretion disk is in a state of neutrino-dominated accretion flow (NDAF) and emits numerous anisotropic MeV neutrinos. Only a small fraction of these neutrinos annihilates in the space outside of the disk and then launch ultrarelativistic jets that break away from the merger ejecta to power gamma-ray bursts. Mergers and their remnants are generally considered sources of gravitational waves (GWs), neutrinos, and kilonovae. Anisotropic neutrino emission and anisotropic high-velocity material outflows from central BH-NDAF systems can also trigger strong GWs and luminous disk-outflow-driven (DOD) kilonovae, respectively. In this paper, the anisotropic multimessenger signals from NDAFs with outflows, including DOD kilonovae, MeV neutrinos, and GWs, are presented. As the results, the typical AB magnitude of the DOD kilonovae is lower than that of AT 2017gfo at the same distance, and it decreases with increasing viewing angles and its anisotropy is not sensitive to the outflow mass distribution but mainly determined by the velocity distribution. Since neutrinos with $\gtrsim 10~\rm MeV$ are mainly produced in the inner region of the disk, they will be dramatically deflected to a large viewing angle by relativity effects. Moreover, the strains of GWs induced by anisotropic neutrinos increase with increasing viewing angles. The accumulation of multimessenger detection of the BH-NS/BNS mergers with different viewing angles might further verify the existence of NDAFs with outflows.

preprint2020arXiv

Gas flows in an Active galactic nuclei. I. Two-phase gas inflow

We perform two-dimensional hydrodynamic simulations to study the non-rotating gas flow in the region from 1-1000 parsec. This region connects the galaxy scale and the accretion disk scale. The gas is irradiated by the photons generated by the accretion disk system near the central black hole with $10^8$ solar mass. We assume that the luminosity of the central accretion system is constant and equals to $0.3$ Eddington luminosity. Gas with density higher than $10^{-24} {\rm g \ cm^{-3}}$ is found to be thermally unstable. Two phase, cold and hot, gas flow will form due to thermal instability. We calculated the ratio of cold gas accretion rate to hot gas accretion rate. This ratio is a function of gas density and generally increases with the increase of gas density. The maximum value of this ratio is $16$, when gas density is $10^{-21} {\rm g \ cm^{-3}} $. Gas with density lower than $10^{-24} {\rm g \ cm^{-3}}$ is found to be thermally stable and the gas flow is one-phase. The applications of the results are briefly discussed.

preprint2020arXiv

Two-temperature radiative hot accretion flow around neutron stars

Numerical simulations of radiative two-temperature hot accretion flows (HAFs) around Neutron stars (NSs) are performed. We assume that all of the energy carried by the HAF around a NS will be thermalized and radiated out at the surface of the NS. The thermal photons will propagate outwards radially and cool the HAF vis Comptonization. We define $\dot m$ as the mass accretion rate at the surface of the central object in unit of Eddington accretion rate ($\dot M_{\rm Edd}=10L_{\rm Edd}/c^2$, with $L_{\rm Edd}$ and $c$ being Eddington luminosity and speed of light, respectively). When $\dot m$ is lower than $\sim 10^{-4}$, the cooling of the HAF is not important and outflows are very strong. When $\dot m > \sim 10^{-3}$, cooling becomes important and outflows are significantly weak. In the range $10^{-4} < \dot m < 10^{-3}$, the HAFs transients from a strong outflow phase to a very weak outflow phase with increase of $\dot m$. The properties of the HAF around a NS are also compared to those of the HAF around a BH. We find that with a similar $\dot m$, the dynamical properties of the HAF around a NS are quite similar as those of the HAF around a BH. However, the emitted spectrum of a HAF around a NS can be quite different from that of a HAF around a BH due to the presence of a thermal soft X-ray component coming from the surface of the NS.

preprint2019arXiv

Two-dimensional Inflow-outflow Solution of Supercritical Accretion Flow

We present the two-dimensional inflow-outflow solutions of radiation hydrodynamic (RHD) equations of supercritical accretion flows. Compared with prior studies, we include all components of the viscous stress tensor. We assume steady state flow and use self-similar solutions in the radial direction to solve the equations in $ r-θ$ domain of the spherical coordinates. The set of differential equations have been integrated from the rotation axis to the equatorial plane. We find that the self-similarity assumption requires that the radial profile of density is described by $ ρ(r) \propto r^{-0.5} $. Correspondingly, the radial profile of the mass inflow rate decreases with decreasing radii as $ \dot{M}_\mathrm{in} \propto r $. Inflow-outflow structure has been found in our solution. In the region $ θ> 65^{\circ} $ there exist inflow while above that flow moves outward and outflow could launch. The driving forces of the outflow are analyzed and found that the radiation force is dominant and push the gas particles outwards with poloidal velocity $ \sim 0.25 c $. The properties of outflow are also studied. The results show that the mass flux weighted angular momentum of the inflow is lower than that of outflow, thus the angular momentum of the flow can be transported by the outflow. We also analyze the convective stability of the supercritical disk and find that in the absence of the magnetic field, the flow is convectively unstable. Our analytical results are fully consistent with the previous numerical simulations of the supercritical accretion flow.

preprint2016arXiv

Effects of anisotropic thermal conduction on wind properties in hot accretion flow

Previous works have clearly shown the existence of winds from black hole hot accretion flow and investigated their detailed properties. In extremely low accretion rate systems, the collisional mean-free path of electrons is large compared with the length-scale of the system, thus thermal conduction is dynamically important. When the magnetic field is present, the thermal conduction is anisotropic and energy transport is along magnetic field lines. In this paper, we study the effects of anisotropic thermal conduction on the wind production in hot accretion flows by performing two-dimensional magnetohydrodynamic simulations. We find that thermal conduction has only moderate effects on the mass flux of wind. But the energy flux of wind can be increased by a factor of $\sim 10$ due to the increase of wind velocity when thermal conduction is included. The increase of wind velocity is because of the increase of driving forces (e.g. gas pressure gradient force and centrifugal force) when thermal conduction is included. This result demonstrates that thermal conduction plays an important role in determining the properties of wind.

preprint2016arXiv

Hydrodynamical numerical simulation of wind production from black hole hot accretion flows at very large radii

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 σ^2 \ln (r)$, where $σ$ 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}/σ^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.

preprint2016arXiv

Magneto-hydrodynamical Numerical simulation of wind production from black hole hot accretion flows at very large radii

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 σ^2 \ln(r)$, with $σ$ 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}/σ^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.

preprint2016arXiv

Two-dimensional inflow-wind solution of black hole accretion with an evenly symmetric magnetic field

We solve the two-dimensional magnetohydrodynamic (MHD) equations of black hole accretion with the presence of magnetic field. The field includes a turbulent component, whose role is represented by the viscosity, and a large-scale ordered component. The latter is further assumed to be evenly symmetric with the equatorial plane. The equations are solved in the $r-θ$ plane of a spherical coordinate by assuming time-steady and radially self-similar. An inflow-wind solution is found. Around the equatorial plane, the gas is inflowing; while above and below the equatorial plane at a certain critical $θ$ angle, $θ\sim 47^{\circ}$, the inflow changes its direction of radial motion and becomes wind. The driving forces are analyzed and found to be the centrifugal force and the gradient of gas and magnetic pressure. The properties of wind are also calculated. The specific angular momentum of wind is found to be significantly larger than that of inflow, thus wind can transfer angular momentum outward. These analytical results are compared to those obtained by the trajectory analysis based on MHD numerical simulation data and good agreements are found.

preprint2014arXiv

Does the circularization radius exist or not for low angular momentum accretion?

If the specific angular momentum of accretion gas at large radius is small compared to the local Keplerian value, one usually believes that there exists a "circularization radius" beyond which the angular momentum of accretion flow is almost a constant while within which a disk is formed and the angular momentum roughly follows the Keplerian distribution. In this paper, we perform numerical simulations to study whether the picture above is correct in the context of hot accretion flow. We find that for a steady accretion flow, the "circularization radius" does not exist and the angular momentum profile will be smooth throughout the flow. However, for transient accretion systems, such as the tidal disruption of a star by a black hole, a "turning point" should exist in the radial profile of the angular momentum, which is conceptually similar to the "circularization radius". At this radius, the viscous timescale equals the life time of the accretion event. The specific angular momentum is close to Keplerian within this radius, while beyond this radius the angular momentum is roughly constant.

preprint2014arXiv

Two dimensional numerical simulations of Supercritical Accretion Flows revisited

We study the dynamics of super-Eddington accretion flows by performing two-dimensional radiation-hydrodynamic simulations. Compared with previous works, in this paper we include the $T_{θϕ}$ component of the viscous stress and consider various values of the viscous parameter $α$. We find that when $T_{θϕ}$ is included, the rotational speed of the high-latitude flow decreases, while the density increases and decreases at the high and low latitudes, respectively. We calculate the radial profiles of inflow and outflow rates. We find that the inflow rate decreases inward, following a power law form of $\dot{M}_{\rm in}\propto r^s$. The value of $s$ depends on the magnitude of $α$ and is within the range of $\sim 0.4-1.0$. Correspondingly, the radial profile of density becomes flatter compared with the case of a constant $\dot{M}(r)$. We find that the density profile can be described by $ρ(r)\propto r^{-p}$, and the value of $p$ is almost same for a wide range of $α$ ranging from $α=0.1$ to $0.005$. The inward decrease of inflow accretion rate is very similar to hot accretion flows, which is attributed to the mass loss in outflows. To study the origin of outflow, we analyze the convective stability of slim disk. We find that depending on the value of $α$, the flow is marginally stable (when $α$ is small) or unstable (when $α$ is large). This is different from the case of hydrodynamical hot accretion flow where radiation is dynamically unimportant and the flow is always convectively unstable. We speculate that the reason for the difference is because radiation can stabilize convection. The origin of outflow is thus likely because of the joint function of convection and radiation, but further investigation is required.

preprint2013arXiv

On the role of initial and boundary conditions in numerical simulations of accretion flows

We study the effects of initial and boundary conditions, taking two-dimensional hydrodynamical numerical simulations of hot accretion flow as an example. The initial conditions considered include a rotating torus, a solution expanded from the one-dimensional global solution of hot accretion flows, injected gas with various angular momentum distributions, and the gas from a large-scale numerical simulation. Special attention is paid to the radial profiles of the mass accretion rate and density. Both can be described by a power-law function, $\dot{M}\propto r^s$ and $ρ\propto r^{-p}$. We find that if the angular momentum is not very low, the value of $s$ is not sensitive to the initial condition and lies within a narrow range, $0.47\la s \la 0.55$. However, the value of $p$ is more sensitive to the initial condition and lies in the range $0.48\la p \la 0.8$. The diversity of the density profile is because different initial conditions give different radial profiles of radial velocity due to the different angular momentum of the initial conditions. When the angular momentum of the accretion flow is very low, the inflow rate is constant with radius. Taking the torus model as an example, we have also investigated the effects of inner and outer boundary conditions by considering the widely adopted "outflow" boundary condition and the "mass flux conservation" condition. We find that the results are not sensitive to these two boundary conditions.

preprint2011arXiv

Magnetothermal and magnetorotational instabilities in hot accretion flows

In a hot, dilute, magnetized accretion flow, the electron mean-free path can be much greater than the Larmor radius, thus thermal conduction is anisotropic and along magnetic field lines. In this case, if the temperature decreases outward, the flow may be subject to a buoyancy instability (the magnetothermal instability, or MTI). The MTI amplifies the magnetic field, and aligns field lines with the radial direction. If the accretion flow is differentially rotating, the magnetorotational instability (MRI) may also be present. Using two-dimensional, time-dependent magnetohydrodynamic simulations, we investigate the interaction between these two instabilities. We use global simulations that span over two orders of magnitude in radius, centered on the region around the Bondi radius where the infall time of gas is longer than the growth time of both the MTI and MRI. Significant amplification of the magnetic field is produced by both instabilities, although we find that the MTI primarily amplifies the radial component, and the MRI primarily the toroidal component, of the field, respectively. Most importantly, we find that if the MTI can amplify the magnetic energy by a factor $F_t$, and the MRI by a factor $F_r$, then when the MTI and MRI are both present, the magnetic energy can be amplified by a factor of $F_t \cdot F_r$. We therefore conclude that amplification of the magnetic energy by the MTI and MRI operates independently. We also find that the MTI contributes to the transport of angular momentum, because radial motions induced by the MTI increase the Maxwell (by amplifying the magnetic field) and Reynolds stresses. Finally, we find that thermal conduction decreases the slope of the radial temperature profile. The increased temperature near the Bondi radius decreases the mass accretion rate.

preprint2010arXiv

On the convective instability of hot radiative accretion flows

How many fraction of gas available at the outer boundary can finally fall onto the black hole is an important question. It determines the observational appearance of accretion flows, and is also related with the evolution of black hole mass and spin. Previous two-dimensional hydrodynamical simulations of hot accretion flows find that the flow is convectively unstable because of its inward increase of entropy. As a result, the mass accretion rate decreases inward, i.e., only a small fraction of accretion gas can fall onto the black hole, while the rest circulates in the convective eddies or lost in convective outflows. Radiation is usually neglected in these simulations. In many cases, however, radiative cooling is important. In the regime of the luminous hot accretion flow (LHAF), radiative cooling is even stronger than the viscous dissipation. In the one dimensional case, this implies that the inward increase of entropy will become slower or the entropy even decreases inward in the case of an LHAF. We therefore expect that convective instability becomes weaker or completely disappears when radiative cooling is important. To examine the validity of this expectation, in this paper we perform two-dimensional hydrodynamical simulations of hot accretion flows with strong radiative cooling. We find that compared to the case of negligible radiation, convection only becomes slightly weaker. Even an LHAF is still strongly convectively unstable, its radial profile of accretion rate correspondingly changes little. We find the reason is that the entropy still increases inward in the two-dimensional case.