Vertical Structure of Radiation-Pressure-Dominated Thin Disks: Link between Vertical Advection and Convective Stability

TL;DR

This paper investigates how vertical advection due to magnetic buoyancy influences the vertical structure and convective stability of radiation-pressure-dominated thin disks, reconciling classical theory with recent simulation results.

Contribution

It demonstrates that vertical advection significantly contributes to energy transport and stabilizes the disk against convection, providing new insights into disk structure.

Findings

Vertical advection enhances energy transport in thin disks.

Vertical advection stabilizes the disk against convective instability.

The study links energy transport mechanisms with convective stability.

Abstract

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

Equations

$$P_{\mathrm{gas}}=\frac{\rho k_{\mathrm{B}}T}{\mu m_{\mathrm{p}}},$$

$$P_{\mathrm{rad}}=\frac{1}{3}a{T}^4,$$

$$\frac{d{P}_{\mathrm{rad}}}{dz}=-\frac{{\kappa }_{\mathrm{es}}\rho F_z}{c},$$

$$\frac{d{P}_{\mathrm{gas}}}{dz}=-\rho {\Omega }_{\mathrm{K}}^{2}z+\frac{{\kappa }_{\mathrm{es}}\rho F_z}{c},$$

$$\frac{d{F}_z}{dz}=-{\tau }_{r\varphi }{\Omega }_{\mathrm{K}}{g}_{\ast }-\frac{d{F}_{\mathrm{adv}}}{dz},$$

$$S=S_{\mathrm{gas}}+S_{\mathrm{rad}}=\frac{3}{2}\frac{k_{\mathrm{B}}}{\mu m_{\mathrm{p}}}\ln \left(\frac{P_{\mathrm{gas}}}{{\rho }^{5/3}}\right)+\frac{4}{3}\frac{a{T}^3}{\rho }+const.$$

$$\frac{dS}{dz}>0(z>0),$$

$${\nabla }_{\mathrm{rad}}<{\nabla }_{\mathrm{ad}},$$

$${\nabla }_{\mathrm{rad}}\equiv {\left(\frac{\partial \ln T}{\partial \ln P}\right)}_{\mathrm{rad}},{\nabla }_{\mathrm{ad}}\equiv {\left(\frac{\partial \ln T}{\partial \ln P}\right)}_S.$$

$${\nabla }_{\mathrm{ad}}=\frac{4-3\beta }{1.5{\beta }^2+12\beta (1-\beta )+(4-3\beta {)}^2},$$

$$\frac{d\ln T}{d\ln P}=\min ({\nabla }_{\mathrm{rad}},\lambda {\nabla }_{\mathrm{ad}}),$$

$$F_z=0,$$

$$F_{\mathrm{adv}}=0.$$

$$-{\kappa }_{\mathrm{es}}{\rho }^2{\left(\frac{d\rho }{dz}\right)}^{-1}=1,$$

$$F_z=2\sigma T^4,$$

$$F_z+F_{\mathrm{adv}}=\frac{1}{4\pi }\dot{M}{\Omega }_{\mathrm{K}}^2{f}_{\ast }{g}_{\ast },$$

$$H=\frac{3}{8\pi }\frac{{\kappa }_{\mathrm{es}}\dot{M}}{c},$$

Full text

accretion, accretion disks — black hole physics — convection — instabilities

1 Introduction

The standard thin accretion disk model under the alpha description was constructed by Shakura and Sunyaev (Shakura & Sunyaev, 1973), which has been successfully applied to X-ray binaries and active galactic nuclei (for a review, see Frank et al., 2002; Kato et al., 2008). However, when we compare the theoretical results with observations and simulations, there exist some basic problems such as thermal stability and convective stability. Recent simulations (Hirose et al., 2009; Jiang et al., 2013) on thin disks based on a shearing box showed that, for high mass accretion rates around $0.1\dot{M}_{\rm Edd}$, where $\dot{M}_{\rm Edd}$ is the Eddington accretion rate, the inner disk is radiation pressure dominated, and the vertical advection can be an efficient process for energy transport, which is probably related to the magnetic buoyancy (Hirose et al., 2009; Jiang et al., 2013). In addition, the local convective stability criterion, $dS/dz>0$, where $S$ is the entropy, is well satisfied according to the simulation results. On the other hand, in the classic theory of thin disks, the vertical energy transport is completely dominated by the diffusion process. Moreover, when the radiation pressure dominates over the gas pressure, the local convective stability criterion $dS/dz>0$ is not satisfied and therefore the disk may suffer convective instability (Sa̧dowski et al., 2011). Then a question arises that, is there any link between the vertical advection and the convective stability?

