恒星结构和演化-学习记录3-第三章-物态方程2

第三章：物态方程（2）

多方球与白矮星

压强随距球心距离的变化

  多方球：物态变化是一个多方过程

  白矮星随着时间的演化，温度会变得比较低，密度会变得非常高。中间没有能源，因此白矮星的情形相对简单

  物态为电子简并，因此压强（电子简并压）几乎与温度没有关系：$P\sim P_e(\rho )$，而费米能与密度有关

  在非相对论性的情况下：
$$\begin{gathered} P_e=\frac{1}{20}{\left(\frac{3}{\pi }\right)}^{\frac{2}{3}}{\left(\frac{h}{m_{e}c}\right)}^2{\left(\frac{1}{{\mu }_e{m}_A}\right)}^{\frac{5}{3}}{m}_e{c}^2{\rho }^{\frac{5}{3}} \\ {x}_F={\left(\frac{\rho }{9.7\times 10^5{\mu }_{e}gc{m}^{-3}}\right)}^{\frac{1}{3}}\ll 1 \end{gathered}$$
  在相对论性的情况下：
$$\begin{gathered} P_e=\frac{1}{8}{\left(\frac{3}{\pi }\right)}^{\frac{1}{3}}\left(\frac{h}{m_{e}c}\right){\left(\frac{1}{{\mu }_e{m}_A}\right)}^{\frac{4}{3}}{m}_e{c}^2{\rho }^{\frac{4}{3}} \\ {x}_F\gg 1 \end{gathered}$$
  再次注意到压强只和密度有关

  因为离子的压强比电子的压强要低得多，所以只考虑电子压强（因为在简并的情况下，费米能比热动能大得多）

  可以将上述的压强表达式写为$P=K_{\gamma }{\rho }^{\gamma }=K_n{\rho }^{1+\frac{1}{n}}$形式，其中，$n=\frac{1}{\gamma -1}$称为多方指数，由此，可以由上述的表达式得知$K_n$系数的值：
$$\begin{gathered} K_{\frac{5}{3}}=1.0036\times 10^{13}{\mu }_e^{-5/3}(cgs) \\ {K}_{\frac{4}{3}}=1.2435\times 10^{15}{\mu }_e^{-4/3}(cgs) \end{gathered}$$
  现在考虑稳态球对称的前两个方程：
$$\begin{gathered} \frac{dP}{dr}=-\frac{Gm}{r^2}\rho \\ \frac{dm}{dr}=4\pi r^2\rho \end{gathered}$$
  结合后得到：
$$\frac{1}{r^2}\frac{d}{dr}\left(\frac{r^2}{\rho }\frac{dP}{dr}\right)=-4\pi G\rho$$
  将$P=K_n{\rho }^{1+\frac{1}{n}}$代入到上式中，即可得到一个关于$\rho$的微分方程

微分方程求解

  现在求解这个微分方程。首先，将$\rho$写作$\rho \equiv {\rho }_c{\theta }^n$。其中${\rho }_c$为中心密度，$\theta$是$r$的函数，无量纲。由此，$P=K{\rho }^{1+\frac{1}{n}}=K{\rho }_c^{1+\frac{1}{n}}{\theta }^{1+n}$

  将这个形式的$P$，代入到微分方程中，得到：
$$\left[\frac{(n+1)K{\rho }_c^{\frac{1}{n}-1}}{4\pi G}\right]\frac{1}{r^2}\frac{d}{dr}\left(r^2\frac{d\theta }{dr}\right)=-{\theta }^n$$
  由于$\left[\frac{(n+1)K{\rho }_c^{\frac{1}{n}-1}}{4\pi G}\right]$的量纲为长度的平方，因此，可以记$r_n={\left[\frac{(n+1)K{\rho }_c^{\frac{1}{n}-1}}{4\pi G}\right]}^{1/2}$，$r=\xi r_n$，方程最后写为：

$$\frac{1}{{\xi }^2}\frac{d}{d\xi }\left({\xi }^2\frac{d\theta }{d\xi }\right)=-{\theta }^n$$
  称作Lane-Rmden方程

  对于二阶常微分方程，需要两个边界条件。在这里的条件为，在中心$\xi =0$处：
$$\begin{gathered} \rho ={\rho }_c,\theta (0)=1 \\ \frac{dP}{dr}=0,{\theta }^{\prime }(0)=0 \end{gathered}$$
  对方程从里向外积分，积到$\theta ({\xi }_1)=0$处，就积到了表面

  当$n=0,1,5$时，积分可以得到解析解：
$$\begin{gathered} {\theta }_0(\xi )=1-\frac{{\xi }^2}{6} \\ {\theta }_1(\xi )=\frac{\sin \xi }{\xi } \\ {\theta }_5(\xi )={\left[1+\frac{{\xi }^2}{3}\right]}^{-\frac{1}{2}} \end{gathered}$$
  而对于其它的$n$值，只能求数值解。对于相同半径的恒星，$n$值越大，质量越聚集在中心（注意归一化）

