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

第三章：物态方程

  在统计力学中，平衡体系，单位相空间体积内的粒子数由下述分布表示：
n ( p ) = 1 h 3 g j exp ⁡ { [ ϵ j ( p ⃗ ) − μ ] / k T } ∓ 1 n(p)=\frac{1}{h^3}\frac{g_j}{\exp\{[\epsilon_j(\vec p)-\mu]/kT\}\mp 1} n(p)=h31​exp{[ϵj​(p ​)−μ]/kT}∓1gj​​
  其中，当粒子为玻色子时，符号取负；粒子为费米子时，符号取正。$g_j$为能级简并度，$ϵ(\vec{p})$为粒子动能，$\mu ={\left(\frac{\partial E}{\partial N}\right)}_{S,V}$为化学势，是在不对外做功和熵不变时，每增加一个粒子的能量

物态方程

粒子数密度

  由粒子分布函数，可以积分得到粒子数密度（积分结果是单位体积内的粒子数）：
$$n={\int }_0^{\infty }n(p)4\pi p^{2}dp$$
  积分后得到$n=n(\mu ,T)$，或者写作$\mu =\mu (n,T)$，因此密度为$\rho =nm$

压强

  而对于压强：
$$\begin{array}{rl}P& =\iiint n(p)v2p{\cos }^{2}i{p}^{2}d\Omega dp\\ & =\frac{1}{3}{\int }_0^{\infty }n(p)vp4\pi p^{2}dp\end{array}$$
  右侧第一行是对于任意方向的动量进行的积分，第二行则是只对一个特定的方向进行的积分

  积分后的$P=P(\mu ,T)$，如果将$\mu$用$\rho$来表示，则$P=P(\rho ,T)$，是熟悉的物态方程

  粒子的能量与动量的关系为：
$$ϵ(p)=(p^2{c}^2+m^2{c}^4{)}^{1/2}-m{c}^2$$
  速度与能量的关系为：
$$v=\frac{dϵ}{dp}$$

内能密度

  内能密度$u$为：
$$u={\int }_0^{\infty }n(p)ϵ(p)4\pi p^{2}dp$$
  对于理想气体而言，在非相对论极限下，$ϵ(p)=\frac{p^2}{2m}$，而在相对论极限下，$ϵ(p)=pc$

  此积分结果的物理意义：单位体积的内能。积分后，内能$u=u(n,T)=u(\rho ,T)$

  而如果是对于非理想气体而言，还需要考虑内部能级的贡献

  对于其他的热力学量，可以通过热力学关系给出

热平衡辐射场

压强、内能与熵

  辐射场实际是光子场

  光子自旋为1，但质量为0，因此就只有两个偏振态，简并度$g=2$，化学势$\mu =0$（因为光子数不守恒），能量$ϵ=pc$（速度为光速），将这些值代入到前面的$n(p)$中后，得到：
$$n(p)=\frac{2}{h^3}\frac{1}{e^{pc/kT}-1}$$
  再根据前面的$P$和$u$的公式，积分得到：
$$\begin{array}{rl}& P=\frac{8{\pi }^5}{15}\frac{k^4}{c^3{h}^3}\frac{T^4}{3}=\frac{1}{3}a{T}^4\\ & u=a{T}^4\end{array}$$
  其中，$a=\frac{8{\pi }^5}{15}\frac{k^4}{c^3{h}^3}=7.566\times 10^{-15}ergc{m}^{-3}{K}^{-4}$

  由积分得到的$P$和$u$，可以将它们之间的关系写作：
$$P=(\gamma -1)u$$
  其中，$\gamma =\frac{4}{3}$

  通过热力学关系，还可以给出辐射场的熵：
$$\begin{array}{rl}Td{S}_r& =dU+PdV\\ Td{S}_r& =d(a{T}^{4}V)+\frac{1}{3}a{T}^{4}dV\\ d{S}_r& =\frac{4}{3}a{T}^{3}dV+Vd\left(\frac{4}{3}a{T}^3\right)\\ S_r& =\frac{4}{3}a{T}^{3}V\end{array}$$
  因为光子的化学势$\mu =0$，所以上述的积分常数为0（？？？？？？？？？？？？？？？）

