原子核结构基本图像
原子由核和电子构成,其中电子轨道的尺度为 1 0 − 8 c m 10^{-8}cm 10−8cm,核的尺度为 1 0 − 13 c m 10^{-13}cm 10−13cm。原子核由质子和中子组成,尺度为 1 0 − 5 10^{-5} 10−5原子半径。质子和中子的相互作用是强作用力,通过交换质量为 m m m的 π \pi π介子产生(按现有的理论,所有的相互作用力都是通过交换粒子产生的)。质子和中子这两个核子都由夸克组成,夸克的电荷是分数电荷
介子有三类:
π
+
\pi^+
π+,
π
0
\pi^0
π0和
π
−
\pi^-
π−,这些介子是会发生衰变的(介子由两个夸克组成):
u
d
ˉ
:
π
+
→
μ
+
ν
μ
(
u
u
ˉ
−
d
d
ˉ
)
/
2
:
π
0
→
2
γ
d
u
ˉ
:
π
−
→
μ
ν
ˉ
μ
u\bar d:\ \pi^+ \rightarrow \mu^+ \nu_\mu\\ (u\bar u-d\bar d)/\sqrt{2}:\ \pi^0\rightarrow 2\gamma\\ d\bar u:\ \pi^-\rightarrow \mu\bar \nu_\mu
udˉ: π+→μ+νμ(uuˉ−ddˉ)/2: π0→2γduˉ: π−→μνˉμ
两个核子的作用由Yukawa势描述:
V
12
(
r
)
=
−
g
2
exp
[
−
(
m
c
h
)
r
]
r
V_{12}(r)=-g^2\frac{\exp\left[-\left(\frac{mc}{h}\right)r\right]}{r}
V12(r)=−g2rexp[−(hmc)r]
其中,
m
c
2
≈
164
M
e
V
mc^2\approx 164 MeV
mc2≈164MeV,耦合常数在强相互作用下,
g
2
/
ℏ
c
≈
10
g^2/\hbar c\approx 10
g2/ℏc≈10。与
e
2
/
ℏ
c
=
1
/
137
e^2/\hbar c=1/137
e2/ℏc=1/137相对比,核作用比电磁作用强
1
0
3
10^3
103倍
在距离非常近的时候,核子之间的作用非常强(呈现出很深的势阱),但距离在超过1 f m fm fm后,强相互作用迅速减弱
核子也是费米子,在势场中按能级分布。质子和质子之间有库仑作用,核子之间有弱相互作用。质子和中子看做是核子的不同态(同位旋不同),质子和中子之间可以通过弱相互作用转换
经验的核质量
将质子和中子散开后所差的质量就是核的质量,经验上,核质量为:
M
(
A
,
Z
)
=
(
A
−
Z
)
m
n
+
Z
m
p
−
a
V
A
+
a
s
A
2
3
+
a
c
Z
2
A
−
1
3
+
a
a
(
A
−
2
Z
)
2
4
A
+
δ
A
1
/
2
M(A,Z)=(A-Z)m_n+Zm_p-a_V A+a_s A^\frac{2}{3}+a_c Z^2A^{-\frac{1}{3}}+a_a\frac{(A-2Z)^2}{4A}+\frac{\delta}{A^{1/2}}
M(A,Z)=(A−Z)mn+Zmp−aVA+asA32+acZ2A−31+aa4A(A−2Z)2+A1/2δ
其中,前两项为组份的贡献,第三项为短程力作用
∝
\propto
∝核子数,第四项为表面项,第五项为库仑作用项,第六项为不对称性能量项,最后一项为成对项,其中:
δ
=
{
+
δ
p
,
o
d
d
−
o
d
d
0
,
o
d
d
−
e
v
e
n
−
δ
p
,
e
v
e
n
−
e
v
e
n
\delta=\left\{ \begin{aligned} &+\delta_p,\ odd-odd\\ &0,\ odd-even\\ &-\delta_p,\ even-even \end{aligned} \right.
δ=⎩⎪⎨⎪⎧+δp, odd−odd0, odd−even−δp, even−even
能级结构求解
方势阱
对于方势阱:
V
(
r
)
=
{
−
V
0
,
r
<
r
n
u
c
l
0
,
r
>
r
n
u
c
l
V(r)=\left\{ \begin{aligned} &-V_0,\ r<r_{nucl}\\ & 0,\ r> r_{nucl} \end{aligned} \right.
V(r)={−V0, r<rnucl0, r>rnucl
可以求出能量本征值的解析解:
E
n
,
l
=
−
∣
V
n
u
c
∣
+
h
2
2
M
R
2
[
π
2
(
n
+
l
2
)
2
−
l
(
l
+
1
)
]
E_{n,l}=-|V_{nuc}|+\frac{h^2}{2MR^2}\left[\pi^2\left(n+\frac{l}{2}\right)^2-l(l+1)\right]
En,l=−∣Vnuc∣+2MR2h2[π2(n+2l)2−l(l+1)]
Woods-Saxon势能
Woods-Saxon势能为:
V
(
r
)
=
−
V
0
1
+
exp
(
−
r
−
R
a
)
V(r)=-\frac{V_0}{1+\exp\left(-\frac{r-R}{a}\right)}
V(r)=−1+exp(−ar−R)V0
其中,
R
=
r
0
A
1
/
3
R=r_0 A^{1/3}
R=r0A1/3,
r
0
=
1.25
f
m
r_0=1.25\ fm
r0=1.25 fm,
A
A
A为原子质量数。典型情况下,
V
0
≈
50
M
e
V
,
a
=
0.5
f
m
V_0\approx 50\ MeV,\ a=0.5\ fm
V0≈50 MeV, a=0.5 fm
求出能级:
E
=
−
δ
2
R
0
2
h
2
2
m
[
β
2
γ
−
(
n
+
1
2
)
]
2
+
V
0
+
h
2
2
m
D
0
l
(
l
+
1
)
E=-\frac{\delta^2}{R_0^2}\frac{h^2}{2m}\left[\frac{\beta}{2\sqrt{\gamma}}-\left(n+\frac{1}{2}\right)\right]^2+V_0+\frac{h^2}{2m}D_0 l(l+1)
E=−R02δ22mh2[2γβ−(n+21)]2+V0+2mh2D0l(l+1)
其中:
δ
=
R
0
a
,
β
=
β
0
δ
2
,
γ
=
γ
0
δ
2
β
0
=
−
2
m
R
0
2
ℏ
2
[
−
V
0
+
h
2
2
m
D
1
l
(
l
+
1
)
]
γ
0
=
2
m
R
0
2
ℏ
2
[
V
0
+
h
2
2
m
D
2
l
(
l
+
1
)
]
D
0
=
1
+
3
δ
+
3
δ
2
D
1
=
−
4
δ
−
6
δ
2
D
2
=
1
δ
+
3
δ
2
\begin{aligned} &\delta=\frac{R_0}{a},\ \beta=\frac{\beta_0}{\delta^2},\ \gamma=\frac{\gamma_0}{\delta^2}\\ &\beta_0=-\frac{2mR_0^2}{\hbar^2}\left[-V_0+\frac{h^2}{2m}D_1l(l+1)\right]\\ & \gamma_0=\frac{2mR_0^2}{\hbar^2}\left[V_0+\frac{h^2}{2m}D_2l(l+1)\right]\\ & D_0=1+\frac{3}{\delta}+\frac{3}{\delta^2}\\ & D_1=-\frac{4}{\delta}-\frac{6}{\delta^2}\\ & D_2=\frac{1}{\delta}+\frac{3}{\delta^2} \end{aligned}
δ=aR0, β=δ2β0, γ=δ2γ0β0=−ℏ22mR02[−V0+2mh2D1l(l+1)]γ0=ℏ22mR02[V0+2mh2D2l(l+1)]D0=1+δ3+δ23D1=−δ4−δ26D2=δ1+δ23
LS coupling
产能率及元素变化与核反应速率的关系
现在考虑下述核反应过程:
I
(
A
i
,
Z
i
)
+
J
(
A
j
,
Z
j
)
⇔
K
(
A
k
,
Z
k
)
+
L
(
A
l
,
Z
l
)
I(A_i, Z_i)+J(A_j,Z_j)\Leftrightarrow K(A_k,Z_k)+L(A_l,Z_l)
I(Ai,Zi)+J(Aj,Zj)⇔K(Ak,Zk)+L(Al,Zl)
反应过程中,核子数守恒:
A
i
+
A
j
=
A
k
+
A
l
A_i+A_j = A_k+A_l
Ai+Aj=Ak+Al
电荷守恒:
Z
i
+
Z
j
=
Z
k
+
Z
l
Z_i+Z_j = Z_k+Z_l
Zi+Zj=Zk+Zl
如果在上述反应中,物质
i
i
i、
j
j
j的反应截面为
σ
\sigma
σ,粒子相对速度为
v
v
v,单位体积、单位时间的粒子发生反应的次数为
n
i
n
j
v
σ
n_i n_j v \sigma
ninjvσ,则反应速率为:
R
i
j
k
=
⟨
v
σ
⟩
R_{ijk}=\langle v\sigma \rangle
Rijk=⟨vσ⟩
粒子
i
i
i的数密度变化率为:
n
˙
i
=
−
n
i
∑
j
,
k
(
1
+
δ
i
j
)
n
j
(
1
+
δ
i
j
)
R
i
j
k
+
∑
k
,
l
n
l
n
k
(
1
+
δ
l
k
)
R
l
k
i
\dot n_i=-n_i \sum_{j,k}(1+\delta_{ij})\frac{n_j}{(1+\delta_{ij})}R_{ijk}+\sum_{k,l}\frac{n_l n_k}{(1+\delta_{lk})}R_{lki}
n˙i=−nij,k∑(1+δij)(1+δij)njRijk+k,l∑(1+δlk)nlnkRlki
δ
\delta
δ项的存在是因为同种粒子在反应时,是一半和另一半相互消耗的
由于原子核的数目为
n
i
=
ρ
X
i
m
A
A
i
n_i=\frac{\rho X_i}{m_A A_i}
ni=mAAiρXi,因此,粒子数密度的变化率还可以写成
X
X
X的变化:
X
˙
i
A
i
=
[
−
X
i
A
i
∑
j
,
k
R
i
j
k
X
j
A
j
+
∑
k
,
l
R
l
k
i
(
1
+
δ
l
k
)
X
l
A
l
X
k
A
k
]
ρ
m
A
\frac{\dot X_i}{A_i}=\left[-\frac{X_i}{A_i}\sum_{j,k}R_{ijk}\frac{X_j}{A_j}+\sum_{k,l}\frac{R_{lki}}{(1+\delta_{lk})}\frac{X_l}{A_l}\frac{X_k}{A_k}\right]\frac{\rho}{m_A}
AiX˙i=⎣⎡−AiXij,k∑RijkAjXj+k,l∑(1+δlk)RlkiAlXlAkXk⎦⎤mAρ
由核物理得知,
R
i
j
k
=
R
i
j
k
(
T
)
R_{ijk}=R_{ijk}(T)
Rijk=Rijk(T),因此,可以将
X
˙
i
\dot X_i
X˙i写为下述三个变量表示的形式:
X
˙
i
=
X
˙
(
ρ
,
T
,
X
⃗
)
\dot X_i = \dot X(\rho, T, \vec X)
X˙i=X˙(ρ,T,X)
反应的产能率为:
Q
i
j
k
=
Δ
M
c
2
=
(
M
i
+
M
j
−
M
k
−
M
l
)
c
2
Q_{ijk}=\Delta M c^2=(M_i + M_j - M_k - M_l)c^2
Qijk=ΔMc2=(Mi+Mj−Mk−Ml)c2
一般情况下,应该去掉反应产生的中微子的能量,因为中微子几乎是自由逃逸的。对所有的反应求和,再除以密度:
q
=
ρ
m
A
2
∑
i
j
k
1
1
+
δ
i
j
X
i
A
i
X
j
A
j
R
i
j
k
Q
i
j
k
=
q
(
ρ
,
T
,
X
⃗
)
q=\frac{\rho}{m_A^2}\sum_{ijk}\frac{1}{1+\delta_{ij}}\frac{X_i}{A_i}\frac{X_j}{A_j}R_{ijk}Q_{ijk}=q(\rho,T,\vec X)
q=mA2ρijk∑1+δij1AiXiAjXjRijkQijk=q(ρ,T,X)
核反应的速率
由于短程力的特点,只有当离核足够近的时候才会发生反应。核带来的势能为:
V
(
r
)
=
−
g
2
exp
(
−
m
c
h
r
)
r
+
Z
1
Z
2
e
2
r
V(r)=-g^2\frac{\exp(-\frac{mc}{h}r)}{r}+\frac{Z_1Z_2 e^2}{r}
V(r)=−g2rexp(−hmcr)+rZ1Z2e2
大致图像如下:
其中,
m
c
2
≈
164
M
e
V
mc^2\approx 164\ MeV
mc2≈164 MeV,
r
0
≈
1.44
×
1
0
−
13
A
1
3
c
m
r_0\approx 1.44\times 10^{-13}A^\frac{1}{3}\ cm
r0≈1.44×10−13A31 cm,
r
0
r_0
r0处的库仑势为:
R
C
=
Z
1
Z
2
e
2
r
0
∼
Z
1
Z
2
M
e
V
R_C=\frac{Z_1 Z_2 e^2}{r_0}\sim Z_1 Z_2\ MeV
RC=r0Z1Z2e2∼Z1Z2 MeV
而粒子的动能约为
E
K
≈
T
1.2
×
1
0
7
K
k
e
V
E_K\approx \dfrac{T}{1.2\times 10^7\ K}\ keV
EK≈1.2×107 KT keV,这个量级是远小于
E
C
E_C
EC的。在经典物理中,核碰撞的概率
∼
e
−
E
C
/
k
T
\sim e^{-E_C/kT}
∼e−EC/kT非常地小,但在量子力学的理论下,却有足够的概率能够发生碰撞
记动能与势能相交的那个点(势转折点)的位置为 r t r_t rt(图中为 r C r_C rC),在经典物理中,粒子到达那个位置后就会被推回去,而在量子力学中,则可能会穿透
定义 χ 2 = 2 m E 0 ℏ 2 \chi^2=\frac{2mE_0}{\hbar ^2} χ2=ℏ22mE0,此为粒子在很远的地方的波矢。 r t = Z 1 Z 2 e 2 E 0 r_t=\frac{Z_1Z_2 e^2}{E_0} rt=E0Z1Z2e2,且 r t ≫ r N r_t\gg r_N rt≫rN
根据动能和势能的总和,求解径向部分的不含时薛定谔方程:
−
ℏ
2
2
μ
1
r
2
1
d
r
(
r
2
d
ψ
d
r
)
+
(
E
−
q
1
q
2
r
)
ψ
=
0
-\frac{\hbar ^2}{2\mu}\frac{1}{r^2}\frac{1}{dr}\left(r^2\frac{d\psi}{dr}\right)+\left(E-\frac{q_1q_2}{r}\right)\psi = 0
−2μℏ2r21dr1(r2drdψ)+(E−rq1q2)ψ=0
前面提到的势转折点为
r
t
=
q
1
q
2
ϵ
=
2
q
1
q
2
m
v
2
r_t=\frac{q_1q_2}{\epsilon}=\frac{2q_1q_2}{mv^2}
rt=ϵq1q2=mv22q1q2
在 r ≫ r t r\gg r_t r≫rt处,得到的解 ψ ∝ e − i k r r \psi\propto \frac{e^{-ikr}}{r} ψ∝re−ikr是一个球面波函数
而在
r
<
r
t
r<r_t
r<rt处,得到的波函数应该是指数衰减的,设
ψ
=
ϕ
/
r
\psi=\phi/r
ψ=ϕ/r,则薛定谔方程化为:
−
ℏ
2
2
μ
d
2
ϕ
d
r
2
+
(
E
−
q
1
q
2
r
)
ϕ
=
0
-\frac{\hbar^2}{2\mu}\frac{d^2\phi}{dr^2}+\left(E-\frac{q_1q_2}{r}\right)\phi=0
−2μℏ2dr2d2ϕ+(E−rq1q2)ϕ=0
如果势能项变化足够慢,振荡但幅度与波长变小,WKB近似:
ϕ
=
A
(
r
)
e
i
B
(
r
)
\phi= A(r)e^{iB(r)}
ϕ=A(r)eiB(r)
代入方程
→
\rightarrow
→实部与虚部都满足,忽略高阶项导数,最后得到的解:
ϕ
=
C
p
(
r
)
e
−
i
ℏ
∫
0
r
p
(
r
′
)
d
r
′
,
p
(
r
′
)
=
2
μ
(
E
−
q
1
q
2
r
)
\phi=\frac{C}{\sqrt{p(r)}}e^{-\frac{i}{\hbar}\int^r_0 p(r')dr'},\ p(r')=\sqrt{2\mu\left(E-\frac{q_1q_2}{r}\right)}
ϕ=p(r)Ce−ℏi∫0rp(r′)dr′, p(r′)=2μ(E−rq1q2)
WKB近似:将
ϕ
\phi
ϕ写作
A
(
r
)
e
i
B
(
r
)
A(r)e^{iB(r)}
A(r)eiB(r)的形式(
A
(
r
)
A(r)
A(r)和
B
(
r
)
B(r)
B(r)都是实函数),由此,方程既有实数部分又有虚数部分
利用WKB近似后,得到:
ϕ
(
r
)
=
{
exp
(
i
∫
r
t
t
K
d
r
)
,
∞
>
r
≥
r
t
exp
(
−
∫
r
t
r
K
d
r
)
,
r
≤
r
t
\phi(r)=\left\{\begin{aligned} &\exp(i\int^t_{r_t}Kdr),\ \infty > r \ge r_t\\ &\exp(-\int^r_{r_t}Kdr),\ r\le r_t \end{aligned}\right.
ϕ(r)=⎩⎪⎪⎪⎨⎪⎪⎪⎧exp(i∫rttKdr), ∞>r≥rtexp(−∫rtrKdr), r≤rt
其中,
K
≡
2
μ
ℏ
2
∣
E
−
q
1
q
2
r
∣
K\equiv\sqrt{\frac{2\mu}{\hbar^2}\left|E-\frac{q_1q_2}{r}\right|}
K≡ℏ22μ∣∣E−rq1q2∣∣
则库仑穿透概率为:
P
=
4
π
r
2
∣
ψ
∣
2
(
r
=
0
)
4
π
r
2
∣
ψ
∣
2
(
r
=
∞
)
=
∣
ϕ
(
0
)
∣
2
∣
ϕ
(
∞
)
∣
2
=
exp
(
−
2
χ
∫
0
r
t
r
t
r
−
1
d
r
)
P=\frac{4\pi r^2|\psi|^2(r=0)}{4\pi r^2 |\psi|^2(r=\infty)}=\frac{|\phi(0)|^2}{|\phi(\infty)|^2}=\exp(-2\chi\int^{r_t}_0\sqrt{\frac{r_t}{r}-1}dr)
P=4πr2∣ψ∣2(r=∞)4πr2∣ψ∣2(r=0)=∣ϕ(∞)∣2∣ϕ(0)∣2=exp(−2χ∫0rtrrt−1dr)
其中,
χ
=
2
μ
E
ℏ
2
\chi=\sqrt{\frac{2\mu E}{\hbar ^2}}
χ=ℏ22μE
从0开始积分会有极点,但因为
r
0
≪
r
t
r_0\ll r_t
r0≪rt,所以从
r
0
r_0
r0开始与从0开始几乎无异,积分后得到:
P
=
exp
(
−
2
α
χ
r
t
)
=
exp
(
−
4
γ
q
1
q
2
ℏ
ν
)
P=\exp(-2\alpha\chi r_t)=\exp(-4\gamma\frac{q_1q_2}{\hbar \nu})
P=exp(−2αχrt)=exp(−4γℏνq1q2)
其中,
γ
=
∫
0
1
1
−
x
x
d
x
=
π
2
\gamma=\int^1_0\sqrt{\frac{1-x}{x}}dx=\frac{\pi}{2}
γ=∫01x1−xdx=2π
上述形式也可以写为:
P
=
e
−
2
π
η
P=e^{-2\pi\eta}
P=e−2πη
其中,
η
=
(
μ
/
2
)
1
2
Z
1
Z
2
e
2
ℏ
E
1
/
2
\eta=(\mu/2)^\frac{1}{2}\frac{Z_1Z_2e^2}{\hbar E^{1/2}}
η=(μ/2)21ℏE1/2Z1Z2e2
在非共振反应(没有角动量)下,碰撞参数为德布罗意波长: b = λ = h p b=\lambda=\frac{h}{p} b=λ=ph,碰撞截面为: σ = π λ 2 ∝ 1 E \sigma=\pi\lambda^2\propto \frac{1}{E} σ=πλ2∝E1,非共振反应截面为: σ ( E ) = S 0 E − 1 e − 2 π η \sigma(E)=S_0 E^{-1}e^{-2\pi\eta} σ(E)=S0E−1e−2πη,其中 S 0 S_0 S0需要由实验给出
在热核反应下,粒子能量热分布为:
f
(
E
)
d
E
=
2
π
E
1
2
(
k
T
)
3
2
e
−
E
k
T
d
E
f(E)dE=\frac{2}{\sqrt{\pi}}\frac{E^\frac{1}{2}}{(kT)^\frac{3}{2}}e^{-\frac{E}{kT}}dE
f(E)dE=π2(kT)23E21e−kTEdE
其中,
v
=
2
E
μ
v=\sqrt{\frac{2E}{\mu}}
v=μ2E,
μ
=
m
i
m
j
m
i
+
m
j
\mu=\frac{m_im_j}{m_i+m_j}
μ=mi+mjmimj
核作用截面的理论计算
库仑屏蔽修正
在高密度低温的情形下,原子核附近更有可能有电子分布,从而导致库仑屏蔽效应。对于低质量恒星而言,中心温度低密度高,此时就需要考虑此情形
H燃烧阶段
pp过程
在恒星的氢燃烧阶段,第一种反应是pp过程:
1
H
+
1
H
→
2
H
+
e
+
+
ν
e
^1H+\ ^1H\rightarrow \ ^2H+e^++\nu_e
1H+ 1H→ 2H+e++νe
pp1:
2
H
+
1
H
→
3
H
e
+
γ
3
H
e
+
3
H
e
→
4
H
e
+
1
H
+
1
H
^2H+\ ^1H\rightarrow\ ^3 He+\gamma\\ \ ^3 He+\ ^3He\rightarrow\ ^4He+\ ^1H+\ ^1H
2H+ 1H→ 3He+γ 3He+ 3He→ 4He+ 1H+ 1H
pp2:
3
H
e
+
4
H
e
→
7
B
e
+
γ
7
B
e
+
e
−
→
7
L
i
+
ν
e
+
γ
7
L
i
+
1
H
→
4
H
e
+
4
H
e
^3He+\ ^4He\rightarrow\ ^7Be+\gamma\\ \ ^7Be+ e^-\rightarrow\ ^7Li+\nu_e+\gamma\\ \ ^7Li+\ ^1H\rightarrow\ ^4He+\ ^4He
3He+ 4He→ 7Be+γ 7Be+e−→ 7Li+νe+γ 7Li+ 1H→ 4He+ 4He
pp3:
7
B
e
+
1
H
→
8
B
+
γ
8
B
→
8
B
e
+
e
+
+
ν
e
8
B
e
→
4
H
e
+
4
H
e
^7Be+\ ^1H\rightarrow\ ^8B+\gamma\\ \ ^8B\rightarrow\ ^8Be+e^+ +\nu_e\\ \ ^8Be\rightarrow \ ^4He+\ ^4He
7Be+ 1H→ 8B+γ 8B→ 8Be+e++νe 8Be→ 4He+ 4He
此反应的产能为26.73
M
e
V
MeV
MeV,扣除中微子后约为26
M
e
V
MeV
MeV
pp1的反应速度最慢,弱作用,决定了整个氢燃烧过程的速度。温度越高,pp2和pp3越重要。最慢的过程: 1 H + 1 H → 2 H + e + + ν e ^1H+\ ^1H\rightarrow \ ^2H+e^++\nu_e 1H+ 1H→ 2H+e++νe决定了pp过程和反应速度
pp反应在单位体积、单位时间内的速率:
r
p
p
=
5.34
×
1
0
7
T
7
−
2
3
X
2
ρ
2
e
−
15.689
T
7
1
/
3
c
m
−
3
s
−
1
r_{pp}=5.34\times 10^7 T_7^{-\frac{2}{3}}X^2\rho^2 e^{-\frac{15.689}{T_7^{1/3}}}\ cm^{-3}s^{-1}
rpp=5.34×107T7−32X2ρ2e−T71/315.689 cm−3s−1
单位质量、单位时间的产能率:
q
e
f
f
(
p
p
−
c
h
a
i
n
s
)
≈
1.1
×
1
0
3
ρ
X
2
T
7
2
/
3
e
−
15.689
T
7
1
/
3
e
r
g
g
−
1
s
−
1
q_{eff}(pp-chains)\approx 1.1\times 10^3 \frac{\rho X^2}{T_7^{2/3}}e^{-\frac{15.689}{T_7^{1/3}}}\ erg\ g^{-1}\ s^{-1}
qeff(pp−chains)≈1.1×103T72/3ρX2e−T71/315.689 erg g−1 s−1
幂律指数:
ν
p
p
=
5.25
T
7
1
/
3
−
2
3
\nu_{pp}=\frac{5.25}{T_7^{1/3}}-\frac{2}{3}
νpp=T71/35.25−32
氢燃烧的时标为: τ p = − n p d n p d t = n p 2 r p p \tau_p=-\frac{n_p}{\frac{dn_p}{dt}}=\frac{n_p}{2r_{pp}} τp=−dtdnpnp=2rppnp
在太阳中心,密度约为 100 g c m − 3 100\ g\ cm^{-3} 100 g cm−3,温度为 15 M K 15\ MK 15 MK,燃烧时标为 6 × 1 0 9 y r 6\times 10^9\ yr 6×109 yr
在pp反应中,有一些过程非常快,这些快速反应的元素在反应开始后不久就达到了反应平衡,从而数密度几乎不随时间变化,数密度取决于反应速率
CNO循环
CN循环
CN循环的过程:
12
C
+
1
H
→
13
N
+
γ
13
N
→
13
C
+
e
+
+
ν
e
13
C
+
1
H
→
14
N
+
γ
14
N
+
1
H
→
15
O
+
γ
15
O
→
15
N
+
e
+
+
ν
e
15
N
+
1
H
→
12
C
+
4
H
e
^{12}C+\ ^1H\rightarrow\ ^{13}N+\gamma\\ \ ^{13}N\rightarrow\ ^{13}C+e^++\nu_e\\ \ ^{13}C+\ ^1H\rightarrow\ ^{14}N+\gamma\\ \ ^{14}N+\ ^1H\rightarrow\ ^{15}O+\gamma\\ \ ^{15}O\rightarrow\ ^{15}N+ e^+ + \nu_e\\ \ ^{15}N+\ ^1H\rightarrow\ ^{12}C+\ ^4He
12C+ 1H→ 13N+γ 13N→ 13C+e++νe 13C+ 1H→ 14N+γ 14N+ 1H→ 15O+γ 15O→ 15N+e++νe 15N+ 1H→ 12C+ 4He
上述反应链可以简写为:
12
C
(
p
,
γ
)
13
N
(
,
e
+
ν
e
)
13
C
(
p
,
γ
)
14
N
(
p
,
γ
)
15
O
(
,
e
+
ν
e
)
15
N
(
p
,
α
)
12
C
^{12}C(p,\gamma)^{13}N(,e^+\nu_e)^{13}C(p,\gamma)^{14}N(p,\gamma)^{15}O(,e^+\nu_e)^{15}N(p,\alpha)^{12}C
12C(p,γ)13N(,e+νe)13C(p,γ)14N(p,γ)15O(,e+νe)15N(p,α)12C
整条CN反应链的结果为:
4
1
H
→
4
H
e
+
2
e
+
+
2
ν
e
+
3
γ
4^1H\rightarrow\ ^4He+2e^+ + 2\nu_e+3\gamma
41H→ 4He+2e++2νe+3γ
完整的CNO循环
在CN循环中,中间产物
15
N
^{15}N
15N还有另外的反应途径,由此,整个CNO循环为:
12
C
(
p
,
γ
)
13
N
(
,
e
+
ν
e
)
13
C
(
p
,
γ
)
14
N
(
p
,
γ
)
15
O
(
,
e
+
ν
e
)
15
N
(
p
,
α
)
12
C
15
N
(
p
,
γ
)
16
O
(
p
,
γ
)
17
F
(
,
e
+
ν
e
)
17
O
(
p
,
α
)
14
N
17
O
(
p
,
γ
)
18
F
(
e
−
,
ν
e
)
18
O
(
p
,
γ
)
19
F
(
p
,
α
)
16
O
^{12}C(p,\gamma)^{13}N(,e^+\nu_e)^{13}C(p,\gamma)^{14}N(p,\gamma)^{15}O(,e^+\nu_e)^{15}N(p,\alpha)^{12}C\\ \ ^{15}N(p,\gamma)^{16}O(p,\gamma)^{17}F(,e^+\nu_e)^{17}O(p,\alpha)^{14}N\\ \ ^{17}O(p,\gamma)^{18}F(e^-,\nu_e)^{18}O(p,\gamma)^{19}F(p,\alpha)^{16}O
12C(p,γ)13N(,e+νe)13C(p,γ)14N(p,γ)15O(,e+νe)15N(p,α)12C 15N(p,γ)16O(p,γ)17F(,e+νe)17O(p,α)14N 17O(p,γ)18F(e−,νe)18O(p,γ)19F(p,α)16O
在CNO循环中,CNO数目不变,每次循环的结果都是
4
1
H
→
4
H
e
4^1H\rightarrow\ ^4He
41H→ 4He
CNO循环开始时只要有其中一种元素就可以,平衡时得出三者的比例,这个比例与温度和密度有关
在低温情况下,反应速率由俘获速率决定。而在高温情况下,反应速率由逆 β \beta β衰变决定。温度极高时,反应速率与温度无关
CN循环的反应速度最快,因此CN循环这一部分最先达到反应平衡。CN循环中的
15
N
(
p
,
α
)
12
C
^{15}N(p,\alpha)^{12}C
15N(p,α)12C是强相互作用,而其他的反应则是电磁作用,因此
15
N
^{15}N
15N最容易达到平衡。在低温状况下,
β
\beta
β衰变过程相对来说较快,可以看成是瞬间发生的,由此可以写出CN平衡方程:
d
n
15
d
t
=
n
p
(
n
14
R
14
−
n
15
R
15
)
=
0
d
n
14
d
t
=
n
p
(
n
13
R
13
−
n
14
R
14
)
d
n
13
d
t
=
n
p
(
n
12
R
12
−
n
13
R
13
)
d
n
12
d
t
=
n
p
(
n
15
R
15
−
n
12
R
12
)
\begin{aligned} &\frac{dn_{15}}{dt}=n_p(n_{14}R_{14}-n_{15}R_{15})=0\\ &\frac{dn_{14}}{dt}=n_p(n_{13}R_{13}-n_{14}R_{14})\\ &\frac{dn_{13}}{dt}=n_p(n_{12}R_{12}-n_{13}R_{13})\\ &\frac{dn_{12}}{dt}=n_p(n_{15}R_{15}-n_{12}R_{12}) \end{aligned}
dtdn15=np(n14R14−n15R15)=0dtdn14=np(n13R13−n14R14)dtdn13=np(n12R12−n13R13)dtdn12=np(n15R15−n12R12)
以及:
d
d
t
(
n
12
+
n
13
+
n
14
)
=
0
\frac{d}{dt}(n_{12}+n_{13}+n_{14})=0
dtd(n12+n13+n14)=0
其中,反应速率系数为:
R
12
=
3.39
×
1
0
−
17
T
9
−
2
/
3
exp
(
−
13.690
T
9
−
1
/
3
)
R
13
=
1.33
×
1
0
−
16
T
9
−
2
/
3
exp
(
−
13.717
T
9
−
1
/
3
)
R
14
=
8.44
×
1
0
−
17
T
9
−
2
/
3
exp
(
−
15.228
T
9
−
1
/
3
)
R
15
=
1.36
×
1
0
−
12
T
9
−
2
/
3
exp
(
−
15.251
T
9
−
1
/
3
)
\begin{aligned} R_{12}&= 3.39\times 10^{-17}T_9^{-2/3}\exp (-13.690T_9^{-1/3})\\ R_{13}&= 1.33\times 10^{-16}T_9^{-2/3}\exp (-13.717T_9^{-1/3})\\ R_{14}&= 8.44\times 10^{-17}T_9^{-2/3}\exp (-15.228T_9^{-1/3})\\ R_{15}&= 1.36\times 10^{-12}T_9^{-2/3}\exp (-15.251T_9^{-1/3}) \end{aligned}
R12R13R14R15=3.39×10−17T9−2/3exp(−13.690T9−1/3)=1.33×10−16T9−2/3exp(−13.717T9−1/3)=8.44×10−17T9−2/3exp(−15.228T9−1/3)=1.36×10−12T9−2/3exp(−15.251T9−1/3)
CNO循环的主要特点在于:
12
C
+
13
C
+
14
N
^{12}C+\ ^{13}C+\ ^{14}N
12C+ 13C+ 14N的总数目不变。
14
N
^{14}N
14N反应最慢,在反应开始时增加数目,为主要合成途径。平衡时数密度与反应速率系数成反比
对于CN循环而言,单位质量的产能率为:
q
C
N
O
≈
4.4
×
1
0
25
ρ
X
Z
T
9
2
/
3
e
−
15.228
T
9
1
/
3
e
r
g
g
−
1
s
−
1
q_{CNO}\approx 4.4\times 10^{25}\frac{\rho XZ}{T_9^{2/3}}e^{-\frac{15.228}{T_9^{1/3}}}\ erg\ g^{-1}\ s^{-1}
qCNO≈4.4×1025T92/3ρXZe−T91/315.228 erg g−1 s−1
其中,
T
9
=
T
1
0
9
K
T_9=\frac{T}{10^9 K}
T9=109KT
幂律指数随温度变化。当 T = 2 × 1 0 7 K T=2\times 10^7 K T=2×107K时, ν ≈ 18 \nu\approx 18 ν≈18
He燃烧阶段
He燃烧阶段相应的反应为:
4
H
e
+
4
H
e
→
8
B
e
^4He+\ ^4He\rightarrow\ ^8Be
4He+ 4He→ 8Be
8
B
e
^8Be
8Be的能量比He高出
92
k
e
V
92\ keV
92 keV,反应吸热,因此Be不稳定(Be的寿命只有
τ
=
2.6
×
1
0
−
16
s
\tau=2.6\times 10^{-16}\ s
τ=2.6×10−16 s),容易分裂成两个He:
α
+
α
↔
8
B
e
\alpha+\alpha\leftrightarrow\ ^8Be
α+α↔ 8Be
Be的密度由上述反应的平衡给出,类似于Saha方程:
n
α
2
n
(
8
B
e
)
=
(
π
m
α
k
T
h
2
)
3
2
e
−
Q
k
T
Q
=
−
91.78
k
e
V
\frac{n_\alpha^2}{n(^8Be)}=\left(\frac{\pi m_\alpha kT}{h^2}\right)^\frac{3}{2}e^{-\frac{Q}{kT}}\\ Q=-91.78\ keV
n(8Be)nα2=(h2πmαkT)23e−kTQQ=−91.78 keV
为了保证一定的Be丰度,温度必须接近于
92
k
e
V
92\ keV
92 keV
在
ρ
∼
1
0
6
g
c
m
−
3
\rho\sim 10^6\ g\ cm^{-3}
ρ∼106 g cm−3,
n
α
≈
1.5
×
1
0
29
c
m
−
3
n_{\alpha}\approx 1.5\times 10^{29}\ cm^{-3}
nα≈1.5×1029 cm−3,
T
≈
1
0
8
K
T\approx 10^8\ K
T≈108 K时:
n
(
8
B
e
)
n
α
∼
7
×
1
0
−
9
\frac{n(^8Be)}{n_\alpha}\sim 7\times 10^{-9}
nαn(8Be)∼7×10−9
此外,少量的Be会俘获
α
\alpha
α粒子,从而到达一个共振能级
C
∗
C^\ast
C∗,随后它们中的多数将会回到Be+
α
\alpha
α:
8
B
e
+
α
↔
C
∗
^8Be+\alpha \leftrightarrow C^\ast
8Be+α↔C∗
由核的Saha方程,给出
C
∗
C^\ast
C∗的密度:
n
(
8
B
e
)
n
(
4
H
e
)
n
(
12
C
∗
)
=
f
α
B
e
(
2
π
μ
α
B
e
k
T
h
2
)
3
/
2
e
−
χ
k
T
\frac{n(^8Be)n(^4He)}{n(^{12}C^\ast)}=f_{\alpha Be}\left(\frac{2\pi \mu_{\alpha Be}kT}{h^2}\right)^{3/2}e^{-\frac{\chi}{kT}}
n(12C∗)n(8Be)n(4He)=fαBe(h22πμαBekT)3/2e−kTχ
其中,
f
f
f为电子屏蔽因子
3
α
\alpha
α反应的速率为
n
(
12
C
∗
)
τ
\dfrac{n(^{12}C^\ast)}{\tau}
τn(12C∗),
Q
=
7.275
M
e
V
Q=7.275\ MeV
Q=7.275 MeV,单位质量产能率为:
q
=
n
(
12
C
∗
)
τ
ρ
Q
=
5
×
1
0
8
ρ
2
Y
3
(
T
1
0
9
K
)
−
3
e
−
4.4027
T
9
e
r
g
g
−
1
s
−
1
q=\frac{n(^{12}C^\ast)}{\tau\rho}Q=5\times 10^8 \rho^2 Y^3\left(\frac{T}{10^9\ K}\right)^{-3}e^{-\frac{4.4027}{T_9}}\ erg\ g^{-1}\ s^{-1}
q=τρn(12C∗)Q=5×108ρ2Y3(109 KT)−3e−T94.4027 erg g−1 s−1
幂律指数
ν
3
α
≈
4.4
T
9
−
3
\nu_{3\alpha}\approx \dfrac{4.4}{T_9}-3
ν3α≈T94.4−3。当
T
=
1
0
8
K
T=10^8\ K
T=108 K时,指数约为40
He燃烧原则上可以通过
α
\alpha
α俘获,产生O和Ne:
12
C
+
α
→
16
O
+
γ
,
Q
=
7.162
M
e
V
16
O
+
α
→
20
N
e
+
γ
,
Q
=
4.734
M
e
V
20
N
e
+
α
→
24
M
g
+
γ
,
Q
=
9.31
M
e
V
\begin{aligned} ^{12}C+\alpha\rightarrow\ ^{16}O+\gamma,\ Q=7.162\ MeV\\ ^{16}O+\alpha\rightarrow\ ^{20}Ne+\gamma,\ Q=4.734\ MeV\\ ^{20}Ne+\alpha\rightarrow\ ^{24}Mg+\gamma,\ Q=9.31\ MeV \end{aligned}
12C+α→ 16O+γ, Q=7.162 MeV16O+α→ 20Ne+γ, Q=4.734 MeV20Ne+α→ 24Mg+γ, Q=9.31 MeV
需要高温条件,才能克服库仑势垒。在一般的温度条件下,这样的反应比例不高
C、O燃烧阶段
C、O的燃烧对温度的要求更高。其中,C的燃烧温度高于 5 × 1 0 8 K 5\times 10^8\ K 5×108 K,O的燃烧温度高于 1 0 9 K 10^9 K 109K
C燃烧
在此阶段,主要的反应为:
12
C
+
12
C
→
24
M
g
+
γ
,
13.930
M
e
V
12
C
+
12
C
→
23
N
a
+
p
,
2.238
M
e
V
12
C
+
12
C
→
20
N
e
+
α
,
4.616
M
e
V
12
C
+
12
C
→
23
M
g
+
n
,
−
2.605
M
e
V
12
C
+
12
C
→
16
O
+
2
α
,
−
0.114
M
e
V
\begin{aligned} ^{12}C+\ ^{12}C&\rightarrow\ ^{24}Mg+\gamma,\ 13.930\ MeV\\ ^{12}C+\ ^{12}C&\rightarrow\ ^{23}Na+p,\ 2.238\ MeV\\ ^{12}C+\ ^{12}C&\rightarrow\ ^{20}Ne+\alpha,\ 4.616\ MeV\\ ^{12}C+\ ^{12}C&\rightarrow\ ^{23}Mg+n,\ -2.605\ MeV\\ ^{12}C+\ ^{12}C&\rightarrow\ ^{16}O+2\alpha,\ -0.114\ MeV \end{aligned}
12C+ 12C12C+ 12C12C+ 12C12C+ 12C12C+ 12C→ 24Mg+γ, 13.930 MeV→ 23Na+p, 2.238 MeV→ 20Ne+α, 4.616 MeV→ 23Mg+n, −2.605 MeV→ 16O+2α, −0.114 MeV
其中,反应1是最主要的,反应2/3等概率次之,反应4/5则是最低的
此反应的反应速率为:
R
12
,
12
=
2
×
1
0
4
f
12
,
12
T
9
−
2
/
3
exp
(
−
36.55
(
1
+
0.1
T
9
)
1
/
3
T
9
1
/
3
)
R_{12,12}=2\times 10^4 f_{12,12}T_9^{-2/3}\exp(-\frac{36.55(1+0.1T_9)^{1/3}}{T_9^{1/3}})
R12,12=2×104f12,12T9−2/3exp(−T91/336.55(1+0.1T9)1/3)
单位质量产能率为:
q
=
2.6
×
1
0
40
ρ
X
12
2
R
12
,
12
e
r
g
g
−
1
s
−
1
q=2.6\times 10^{40}\rho X_{12}^2 R_{12,12}\ erg\ g^{-1}\ s^{-1}
q=2.6×1040ρX122R12,12 erg g−1 s−1
在C燃烧过程中,也有产生中子的反应:
12
C
(
p
,
γ
)
13
N
(
,
e
+
ν
e
)
13
C
(
α
,
n
)
16
O
^{12}C(p,\gamma)^{13}N(,e^+\nu_e)^{13}C(\alpha,n)^{16}O
12C(p,γ)13N(,e+νe)13C(α,n)16O
O燃烧
O的燃烧与C的燃烧类似,但需要更高的温度:
T
=
1
0
9
K
T=10^9\ K
T=109 K:
16
O
+
16
O
→
32
S
+
γ
,
16.539
M
e
V
16
O
+
16
O
→
31
P
+
p
,
7.676
M
e
V
16
O
+
16
O
→
28
S
i
+
α
,
9.593
M
e
V
16
O
+
16
O
→
31
S
+
n
,
1.459
M
e
V
16
O
+
16
O
→
24
M
g
+
2
α
,
−
0.393
M
e
V
\begin{aligned} ^{16}O+\ ^{16}O&\rightarrow\ ^{32}S+\gamma,\ 16.539\ MeV\\ ^{16}O+\ ^{16}O&\rightarrow\ ^{31}P+p,\ 7.676\ MeV\\ ^{16}O+\ ^{16}O&\rightarrow\ ^{28}Si+\alpha,\ 9.593\ MeV\\ ^{16}O+\ ^{16}O&\rightarrow\ ^{31}S+n,\ 1.459\ MeV\\ ^{16}O+\ ^{16}O&\rightarrow\ ^{24}Mg+2\alpha,\ -0.393\ MeV \end{aligned}
16O+ 16O16O+ 16O16O+ 16O16O+ 16O16O+ 16O→ 32S+γ, 16.539 MeV→ 31P+p, 7.676 MeV→ 28Si+α, 9.593 MeV→ 31S+n, 1.459 MeV→ 24Mg+2α, −0.393 MeV
主要产物为
28
S
i
^{28}Si
28Si,反应速率为:
R
16
,
16
=
7.94
×
1
0
17
f
16
,
16
T
9
−
2
/
3
exp
(
−
59.01
(
1
+
0.14
T
9
)
1
/
3
T
9
1
/
3
)
R_{16,16}=7.94\times 10^{17}f_{16,16} T^{-2/3}_9\exp(-\frac{59.01(1+0.14T_9)^{1/3}}{T_9^{1/3}})
R16,16=7.94×1017f16,16T9−2/3exp(−T91/359.01(1+0.14T9)1/3)
单位质量产能率为:
q
=
2
×
1
0
40
ρ
X
16
2
R
16
,
16
e
r
g
g
−
1
s
−
1
q=2\times 10^{40}\rho X^2_{16}R_{16,16}\ erg\ g^{-1}\ s^{-1}
q=2×1040ρX162R16,16 erg g−1 s−1
Si燃烧
在温度较高时,光子的能量高,足以分解粒子:
20
N
e
+
γ
↔
16
O
+
α
^{20}Ne+\gamma\leftrightarrow\ ^{16}O+\alpha
20Ne+γ↔ 16O+α
通过光分解产生高能的
α
\alpha
α粒子具有较高能量,可以被其他粒子俘获:
20
N
e
+
α
→
24
M
g
+
γ
^{20}Ne+\alpha\rightarrow\ ^{24}Mg+\gamma
20Ne+α→ 24Mg+γ
Si燃烧产生Fe族元素,直接反应温度要求很高(Si离解的温度大于
3
×
1
0
9
K
3\times 10^9\ K
3×109 K),Si离解产生的轻粒子被Si吸收,从而合成更重的元素。此过程中的反应速度基本由光离解速度决定
光离解过程:
X
+
Y
↔
Z
+
γ
X+Y\leftrightarrow Z+\gamma
X+Y↔Z+γ
在反应达到平衡时:
n
Z
R
γ
=
n
X
n
Y
R
X
Y
n_ZR_\gamma=n_Xn_YR_{XY}
nZRγ=nXnYRXY
在平衡后,Saha方程:
n
X
n
Y
n
Z
=
g
X
g
Y
g
Z
(
2
π
μ
k
T
h
2
)
3
/
2
e
−
Q
k
T
\frac{n_Xn_Y}{n_Z}=\frac{g_Xg_Y}{g_Z}\left(\frac{2\pi\mu kT}{h^2}\right)^{3/2}e^{-\frac{Q}{kT}}
nZnXnY=gZgXgY(h22πμkT)3/2e−kTQ
由上述两式得到光离解速率系数:
R
γ
=
g
X
g
Y
g
Z
(
2
π
μ
k
T
h
2
)
3
/
2
e
−
Q
k
T
R
X
Y
=
5.943
×
1
0
33
g
X
g
Y
g
Z
(
A
T
9
)
3
/
2
exp
(
−
11.605
Q
T
9
)
R
X
Y
\begin{aligned} R_\gamma&=\frac{g_Xg_Y}{g_Z}\left(\frac{2\pi\mu kT}{h^2}\right)^{3/2}e^{-\frac{Q}{kT}}R_{XY}\\ &=5.943\times 10^{33}\frac{g_Xg_Y}{g_Z}(AT_9)^{3/2}\exp(-\frac{11.605Q}{T_9})R_{XY} \end{aligned}
Rγ=gZgXgY(h22πμkT)3/2e−kTQRXY=5.943×1033gZgXgY(AT9)3/2exp(−T911.605Q)RXY
离解出的小粒子被大的核俘获,形成更大的核,且这个反应是可逆的:
在Si燃烧反应中,就有很多元素参与多种反应,例如质子俘获和中子俘获