In the present work, we will revisit the vertical structure and energy transport of radiation-pressure-dominated thin disks by including the possible vertical advection process. In addition, based on the simulation results, we will take into account the local convective stability criterion to modify the set of equations. The rest part is organized as follows. Equations and boundary conditions are described in Section 2. Numerical results are shown in Section 3. Conclusions and discussion are made in Section 4.

2 Equations and boundary conditions

The set of equations for the vertical structure of thin disks is based on the alpha-stress assumption. The gas pressure $P_{\rm gas}$ and the radiation pressure $P_{\rm rad}$ and their derivates take the following expressions:

[TABLE]

[TABLE]

[TABLE]

[TABLE]

where $\rho$ is the mass density, $T$ is the temperature, $k_{\rm B}$ is the Boltzmann constant, $m_{\rm p}$ is the mass of a proton. The Keplerian angular velocity $\Omega_{\rm K}$ is written as $\Omega_{\rm K}=\sqrt{GM_{\rm BH}/R(R-R_{\rm g})^{2}}$ under the well-known Paczyński-Wiita potential $\Psi_{\rm PW}=-GM_{\rm BH}/(R-R_{\rm g})$ (Paczyński & Wiita, 1980), where the gravitational radius is defined as $R_{\rm g}\equiv 2GM_{\rm BH}/{c^{2}}$.

In the standard disk model, the angular velocity is assumed to be Keplerian, i.e., $\Omega=\Omega_{\rm K}$. Moreover, the local energy balance at each radius is that the viscous heating rate equals the radiative cooling rate, and the vertical energy transport is dominated by the photon diffusion process. Here, following the spirit of simulation results, the vertical advection may have an additional contribution to the energy transport. Thus, the vertical flux $F_{z}$ due to the diffusion is expressed as

[TABLE]

where $F_{\rm adv}$ is the vertical flux owing to the vertical advection process, and the parameter $g_{\ast}$ takes the form $g_{\ast}\equiv-d\ln\Omega_{\rm K}/d\ln R=3/2+1/(R/R_{\rm g}-1)$. The $r\phi$-component of the shear stress $\tau_{r\phi}$ is assumed to be proportional to the total pressure, i.e., $\tau_{r\phi}=-\alpha(P_{\rm gas}+P_{\rm rad})$, where $\alpha$ is a constant parameter.

The total entropy $S$ is the sum of gas and radiation, which is written as

[TABLE]

It is well-known that the local convective stability criterion can be expressed as

[TABLE]

which is equivalent to the relationship between the radiative gradient $\nabla_{\rm rad}$ and the adiabatic gradient $\nabla_{\rm ad}$:

[TABLE]

where

[TABLE]

With the above expression of the total entropy $S$, we can derive the explicit form of $\nabla_{\rm ad}$ as

[TABLE]

where $\beta$ is defined as $\beta\equiv P_{\rm gas}/(P_{\rm gas}+P_{\rm rad})$.

As mentioned in the first section, simulations have not shown convective instability, which indicates that the relationship $\nabla_{\rm rad}\leqslant\nabla_{\rm ad}$ may be well satisfied in the disk. Following this spirit, we assume the thermodynamical gradient to be

[TABLE]

where $\lambda\leqslant 1$ can guarantee the disk to be convectively stable.

The system consists of six equations, Equations (1-6), for the six unknown variables $\rho$, $T$, $P_{\rm gas}$, $P_{\rm rad}$, $F_{z}$, and $F_{\rm adv}$. There are four first-order differential equations in this system. In addition, the position of photosphere $H_{\rm phot}$ is unknown. Thus, totally five boundary conditions are required to solve the system between the equatorial plane ($z=0$) and the photosphere ($z=H_{\rm phot}$).

