Black Holes in Binary Systems. Observational Appearance

Shakura, N. I., & Sunyaev, R. A. (1973). Black holes in binary systems. Observational appearance. Astronomy & Astrophysics, 24, 337–355.

Mechanisms of Angular Momentum Transfer

Magnetic field, turbulence, molecular and radiative viscosity.

The tangential stresses due to the turbulence and magnetic field

wrϕρvs2(vtvs)+ρvs2(H24πρvs2)=αρvs2-w_{r\phi}\sim\rho v_s^2\left(\frac{v_t}{v_s}\right)+\rho v_s^2\left(\frac{H^2}{4\pi\rho v_s^2}\right)=\alpha\rho v_s^2

can be characterized by only one parameter α\alpha.

The Structure of the Disk

The energy flux, radiated from surface unit of the disk in unit of the time

Q=12WrϕRdωdR=38πM˙GMR3[1(R0R)1/2](2.6)\begin{aligned} Q&=\frac{1}{2}W_{r\phi} R\frac{\mathrm{d}\omega}{\mathrm{d}R} \\ &=\frac{3}{8\pi}\dot{M}\frac{GM}{R^3}\left[1-\left(\frac{R_0}{R}\right)^{1/2}\right] \end{aligned}\tag{2.6}

Define parameters

m=MM,m˙=M˙M˙cr=M˙3×108Myr(MM)m=\frac{M}{M_\odot},\quad\dot{m}=\frac{\dot{M}}{\dot{M}_\text{cr}}=\frac{\dot{M}}{3\times10^{-8}\dfrac{M_\odot}{\text{yr}}}\left(\frac{M_\odot}{M}\right)

r=R3Rg=16Rc2GM=MMR9 kmr=\frac{R}{3R_g}=\frac{1}{6}\frac{Rc^2}{GM}=\frac{M_\odot}{M}\frac{R}{9\text{ km}}

Radiation Spectrum of the Disk

Local Radiation Spectrum

  1. In the outer r>800α4/57m46/57m˙37/57r>800\alpha^{4/57}m^{-46/57}\dot{m}^{37/57} regions
    free-free processes give the main contribution to the opacity a planckian spectrum of radiation

    F(x)=B(x)=2πhc2(kTh)3x3ex1,where x=hνkT(3.1)F(x)=B(x)=\frac{2\pi h}{c^2}\left(\frac{kT}{h}\right)^3\frac{x^3}{e^x-1},\quad\text{where }x=\frac{h\nu}{kT}\tag{3.1}

    The corresponding flux of energy is

    Q=F(x)dx=c4bTs4ergcm2 sQ=\int F(x)\mathrm{d}x=\frac{c}{4}bT_s^4\,\frac{\text{erg}}{\text{cm}^2\text{ s}}

  2. In the intermediate region
    800α4/57m46/57m˙37/57>r>25α2/9m˙2/3800 \alpha^{4/57}m^{-46/57}\dot{m}^{37/57}>r>25\alpha^{2/9}\dot{m}^{2/3}
    where Thomoson scattering dominates the opacity

    1. In the case of a homogeneous medium with a sharp boundary

      F(x)=3ϰ(x)nσTmPB(x)const nT5/4x3/2ex(1ex)1/2(3.2)F(x)=\sqrt{\frac{3\varkappa(x)n}{\sigma_Tm_P}}B(x)\sim\text{const }\sqrt{n}T^{5/4}\frac{x^{3/2}e^{-x}}{(1-e^{-x})^{1/2}}\tag{3.2}

      Q=1.8×104nT2.25ergcm2 sQ=1.8\times10^{-4}\sqrt{n}T^{2.25}\,\frac{\text{erg}}{\text{cm}^2\text{ s}}

    2. In the case of exponential varying atmosphere n=n(z1)ez/H0n=n(z_1)e^{-z/H_0}

      F(x)=(3ϰ(x)σT2mP2H0)1/3B(x)const H01/3T11/6x2ex(1ex)1/3(3.3)\begin{aligned} F(x)&=\left(\frac{3\varkappa(x)}{\sigma_T^2m_P^2H_0}\right)^{1/3}B(x) \\ &\sim\text{const }H_0^{-1/3}T^{11/6}x^2\frac{e^{-x}}{(1-e^{-x})^{1/3}} \end{aligned}\tag{3.3}

      Q=1.3×104H01/3T17/6ergcm2 sT2.5Q=1.3\times10^4H_0^{-1/3}T^{17/6}\,\frac{\text{erg}}{\text{cm}^2\text{ s}}\sim T^{2.5}

    the coefficient of free-free absorption

    ϰ(x)=4.1×1023(1ex)T7/2x3cm5\varkappa(x)=\frac{4.1\times10^{-23}(1-e^{-x})}{T^{7/2}x^3}\text{cm}^5

  3. In the inner part of the disk r<25α2/9m˙2/3r<25\alpha^{2/9}\dot{m}^{2/3}
    the processes of comptonization effect strongly the shape of the emitted spectrum

    F(x)x3ex,Q=cd(r)4T4,d(r)b.(3.4)F(x)\sim x^3e^{-x},\quad Q=\frac{cd(r)}{4}T^4,\quad d(r)\ll b.\tag{3.4}

