大光深下的辐射转移
由之前的结论,能流密度表示为:
F
⃗
=
−
c
3
κ
ρ
∇
U
\vec F=-\frac{c}{3\kappa\rho}\nabla U
F=−3κρc∇U
其中
κ
\kappa
κ为单位质量辐射跟物质作用的截面。这个作用和频率有关,不同的频率对应的能流不一样:
F
⃗
ν
=
−
c
3
κ
ν
ρ
∇
U
ν
\vec F_\nu = -\frac{c}{3\kappa_\nu\rho}\nabla U_\nu
Fν=−3κνρc∇Uν
上述的方程也可以从平衡态辐射转移方程在
τ
ν
≫
1
\tau_\nu\gg 1
τν≫1时近似求解得到。平衡态辐射转移方程为:
μ
d
I
ν
d
τ
ν
=
I
ν
−
S
ν
\mu\frac{d I_\nu}{d\tau_\nu}=I_\nu-S_\nu
μdτνdIν=Iν−Sν
其中
I
ν
I_\nu
Iν为单位面积单位时间,沿着单位立体角的辐射强度;
μ
=
cos
θ
\mu=\cos\theta
μ=cosθ。在光深
τ
μ
≫
1
\tau_\mu\gg 1
τμ≫1时,辐射强度几乎是各向同性的,此时把辐射场按
μ
\mu
μ展开:
I
ν
(
μ
)
=
I
0
+
I
1
μ
+
I
2
μ
2
+
…
I_\nu(\mu)=I_0+I_1\mu+I_2\mu^2+\dots
Iν(μ)=I0+I1μ+I2μ2+…
在
τ
ν
\tau_\nu
τν很大时,
S
ν
=
B
ν
S_\nu=B_\nu
Sν=Bν,将展开后的辐射强度代回到辐射转移方程中,得到:
∑
n
μ
n
+
1
d
I
n
d
τ
ν
=
∑
n
μ
n
I
n
−
B
ν
\sum_n\mu^{n+1}\frac{dI_n}{d\tau_\nu}=\sum_n\mu^n I_n-B_\nu
n∑μn+1dτνdIn=n∑μnIn−Bν
从而求出:
I
0
=
B
ν
,
I
1
=
d
B
ν
d
τ
ν
,
I
2
=
d
2
B
ν
d
τ
ν
2
,
…
I_0=B_\nu,\ I_1=\frac{dB_\nu}{d\tau_\nu},\ I_2=\frac{d^2 B_\nu}{d\tau_\nu^2},\ \dots
I0=Bν, I1=dτνdBν, I2=dτν2d2Bν, …
辐射能流与辐射强度的关系(辐射强度保留到最低的非零阶)(第零阶怎么没了?最后一步怎么变过去的?):
F
ν
,
z
=
∬
I
ν
μ
d
Ω
≈
2
π
∫
−
1
1
d
B
ν
d
τ
ν
μ
2
d
μ
=
4
π
3
d
B
ν
d
τ
ν
=
−
c
3
κ
ν
ρ
d
U
ν
d
z
F_{\nu,z}=\iint I_\nu\mu d\Omega\approx 2\pi\int_{-1}^{1}\frac{dB_\nu}{d\tau_\nu}\mu^2 d\mu=\frac{4\pi}{3}\frac{dB_\nu}{d\tau_\nu}=-\frac{c}{3\kappa_\nu\rho}\frac{dU_\nu}{dz}
Fν,z=∬IνμdΩ≈2π∫−11dτνdBνμ2dμ=34πdτνdBν=−3κνρcdzdUν
对所有的频率进行积分:
F
⃗
=
−
∫
0
∞
c
3
κ
ν
ρ
∇
U
ν
d
ν
=
−
c
3
ρ
∫
0
∞
∇
U
ν
κ
ν
d
ν
=
−
4
π
3
ρ
(
∫
0
∞
1
κ
ν
d
B
ν
d
T
d
ν
)
∇
T
\begin{aligned} \vec F&=-\int^\infty_0\frac{c}{3\kappa_\nu\rho}\nabla U_\nu d\nu\\ &=-\frac{c}{3\rho}\int^\infty_0\frac{\nabla U_\nu}{\kappa_\nu}d\nu\\ &=-\frac{4\pi}{3\rho}\left(\int^\infty_0\frac{1}{\kappa_\nu}\frac{dB_\nu}{dT}d\nu\right)\nabla T \end{aligned}
F=−∫0∞3κνρc∇Uνdν=−3ρc∫0∞κν∇Uνdν=−3ρ4π(∫0∞κν1dTdBνdν)∇T
其中:
∫
0
∞
d
B
ν
d
T
d
ν
=
d
d
T
∫
0
∞
B
ν
d
ν
=
a
c
π
T
3
\int^\infty_0\frac{dB_\nu}{dT}d\nu=\frac{d}{dT}\int^\infty_0 B_\nu d\nu=\frac{ac}{\pi}T^3
∫0∞dTdBνdν=dTd∫0∞Bνdν=πacT3
由此引入Rosseland平均值:
1
κ
ˉ
≡
∫
0
∞
1
κ
ν
d
B
ν
d
T
d
ν
∫
0
∞
d
B
ν
d
T
d
ν
=
π
a
c
T
3
∫
0
∞
1
κ
ν
d
B
ν
d
T
d
ν
\frac{1}{\bar \kappa}\equiv \frac{\int^\infty_0\frac{1}{\kappa_\nu}\frac{dB_\nu}{dT}d\nu}{\int^\infty_0\frac{dB_\nu}{dT}d\nu}=\frac{\pi}{acT^3}\int^\infty_0\frac{1}{\kappa_\nu}\frac{dB_\nu}{dT}d\nu
κˉ1≡∫0∞dTdBνdν∫0∞κν1dTdBνdν=acT3π∫0∞κν1dTdBνdν
从而得到温度梯度和辐射能流的关系:
F
=
−
4
a
c
T
3
3
ρ
κ
ˉ
d
T
d
r
F=-\frac{4acT^3}{3\rho\bar\kappa}\frac{dT}{dr}
F=−3ρκˉ4acT3drdT
如果
κ
ν
∝
ν
−
n
\kappa_\nu\propto \nu^{-n}
κν∝ν−n,则有
κ
ˉ
∝
T
−
n
\bar \kappa\propto T^{-n}
κˉ∝T−n。如果
κ
ν
=
a
T
α
ν
β
\kappa_\nu=aT^\alpha\nu^\beta
κν=aTανβ,则Rosseland平均为:
κ
R
∝
T
α
+
β
\kappa_R\propto T^{\alpha+\beta}
κR∝Tα+β
自由-自由吸收
自由-自由吸收:光子把能量给了自由电子
物质在吸收能量后,从初始态
i
i
i到终态
f
f
f。这之间辐射场的吸收从量子的角度来看,辐射可以看做电磁波(平面波),电场强度为
E
⃗
\vec E
E,辐射又有振荡。在辐射场的作用下,粒子发生了状态变化。原子(或离子)的这个过程可以用微扰的方式写出,表示为原有的哈密顿加上电场的微扰:
H
+
H
i
n
t
H+H^{in t}
H+Hint
计算跃迁概率:
⟨
i
∣
H
i
n
t
∣
f
⟩
\langle i|H^{int}|f\rangle
⟨i∣Hint∣f⟩
与电场的相互作用:
H
i
n
t
∼
∑
i
q
i
E
⃗
⋅
r
⃗
=
E
⃗
⋅
d
⃗
H^{in t}\sim \sum_i q_i \vec E\cdot \vec r=\vec E\cdot \vec d
Hint∼i∑qiE⋅r=E⋅d
光子的流量可以大致表示为
n
p
h
c
=
c
V
n_{ph}c=\dfrac{c}{V}
nphc=Vc,其中
n
p
h
n_{ph}
nph为光子的数密度,
V
V
V为每个光子的体积
末态可以是自由态,也可以是束缚态。如果研究对象是自由电子,则初态和末态就都是自由态。自由态有两种情况,一种是存在这个电子,这时拿光子能量 h ν h\nu hν和电子的静能 m e c 2 m_ec^2 mec2相比较。对于正常的恒星而言, k T ≪ m e c 2 kT\ll m_e c^2 kT≪mec2,此时光子的能量远小于电子的静能,因此只考虑自由电子散射
在自由电子散射中,(每个电子的)散射截面是一个常数。有了电子散射截面后,就可以计算出
κ
\kappa
κ:
σ
T
=
6.65
×
1
0
−
24
c
m
−
2
κ
=
σ
T
μ
e
m
H
=
0.4
μ
e
c
m
2
g
−
1
≈
1
+
X
2
0.4
c
m
2
g
−
1
\sigma_T=6.65\times 10^{-24}cm^{-2}\\ \kappa=\frac{\sigma_T}{\mu_e m_H}=\frac{0.4}{\mu_e}cm^2g^{-1}\approx\frac{1+X}{2}0.4 cm^2 g^{-1}
σT=6.65×10−24cm−2κ=μemHσT=μe0.4cm2g−1≈21+X0.4cm2g−1
在没有原子核的情况下,自由电子不能吸收光子(因为自由电子吸收光子后不可能满足能量动量守恒)
自由-自由吸收中的发射系数:
ϵ
f
f
(
ν
)
=
4
π
j
f
f
=
2
5
π
e
6
3
m
e
c
3
(
2
π
3
m
e
k
T
)
1
/
2
g
f
f
n
e
n
I
Z
2
e
−
h
ν
k
T
\epsilon_{ff}(\nu)=4\pi j_{ff}=\frac{2^5\pi e^6}{3m_ec^3}\left(\frac{2\pi}{3m_ekT}\right)^{1/2}g_{ff}n_en_I Z^2 e^{-\frac{h\nu}{kT}}
ϵff(ν)=4πjff=3mec325πe6(3mekT2π)1/2gffnenIZ2e−kThν
发射和吸收平衡时,根据Kirchoff定律,有
S
ν
=
j
ν
a
ν
=
B
ν
(
T
)
S_\nu=\dfrac{j_\nu}{a_\nu}=B_\nu(T)
Sν=aνjν=Bν(T)。其中,
j
ν
j_\nu
jν为体发射系数,从而得到吸收系数
a
ν
a_\nu
aν:
a
f
f
(
ν
)
=
ρ
κ
ν
=
4
e
6
3
m
e
h
c
(
2
π
3
m
e
k
T
)
1
/
2
g
f
f
n
e
n
I
Z
2
ν
−
3
{
1
−
e
−
h
ν
k
T
}
a_{ff}(\nu)=\rho \kappa_\nu=\frac{4e^6}{3m_e hc}\left(\frac{2\pi}{3m_e kT}\right)^{1/2}g_{ff}n_en_I Z^2 \nu^{-3}\left\{1-e^{-\frac{h\nu}{kT}}\right\}
aff(ν)=ρκν=3mehc4e6(3mekT2π)1/2gffnenIZ2ν−3{1−e−kThν}
电子和离子的数密度分别为
n
e
=
ρ
μ
e
m
H
n_e=\dfrac{\rho}{\mu_e m_H}
ne=μemHρ和
n
I
=
ρ
μ
I
m
H
n_I=\dfrac{\rho}{\mu_I m_H}
nI=μImHρ,由此得到:
κ
ν
=
C
(
Z
2
μ
e
μ
I
)
g
f
f
ρ
T
−
1
/
2
ν
−
3
{
1
−
e
−
h
ν
k
T
}
C
=
4
e
6
3
m
e
h
c
N
A
2
(
2
π
3
m
e
k
)
1
/
2
=
1.34
×
1
0
56
(
c
g
s
)
\kappa_\nu=C\left(\frac{Z^2}{\mu_e\mu_I}\right)g_{ff}\rho T^{-1/2}\nu^{-3}\left\{1-e^{-\frac{h\nu}{kT}}\right\}\\ C=\frac{4e^6}{3m_e hc}N_A^2\left(\frac{2\pi}{3m_e k}\right)^{1/2}=1.34\times 10^{56} (cgs)
κν=C(μeμIZ2)gffρT−1/2ν−3{1−e−kThν}C=3mehc4e6NA2(3mek2π)1/2=1.34×1056(cgs)
最后得到Rosseland平均:
κ
f
f
∼
1
0
23
Z
2
μ
e
μ
I
ρ
T
−
7
/
2
c
m
2
g
−
1
κ
f
f
≈
4
×
1
0
22
(
X
+
Y
)
(
1
+
X
)
ρ
T
−
3.5
c
m
2
g
−
1
\kappa_{ff}\sim 10^{23}\frac{Z^2}{\mu_e\mu_I}\rho T^{-7/2}\ cm^2 g^{-1}\\ \kappa_{ff}\approx 4\times 10^{22}(X+Y)(1+X)\rho T^{-3.5}\ cm^2 g^{-1}
κff∼1023μeμIZ2ρT−7/2 cm2g−1κff≈4×1022(X+Y)(1+X)ρT−3.5 cm2g−1
(怎么吸收?怎么发射?)
如果密度高,则还需要考虑电子对原子核的屏蔽效应
束缚-自由跃迁
电子开始时处于束缚态,结束时处于自由态。从束缚态到自由态的过程,光电吸收
对于类氢原子,态
(
n
,
l
)
(n,l)
(n,l)的吸收截面为:
σ
b
f
(
Z
,
n
,
l
)
=
2
9
π
7
m
e
e
10
Z
4
3
3
c
h
6
n
5
g
b
f
(
2
π
ν
)
3
=
64
π
4
m
e
e
10
3
3
c
h
6
Z
4
n
5
g
b
f
ν
3
\sigma_{bf}(Z,n,l)=\frac{2^9\pi^7 m_e e^{10}Z^4}{3\sqrt{3} ch^6 n^5}\frac{g_{bf}}{(2\pi\nu)^3}=\frac{64\pi^4 m_e e^{10}}{3\sqrt{3}ch^6}\frac{Z^4}{n^5}\frac{g_{bf}}{\nu^3}
σbf(Z,n,l)=33ch6n529π7mee10Z4(2πν)3gbf=33ch664π4mee10n5Z4ν3gbf
从能量守恒的角度来看,只有光子
h
ν
>
χ
n
=
2
π
2
m
e
e
4
h
2
Z
2
n
2
h\nu>\chi_n=\dfrac{2\pi^2 m_e e^4}{h^2}\dfrac{Z^2}{n^2}
hν>χn=h22π2mee4n2Z2时,才会发生相互作用
类氢原子散射截面的波尔近似:
σ
b
f
(
ν
)
=
8
π
3
r
0
2
1
a
3
32
Z
2
(
χ
K
h
ν
)
7
/
2
(
E
h
ν
)
3
/
2
\sigma_{bf}(\nu)=\frac{8\pi}{3}r_0^2\frac{1}{a^3}\frac{32}{Z^2}\left(\frac{\chi_K}{h\nu}\right)^{7/2}\left(\frac{E}{h\nu}\right)^{3/2}
σbf(ν)=38πr02a31Z232(hνχK)7/2(hνE)3/2
其中:
r
0
=
e
2
m
e
c
2
=
α
2
a
0
,
α
=
e
2
ℏ
c
,
E
=
h
ν
−
χ
K
r_0=\frac{e^2}{m_e c^2}=\alpha^2 a_0,\ \alpha=\frac{e^2}{\hbar c},\ E=h\nu-\chi_K
r0=mec2e2=α2a0, α=ℏce2, E=hν−χK
电离能为:
χ
K
=
Z
e
2
2
a
=
Z
2
e
2
2
K
2
a
0
\chi_K=\frac{Ze^2}{2a}=\frac{Z^2 e^2}{2K^2 a_0}
χK=2aZe2=2K2a0Z2e2
对于
i
i
i-电离态而言:
κ
ν
=
∑
Z
∑
n
σ
b
f
(
Z
,
n
,
l
)
(
n
i
,
n
n
i
)
Z
(
n
i
∑
i
n
i
)
Z
ρ
A
Z
m
H
\kappa_\nu=\sum_Z\sum_n\sigma_{bf}(Z,n,l)\left(\frac{n_{i,n}}{n_i}\right)_Z\left(\frac{n_i}{\sum_i n_i}\right)_Z\frac{\rho}{A_Z m_H}
κν=Z∑n∑σbf(Z,n,l)(nini,n)Z(∑inini)ZAZmHρ
对于高电离物质,只考虑类氢原子和完全电离物质这两种电离状态,利用Saha定理和玻尔兹曼分布,有:
n
i
,
n
n
i
=
g
n
Z
e
−
ϵ
i
,
n
/
k
T
n
i
+
1
n
e
n
i
=
2
Z
i
+
1
(
T
)
Z
i
(
T
)
(
2
π
m
e
k
T
h
2
)
3
/
2
e
−
χ
i
,
1
/
K
T
\frac{n_{i,n}}{n_i}=\frac{g_n}{Z}e^{-\epsilon_{i,n}/kT}\\ \frac{n_{i+1}n_e}{n_i}=\frac{2Z_{i+1}(T)}{Z_i(T)}\left(\frac{2\pi m_ekT}{h^2}\right)^{3/2}e^{-\chi_{i,1}/KT}
nini,n=Zgne−ϵi,n/kTnini+1ne=Zi(T)2Zi+1(T)(h22πmekT)3/2e−χi,1/KT
最后得到吸收系数为:
κ
b
f
=
a
ν
ρ
≈
∑
z
,
n
C
(
χ
i
,
n
k
T
e
χ
i
,
n
/
k
T
)
(
Z
2
χ
Z
μ
e
A
Z
)
ρ
T
−
1
/
2
ν
−
3
\kappa_{bf}=\frac{a_\nu}{\rho}\approx \sum_{z,n}C\left(\frac{\chi_{i,n}}{kT}e^{\chi_{i,n}/kT}\right)\left(\frac{Z^2\chi_Z}{\mu_e A_Z}\right)\rho T^{-1/2}\nu^{-3}
κbf=ρaν≈z,n∑C(kTχi,neχi,n/kT)(μeAZZ2χZ)ρT−1/2ν−3
当
k
T
≪
χ
i
,
1
kT\ll \chi_{i,1}
kT≪χi,1时,
κ
ˉ
b
f
∝
ρ
T
−
7
/
2
\bar\kappa_{bf}\propto \rho T^{-7/2}
κˉbf∝ρT−7/2
当 k T ≫ χ i , 1 kT\gg \chi_{i,1} kT≫χi,1时, κ ˉ b f ∝ ρ T − 9 / 2 \bar\kappa_{bf}\propto \rho T^{-9/2} κˉbf∝ρT−9/2
低温时,在恒星内部,数值比自由-自由吸收大三个量级
κ
ˉ
b
f
≈
4
×
1
0
25
Z
(
1
+
X
)
ρ
T
−
7
/
2
c
m
2
g
−
1
\bar \kappa_{bf}\approx 4\times 10^{25} Z(1+X)\rho T^{-7/2}\ cm^2g^{-1}
κˉbf≈4×1025Z(1+X)ρT−7/2 cm2g−1
束缚-束缚跃迁
在高温( T > 1 0 6 K T> 10^6 K T>106K)时,几乎都是电离物质,此时束缚-束缚跃迁不重要。但在低温情况下,束缚-束缚跃迁重要。在恒星的内部,压强比较大,原子之间距离更近,缩短了寿命,导致谱线变得更宽,这个现象名为压力致宽
在束缚-束缚跃迁中,共有三种过程:
(1)吸收: h ν + i → j h\nu+i\rightarrow j hν+i→j,速率为 B i j J ˉ ν B_{ij}\bar J_\nu BijJˉν
(2)自发跃迁: j → h ν + i j\rightarrow h\nu +i j→hν+i,速率为 A j i A_{ji} Aji
(3)受激跃迁: j + h ν → 2 h ν + i j+h\nu \rightarrow 2h\nu +i j+hν→2hν+i,速率为 B j i J ˉ ν B_{ji}\bar J_\nu BjiJˉν
其中, J ˉ ν = ∫ 0 ∞ J ν ϕ ( ν ) d ν \bar J_\nu=\int^\infty_0 J_\nu \phi(\nu)d\nu Jˉν=∫0∞Jνϕ(ν)dν, J ν ≡ 1 4 π ∫ I ν d Ω J_\nu\equiv \dfrac{1}{4\pi}\int I_\nu d\Omega Jν≡4π1∫IνdΩ, ∫ 0 ∞ ϕ ( ν ) d ν = 1 \int^\infty_0\phi(\nu)d\nu=1 ∫0∞ϕ(ν)dν=1
线致宽的来源:自发辐射,热运动多普勒致宽和碰撞/压强致宽
爱因斯坦系数间关系
对于完全热平衡的气体,在低密度下的能级分布不变:
n
i
B
i
j
J
ˉ
ν
=
n
j
A
j
i
+
n
j
B
j
i
J
ˉ
ν
n_i B_{ij}\bar J_\nu=n_j A_{ji}+n_j B_{ji}\bar J_\nu
niBijJˉν=njAji+njBjiJˉν
即:
J
ˉ
ν
=
A
j
i
/
B
j
i
(
n
i
/
n
j
)
(
B
i
j
/
B
j
i
)
−
1
\bar J_\nu = \frac{A_{ji}/B_{ji}}{(n_i/n_j)(B_{ij}/B_{ji})-1}
Jˉν=(ni/nj)(Bij/Bji)−1Aji/Bji
在热平衡下:
J
ˉ
ν
≈
B
ν
=
2
h
ν
3
/
c
2
e
−
Δ
E
/
k
T
−
1
n
j
=
n
i
g
j
g
i
e
−
Δ
E
k
T
\bar J_\nu\approx B_\nu=\frac{2h\nu^3/c^2}{e^{-\Delta E/kT}-1}\\ n_j=n_i\frac{g_j}{g_i}e^{-\frac{\Delta E}{kT}}
Jˉν≈Bν=e−ΔE/kT−12hν3/c2nj=nigigje−kTΔE
由此得到爱因斯坦系数之间的关系:
g
i
B
i
j
=
g
j
B
j
i
A
j
i
B
j
i
=
2
h
ν
3
c
2
g_i B_{ij}=g_j B_{ji}\\ \frac{A_{ji}}{B_{ji}}=\frac{2h\nu^3}{c^2}
giBij=gjBjiBjiAji=c22hν3
这一关系无论是否平衡态都成立,因为吸收和发射系数都是物质性质
爱因斯坦系数与吸收系数的关系
在体积元
Δ
V
\Delta V
ΔV中,单位时间内,原子从
i
i
i跃迁
j
j
j这个过程中吸收的辐射能量为:
(
n
1
B
12
−
n
2
B
21
)
I
ν
d
Ω
4
π
ϕ
(
ν
)
d
ν
h
ν
d
t
d
V
=
I
ν
α
ν
d
s
d
t
d
Ω
d
A
d
ν
(n_1B_{12}-n_2B_{21})\frac{I_\nu d\Omega}{4\pi}\phi(\nu)d\nu h\nu dtdV=I_\nu \alpha_\nu dsdtd\Omega dAd\nu
(n1B12−n2B21)4πIνdΩϕ(ν)dνhνdtdV=IνανdsdtdΩdAdν
得到:
α
ν
=
(
n
1
B
12
−
n
2
B
21
)
ϕ
(
ν
)
h
ν
4
π
\alpha_\nu=\frac{(n_1B_{12}-n_2B_{21})\phi(\nu)h\nu}{4\pi}
αν=4π(n1B12−n2B21)ϕ(ν)hν
吸收线的轮廓
吸收线的轮廓 ϕ ( ν ) \phi(\nu) ϕ(ν)与具体的物理过程有关:
(a)自然阻尼情形,能级的寿命为 τ = 1 A j i \tau=\dfrac{1}{A_{ji}} τ=Aji1, Δ E τ ≥ h 2 π \Delta E\tau\ge\dfrac{h}{2\pi} ΔEτ≥2πh, Δ ν 1 / 2 = A j i 2 π \Delta\nu_{1/2}=\dfrac{A_{ji}}{2\pi} Δν1/2=2πAji
此时这个轮廓为洛伦兹曲线:
ϕ
(
ν
)
=
1
π
Γ
/
4
π
(
ν
−
ν
0
)
2
+
(
Γ
/
4
π
)
2
Γ
j
=
∑
i
A
j
i
\phi(\nu)=\frac{1}{\pi}\frac{\Gamma/4\pi}{(\nu-\nu_0)^2+(\Gamma/4\pi)^2}\\ \Gamma_j=\sum_i A_{ji}
ϕ(ν)=π1(ν−ν0)2+(Γ/4π)2Γ/4πΓj=i∑Aji
(b)热运动致宽情形:
ϕ
(
ν
)
=
1
π
Δ
ν
D
e
−
(
ν
−
ν
0
Δ
ν
D
)
2
Δ
ν
D
=
ν
0
(
v
T
/
c
)
,
v
T
=
(
2
k
T
m
)
1
/
2
\phi(\nu)=\frac{1}{\sqrt{\pi}\Delta \nu_D}e^{-\left(\frac{\nu-\nu_0}{\Delta\nu_D}\right)^2}\\ \Delta\nu_D=\nu_0(v_T/c),\ v_T=\left(\frac{2kT}{m}\right)^{1/2}
ϕ(ν)=πΔνD1e−(ΔνDν−ν0)2ΔνD=ν0(vT/c), vT=(m2kT)1/2
(c)碰撞致宽情形:
斯塔克效应扰动,取决于具体粒子以及运动速度(作用时间与原子能级寿命的比)。 Δ E = C r − n ( t ) \Delta E=C r^{-n}(t) ΔE=Cr−n(t),不同类型斯塔克效应的 n n n值不同
I. 使用impact近似时,作用时间远短于原子的能级寿命:
W
(
t
)
=
e
−
t
/
τ
d
t
/
τ
W(t)=e^{-t/\tau}dt/\tau
W(t)=e−t/τdt/τ
此时进行傅里叶变换可以得到洛伦兹轮廓
II. 慢运动近似,分布取决于距离
r
r
r和相互场量
F
F
F的分布
W
(
F
)
\mathcal W(F)
W(F):
ϕ
(
ν
)
d
ν
=
W
(
F
)
∂
F
∂
ν
d
ν
\phi(\nu)d\nu=\mathcal W(F)\frac{\partial F}{\partial \nu}d\nu
ϕ(ν)dν=W(F)∂ν∂Fdν
当多种致宽机制存在时,情况更为复杂,这时需要使用卷积的方式进行计算:
ϕ
(
Δ
ν
)
=
∫
−
∞
∞
ϕ
1
(
ν
−
ν
1
)
ϕ
2
(
ν
−
ν
0
)
d
(
ν
1
−
ν
0
)
\phi(\Delta \nu)=\int^\infty_{-\infty}\phi_1(\nu-\nu_1)\phi_2(\nu-\nu_0)d(\nu_1-\nu_0)
ϕ(Δν)=∫−∞∞ϕ1(ν−ν1)ϕ2(ν−ν0)d(ν1−ν0)
例如,洛伦兹和热运动致宽机制同时存在时,进行卷积得到:
ϕ
(
Δ
v
)
=
1
π
Δ
ν
D
H
(
a
,
Δ
ν
Δ
ν
D
)
\phi(\Delta v)=\frac{1}{\sqrt{\pi}\Delta \nu_D}H\left(a,\frac{\Delta\nu}{\Delta\nu_D}\right)
ϕ(Δv)=πΔνD1H(a,ΔνDΔν)
其中,
a
=
Γ
4
π
Δ
ν
D
a=\dfrac{\Gamma}{4\pi\Delta \nu_D}
a=4πΔνDΓ,
Δ
ν
=
ν
−
ν
0
\Delta\nu=\nu-\nu_0
Δν=ν−ν0,
y
=
Δ
ν
1
Δ
ν
D
y=\dfrac{\Delta\nu_1}{\Delta \nu_D}
y=ΔνDΔν1
H − ^- −吸收
H
−
H^-
H−:一个质子外面有两个电子。正常的氢原子,加上一个电子
H
+
e
−
⇄
H
−
+
γ
H+e^-\rightleftarrows H^- + \gamma
H+e−⇄H−+γ
韧致辐射,碰撞产生光子——自由-自由吸收:
H
−
+
e
−
(
v
)
+
h
ν
⇄
H
−
+
e
−
(
v
)
H^-+e^-(v)+h\nu\rightleftarrows H^-+e^-(v)
H−+e−(v)+hν⇄H−+e−(v)
光致电离,相当于将
H
−
H^-
H−离解——束缚-自由跃迁:
H
−
+
h
ν
⇄
H
+
e
−
(
v
)
H^-+h\nu\rightleftarrows H+e^-(v)
H−+hν⇄H+e−(v)
利用Saha方程,计算
H
−
H^-
H−离子的数量:
n
H
−
n
H
0
=
n
e
2
Z
H
−
Z
H
0
(
h
2
2
π
m
e
k
T
)
3
2
e
χ
k
T
\frac{n_{H^-}}{n_{H_0}}=\frac{n_e}{2}\frac{Z_{H^-}}{Z_{H_0}}\left(\frac{h^2}{2\pi m_e kT}\right)^\frac{3}{2}e^\frac{\chi}{kT}
nH0nH−=2neZH0ZH−(2πmekTh2)23ekTχ
其中,
Z
H
−
=
1
Z_{H^-}=1
ZH−=1,
Z
H
0
=
2
Z_{H_0}=2
ZH0=2,
χ
=
0.754
e
V
\chi=0.754eV
χ=0.754eV
低温情况下,主要的 H − H^- H−来自于金属。而如果缺乏金属的话(贫金属星),则主要来自于 H H H的电离。因为 H H H的电离势高,所以在低温时,电子密度随温度上升很快
在纯氢气体中,利用含有电离度
x
x
x的方程进行电离度的计算:
x
2
1
−
x
2
=
1
n
(
2
π
m
e
k
T
h
2
)
3
2
e
−
χ
I
k
T
\frac{x^2}{1-x^2}=\frac{1}{n}\left(\frac{2\pi m_e kT}{h^2}\right)^\frac{3}{2}e^{-\frac{\chi_I}{kT}}
1−x2x2=n1(h22πmekT)23e−kTχI
对于太阳大气而言,
T
∼
6600
K
T\sim 6600K
T∼6600K,
n
∼
1
0
18
c
m
−
3
n\sim 10^{18}cm^{-3}
n∼1018cm−3,电离度的估计值为
n
H
−
n
H
0
∼
1.3
×
1
0
−
7
\dfrac{n_{H^-}}{n_{H_0}}\sim 1.3\times 10^{-7}
nH0nH−∼1.3×10−7
在6600 K K K的大气中,主要辐射波长大于3650KaTeX parse error: Can't use function '\r' in math mode at position 1: \̲r̲ ̲A,氢原子主要吸收 n = 3 n=3 n=3的光电吸收,数目占比 n 3 n H 0 ∼ 2.5 × 1 0 − 9 \dfrac{n_3}{n_{H_0}}\sim 2.5\times 10^{-9} nH0n3∼2.5×10−9
在富金属物质中,因为金属价电子结合能低,所以容易电离,从而使得电子主要由金属提供,电子数密度对温度不敏感,这样的现象在低温恒星中更明显:
κ
H
−
≈
2.5
×
1
0
−
31
(
Z
0.02
)
ρ
1
2
T
9
c
m
2
g
−
1
\kappa_{H^-}\approx 2.5\times 10^{-31}\left(\frac{Z}{0.02}\right)\rho^\frac{1}{2}T^9\ cm^2g^{-1}
κH−≈2.5×10−31(0.02Z)ρ21T9 cm2g−1
在此情形下,需要对每一种重要的金属都进行考虑:
n
j
=
ρ
μ
I
m
u
Y
j
Y
i
+
Y
e
Y
i
=
Z
i
(
ρ
,
T
)
=
2
g
i
+
g
i
μ
I
m
u
ρ
(
2
π
m
e
k
T
h
2
)
3
2
e
−
χ
i
k
T
n_j=\frac{\rho}{\mu_I m_u}Y_j\\ \frac{Y_i+Y_e}{Y_i}=Z_i(\rho,T)=\frac{2g_{i^+}}{g_i}\frac{\mu_I m_u}{\rho}\left(\frac{2\pi m_e kT}{h^2}\right)^\frac{3}{2}e^{-\frac{\chi_i}{kT}}
nj=μImuρYjYiYi+Ye=Zi(ρ,T)=gi2gi+ρμImu(h22πmekT)23e−kTχi
其中,
Y
j
=
X
j
A
j
Y_j=\dfrac{X_j}{A_j}
Yj=AjXj为粒子数丰度
粒子数: Y i + + Y i = Y i , 0 Y_{i^+}+Y_i=Y_{i,0} Yi++Yi=Yi,0
电荷数: Y e + Y H − = ∑ Y i + + Y H + Y_e+Y_{H^-}=\sum Y_{i^+}+Y_{H^+} Ye+YH−=∑Yi++YH+
辐射极限:Eddington光度
对于由辐射压主导的恒星:
P
r
a
d
=
1
3
a
T
4
<
P
P_{rad}=\frac{1}{3}aT^4<P
Prad=31aT4<P
其中,
P
P
P为总的压强,辐射压总是会比总压强要小
在恒星大气中,能量输运以辐射传导为主:
F
=
L
4
π
R
2
=
−
4
a
c
3
κ
ρ
T
3
d
T
d
r
F=\frac{L}{4\pi R^2}=-\frac{4ac}{3\kappa\rho}T^3\frac{dT}{dr}
F=4πR2L=−3κρ4acT3drdT
上式可以改写为:
L
=
−
16
π
a
c
3
κ
ρ
R
2
T
3
d
T
d
r
=
−
4
π
c
κ
ρ
R
2
d
P
r
a
d
d
r
<
−
4
π
c
κ
ρ
R
2
d
P
d
r
\begin{aligned} L&=-\frac{16\pi ac}{3\kappa\rho}R^2T^3\frac{dT}{dr}\\ &=-\frac{4\pi c}{\kappa\rho}R^2\frac{dP_{rad}}{dr}\\ &< -\frac{4\pi c}{\kappa\rho}R^2\frac{dP}{dr} \end{aligned}
L=−3κρ16πacR2T3drdT=−κρ4πcR2drdPrad<−κρ4πcR2drdP
其中,用到了关系式
d
P
r
a
d
d
T
=
4
3
a
T
3
\dfrac{dP_{rad}}{dT}=\dfrac{4}{3}aT^3
dTdPrad=34aT3
由流体静力平衡方程:
d
P
d
r
=
−
G
M
R
2
ρ
\frac{dP}{dr}=-\frac{GM}{R^2}\rho
drdP=−R2GMρ
有:
L
<
L
E
d
d
=
4
π
c
G
M
κ
L<L_{Edd}=\frac{4\pi c GM}{\kappa}
L<LEdd=κ4πcGM
此即Eddington极限光度,这是辐射压与引力平衡时的光度,如果超出了这个极限,则不能达到平衡
通过恒星的质量与光度的关系,再与Eddington极限相比较,即可得到最大稳定质量
对于主序星而言,
L
∝
M
3.5
L\propto M^{3.5}
L∝M3.5,而主序星是以电子散射为主,因此
κ
T
≈
0.2
(
1
+
X
)
c
m
2
g
−
1
\kappa_T\approx 0.2(1+X)\ cm^2g^{-1}
κT≈0.2(1+X) cm2g−1,代入到Eddington极限光度表达式中,得到:
(
M
M
⊙
)
2.5
≤
4
π
G
c
0.2
(
1
+
X
)
M
⊙
L
⊙
\left(\frac{M}{M_\odot}\right)^{2.5}\le \frac{4\pi Gc}{0.2(1+X)}\frac{M_\odot}{L_\odot}
(M⊙M)2.5≤0.2(1+X)4πGcL⊙M⊙
这样的质量算出来是
60
M
⊙
60M_\odot
60M⊙,当质量大于
60
M
⊙
60M_\odot
60M⊙时,一定会有很强的风。星风在大质量恒星演化过程和恒星演化的晚期起着重要的作用
对于AGB阶段的恒星,温度比较低,不透明度大,因此其Eddington光度远低于上述值
热不稳定性与对流
当物质沿引力方向温度梯度较高时,可能会发展对流不稳定性。对流在恒星中起到能量输运和物质混合的作用,在转动时会有其他动力学效应
化学均匀情形
假定在 r r r处有一个密度为 ρ \rho ρ,压强为 P P P,温度为 T T T的泡,在密度稍低于周围的密度 ρ ∗ \rho^\ast ρ∗时,受浮力作用上升到 r + d r r+dr r+dr处。如果在 r + d r r+dr r+dr处密度反差比在 r r r处更大,则浮力会变大,导致继续上升,此时对流不稳定。由此,对流稳定的条件为: ρ ( r + d r ) > ρ ∗ ( r + d r ) \rho(r+dr)>\rho^\ast(r+dr) ρ(r+dr)>ρ∗(r+dr)
为了简化讨论,作出下述三个假设:
(1)物质元在浮力作用下上升很快,近似于是绝热过程:
ρ
∗
(
r
+
d
r
)
=
ρ
(
r
)
+
d
ρ
∗
=
ρ
(
r
)
+
(
∂
ρ
∂
P
)
T
d
P
∗
+
(
∂
ρ
∂
T
)
P
d
T
∗
+
(
∂
ρ
∂
μ
)
P
d
μ
∗
\begin{aligned} \rho^\ast(r+dr)&=\rho(r)+d\rho^\ast\\ &=\rho(r)+\left(\frac{\partial \rho}{\partial P}\right)_TdP^\ast+\left(\frac{\partial \rho}{\partial T}\right)_PdT^\ast+\left(\frac{\partial \rho}{\partial \mu}\right)_Pd\mu^\ast \end{aligned}
ρ∗(r+dr)=ρ(r)+dρ∗=ρ(r)+(∂P∂ρ)TdP∗+(∂T∂ρ)PdT∗+(∂μ∂ρ)Pdμ∗
ρ ( r + d r ) = ρ ( r ) + d ρ = ρ ( r ) + ( ∂ ρ ∂ P ) T d P + ( ∂ ρ ∂ T ) P d T + ( ∂ ρ ∂ μ ) P d μ \begin{aligned} \rho(r+dr)&=\rho(r)+d\rho\\ &=\rho(r)+\left(\frac{\partial \rho}{\partial P}\right)_TdP+\left(\frac{\partial \rho}{\partial T}\right)_PdT+\left(\frac{\partial \rho}{\partial \mu}\right)_Pd\mu \end{aligned} ρ(r+dr)=ρ(r)+dρ=ρ(r)+(∂P∂ρ)TdP+(∂T∂ρ)PdT+(∂μ∂ρ)Pdμ
(2)化学均匀, d μ = d μ ∗ = 0 d\mu=d\mu^\ast=0 dμ=dμ∗=0
(3)上升的过程中和外界近似压力平衡(速度远小于声速 c s c_s cs): P ≈ P ∗ P\approx P^\ast P≈P∗
因为
(
∂
ρ
∂
T
)
P
<
0
\left(\dfrac{\partial \rho}{\partial T}\right)_P<0
(∂T∂ρ)P<0,所以稳定条件为:
(
d
T
d
r
)
∗
>
(
d
T
d
r
)
a
d
或
∣
d
T
d
r
∣
∗
<
∣
d
T
d
r
∣
a
d
\left(\frac{dT}{dr}\right)_\ast>\left(\frac{dT}{dr}\right)_{ad}\ 或\ \left|\frac{dT}{dr}\right|_\ast<\left|\frac{dT}{dr}\right|_{ad}
(drdT)∗>(drdT)ad 或 ∣∣∣∣drdT∣∣∣∣∗<∣∣∣∣drdT∣∣∣∣ad
上式为Schwarzschild判据,在温度梯度高于绝热温度梯度时则会导致对流不稳定
Schwarzschild判据可以写成其他的形式,例如
ρ
=
ρ
(
s
,
P
)
\rho=\rho(s,P)
ρ=ρ(s,P):
d
ρ
=
(
∂
ρ
∂
s
)
P
d
s
+
(
∂
ρ
∂
P
)
s
d
P
d
ρ
∗
=
(
∂
ρ
∂
s
)
P
d
s
∗
+
(
∂
ρ
∂
P
)
s
d
P
∗
\begin{aligned} d\rho&=\left(\frac{\partial \rho}{\partial s}\right)_P ds + \left(\frac{\partial \rho}{\partial P}\right)_s dP\\ d\rho_\ast &= \left(\frac{\partial\rho}{\partial s}\right)_P ds_\ast + \left(\frac{\partial\rho}{\partial P}\right)_s dP_\ast \end{aligned}
dρdρ∗=(∂s∂ρ)Pds+(∂P∂ρ)sdP=(∂s∂ρ)Pds∗+(∂P∂ρ)sdP∗
根据
(
∂
ρ
∂
s
)
P
<
0
\left(\dfrac{\partial\rho}{\partial s}\right)_P<0
(∂s∂ρ)P<0和
d
s
=
0
ds=0
ds=0,得到稳定条件:
(
d
s
d
r
)
∗
>
0
\left(\frac{ds}{dr}\right)_\ast > 0
(drds)∗>0
由此得到结论:低熵物质在高熵物质上不稳定,引力方向和熵梯度方向相反时稳定
由绝热方程:
d
P
P
+
Γ
2
1
−
Γ
2
d
T
T
=
0
\frac{dP}{P}+\frac{\Gamma_2}{1-\Gamma_2}\frac{dT}{T}=0
PdP+1−Γ2Γ2TdT=0
得到稳定条件(
d
P
d
r
\dfrac{dP}{dr}
drdP为负):
(
d
T
d
r
)
>
(
1
−
1
Γ
2
)
T
P
(
d
P
d
r
)
∇
≡
(
d
ln
T
d
ln
P
)
<
Γ
2
−
1
Γ
2
\begin{aligned} &\left(\frac{dT}{dr}\right)>\left(1-\frac{1}{\Gamma_2}\right)\frac{T}{P}\left(\frac{dP}{dr}\right)\\ &\nabla\equiv\left(\frac{d\ln T}{d\ln P}\right)<\frac{\Gamma_2-1}{\Gamma_2} \end{aligned}
(drdT)>(1−Γ21)PT(drdP)∇≡(dlnPdlnT)<Γ2Γ2−1
则不稳定条件为:
∇
≡
(
d
ln
T
d
ln
P
)
>
Γ
2
−
1
Γ
2
\nabla\equiv\left(\frac{d\ln T}{d\ln P}\right)>\frac{\Gamma_2-1}{\Gamma_2}
∇≡(dlnPdlnT)>Γ2Γ2−1
假设恒星以辐射输运为主,光度为
L
L
L,则有:
l
(
r
)
4
π
r
2
=
F
(
r
)
=
−
4
a
c
3
κ
R
ρ
T
3
d
T
d
r
d
P
d
r
=
−
ρ
G
m
r
2
\frac{l(r)}{4\pi r^2}=F(r)=-\frac{4ac}{3\kappa_R\rho}T^3\frac{dT}{dr}\\ \frac{dP}{dr}=-\rho\frac{Gm}{r^2}
4πr2l(r)=F(r)=−3κRρ4acT3drdTdrdP=−ρr2Gm
辐射的温度梯度为:
∇
r
a
d
≡
(
d
ln
T
d
ln
P
)
r
a
d
=
3
κ
P
l
16
π
a
c
G
m
T
4
\nabla_{rad}\equiv\left(\frac{d\ln T}{d\ln P}\right)_{rad}=\frac{3\kappa Pl}{16\pi acGmT^4}
∇rad≡(dlnPdlnT)rad=16πacGmT43κPl
对流稳定的条件为
∇
r
a
d
≤
∇
a
d
\nabla_{rad}\le\nabla_{ad}
∇rad≤∇ad,于是:
l
r
≤
16
π
a
c
G
3
κ
R
(
1
−
1
Γ
2
)
T
4
P
m
(
r
)
l_r\le \frac{16\pi acG}{3\kappa_R}\left(1-\frac{1}{\Gamma_2}\right)\frac{T^4}{P}m(r)
lr≤3κR16πacG(1−Γ21)PT4m(r)
因此,辐射运输极限
l
m
a
x
∝
1
κ
R
(
1
−
1
Γ
2
)
l_{max}\propto \dfrac{1}{\kappa_R}\left(1-\dfrac{1}{\Gamma_2}\right)
lmax∝κR1(1−Γ21)
在 κ \kappa κ很大, Γ \Gamma Γ接近1,光度很大和部分电离区,容易发生对流不稳定
化学均匀对流不稳定的条件为:
∇
r
a
d
>
∇
a
d
;
∇
r
a
d
≡
(
d
ln
T
d
ln
P
)
r
a
d
;
∇
a
d
≡
(
d
ln
T
d
ln
P
)
a
d
\nabla_{rad}>\nabla_{ad};\ \nabla_{rad}\equiv\left(\frac{d\ln T}{d\ln P}\right)_{rad};\ \nabla_{ad}\equiv\left(\frac{d\ln T}{d\ln P}\right)_{ad}
∇rad>∇ad; ∇rad≡(dlnPdlnT)rad; ∇ad≡(dlnPdlnT)ad
化学不均匀情形
恒星中由于核合成的原因,使得在中心存在化学不均匀。化学不均匀情形下,对流不稳定条件为:
∇
r
a
d
>
∇
a
d
+
φ
δ
∇
μ
\nabla_{rad}>\nabla_{ad}+\frac{\varphi}{\delta}\nabla_\mu
∇rad>∇ad+δφ∇μ
其中:
δ
=
−
(
∂
ln
ρ
∂
ln
T
)
P
,
φ
=
(
∂
ln
ρ
∂
ln
μ
)
P
,
∇
μ
=
(
d
ln
μ
d
ln
P
)
T
\delta=-\left(\frac{\partial \ln \rho}{\partial\ln T}\right)_P,\ \varphi=\left(\frac{\partial \ln \rho}{\partial \ln \mu}\right)_P,\ \nabla_\mu =\left(\frac{d\ln\mu}{d\ln P}\right)_T
δ=−(∂lnT∂lnρ)P, φ=(∂lnμ∂lnρ)P, ∇μ=(dlnPdlnμ)T
恒星核合成的化学分层,导致
∇
μ
>
0
,
φ
δ
∇
μ
>
0
\nabla_\mu >0,\ \dfrac{\varphi}{\delta}\nabla_\mu>0
∇μ>0, δφ∇μ>0,因此这样的化学分层起到了稳定的作用
在化学不均匀的情况下,存在半对流区:
∇
a
d
<
∇
r
a
d
<
∇
a
d
+
φ
δ
∇
μ
\nabla_{ad}<\nabla_{rad}<\nabla_{ad}+\frac{\varphi}{\delta}\nabla_\mu
∇ad<∇rad<∇ad+δφ∇μ
在半对流区,对流元在平衡位置附近震荡,向外运动时由于高于背景温度而放热,向内运动时则会吸热,热能转换为机械能——震动不稳定性