On the equatorial plane there exist two natural boundary conditions:

[TABLE]

[TABLE]

At the photosphere, the other three boundary conditions can be derived:

[TABLE]

[TABLE]

[TABLE]

where Equation (9) can be regarded as the definition of the photosphere position. The five boundary conditions (7-11) together with Equations (1-6) enable us to solve the system and derive the vertical structure.

3 Numerical results

Following Hirose et al. (2009), we define a vertical thickness $H$ as

[TABLE]

which is used as a scale of the height. For simplicity, we choose $\lambda=1$ in Equation (6) for our numerical calculations. The other parameters are $M_{\rm BH}=10M_{\sun}$, $\dot{M}=0.1\dot{M}_{\rm Edd}$, and $\alpha=0.01$.

In order to directly compare the structure including vertical advection process with that in the classic picture, we made numerical calculations for both of these two models. Figure 1 shows the vertical density profile, where the solid line corresponds to the results including vertical advection, and the dashed line corresponds to the results of the classic model. The solid line has a peak at $z=0$ whereas the dashed line has two peaks at $z/H\approx\pm 0.3$. The peculiar shape of the dashed line, i.e., increasing $\rho$ with $z$ in the range $0<z/H<0.3$, indicates that the disk suffers convective instability. The well-known local convective stability criterion, $dS/dz>0$ for $z>0$ (or $dS/dz<0$ for $z<0$) may work well in geometrically thin disks, where the radial velocity is low and therefore the advection effects may be negligible. Even though the entropy profile is not plotted in Figure 1, it is obvious that the entropy $S$ will decrease with increasing $z$ in the range $0<z/H<0.3$ due to the peculiar profile of density. Thus, the local convective stability criterion is not satisfied and therefore the disk is likely to be convectively unstable, as previously investigated by Sa̧dowski et al. (2011). On the contrary, the profile of the density in the case including vertical advection (solid line) is similar to that of simulations (Hirose et al., 2009), which shows a continuously decreasing density with increasing $z$ for $z>0$. Thus, a convectively stable disk is quite possible.

Figure 2 provides the vertical profiles of flux, where the black dashed line corresponds to the diffusive flux $F_{z}$ in the classic thin disk model, and the two solid lines correspond to the case including the vertical advection, where the black and red lines show the variations of the diffusive flux $F_{z}$ and the advection flux $F_{\rm adv}$, respectively. It is seen that the vertical advection has significant contribution to the total radiation flux, in particular for the region near the equatorial plane ($0<z/H<0.1$). In some regions on the right part ($z>0$), the red line is even higher than the black solid line, which means that $F_{\rm adv}$ can dominate over $F_{z}$. In addition, the profile of $F_{\rm adv}$ in Figure 2 is quite similar to that in simulations (e.g., Hirose et al., 2009). Since our numerical calculation is based on a convectively stable disk, the results may indicate a link between the energy transport due to vertical advection and the convectively stable property. The physical reason for the link is probably that, the additional energy transport decreases the entropy in the region near the equatorial plane and therefore the entropy can keep to increase with increasing vertical height, which satisfies the local convective stability criterion.

Taking $\dot{M}=0.1\dot{M}_{\rm Edd}$ as a typical accretion rate, we also investigate the strength of vertical advection for different radii. Obviously, the ratio of the gas pressure to the total pressure $\beta$ will increase with increasing radius for a fixed $\dot{M}$. In the outer part where gas pressure is dominant, the vertical advection may be negligible. In other words, the diffusion process is the dominant energy transport and the disk is well convectively stable. However, in the inner part where the radiation pressure dominates over the gas pressure, the diffusion and the vertical advection may both be of importance for the energy transport.