Distribution of the Surface Temperature

The average energy of each photon

xˉ=FνdνFνhνdν=hνkT\bar{x}=\frac{\int F_\nu\mathrm{d}\nu}{\int \dfrac{F_\nu}{h\nu}\mathrm{d}\nu}=\frac{\langle h\nu\rangle}{kT}

  1. In the outer regions xˉ=2.7\bar{x}=2.7

    Ts=3×107m1/4m˙1/4r3/4(1r1/2)1/4 K(3.5)T_s=3\times10^7 m^{-1/4}\dot{m}^{1/4}r^{-3/4}(1-r^{-1/2})^{1/4}\text{ K}\tag{3.5}

    Ts=[3GMM˙8πσR3(1R0R)]1/4T_s=\left[\frac{3GM\dot{M}}{8\pi\sigma R^3}\left(1-\sqrt{\frac{R_0}{R}}\right)\right]^{1/4}

  2. In the intermedian region
    1. according to (3.3) (exponential) xˉ=1.66\bar{x}=1.66

      Ts=108α1/75m˙28/75m19/75r141/150(1r1/2)28/75 K(3.6)T_s=10^8\alpha^{1/75}\dot{m}^{28/75}m^{-19/75}r^{-141/150}(1-r^{-1/2})^{28/75}\text{ K}\tag{3.6}

    2. according to (3.2) (homogeneous) xˉ=1.2\bar{x}=1.2

      Ts=1.4×109α2/9m˙8/9m2/9r5/3(1r1/2)8/9 K(3.7)T_s=1.4\times10^9\alpha^{2/9}\dot{m}^{8/9}m^{-2/9}r^{-5/3}(1-r^{-1/2})^{8/9}\text{ K}\tag{3.7}

  3. In the inner region xˉ=3\bar{x}=3
    1. if y=kT(z1)mec2τT2(z1)>1y=\dfrac{kT(z_1)}{m_ec^2}\tau_T^2(z_1)>1

      Ts=1.4×109A2/9α2/9m˙8/9m2/9r5/3(1r1/2)8/95×108α1/5m˙4/5m1/5r3/2(1r1/2)4/5 K(3.8)\begin{aligned} T_s&=1.4\times10^9 A^{-2/9}\alpha^{2/9}\dot{m}^{8/9}m^{-2/9}r^{-5/3}(1-r^{-1/2})^{8/9} \\ &\simeq 5\times10^8\alpha^{1/5}\dot{m}^{4/5}m^{-1/5}r^{-3/2}(1-r^{-1/2})^{4/5}\text{ K} \end{aligned}\tag{3.8}

    2. if τ(u0)=τffτs<1\tau^*(u_0)=\sqrt{\tau_{ff}\tau_s}<1

      Ts=1014α12/5m˙24/5r36/5(1r1/2)24/5(3.9)T_s=10^{14}\alpha^{12/5}\dot{m}^{24/5}r^{-36/5}(1-r^{-1/2})^{24/5}\tag{3.9}

Integral Spectrum of the Outgoing Disk Radiation

Integrating the local spectrum

Iν=4πR0R1Fν[Ts(R)]RdRI_\nu=4\pi \int_{R_0}^{R_1}F_\nu[T_s(R)]R\mathrm{d}R

where R1R_1 is the external boundary of the disk.

For a local planckian spectrum and dependence Ts(R)T_s(R) (3.5), we obtain for νkT0h\nu\ll\dfrac{kT_0}{h}

Iν=16π2R02hc2(kT0h)8/3ν1/3I_\nu=\frac{16\pi^2 R_0^2 h}{c^2}\left(\frac{kT_0}{h}\right)^{8/3}\nu^{1/3}

The spectral index of radiation γ=dlnIνdlnν=13\gamma=\dfrac{\mathrm{d}\ln I_\nu}{\mathrm{d}\ln \nu}=\dfrac{1}{3}

For r1r\gg1,
γ=0.07\gamma=0.07 for spectrum (3.3) and temperature (3.6), (intermediate, exponential)
γ=0.04\gamma=0.04 for spectrum (3.2) and temperature (3.7), (intermediate, homogeneous)
γ=1/3\gamma=-1/3 for spectrum (3.4) and temperature (3.8), (inner)
γ=1\gamma=-1 for spectrum (3.4) and temperature (3.9). (inner)


Black Holes in Binary Systems. Observational Appearance
http://jingliangwei.github.io/blog-hexo/2026/07/03/Black-Holes-in-Binary-Systems-Observational-Appearance/
Author
Arwell
Posted on
July 3, 2026
Licensed under