$\theta$的物理意义

  将基本方程$\frac{dP(r)}{dr}=-\rho \frac{d\Phi }{dr}$写成$\frac{d\Phi }{dr}=-\frac{1}{\rho }\frac{dP(r)}{dr}$

  代入$P=K{\rho }^{1+\frac{1}{n}}$，得到：
$$\frac{d\Phi }{dr}=-K\frac{n+1}{n}{\rho }^{\frac{1}{n}-1}\frac{d\rho }{dr}=-(n+1)K\frac{d{\rho }^{\frac{1}{n}}}{dr}$$
  进一步得到：
$$\Phi =-(n+1)K{\rho }^{\frac{1}{n}}+C=-(n+1)K{\rho }_c^{\frac{1}{n}}\theta +C={\Phi }_c\theta$$
  积分常数$C$取0，即表面引力势取零，得到关系：
$$\theta =\frac{\Phi }{{\Phi }_c},\frac{P}{\rho }=-\frac{\Phi }{n+1}$$
  由此得知$\theta$的物理意义

其他物理量

  在得到多方球的解后，可以分别求出物理量$\frac{M(r)}{M}$，$\frac{\rho (r)}{{\rho }_c}$，$\frac{P(r)}{P_c}$和$\frac{g(r)}{g(R)}$

  首先：
$$m(r)={\int }_0^{r}4\pi r^2\rho dr=4\pi {\rho }_c{\int }_0^r{\theta }^n{r}^{2}dr=4\pi {\rho }_c\frac{r^3}{{\xi }^3}{\int }_0^{\xi }{\theta }^n{\xi }^{2}d\xi$$
  由前面提到的Lana-Emden方程，得到：
$$\begin{gathered} m(r)=4\pi {\rho }_c{r}^3\left(-\frac{1}{\xi }\frac{d\theta }{d\xi }\right) \\ M(r)=4\pi {\rho }_c{R}^3{\left(-\frac{1}{\xi }\frac{d\theta }{d\xi }\right)}_{\xi ={\xi }_1} \\ \bar{\rho }=3{\rho }_c{\left(-\frac{1}{\xi }\frac{d\theta }{d\xi }\right)}_{\xi ={\xi }_1} \end{gathered}$$
  对$\frac{d\Phi }{dr}=-\frac{Gm(r)}{r^2}$积分：
$$\begin{array}{rl}\Phi (r)& =-4\pi G{\rho }_c{\int }_r^R{r}^{\prime }\left(-\frac{1}{{\xi }^{\prime }}\frac{d\theta }{d{\xi }^{\prime }}\right)d{r}^{\prime }\\ & =-4\pi G{\rho }_c{r}_n^2{\int }_{\xi }^{{\xi }_1}{\xi }^{\prime }\left(-\frac{1}{{\xi }^{\prime }}\frac{d\theta }{d{\xi }^{\prime }}\right)d{\xi }^{\prime }\\ & =-4\pi G{\rho }_c{r}_n^2\theta (\xi )\end{array}$$

$${\Phi }_c=-4\pi G\frac{M}{4\pi R^3}\frac{1}{{\left(-\frac{1}{\xi }\frac{d\theta }{d\xi }\right)}_{{\xi }_1}}\frac{R^2}{{\xi }_1^2}=-\frac{GM}{R}\frac{1}{{\left(-\xi \frac{d\theta }{d\xi }\right)}_{{\xi }_1}}$$

  中心压强：
$$\begin{array}{rl}{P}_c& =-\frac{{\Phi }_c}{n+1}{\rho }_c\\ & =\frac{GM}{R}\frac{1}{{\left(-\xi \frac{d\theta }{d\xi }\right)}_{{\xi }_1}}\frac{M}{4\pi R^3}\frac{1}{{\left(-\frac{1}{\xi }\frac{d\theta }{d\xi }\right)}_{{\xi }_1}}\\ & =\frac{G{M}^2}{4\pi R^4}\frac{1}{{\left(\frac{d\theta }{d\xi }\right)}_{{\xi }_1}^2}\end{array}$$

应用到白矮星

  在非相对论性情形下：
$$R={\xi }_1{r}_n=3.65r_n=3.65{\left[\frac{(n+1)K{\rho }_c^{\frac{1}{n}-1}}{4\pi G}\right]}^{\frac{1}{2}}$$
  将$n=\frac{3}{2}$，$\frac{{\rho }_c}{\bar{\rho }}=5.99$及$\bar{\rho }=\frac{3}{4\pi }\frac{M}{R^3}$代入上式，得到质量与半径的关系：