Figure 3 shows the profiles of diffusive flux (black) and advective flux (red) at the three radii $R=10R_{\rm g}$ (solid lines), $100R_{\rm g}$ (dashed lines), and $1000R_{\rm g}$ (dotted line). Here, we take the height of photosphere $H_{\rm phot}$ as the length unit. The red solid line shows that, at $R=10R_{\rm g}$, the advective flux $F_{\rm adv}$ covers the range from the equatorial plane ($z=0$) to the photosphere ($z=H_{\rm phot}$). At a larger radius $R=100R_{\rm g}$, the red dashed line shows a smaller vertical range for $F_{\rm adv}$ ($\lesssim 2H_{\rm phot}/3$). Moreover, for a sufficiently large radius $R=1000R_{\rm g}$, Figure 3 shows that the advective flux disappears and the diffusion is the only mechanism for the vertical energy transport, as shown by the dotted line. Thus, for thin disks, the two conditions for the occurrence of vertical advection are high accretion rates $\dot{M}\gtrsim 0.1M_{\sun}~{}{\rm s}^{-1}$ and small radii $R\lesssim 100R_{\rm g}$.

4 Conclusions and discussion

In this work, we have revisited the vertical structure of radiation-pressure-dominated thin disks by taking into account the role of vertical advection process. Our study has shown that the vertical advection not only has significant contribution to the energy transport, but also plays an important role to make the disk convectively stable. Thus, a link may exist between the vertical advective energy transport and the convectively stable property. The physical reason for the link is probably that, the additional energy transport decreases the entropy in the region near the equatorial plane and therefore the entropy can keep to increase with increasing vertical height, which satisfies the local convective stability criterion. Our work may be helpful to understand the discrepancy between the classic theory and simulation results.

We would point out that, the detailed study of convective stability may require global rather than local stability analyses. For example, Abramowicz et al. (1993) demonstrated that when the viscosity is taken fully into account, stability analyses cannot be discussed within the framework of a local analysis, and a fully global treatment is required. Moreover, the global stability analyses of vertical convection of a thin gaseous disk were performed by several works (e.g., Ruden et al., 1988). On the other hand, it is known that the advection may play an important role in stabilizing the disk against dynamic, thermal, and viscous perturbations. For instance, the geometrically thick disk without radial motion (Abramowicz et al., 1980; Paczyński & Wiita, 1980) may suffer the Papaloizou-Pringle instability (Papaloizou & Pringle, 1984), which is a dynamic instability based on the acoustic perturbations propagating between two boundaries in a differential rotating system. Later, by including the advection terms, Blaes (1987) found that all the unstable modes for the purely rotating flow are quickly stabilized by the advection process. Furthermore, Abramowicz et al. (1988) proposed the well-known slim disk model (or named as the optically thick, advection-dominated accretion disk) for super-Eddington accretion systems. The slim disk was found to be dynamically stable, thermally stable, and viscously stable, which may be related to the radial advection process. In the present work, we focus on the stability of geometrically thin disks, where the radial velocity ($v_{R}\approx-\alpha(H/R)^{2}v_{\rm K}$, where $v_{\rm K}$ is the Keplerian velocity) is quite low, and therefore the effects of advection may also be quite weak.

In a previous work, Gu (2012) showed that, for high mass accretion rates $\dot{M}\gtrsim\dot{M}_{\rm Edd}$ where the radiation pressure completely dominates over the gas pressure, the disk is likely to be convectively stable without including the energy transport through the vertical advection. The physics of the stable property is probably related to the radial advection effects and the geometrically thick structure. On the other hand, in recent years many global simulation works have been done on the super-Eddington accretion flows (e.g., Ohsuga et al., 2005; Jiang et al., 2014; Sa̧dowski & Narayan, 2015). The simulations of Jiang et al. (2014) revealed the importance of the vertical advection, which can essentially enhance the radiative efficiency. In addition, outflows may play another important role in such flows (Jiang et al., 2014; Sa̧dowski & Narayan, 2015). In our opinion, the theory of super-Eddington accretion flows is worth further investigation.

The authors would thank Yan-Fei Jiang for providing the vertical profile of entropy in simulations, and thank the referee for helpful comments that improved the paper. This work was supported by the National Basic Research Program of China (973 Program) under grants 2014CB845800, the National Natural Science Foundation of China under grants 11573023, 11333004, and 11222328, and the CAS Open Research Program of Key Laboratory for the Structure and Evolution of Celestial Objects under grant OP201503.