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
− w r ϕ ∼ ρ v s 2 ( v t v s ) + ρ v s 2 ( H 2 4 π ρ v s 2 ) = α ρ v s 2 -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
− w r ϕ ∼ ρ v s 2 ( v s v t ) + ρ v s 2 ( 4 π ρ v s 2 H 2 ) = α ρ 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 = 1 2 W r ϕ R d ω d R = 3 8 π M ˙ G M R 3 [ 1 − ( R 0 R ) 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}
Q = 2 1 W r ϕ R d R d ω = 8 π 3 M ˙ R 3 GM [ 1 − ( R R 0 ) 1/2 ] ( 2.6 )
Define parameters
m = M M ⊙ , m ˙ = M ˙ M ˙ cr = M ˙ 3 × 10 − 8 M ⊙ yr ( M ⊙ M ) 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)
m = M ⊙ M , m ˙ = M ˙ cr M ˙ = 3 × 1 0 − 8 yr M ⊙ M ˙ ( M M ⊙ )
r = R 3 R g = 1 6 R c 2 G M = M ⊙ M R 9 km r=\frac{R}{3R_g}=\frac{1}{6}\frac{Rc^2}{GM}=\frac{M_\odot}{M}\frac{R}{9\text{ km}}
r = 3 R g R = 6 1 GM R c 2 = M M ⊙ 9 km R
Radiation Spectrum of the Disk
Local Radiation Spectrum
In the outer r > 800 α 4 / 57 m − 46 / 57 m ˙ 37 / 57 r>800\alpha^{4/57}m^{-46/57}\dot{m}^{37/57} r > 800 α 4/57 m − 46/57 m ˙ 37/57 regions
free-free processes give the main contribution to the opacity a planckian spectrum of radiation
F ( x ) = B ( x ) = 2 π h c 2 ( k T h ) 3 x 3 e x − 1 , where x = h ν k T (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}
F ( x ) = B ( x ) = c 2 2 π h ( h k T ) 3 e x − 1 x 3 , where x = k T h ν ( 3.1 )
The corresponding flux of energy is
Q = ∫ F ( x ) d x = c 4 b T s 4 erg cm 2 s Q=\int F(x)\mathrm{d}x=\frac{c}{4}bT_s^4\,\frac{\text{erg}}{\text{cm}^2\text{ s}}
Q = ∫ F ( x ) d x = 4 c b T s 4 cm 2 s erg
In the intermediate region
800 α 4 / 57 m − 46 / 57 m ˙ 37 / 57 > r > 25 α 2 / 9 m ˙ 2 / 3 800 \alpha^{4/57}m^{-46/57}\dot{m}^{37/57}>r>25\alpha^{2/9}\dot{m}^{2/3} 800 α 4/57 m − 46/57 m ˙ 37/57 > r > 25 α 2/9 m ˙ 2/3
where Thomoson scattering dominates the opacity
In the case of a homogeneous medium with a sharp boundaryF ( x ) = 3 ϰ ( x ) n σ T m P B ( x ) ∼ const n T 5 / 4 x 3 / 2 e − x ( 1 − e − x ) 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}
F ( x ) = σ T m P 3 ϰ ( x ) n B ( x ) ∼ const n T 5/4 ( 1 − e − x ) 1/2 x 3/2 e − x ( 3.2 )
Q = 1.8 × 10 − 4 n T 2.25 erg cm 2 s Q=1.8\times10^{-4}\sqrt{n}T^{2.25}\,\frac{\text{erg}}{\text{cm}^2\text{ s}}
Q = 1.8 × 1 0 − 4 n T 2.25 cm 2 s erg
In the case of exponential varying atmosphere n = n ( z 1 ) e − z / H 0 n=n(z_1)e^{-z/H_0} n = n ( z 1 ) e − z / H 0 F ( x ) = ( 3 ϰ ( x ) σ T 2 m P 2 H 0 ) 1 / 3 B ( x ) ∼ const H 0 − 1 / 3 T 11 / 6 x 2 e − x ( 1 − e − x ) 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}
F ( x ) = ( σ T 2 m P 2 H 0 3 ϰ ( x ) ) 1/3 B ( x ) ∼ const H 0 − 1/3 T 11/6 x 2 ( 1 − e − x ) 1/3 e − x ( 3.3 )
Q = 1.3 × 10 4 H 0 − 1 / 3 T 17 / 6 erg cm 2 s ∼ T 2.5 Q=1.3\times10^4H_0^{-1/3}T^{17/6}\,\frac{\text{erg}}{\text{cm}^2\text{ s}}\sim T^{2.5}
Q = 1.3 × 1 0 4 H 0 − 1/3 T 17/6 cm 2 s erg ∼ T 2.5
the coefficient of free-free absorption
ϰ ( x ) = 4.1 × 10 − 23 ( 1 − e − x ) T 7 / 2 x 3 cm 5 \varkappa(x)=\frac{4.1\times10^{-23}(1-e^{-x})}{T^{7/2}x^3}\text{cm}^5
ϰ ( x ) = T 7/2 x 3 4.1 × 1 0 − 23 ( 1 − e − x ) cm 5
In the inner part of the disk r < 25 α 2 / 9 m ˙ 2 / 3 r<25\alpha^{2/9}\dot{m}^{2/3} r < 25 α 2/9 m ˙ 2/3
the processes of comptonization effect strongly the shape of the emitted spectrum
F ( x ) ∼ x 3 e − x , Q = c d ( r ) 4 T 4 , 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}
F ( x ) ∼ x 3 e − x , Q = 4 c d ( r ) T 4 , d ( r ) ≪ b . ( 3.4 )
Distribution of the Surface Temperature
The average energy of each photon
x ˉ = ∫ F ν d ν ∫ F ν h ν d ν = ⟨ h ν ⟩ k T \bar{x}=\frac{\int F_\nu\mathrm{d}\nu}{\int \dfrac{F_\nu}{h\nu}\mathrm{d}\nu}=\frac{\langle h\nu\rangle}{kT}
x ˉ = ∫ h ν F ν d ν ∫ F ν d ν = k T ⟨ h ν ⟩
In the outer regions x ˉ = 2.7 \bar{x}=2.7 x ˉ = 2.7 T s = 3 × 10 7 m − 1 / 4 m ˙ 1 / 4 r − 3 / 4 ( 1 − r − 1 / 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}
T s = 3 × 1 0 7 m − 1/4 m ˙ 1/4 r − 3/4 ( 1 − r − 1/2 ) 1/4 K ( 3.5 )
T s = [ 3 G M M ˙ 8 π σ R 3 ( 1 − R 0 R ) ] 1 / 4 T_s=\left[\frac{3GM\dot{M}}{8\pi\sigma R^3}\left(1-\sqrt{\frac{R_0}{R}}\right)\right]^{1/4}
T s = [ 8 π σ R 3 3 GM M ˙ ( 1 − R R 0 ) ] 1/4
In the intermedian region
according to (3.3) (exponential) x ˉ = 1.66 \bar{x}=1.66 x ˉ = 1.66 T s = 10 8 α 1 / 75 m ˙ 28 / 75 m − 19 / 75 r − 141 / 150 ( 1 − r − 1 / 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}
T s = 1 0 8 α 1/75 m ˙ 28/75 m − 19/75 r − 141/150 ( 1 − r − 1/2 ) 28/75 K ( 3.6 )
according to (3.2) (homogeneous) x ˉ = 1.2 \bar{x}=1.2 x ˉ = 1.2 T s = 1.4 × 10 9 α 2 / 9 m ˙ 8 / 9 m − 2 / 9 r − 5 / 3 ( 1 − r − 1 / 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}
T s = 1.4 × 1 0 9 α 2/9 m ˙ 8/9 m − 2/9 r − 5/3 ( 1 − r − 1/2 ) 8/9 K ( 3.7 )
In the inner region x ˉ = 3 \bar{x}=3 x ˉ = 3
if y = k T ( z 1 ) m e c 2 τ T 2 ( z 1 ) > 1 y=\dfrac{kT(z_1)}{m_ec^2}\tau_T^2(z_1)>1 y = m e c 2 k T ( z 1 ) τ T 2 ( z 1 ) > 1 T s = 1.4 × 10 9 A − 2 / 9 α 2 / 9 m ˙ 8 / 9 m − 2 / 9 r − 5 / 3 ( 1 − r − 1 / 2 ) 8 / 9 ≃ 5 × 10 8 α 1 / 5 m ˙ 4 / 5 m − 1 / 5 r − 3 / 2 ( 1 − r − 1 / 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}
T s = 1.4 × 1 0 9 A − 2/9 α 2/9 m ˙ 8/9 m − 2/9 r − 5/3 ( 1 − r − 1/2 ) 8/9 ≃ 5 × 1 0 8 α 1/5 m ˙ 4/5 m − 1/5 r − 3/2 ( 1 − r − 1/2 ) 4/5 K ( 3.8 )
if τ ∗ ( u 0 ) = τ f f τ s < 1 \tau^*(u_0)=\sqrt{\tau_{ff}\tau_s}<1 τ ∗ ( u 0 ) = τ f f τ s < 1 T s = 10 14 α 12 / 5 m ˙ 24 / 5 r − 36 / 5 ( 1 − r − 1 / 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}
T s = 1 0 14 α 12/5 m ˙ 24/5 r − 36/5 ( 1 − r − 1/2 ) 24/5 ( 3.9 )
Integral Spectrum of the Outgoing Disk Radiation
Integrating the local spectrum
I ν = 4 π ∫ R 0 R 1 F ν [ T s ( R ) ] R d R I_\nu=4\pi \int_{R_0}^{R_1}F_\nu[T_s(R)]R\mathrm{d}R
I ν = 4 π ∫ R 0 R 1 F ν [ T s ( R )] R d R
where R 1 R_1 R 1 is the external boundary of the disk.
For a local planckian spectrum and dependence T s ( R ) T_s(R) T s ( R ) (3.5), we obtain for ν ≪ k T 0 h \nu\ll\dfrac{kT_0}{h} ν ≪ h k T 0
I ν = 16 π 2 R 0 2 h c 2 ( k T 0 h ) 8 / 3 ν 1 / 3 I_\nu=\frac{16\pi^2 R_0^2 h}{c^2}\left(\frac{kT_0}{h}\right)^{8/3}\nu^{1/3}
I ν = c 2 16 π 2 R 0 2 h ( h k T 0 ) 8/3 ν 1/3
The spectral index of radiation γ = d ln I ν d ln ν = 1 3 \gamma=\dfrac{\mathrm{d}\ln I_\nu}{\mathrm{d}\ln \nu}=\dfrac{1}{3} γ = d ln ν d ln I ν = 3 1
For r ≫ 1 r\gg1 r ≫ 1 ,
γ = 0.07 \gamma=0.07 γ = 0.07 for spectrum (3.3) and temperature (3.6), (intermediate, exponential)
γ = 0.04 \gamma=0.04 γ = 0.04 for spectrum (3.2) and temperature (3.7), (intermediate, homogeneous)
γ = − 1 / 3 \gamma=-1/3 γ = − 1/3 for spectrum (3.4) and temperature (3.8), (inner)
γ = − 1 \gamma=-1 γ = − 1 for spectrum (3.4) and temperature (3.9). (inner)