$$\begin{gathered} R=3.65{\left[\frac{5K_{3/2}}{8\pi G}{\left(5.99\frac{3}{4\pi }\frac{M}{R^3}\right)}^{-\frac{1}{3}}\right]}^{\frac{1}{2}} \\ {M}^{\frac{1}{3}}R=\frac{2.353}{G}{K}_{3/2}=3.5404\times 10^{20}{\mu }_e^{-5/3}(cgs)=0.04048{\mu }_e^{-5/3}{M}_{\odot }^{1/3}{R}_{\odot } \end{gathered}$$
  在相对论性情形下：
$$\begin{gathered} R={\xi }_1{r}_n=6.9{\left[\frac{K_3}{\pi G}\right]}^{1/2}{\rho }_c^{-1/3}=6.9{\left[\frac{K_3}{\pi G}\right]}^{1/2}{\left(\frac{4\pi }{3\ast 54.2}\right)}^{1/3}{M}^{-1/3}R \\ M=6.9\times \frac{4}{3\ast 54.2\sqrt{\pi }}{\left[\frac{K_3}{G}\right]}^{\frac{3}{2}}=4.559{\left[\frac{K_3}{G}\right]}^{\frac{3}{2}}=1.161\times \frac{10^{34}g}{{\mu }_e^2}\approx 1.46M_{\odot } \end{gathered}$$
  此即白矮星的钱德拉塞卡质量（临界质量）

多方球的引力能

  利用物态方程和位力定理，可以给出多方球的引力势能

  将物态方程$P=K{\rho }^{1+\frac{1}{n}}$代入到流体静力平衡方程：
$$\frac{1}{\rho }\frac{dP}{dr}=(n+1)\frac{d\left(\frac{P}{\rho }\right)}{dr}=-\frac{Gm}{r^2}$$
  并对$r$进行积分：${\int }_r^{\infty }dr$：
$$\begin{gathered} (n+1)\frac{P}{\rho }=-{\int }_r^{\infty }\frac{GM}{r^{\prime 2}}d{r}^{\prime }+{\int }_R^{\infty }\frac{GM}{r^{\prime 2}}d{r}^{\prime }=-\Phi -\frac{GM}{R} \\ \Omega =\frac{1}{2}{\int }_0^R\Phi dm=-\left(\frac{n+1}{2}{\int }_0^R\frac{P}{\rho }dm+\frac{1}{2}\frac{G{M}^2}{R}\right)=-\frac{n+1}{2}{\int }_0^{R}PdV-\frac{1}{2}\frac{G{M}^2}{R} \end{gathered}$$
  在计算引力势能时，应该取无穷远处势能为0。右侧第一项为引力势，第二项为从$r$积到$R$所出现的项。再对各个位置的引力势进行积分，即可得到引力能$\Omega$

  随后由位力定理：
$${\int }_{V}3PdV=3(\gamma -1)U=-\Omega$$
  得到：
$$-\Omega =\frac{3}{5-n}\frac{G{M}^2}{R}$$
  当$n=5$时，$-\Omega \to \infty$

库仑作用

  在之前的讨论中，没有考虑带电粒子间的库仑作用能。现在假设集体电中性，只考虑近邻的粒子：
$$E_{i,c}\sim \frac{(Z_{i}e{)}^2}{r_i}\sim (Z_{i}e{)}^2{n}_i^{1/3}=9.01\times 10^{-11}{Z}_i^2{\left(\frac{\rho X_i}{100A_i}\right)}^{1/3}erg$$
  相比之下，粒子的热运动动能为：
$$E_k\sim 3kT=6.2\times 10^{-9}(T/1.5\times 10^{7}K)erg$$
  在太阳的中心，$E_k\gg E_{i,c}$。而在低质量主序恒星中，库仑作用能重要

  进行更详细的讨论，考虑电子与离子统计分布。在原子核附近的电荷为$Q$，随后将其余电荷看成空间连续分布，则根据泊松方程有：
$${\nabla }^2\varphi (\vec{r})=-4\pi \sum _j{q}_j{n}_j(\vec{r})-4\pi Q\delta \overrightarrow{(}r)$$
  其中$j$代表电子或原子核

  而根据热平衡，有：
$$n_j(\vec{r})=\overline{n_j}\exp \left(-\frac{q_j\varphi (\vec{r})}{kT}\right)$$
  整体中性：
$$\sum _j{q}_j\overline{n_j}=0$$
  由弱场近似，$有q_j\varphi (\vec{r})\ll kT$，因此：
$$\exp \left(-\frac{q_j\varphi (\vec{r})}{kT}\right)\approx 1-\frac{q_j\varphi (\vec{r})}{kT}$$
  将$n_j(\vec{r})$代入到泊松方程中：