黑体辐射强度

  根据定义，单位时间、单位面积、单位立体角、单位频率间隔的黑体的辐射强度为：
$$I_{\nu }dtdAd\Omega d\nu =h\nu n(p)p^{2}dpd\Omega cdtdA$$
  （为什么又有单位面积$dA$，又有单位立体角$d\Omega$？这不是矛盾了吗？还是说这是在相空间里面的，所以才会同时出现这两个？）——不是在相空间才会同时出现两个。这两个都是在正常三维空间里面的，$dA$是接收的面积，$d\Omega$是源在观测处看来的立体角大小

  由此得到：
$$I_{\nu }=ch\nu np{p}^2\frac{dp}{d\nu }=\frac{2h{\nu }^3}{c^2}\frac{1}{e^{h\nu /kT}-1}$$
  此即普朗克方程

理想气体

粒子数密度

  理想气体的分布是一种特殊情况下的分布

  对于非相对论性理想气体而言，粒子能量${ϵ}_i={ϵ}_0$，$ϵ(p)=\frac{p^2}{2m}$，$v=\frac{p}{m}$

  因此，对$p$积分即可得到粒子数密度：
$$\begin{array}{rl}n& =\frac{4\pi }{h^3}g{\int }_0^{\infty }{p}^2{e}^{\frac{\mu }{kT}}{e}^{-\frac{{ϵ}_0}{kT}}{e}^{-\frac{p^2}{2mkT}}dp\\ & =\frac{4\pi g}{h^3}\frac{\sqrt{\pi }}{4}(2mkT{)}^{3/2}{e}^{(\mu -{ϵ}_0)/kT}\end{array}$$
  如果使用经典统计，则应满足下述条件：
$$\begin{array}{rl}& e^{\frac{\mu -{ϵ}_0}{kT}}=\frac{n{h}^3}{g(2\pi mkT{)}^{\frac{3}{2}}}\ll 1\\ & \lambda =\frac{h}{p},⟨p⟩=\frac{1}{n}{\int }_0^{\infty }pdn=\frac{2}{\pi }\sqrt{2\pi mkT}\\ & ⟨n{\lambda }_m^3⟩=n{h}^3⟨\frac{1}{p^3}⟩\approx \frac{n{h}^3}{(2\pi mkT{)}^{\frac{3}{2}}}\ll 1\\ & n{T}^{-3/2}\ll \frac{(2\pi mk{)}^{3/2}}{h^3}\end{array}$$

压强、内能

  利用前面得到的$P$与$u$的计算方式，可以积分得到：
$$\begin{array}{rl}& P=nkT\\ & u=\frac{3}{2}nkT\end{array}$$
  以及压强与内能之间的关系：
$$P=(\gamma -1)u$$
  其中，$\gamma =\frac{5}{3}$

熵

  根据热力学第一定律（对每个粒子的平均值而言），有：
$$\begin{array}{rl}& Tds=du+PdV=\frac{3}{2}kdT+nkTd\left(\frac{1}{n}\right)\\ & ds=\frac{3}{2}k\frac{dT}{T}-k\frac{dn}{n}\\ & s=k\ln (\frac{T^{3/2}}{n})+s_0\end{array}$$
  常数$s_0$无法由经典给出，需要利用低密度下量子分布与经典分布的化学势相同这一关系给出：
$$\begin{array}{rl}& s_0=\frac{5}{2}k+\frac{3}{2}k\ln \left(\frac{2\pi mk}{h^2}\right)+k\ln 2\\ & s=k\ln \left(\frac{T^{3/2}}{n}\right)/m+\frac{s_0}{m}\end{array}$$

多成分理想气体

  设总粒子数为$n$，其中电子数为$n_e$，离子数为$n_I$。则对于$n_I$而言：
$$n_I=\sum _i{n}_i=\sum _i\frac{\rho }{m_H}\frac{X_i}{A_i}=\frac{\rho }{{\mu }_I{m}_H}$$
  其中，$X_i$为质量丰度，$A_i$为分子量，而对于${\mu }_I$而言：
$$\frac{1}{{\mu }_I}\equiv \sum _i\frac{X_i}{A_i}\approx X+\frac{1}{4}Y+\frac{1-X-Y}{⟨A⟩}$$
  其中$X$为氢占的质量，$Y$为氦占的质量，$⟨A⟩$为其他物质的平均分子量

  对于太阳而言，$X=0.707$，$Y=0.274$，$⟨A⟩\approx 20$，${\mu }_I=1.29$

  接下来看电子$n_e$，假设完全电离，则有：
$$n_e=\sum _i{Z}_i{n}_i=\frac{\rho }{m_H}\sum _i{X}_i\frac{Z_i}{A_i}=\frac{\rho }{m_H{\mu }_e}$$
  引入电子分子权重（$Z$为电子数）：
$$\frac{1}{{\mu }_e}\equiv \sum _i{X}_i\frac{Z_i}{A_i}=X+\frac{1}{2}Y+(1-X-Y)⟨\frac{Z}{A}⟩\approx \frac{1}{2}(1+X)$$
  由此，得到总的气体压强：
$$P=(n_I+n_e)kT=\frac{\rho }{\mu m_H}kT$$
  其中：
$$\frac{1}{\mu }=\frac{1}{{\mu }_e}+\frac{1}{{\mu }_I}$$
  对于太阳而言，$\mu =0.61$

  对于电离物质粒子之间的库仑作用，如果与粒子的热能相比很小，则可以忽略前者的作用，属于理想气体：
$$\begin{array}{rl}& \frac{Z{e}^2}{d}\sim Z{e}^2{n}^{1/3}\ll \frac{3}{2}kT\\ & \rho \ll 85{\left(\frac{T}{10^{6}K}\right)}^{3}gc{m}^{-3}\end{array}$$
  而如果密度远高于临界值，则晶格化（库仑力作用远大于热能，使得粒子固定）将会变得非常重要：
$$\frac{Z{e}^2}{d}\gg \frac{3}{2}kT$$

简并的费米子

  对于自旋半整数的粒子，${ϵ}_0=0$，简并度$g=2s+1=2$，一般的分布函数为Fermi-Dirac分布：
$$F(ϵ)=\frac{1}{\exp \left[\frac{ϵ-(\mu -m{c}^2)}{kT}\right]+1}$$
  其中：
$$\begin{array}{rl}& ϵ(p)=m{c}^2\left[\sqrt{1+{\left(\frac{p}{mc}\right)}^2}-1\right]\\ & v=\frac{\partial ϵ}{\partial p}=\frac{p}{m}{\left[1+{\left(\frac{p}{mc}\right)}^2\right]}^{-\frac{1}{2}}\end{array}$$

粒子数密度

  由此，可以计算粒子数密度：
$$n=\frac{2\times 4\pi }{h^3}{\int }_0^{\infty }\frac{p^{2}dp}{\exp \left[\frac{m{c}^2+ϵ(p)-\mu }{kT}\right]+1}$$
  这样的积分，在一般情况下是很困难的，因此先考虑特殊情况。在完全简并时，有$kT\ll (\mu -m{c}^2)$，此时的分布函数为：