$${\nabla }^2\varphi (\vec{r})=4\pi (\sum _j\frac{q_j^2\overline{n_j}}{kT})\varphi (\vec{r})-4\pi Q\delta (\vec{r})$$
  德拜长度（等离子体屏蔽半径，原子核库仑影响半径）定义为：
$${\lambda }_D\equiv {\left(\frac{kT}{4\pi {\sum }_j{q}_j^2\overline{n_j}}\right)}^{\frac{1}{2}}$$
  由此将泊松方程的解写作：
$$\varphi (\vec{r})=\frac{Q}{r}\exp \left(-\frac{r}{{\lambda }_D}\right)$$
  原子核附近的势为：
$$\varphi (\vec{r})=\frac{Z_{i}e}{r}\exp \left(-\frac{r}{{\lambda }_D}\right)$$
  其他电荷产生的势为：
$${\varphi }_{other}(\vec{r})=\frac{Z_{i}e}{r}\left[\exp \left(-\frac{r}{{\lambda }_D}\right)-1\right]$$
  在离子位置，$\vec{r}\to 0$，电势趋近于$-\frac{Z_{i}e}{{\lambda }_D}$

  从而求得单位体积的平均库仑能：
$$\begin{array}{rl}{u}_{Coulomb}& \sim \frac{1}{2}\sum _i{Z}_{i}e\left(-\frac{Z_{i}e}{{\lambda }_D}\right)n_{0i}\\ & =-\frac{1}{2}\sum _i(Z_{i}e{)}^2{n}_{0i}{\left[\frac{4\pi }{kT}\sum _j(Z_{j}e{)}^2{n}_{0j}\right]}^{\frac{1}{2}}\\ & =-{\left(\frac{\pi }{kT}\right)}^{\frac{1}{2}}{\left[\sum _j(Z_{j}e{)}^2{n}_{0j}\right]}^{\frac{3}{2}}\end{array}$$
  对所有种类的粒子求和：
$$U_{Coulomb}=V{u}_{Coulomb}=\frac{B}{\sqrt{VT}}$$
  利用热力学第一定律，可以给出熵和压强：
$$\begin{gathered} S(V,T)=\frac{3}{2}Nk\ln T+\frac{1}{3}\frac{B}{\sqrt{VT}}\frac{1}{T}+g(V) \\ P=\frac{NkT}{V}+\frac{1}{3}\frac{B}{\sqrt{VT}}\frac{1}{V}=P_{Kinetic}+\frac{1}{3}{u}_{Coulomb} \end{gathered}$$
  其中，$g(V)$是积分常数

  对于纯氢而言：
$$P\sim 2\frac{\rho }{m_H}kT\left[1-0.0155{\left(\frac{\rho }{T^3}\right)}^{\frac{1}{2}}\right]$$
  对于纯氦而言：
$$P\sim \frac{3\rho }{4m_H}kT\left[1-0.038{\left(\frac{\rho }{T^3}\right)}^{\frac{1}{2}}\right]$$
  接下来，估计对原子能级的影响。写出薛定谔方程：
$$\begin{gathered} \hat{H}\psi =E\psi \\ \hat{H}=-\frac{h^2}{2m}{\nabla }^2-e\varphi (\vec{r}) \\ \varphi (r)=\frac{Z_{i}e}{r}\exp \left(-\frac{r}{{\lambda }_D}\right) \end{gathered}$$
  氢原子的波函数基态为（假设影响是小量）：
$$\begin{gathered} {\psi }_0(r)=\frac{1}{\sqrt{\pi a_0^3}}{e}^{-\frac{r}{a_0}} \\ \psi (r)=\frac{1}{\sqrt{\pi a^3}}{e}^{-\frac{r}{a}} \end{gathered}$$
  则能量为：
$$E=⟨\psi |H|\psi ⟩=-\frac{1}{\pi a^3}\int e^{-\frac{r}{a}}H{e}^{-\frac{r}{a}}4\pi r^{2}dr$$
  上式对$a$求极值，得到：
$$\begin{gathered} \frac{a_R}{a_0}\left(1+\frac{3a_R}{2{\lambda }_D}\right)={\left(1+\frac{a_R}{2{\lambda }_D}\right)}^3 \\ E(a_R)=-\frac{e^2}{2a_R}\left(1-\frac{a_R}{2{\lambda }_D}\right){\left(1+\frac{a_R}{2{\lambda }_D}\right)}^{-3} \end{gathered}$$
  其中，$a_R$是极小处的$a$。对于$\frac{a_R}{{\lambda }_D}$很小的情况，$a_R\to a_0$，则有：
$$E(a_R)\to -\frac{e^2}{2a_0}+\frac{e^2}{{\lambda }_D}$$
  电离能为：
$${\chi }_D=-E(a_R)\to \chi -\frac{e^2}{{\lambda }_D}=\chi \left(1-\frac{2a_0}{{\lambda }_D}\right)$$