$$F(ϵ)=\{\begin{array}{rl}& 1,ϵ\le (\mu -m{c}^2)\\ & 0,ϵ>(\mu -m{c}^2)\end{array}$$
  引入：
$$x=\frac{p}{mc},x_F=\frac{p_F}{mc}$$
  由此得到费米能：
$${ϵ}_F=m{c}^2[(1+x_F^2{)}^{\frac{1}{2}}-1]$$
  于是：
$$n=\frac{8\pi }{h^3}{\int }_0^{p_F}{p}^{2}dp=8\pi {\left(\frac{h}{mc}\right)}^{-3}{\int }_0^{x_F}{x}^{2}dx=\frac{8\pi }{3}{\left(\frac{h}{mc}\right)}^{-3}{x}_F^3$$
  由此可以计算电子密度：
$${\mu }_e=\frac{\rho }{n_e{m}_A}$$
  代入$n_e$，得到$\frac{\rho }{{\mu }_e}\approx 9.7\times 10^5{x}_F^{3}gc{m}^{-3}$（$m_A$是什么？）

压强

  于是，即可计算压强：
$$P=\frac{8\pi }{3}\frac{m^4{c}^2}{h^3}{\int }_0^{x_F}\frac{x^{4}dx}{(1+x^2{)}^{\frac{1}{2}}}=Af(x_F)$$
  其中：
$$\begin{array}{rl}& A=\frac{\pi }{3}{\left(\frac{h}{mc}\right)}^{-3}m{c}^2\\ & f(x)=x(2x^2-3)(1+x^2{)}^{\frac{1}{2}}+3{\sinh }^{-1}x\\ & {\sinh }^{-1}x=\ln (x+\sqrt{1+x^2})\end{array}$$

内能

  对应的内能为：
$$\begin{array}{rl}& u=Ag(x_F)\\ & g(x)=8x^3[(1+x^2{)}^{\frac{1}{2}}-1]-f(x)\end{array}$$
  在非相对论情况下：
$$x_F\to 0,f(x)=\frac{8}{5}{x}^5,g(x)=\frac{12}{5}{x}^5$$

$$\begin{array}{rl}{P}_e& =\frac{8\pi }{15}{\left(\frac{h}{m_{e}c}\right)}^{-3}{m}_e{c}^2{x}_F^5\\ & =\frac{1}{20}{\left(\frac{3}{\pi }\right)}^{1/3}{\left(\frac{h}{m_{e}c}\right)}^2{\left(\frac{1}{{\mu }_e{m}_A}\right)}^{5/3}{m}_e{c}^2{\rho }^{5/3}\\ & \approx 1.019\times 10^{18}{\left(\frac{\rho }{{\mu }_e}\right)}^{5/3}\end{array}$$

  （单位：cgs）
$$\begin{gathered} u=\frac{4\pi }{5}{\left(\frac{h}{mc}\right)}^{-3}m{c}^2{x}_F^5\propto {\rho }^{\frac{5}{3}} \\ u=\left(\frac{3}{2}\right)P \\ \gamma =\frac{5}{3} \end{gathered}$$
  在极相对论情况下：
$$\begin{gathered} x_F\to \infty ,f(x)\to 2x^4,g(x)\to 6x^4 \\ {P}_e=\frac{2\pi }{3}{\left(\frac{h}{m_{e}c}\right)}^{-3}{m}_e{c}^2{x}_F^4=\frac{1}{8}{\left(\frac{3}{\pi }\right)}^{1/3}\left(\frac{h}{m_{e}c}\right)m_e{c}^2{n}_e^{4/3}\approx 1.243\times 10^{19}{\left(\frac{\rho }{{\mu }_e}\right)}^{4/3} \end{gathered}$$
  （单位：cgs）
$$U_e=3P_e,\gamma =\frac{4}{3}$$

简并

  由热动能和费米能的大小进行比较，可以得知简并与非简并：

  若$kT\gg E_f$，则情况为非简并

  若$kT\ll E_f$，则情况为简并

  根据$\frac{\rho }{{\mu }_e}=9.7\times 10^5{x}_F^{3}gc{m}^{-3}$，可以得到：

  相对论性简并的条件：$\frac{\rho }{{\mu }_e}\gg 4.6\times 10^6{\left(\frac{T}{10^{10}k}\right)}^3$

  非相对论性非简并：${ϵ}_F=m{c}^2\frac{x_F^2}{2},\frac{\rho }{{\mu }_e}\ll 6.0{\left(\frac{T}{10^{6}k}\right)}^{\frac{3}{2}}gc{m}^{-3}$

  根据前面的积分公式，部分简并的物态-积分形式表示如下：
$$\begin{array}{rl}& \psi =\frac{\mu -m{c}^2}{kT}\\ & n=\frac{8\pi }{h^3}{\int }_0^{\infty }\frac{p^{2}dp}{1+e^{(ϵ/kT-\psi )}}\\ & P_e=\frac{8\pi }{3h^3}{\int }_0^{\infty }\frac{p^{3}v(p)dp}{1+e^{(ϵ/kT-\psi )}}\\ & U_e=\frac{8\pi }{h^3}{\int }_0^{\infty }\frac{p^2ϵ(p)dp}{1+e^{(ϵ/kT-\psi )}}\end{array}$$
  在非相对论情形下，表示如下：
$$\begin{array}{rl}& \eta =\frac{p^2}{2mkT}\\ & n=\frac{4\pi }{h^3}(2mkT{)}^{\frac{3}{2}}{\int }_0^{\infty }\frac{{\eta }^{1/2}d\eta }{1+e^{\eta -\psi }}\\ & P_e=\frac{4\pi }{3m{h}^3}(2mkT{)}^{\frac{5}{2}}{\int }_0^{\infty }\frac{{\eta }^{3/2}d\eta }{1+e^{\eta -\psi }}\\ & u_e=\frac{4\pi }{2m{h}^3}(2mkT{)}^{\frac{5}{2}}{\int }_0^{\infty }\frac{{\eta }^{3/2}d\eta }{1+e^{\eta -\psi }}=\frac{3}{2}P\end{array}$$
  在相对论情形下，表示如下：
$$\begin{array}{rl}& \xi =\frac{pc}{kT}\\ & n=\frac{8\pi }{h^3}{\left(\frac{kT}{c}\right)}^3{\int }_0^{\infty }\frac{{\xi }^{2}d\xi }{1+e^{\xi -\psi }}\\ & P_e=\frac{8\pi c}{3h^3}{\left(\frac{kT}{c}\right)}^4{\int }_0^{\infty }\frac{{\xi }^{3}d\xi }{1+e^{\xi -\psi }}\\ & u_e=\frac{8\pi c}{h^3}{\left(\frac{kT}{c}\right)}^4{\int }_0^{\infty }\frac{{\xi }^{3}d\xi }{1+e^{\xi -\psi }}=3P\end{array}$$
  对于弱简并的情况，由于$\exp (\psi )$是小量，因此可以对上式进行泰勒展开

  弱简并在小质量恒星主序阶段和一些演化晚期非常重要

  以极相对论性为例：
$$\frac{1}{1+e^{\xi -\psi }}=e^{\psi -\xi }[1-e^{\psi -\xi }+e^{2(\psi -\xi )}+\cdots ]$$

$$\begin{array}{rl}v& \equiv n_e\frac{h^3}{8\pi }{\left(\frac{kT}{c}\right)}^3={\int }_0^{\infty }{\xi }^2{e}^{\psi -\xi }(1-e^{\psi -\xi }+\cdots )d\xi \\ & =\Gamma (3)e^{\psi }-\frac{\Gamma (3)}{2^3}{e}^{2\psi }+\cdots \\ & =2e^{\psi }-\frac{1}{4}{e}^{2\psi }+\cdots \end{array}$$

  保留到二阶项后：
$$e^{\psi }\approx 4(1\pm \sqrt{1-v})$$

$$\begin{array}{rl}{P}_e& =\frac{8\pi c}{3h^3}{\left(\frac{kT}{c}\right)}^4{\int }_0^{\infty }{\xi }^3{e}^{\psi -\xi }(1-e^{\psi -\xi }+\cdots )d\xi \\ & =\frac{8\pi c}{3h^3}{\left(\frac{kT}{c}\right)}^4(6e^{\psi }-\frac{3}{8}{e}^{2\psi }+\cdots )\end{array}$$

具有内部结构的气体物态方程

  热平衡时，不同能级的分布服从玻尔兹曼分布，对于两个能级$i$、$j$，满足下述关系：
$$\frac{N_i}{N_j}=\frac{g_i}{g_j}{e}^{\frac{-({ϵ}_i-{ϵ}_j)}{kT}}$$
  取基态能量为0，则有：
$$N_i=N\frac{g_i{e}^{\frac{-{ϵ}_i}{kT}}}{Z}$$
  其中，$Z$为配分函数：
$$Z=\sum _j{g}_j{e}^{\frac{-{ϵ}_j}{kT}}$$
  把每一个粒子的能量加起来，得到内部的能级对内能的贡献：
$$U=N⟨ϵ⟩=\frac{N{\sum }_i{ϵ}_i{g}_i{e}^{-{ϵ}_i/kT}}{Z}$$
  一般不计内部结构对压强的影响

  在恒星内部，密度高，原子间距小，因为原子轨道半径随能级增大的速度很快，所以只要考虑$Z$的有限几个能级就可以

Saha方程——热平衡时电离程度与热力学量关系

  如果考虑到了电子离子的分布的话，就需要考虑电离程度了

  对于被电离出的电子，取原子基态的能量$E=0$，记电离势为${\chi }_I$，则自由电子能量为：
$$\Delta ϵ={\chi }_1+\frac{1}{2}{m}_e{v}^2$$
  则此时的热平衡玻尔兹曼分布写作：
$$\frac{d{n}_0^{+}(v)}{n_0}=\frac{g}{g_0}\exp \left[-\frac{{\chi }_I+\frac{1}{2}{m}_e{v}^2}{kT}\right]$$
  其中，$d{n}_0^{+}(v)$表示离子在基态，电子速度介于$v\sim v+dv$之间的态的粒子数，$g$表示离子在基态，电子自由态的权重$g=g_e{g}_0^{+}$

  $g_e$为电子态的统计权重：
$$g_e=\frac{2d^3\vec{x}{d}^3\overrightarrow{p_e}}{h^3}=\frac{2d^3\overrightarrow{p_e}}{h^3{n}_e}$$
  其中的常数2是因为电子的自旋有两个方向

  由于$d^3\overrightarrow{p_e}=4\pi p_e^{2}d{p}_e=4\pi m_e^3{v}_e^{2}d{v}_e$，因此：
$$\frac{d{n}_0^{+}(v)}{n_0}=\frac{8\pi m_e^3}{h^3}\frac{g_0^{+}}{n_e}\exp \left(-\frac{{\chi }_I+m_e{v}^2/2}{kT}\right)v^{2}dv$$
  对$v$积分，得到：
$$\frac{n_0^{+}{n}_e}{n_0}={\left(\frac{2\pi m_{e}kT}{h^2}\right)}^{3/2}\frac{2g_0^{+}}{g_0}\exp \left[-\frac{{\chi }_I}{kT}\right]$$
  代入到玻尔兹曼分布中，得到：
$$\frac{n_0}{n}=\frac{g}{Z(T)},\frac{n_0^{+}}{n^{+}}=\frac{g_0^{+}}{Z^{+}(T)}$$
  最后，得出Saha方程：
$$\frac{n^{+}{n}_e}{n}=\frac{2Z^{+}(T)}{Z(T)}{\left(\frac{2\pi m_{e}kT}{h^2}\right)}^{3/2}\exp \left[-\frac{{\chi }_I}{kT}\right]$$
  Saha方程用于描述相邻两个电离态$n$与$n^{+}$的粒子数密度关系（更高一级的电离态是什么？）

纯氢气体的电离

  纯氢只能电离为质子和电子，就不能再往后电离了，所以只有两种电离状态。电中性条件为：
$$\begin{gathered} n_e=n^{+} \\ n=n^{+}+n^0 \end{gathered}$$
  其中，$n_0$表示基态，$n^0$表示中性的，$n^{+}$表示电离的。第一条表示纯氢电离后电子和离子数目相等，第二条表示总数目等于电离的加上中性的

  引入电离度的定义$x=\frac{n^{+}}{n}$，由Saha方程得到（配分函数怎么算出来的基态离子:中性氢基态=1:2？）：
$$\begin{gathered} \frac{x^2}{1-x}=\frac{1}{n}{\left(\frac{2\pi m_{e}kT}{h^2}\right)}^{\frac{3}{2}}{e}^{\frac{-{\chi }_I}{kT}} \\ \frac{x^2}{1-x}=\frac{4\times 10^{-9}gc{m}^{-3}}{\rho }{T}^{\frac{3}{2}}{e}^{\frac{-1.578\times 10^5}{T}}\equiv \Lambda \\ x=\frac{-\Lambda \pm \sqrt{{\Lambda }^2+4\Lambda }}{2} \end{gathered}$$
  最后就得到了电离度的表达式

Saha方程的不适用范围

  在太阳冕区等地方，粒子密度低，热平衡假设不成立，因此Saha方程不适用

  而密度高时压强电离，电离程度高了之后，就主要跃迁移到紫外波段了。如果离子之间的距离小于波尔半径，电子就可以自由跑来跑去，和导电一样（金属氢）（这个情况叫做压强电离），Saha方程不适用

  温度低时，离子之间有库仑作用，库仑能比热能大很多

  例如，木星体积大，温度低，压强大，密度高，Saha方程不适用

部分电离物质的热力学性质

压强、内能、热容

  以纯氢为例，设电离度为$x$，则粒子数为$(1+x)n$。其中，$n$为氢原子的粒子数密度。由此得到压强为：
$$P=(1+x)nkT=(1+x)\rho kT/m_A$$
  单位体积内能为：
$$U=\frac{3}{2}(1+x)nkT+xn{\chi }_1$$
  单位质量内能为：
$$ϵ=\frac{3}{2}\frac{P}{\rho }+\frac{x}{m_H}{\chi }_I$$
  热容量为：
$$\begin{array}{rl}{c}_V& ={\left(\frac{\partial ϵ}{\partial T}\right)}_{\rho }\\ & =\frac{3}{2}(1+x)\frac{k}{m_H}+\frac{3}{2}\frac{k}{m_H}{\left(\frac{\partial x}{\partial T}\right)}_{\rho }+\frac{{\chi }_I}{k}{\left(\frac{\partial x}{\partial T}\right)}_{\rho }\\ & =\frac{3}{2}(1+x)\frac{k}{m_H}+\left(\frac{3}{2}\frac{k}{m_H}+\frac{{\chi }_I}{k}\right){\left(\frac{\partial x}{\partial T}\right)}_{\rho }\end{array}$$
  对Saha方程微分:
$${\left[\frac{\partial }{\partial T}\left(\frac{x^2}{1-x}\right)\right]}_n=\frac{1}{n}{\left(\frac{2\pi m_{e}kT}{h^2}\right)}^{\frac{3}{2}}{e}^{-\frac{{\chi }_I}{kT}}\left(\frac{3}{2T}+\frac{{\chi }_I}{k{T}^2}\right)$$
  从而得到：
$${\left(\frac{\partial x}{\partial T}\right)}_n=\frac{(1-x{)}^2}{2x-x^2}\frac{1}{n}{\left(\frac{2\pi m_{e}kT}{h^2}\right)}^{\frac{3}{2}}\left(\frac{3}{2T}+\frac{{\chi }_I}{k{T}^2}\right)e^{-\frac{{\chi }_I}{kT}}$$
  于是，可以得到：
$$c_V=c_V^0\left[1+\frac{2}{3}{\left(\frac{3}{2}+\frac{{\chi }_H}{kT}\right)}^2\frac{x(1-x)}{(1+x)(2-x)}\right]$$
  其中，$c_V^0=\frac{3}{2}(1+x)\frac{k}{m_H}$。类似地，也可以得到：
$$c_P=c_V^0\left[\frac{5}{3}+\frac{1}{3}{\left(\frac{5}{2}+\frac{{\chi }_I}{kT}\right)}^{2}x(1-x)\right]$$

绝热指数

  绝热时，$dQ=0$。定义${\Gamma }_1\equiv {\left(\frac{\partial \ln P}{\partial \ln \rho }\right)}_{ad}$，则有：
$$\frac{dP}{P}+{\Gamma }_1\frac{dV}{V}=0$$
  如果将$dQ$写作$dQ=dϵ+Pd\frac{1}{\rho }=c_{V}dT+\dots d\rho$，定义${\Gamma }_3-1\equiv {\left(\frac{\partial \ln T}{\partial \ln \rho }\right)}_{ad}$，则有：
$$\frac{dT}{T}+({\Gamma }_3-1)\frac{dV}{V}=0$$
  因为$\rho =\rho (P,T)$，所以$d\rho$可以写作$d\rho =\dots dP+c_{P}dT$，定义$\frac{{\Gamma }_2-1}{{\Gamma }_2}\equiv {\left(\frac{\partial \ln T}{\partial \ln P}\right)}_{ad}$，由此即有：
$$\frac{dT}{T}-\frac{{\Gamma }_2-1}{{\Gamma }_2}\frac{dP}{P}=0$$
  由此便得到了三种绝热指数：${\Gamma }_1$为绝热压缩系数的倒数，${\Gamma }_3-1$为绝热膨胀系数的倒数，$\frac{{\Gamma }_2-1}{{\Gamma }_2}$绝热温度梯度。因为有一个物态方程，所以三个指数中只有两个独立：
$$({\Gamma }_3-1)=({\Gamma }_2-1)\frac{{\Gamma }_1}{{\Gamma }_2}$$
  绝热指数与热力学量具有下述关系：
$$\begin{gathered} {\Gamma }_1=-\frac{c_P}{c_V}\frac{V}{P}{\left(\frac{\partial P}{\partial V}\right)}_T \\ {\Gamma }_3-1=\frac{V}{T}\frac{c_P-c_V}{c_V}{\left(\frac{\partial T}{\partial V}\right)}_P \\ \frac{{\Gamma }_2}{{\Gamma }_2-1}=-\frac{c_{P}T}{(c_P-c_V){\left(\frac{\partial T}{\partial P}\right)}_{V}P} \end{gathered}$$
  由此计算得到：
$$\begin{gathered} {\Gamma }_1(x)=\frac{5+{\left(\frac{5}{2}+\frac{{\chi }_I}{kT}\right)}^{2}x(1-x)}{3+\left[\frac{3}{2}+{\left(\frac{3}{2}+\frac{{\chi }_I}{kT}\right)}^2\right]x(1-x)} \\ \frac{{\Gamma }_2(x)}{{\Gamma }_2(x)-1}=\frac{5+2D(x)\left[\frac{{\chi }_I}{kT}+\left(\frac{3}{2}+\frac{{\chi }_I}{kT}\right)\left(\frac{5}{2}+\frac{{\chi }_I}{kT}\right)\right]}{2+D(x)\left(\frac{3}{2+\frac{{\chi }_I}{kT}}\right)} \\ {\Gamma }_3(x)-1=\frac{2+2D(x)\left(\frac{3}{2}+\frac{{\chi }_I}{kT}\right)}{3+2D(x){\left(\frac{3}{2}+\frac{{\chi }_I}{kT}\right)}^2} \\ D(x)=\frac{x(1-x)}{(2-x)(1+x)} \end{gathered}$$
  因为大部分的热量用于电离气体而不是升高温度，所以热容量大

  在部分电离区，绝热指数小，原因在于电离参数对温度很敏感

  部分电离区是恒星振动不稳定性的源

辐射和气体混合热力学性质

压强、内能

  在热平衡下：
$$\begin{gathered} P=nkT+\left(\frac{1}{3}\right)a{T}^4=\frac{\rho }{\mu m_A}kT+\frac{1}{3}a{T}^4 \\ u=\frac{3}{2}\frac{1}{\mu m_A}kT+\frac{a{T}^4}{\rho } \end{gathered}$$
  （$u$为单位质量的内能）

  在准静态变化时，由热力学第一定律，有：
$$\begin{gathered} dQ=du+PdV={\left(\frac{\partial u}{\partial T}\right)}_{V}dT+{\left(\frac{\partial u}{\partial V}\right)}_{T}dV+PdV \\ {\left(\frac{\partial u}{\partial T}\right)}_V=4a{T}^{3}V+\frac{3}{2}\frac{k}{\mu m_A}{\left(\frac{\partial u}{\partial V}\right)}_T=a{T}^4 \end{gathered}$$
  由此得到：
$$dQ=\left(4a{T}^{3}V+\frac{3}{2}\frac{k}{\mu m_A}\right)dT+\left(\frac{4}{3}a{T}^4+\frac{k}{\mu m_A}\frac{T}{V}\right)dV$$

物态方程、绝热系数

  对$P$进行微分，得到：
$$\begin{array}{rl}dP& =\left(\frac{4}{3}a{T}^4+\frac{kT}{\mu m_{A}V}\right)\frac{dT}{T}-\frac{kT}{\mu m_{A}V}\frac{dV}{V}\\ & =(4P_{rad}+P_{gas})\frac{dT}{T}-P_{gas}\frac{dV}{V}\end{array}$$
  引入$\beta$（关于气体贡献压强比例的参数）：
$$\begin{gathered} P_{gas}=\beta P \\ {P}_{rad}=(1-\beta )P \end{gathered}$$
  经过一系列的推导后，得到绝热系数：
$$\begin{array}{rl}{\Gamma }_1& =\frac{32-24\beta -3{\beta }^2}{24-21\beta }\\ {\Gamma }_2& =\frac{32-24\beta -3{\beta }^2}{24-18\beta -3{\beta }^2}\\ {\Gamma }_3& =\frac{32-27\beta }{24-21\beta }\end{array}$$