Rambler's Top100Astronet    
  по текстам   по ключевым словам   в глоссарии   по сайтам   перевод   по каталогу
 

На первую страницу << 3. Уравнение состояния при ... | Оглавление |

4. Вещество при очень больших плотностях, нейтронизация, взаимодействие частиц

При очень больших плотностях в условиях сильною вырождения электронов и нуклонов приближенно вещество можно считать холодным с \( T=S=0 \) [110].

а) Холодная нейтронизация вдоль состояний минимума энергии (СМЭ).

При плотностях \( \rho \lt 8.1\mbox{~г}\cdot \mbox{см}^{-3} \) вещество в СМЭ состоит из \( ^{56}Fe \), ядра которого максимально стабильны6. Когда \( \epsilon _{Fe} \) из (2.21) достигает значения \( \epsilon _{\beta } \) при

$$
%\begin{displaymath}
%\epsilon _{\beta }\left( A,Z\right) =\left( m_{A,Z-1}-m_{A,Z}\right) c^{2}-m_{e}c^{2}=B_{A,Z}-B_{A,Z-1}+\left( m_{n}-m_{p}\right) c^{2}-m_{e}c^{2} ,
%\end{displaymath}
\eqalign{ \epsilon_\beta(A,Z)&=(m_{A,Z-1}-\AZ m)c^2-\mec{}2 \cr &=\AZ B-B_{A,Z-1}+(m_{\mathrm{n}}-m_{\mathrm{p}})c^2-\mec{}2, \cr
}
$$ (4.1)

захват электрона стабильным ядром \( \left( A,Z\right) \) становится энергетически выгодным. Ядро \( \left( A,Z-1\right) \), которое в обычных условиях является \( \beta ^{-} \) радиоактивным, при большой энергии Ферми электронов оказывается устойчивым. Процесс захвата электрона ядром, называемый нейтронизацией, рассчитан впервые в [216].

При нулевой температуре состояние термодинамического равновесия соответствует минимуму полной энергии \( E_{tot} \) как функции \( A \) и \( Z \) при данном числе нуклонов в единице объема. С ростом плотности равновесие смещается в сторону все более переобогащенных нейтронами ядер. При \( \rho \geq \rho _{nd} \) энергия связи последнего нейтрона в ядре \( Q_{n} \) близка к нулю и в равновесии появляются свободные нейтроны. В равновесии плотность, при которой начинают отщепляться нейтроны, равна \( \rho _{nd}=4.3\cdot 10^{11}\mbox{~г}\cdot \mbox{см}^{-3} \) [267].

В отсутствие свободных нейтронов при \( \rho \lt \rho _{nd} \) для энергии единицы объема с учетом энергии покоя имеем [267]

\begin{displaymath}
\hbox{%
\begin{tabular}{l}
$\displaystyle{ E_{tot}\left( A,Z,n_{b}\right) =n_{A,Z}\left( m_{A,Z}c^{2}+W_{L}\right) +E_{e}\left( n_{e}\right) , }$ $\displaystyle{ n_{e}=Zn_{A,Z} , \quad n_{b}=An_{A,Z} . }$ \end{tabular}}
\end{displaymath} (4.2)

Массы стабильных ядер определяются с учетом энергии связи из табл. 4. Для ядер, далеких от области стабильности, экспериментальные измерения массы \( m_{A,Z} \) отсутствуют. В этом случае используются полуэмпирические (теоретические, подправленные имеющимися экспериментальными данными) формулы. Наиболее простая, отражающая все качественные закономерности формула для энергии связи, полученная Вейцзеккером, имеет вид [23]
\begin{displaymath}
B_{A,Z}=15.568-17.226A^{2/3}-0.698\frac{Z^{2}}{A^{1/3}}-23.279\frac{\left( A-2Z\right) ^{2}}{A}+\frac{34}{A^{3/4}}\delta \mbox{~МэВ} \end{displaymath} (4.3)

Здесь \( \delta =1 \) для четных \( A \) и \( Z \), \( \delta =0 \) для нечетных \( A \) и любых \( Z \), \( \delta =-1 \) для четных \( A \) и нечетных \( Z \). В (4.3) учтены статистические и капельные свойства ядер, энергия ядерного взаимодействия с учетом эффектов несимметрии (\( A\neq 2Z \)) и спаривания нуклонов, кулоновская и поверхностная энергии ядер. Расчет равновесного состава при холодной нейтронизации с использованием формулы Вейцзеккера сделан в [77]. В работе [267] для аналогичных расчетов использована более точная, но гораздо более сложная формула Майерса и Святецкого [486]. Энергия холодных электронов дана в (2.22), (2.23), а величина \( W_{L} \) есть энергия электростатического взаимодействия, возникающая из-за наличия точечных положительных зарядов в однородном фоне отрицательных. Минимум \( W_{L} \) достигается для объемно-центрированной кубической решетки и определяется формулой [267]
\begin{displaymath}
W_{L}=-1.819620Z^{2}e^{2}/b , \quad n_{A,Z}b^{3}=2 .
\end{displaymath} (4.4)

Электростатическое взаимодействие уменьшает энергию и давление вещества, так как расстояние между отталкивающимися ядрами в среднем больше, чем между притягивающимися зарядами разного знака [110]. Равновесный состав, соответствующий минимуму \( E_{tot} \) из (4.2) при заданной плотности барионов \( n_{b} \) приведен в табл. 6 из работы [267]. Отметим, что между плотностями \( \rho _{max} \) и \( \rho _{max}\left( 1+\frac{\Delta \rho }{\rho }\right) \) давление почти постоянно (слабо растет за счет давления ядер). Величина \( \rho =E_{tot}/c^{2} \) учитывает кинетическую энергию и энергию взаимодействия. Давление \( P \) к релятивистский показатель адиабаты \( \Gamma \) определяются формулами [267]
\begin{displaymath}
P=\left. n_{b}^{2}\frac{\partial \left( E_{tot}/n_{b}\right) }{\partial n_{b}}\right\vert _{S} , \quad
\Gamma=\left. \frac{n_{b}}{P}\frac{\partial P}{\partial n_{b}}\right\vert _{S}=\frac{\rho +P/c^{2}}{P}\left( \frac{\partial P}{\partial \rho }\right) _{S} .
\end{displaymath} (4.5)


Таблица 6. Равновесный ядерный состав при плотностях, меньших плотности испарения нейтронов (из [267])
Ядро \( B_{n},\mbox{~МэВ} \) Z/A \( \rho _{max},\mbox{~г}\cdot \mbox{см}^{-3} \) \( \mu _{te},\mbox{~МэВ} \) \( \Delta \rho /\rho ,\% \)
\( ^{56}Fe \) 8.7905 0.4643 8.1 (6) 0.95 2.9
\( ^{62}Ni \) 8.7947 0.4516 2.7 (8) 2.6 3.1
\( ^{64}Ni \) 8.7777 0.4375 1.2 (9) 4.2 7.9
\( ^{84}Se \) 8.6797 0.4048 8.2 (9) 7.7 3.5
\( ^{82}Ge \) 8.5964 0.3902 2.2 (10) 10.6 3.8
\( ^{80}Zn \) 8.4675 0.3750 4.8 (10) 13.6 4.1
\( ^{78}Ni \) 8.2873 0.3590 1.6 (11) 20.0 4.6
\( ^{76}Fe \) 7.9967 0.3421 1.8 (11) 20.2 2.2
\( ^{124}Mo \) 7.8577 0.3387 1.9 (11) 20.5 3.1
\( ^{122}Zr \) 7.6705 0.3279 2.7 (11) 22.9 3.3
\( ^{120}Sr \) 7.4522 0.3166 3.7 (11) 25.2 3.5
\( ^{118}Kr \) 7.2002 0.3051 (4.3 (11)) (26.2) ...
Здесь \( B_{n}=B_{A,Z}/A \) - энергия связи на нуклон, \( \rho _{max} \) - максимальная плотность, при которой существует данный нуклид, \( \mu _{te} \) - химический потенциал электронов при этой плотности, \( \Delta \rho /\rho \) - относительное увеличение плотности при переходе к следующему нуклиду. При \( \rho _{max}=4.3\cdot 10^{11}\mbox{~г}\cdot \mbox{см}^{-3} \) - начинается испарение нейтронов.

В таблице 7 даны результаты расчета термодинамических функций при равновесной нейтронизации. При \( \rho \lt 10^{4}\mbox{~г}\cdot \mbox{см}^{-3} \), когда энергия взаимодейстция \( W_{L} \) становится порядка кинетической энергии электронов, значения в табл. 7 взяты из работы [201].


Таблица 7. Уравнение состояния вещества в равновесии, соотетствующем минимуму полной энергии, из [267]
\( \rho ,\mbox{~г}\cdot \mbox{~см}^{-3} \) \( P,\mbox{~дин}\cdot \mbox{см}^{-2} \) \( n_{b},\mbox{~см}^{-3} \) Z A Г
7.86 \( \geq \)1.01 (9) 4.73 (24) 26 56 ...
7.90 1.01 (10) 4.76 (24) 26 56 ...
8.15 1.01 (11) 4.91 (24) 26 56 ...
11.6 1.21 (12) 6.99 (24) 26 56 ...
16.4 1.40 (13) 9.90 (24) 26 56 ...
45.1 1.70 (14) 2.72 (25) 26 56 ...
212 5.82 (15) 1.27 (26) 26 56 ...
1150 1.90 (17) 6.93 (26) 26 56 ...
1.044 (4) 9.744 (18) 6.295 (27) 26 56 1.796
2.622 (4) 4.968 (19) 1.581 (28) 26 56 1.744
6.587 (4) 2.431 (20) 3.972 (28) 26 56 1.706
1.654 (5) 1.151 (21) 9.976 (28) 26 56 1.670
4.156 (5) 5.266 (21) 2.506 (29) 26 56 1.631
1.044 (6) 2.318 (22) 6.294 (29) 26 56 1.586
2.622 (6) 9.755 (22) 1.581 (30) 26 56 1.534
6.588 (6) 3.911 (23) 3.972 (30) 26 56 1.482
8.293 (6) 5.259 (23) 5.000 (30) 28 62 1.471
1.655 (7) 1.435 (24) 9.976 (30) 28 62 1.437
3.302 (7) 3.833 (24) 1.990 (31) 28 62 1.408
6.589 (7) 1.006 (25) 3.972 (31) 28 62 1.386
1.315 (8) 2.604 (25) 7.924 (31) 28 62 1.369
2.624 (8) 6.676 (25) 1.581 (32) 28 62 1.357
3.304 (8) 8.738 (25) 1.990 (32) 28 64 1.355
5.237 (8) 1.629 (26) 3.155 (32) 28 64 1.350
8.301 (8) 3.029 (26) 5.000 (32) 28 64 1.346
1.045 (9) 4.129 (26) 6.294 (32) 28 64 1.344
1.316 (9) 5.036 (26) 7.924 (32) 34 84 1.343
1.657 (9) 6.860 (26) 9.976 (32) 34 84 1.342
2.626 (9) 1.272 (27) 1.581 (33) 34 84 1.340
4.164 (9) 2.356 (27) 2.506 (33) 34 84 1.338
6.601 (9) 4.362 (27) 3.972 (33) 34 84 1.337
8.312 (9) 5.662 (27) 5.000 (33) 32 82 1.336
1.046 (10) 7.702 (27) 6.294 (33) 32 82 1.336
1.318 (10) 1.048 (28) 7.924 (33) 32 82 1.336
1.659 (10) 1.425 (28) 9.976 (33) 32 82 1.335
2.090 (10) 1.938 (28) 1.256 (34) 32 82 1.335
2.631 (10) 2.503 (28) 1.581 (34) 30 80 1.335
3.313 (10) 3.404 (28) 1.990 (34) 30 80 1.335
4.172 (10) 4.628 (28) 2.506 (34) 30 80 1.334
5.254 (10) 5.949 (28) 3.155 (34) 28 78 1.334
6.617 (10) 8.089 (28) 3.972 (34) 28 78 1.334
8.332 (10) 1.100 (29) 5.000 (34) 28 78 1.334
1.049 (11) 1.495 (29) 6.294 (34) 28 78 1.334
1.322 (11) 2.033 (29) 7.924 (34) 28 78 1.334
1.664 (11) 2.597 (29) 9.976 (34) 26 76 1.334
1.844 (11) 2.892 (29) 1.105 (35) 42 124 1.334
2.096 (11) 3.290 (29) 1.256 (35) 40 122 1.334
2.640 (11) 4.473 (29) 1.581 (35) 40 122 1.334
3.325 (11) 5.816 (29) 1.990 (35) 38 120 1.334
4.188 (11) 7.538 (29) 2.506 (35) 36 118 1.334
4.299 (11) 7.805 (29) 2.572 (35) 36 118 1.334

б) Появление свободных нуклонов. СМЭ при субъядерных плотностях с учетом взаимодействия нуклонов.

При \( \rho \gt \rho _{nd} \) все связанные состояния нейтронов в ядрах оказываются заполненными и дальнейший рост плотности приводит к появлению свободных нейтронов. После различных попыток расчета уравнения состояния в этой области (см. обзор [266]), корректный подход к решению проблемы был развит в работе [265]. Этот подход основан на следующих принципах.

  1. Единое описание энергии взаимодействия нуклонов как внутри, так и вне ядер.

  2. Использование выражения для поверхностной энергии ядер, которое учитывает наличие окружающих нейтронов и обращается в нуль при идентичности вещества внутри и снаружи.

  3. Учет энергии кулоновского взаимодействия электронов и ядер решетки, а также протонов внутри ядра.

Если \( n_{A,Z} \) и \( n_{n} \) - концентрации ядер и свободных нейтронов в пространстве вне ядер, \( V_{A,Z} \) - объем ядра, а то выражение для полной энергии единицы объема \( E_{tot} \) запишется в виде


\begin{displaymath}
E_{tot}\left( A,Z,n_{n},n_{A,Z},V_{A,Z}\right) =n_{A,Z}\left( W_{A,Z}+W_{L}\right) +\left( 1-V_{A,Z}n_{A,Z}\right) E_{n}\left( n_{n}\right) +E_{e}\left( n_{e}\right) .
\end{displaymath} (4.6)

Для нахождения равновесного состава и уравнения состояния необходимо минимизировать \( E_{tot} \) относительно аргументов при постоянной концентрации нуклонов \( n_{b} \):


\begin{displaymath}
n_{b}=An_{A,Z}+\left( 1-V_{A,Z}n_{A,Z}\right) n_{n} .
\end{displaymath} (4.7)

Входящие в (4.6) кулоновская энергия решетки \( W_{L} \) и энергия холодных электронов определены в (4.4) и (2.22), (2.23) соответственно. Энергия ядра, представляемого в виде жидкой капли, записывается в виде


\begin{displaymath}
\eqalign{ W_{A,Z}\left( A,Z,V_{A,Z},n_{n}\right) =& \left[ \left( 1-x\right) m_{n}c^{2}+xm_{p}c^{2}+W\left( k,x\right) \right] A+ \cr & +W_{coul}\left( A,Z,V_{A,Z},n_{n}\right) +W_{surf}\left( A,Z,V_{A,Z},n_{n}\right), \quad
x=\frac{Z}{A}, \cr }
\end{displaymath} (4.8)

где \( W\left( k,x\right) \) - энергия на один барион в однородной ядерной материи с концентрацией \( n=k^{3}/1.5\pi ^{2}=A/V_{A,Z} \), \( W_{coul} \) - кулоновская энергия взаимодействия протонов в ядре, \( W_{surf} \) - поверхностная энергия ядра. Формула для энергии нейтронного газа \( E_{n} \) получается аналогично (4.8) при \( x=0 \)


\begin{displaymath}
E_{n}\left( n_{n}\right) =n_{n}W\left( k_{n},0\right) +m_{n}c^{2} , \quad
n_{n}=k_{n}^{3}/1.5\pi ^{2} .
\end{displaymath} (4.9)

Функция \( W\left( k,x\right) \), входящая в (4.8), (4.9), а также \( W_{coul} \) и \( W_{surf} \), входящие в (4.8), вычислены в [266]. Поверхностная энергия в [265] оценивалась довольно грубо из соображений размерности. Более точные выражения для \( W_{surf} \) были найдены в работах [313] методом Томаса-Ферми, [256] вариационным методом и [504548] методом Хартри-Фока. Величина - \( \left[ W\left( k,x\right) A+W_{coul}+W_{surf}\right] \) по физическому смыслу есть энергия связи ядра \( \left( A,Z\right) \) с учетом окружающей плотности нейтронов, которая для нормальных ядер в вакууме аппроксимируется формулой Вейцзеккера (4.3). Минимизация величины \( E_{tot} \) относительно своих аргументов при постоянстве \( n_{b} \) из (4.7) сводится к четырем уравнениям, определяющим равновесие по следующим процессам: 1) сжатие и расширение ядра давлением внешних нейтронов, 2) обмен нейтронами между ядром и нейтронным газом, 3) превращение нейтронов и протонов друг в друга внутри ядра, 4) выделение ядра, минимизирующего полную энергию. Решение этих четырех уравнений определяет свойства равновесного вещества, состоящего из ядер и свободных нейтронов. В табл.8 из [265] приведены результаты расчетов. С увеличением плотности происходит рост массы и объема ядер, пока они не касаются друг друга при \( \rho _{nm}=2.4\cdot 10^{14}\mbox{~г}\cdot \mbox{см}^{-3} \). В [265] предполагается, что при плотности \( \rho _{nn} \) происходит фазовый переход первого рода к однородной ядерной материи.


Таблица 8. Характеристики вещества в минимуме полной энергии при наличии свободных нейтронов (из [265])
\( \rho ,\mbox{~г}\cdot \mbox{см}^{-3} \) \( P,\mbox{~дин}\cdot \mbox{см}^{-3} \) \( n_{b},\mbox{~см}^{-3} \) Z A Г
4.460(11) 7.890(29) 2.,670(35) 40 126 0.40
5.228(11) 8.352(29) 3.126(35) 40 128 0.36
6.610(11) 9.098(29) 3.951(35) 40 130 0.40
7.964(11) 9.831(29) 4.759(35) 41 132 0.46
9.728(11) 1.083(30) 5.812(35) 41 135 0.54
1.196(12) 1.218(30) 7.143(35) 42 137 0.63
1.471(12) 1.399(30) 8.786(35) 42 140 0.73
1.805(12) 1.638(30) 1.077(36) 43 142 0.83
2.202(12) 1.950(30) 1.314(36) 43 146 0.93
2.930(12) 2.592(30) 1.748(36) 44 151 1.06
3.833(12) 3.506(30) 2.287(36) 45 156 1.17
4.933(12) 4.771(30) 2.942(36) 46 163 1.25
6.248(12) 6.481(30) 3.726(36) 48 170 1.31
7.801(12) 8.748(30) 4.650(36) 49 178 1.36
9.611(12) 1.170(31) 5.728(36) 50 186 1.39
1.246(13) 1.695(31) 7.424(36) 52 200 1.43
1.496(13) 2.209(31) 8.907(36) 54 211 1.44
1.778(13) 2.848(31) 1.059(37) 56 223 1.46
2.210(13) 3.931(31) 1.315(37) 58 241 1.47
2.988(13) 6.178(31) 1.777(37) 63 275 1.49
3.767(13) 8.774(31) 2.239(37) 67 311 1.51
5.081(13) 1.386(32) 3.017(37) 74 375 1.53
6.193(13) 1.882(32) 3.675(37) 79 435 1.54
7.732(13) 2.662(32) 4.585(37) 88 529 1.56
9.826(13) 3.897(32) 5.821(37) 100 683 1.60
1.262(14) 5.861(32) 7.468(37) 117 947 1.65
1.586 (14) 8.595 (32) 9.371 (37) 143 1390 1.70
2.004(14) 1.286(33) 1.182(38) 201 2500 1.74
2.520(14) 1.900(33) 1.484(38) ... ... 1.81
2.761(14) 2.242(33) 1.625(38) ... ... 1.82
3.085(14) 2.751(33) 1.814(38) ... ... 1.87
3.433(14) 3.369(33) 2.017(38) ... ... 1.92
3.885(14) 4.286(33) 2.280(38) ... ... 1.97
4.636(14) 6.103(33) 2.715(38) ... ... 2.03
5.094(14) 7.391(33) 2.979(38) ... ... 2.05

Уточнение формулы для поверхностной энергии сказалось только на изменении заряда ядра \( Z\left( \rho \right) \), рост которого замедлился по сравнению с [265]. На рис. 7 из [266] приведены результаты расчета \( Z\left( \rho \right) \) различными авторами. Уменьшение величины \( Z \) при больших плотностях приводит к более плавному изменению \( Z \) в процессе фазового перехода к однородной ядерной материи. Несовпадение зависимостей \( Z\left( \rho \right) \) у авторов [256, 313, 504, 548] отражает как недостаток физических знаний о законах взаимодействий между нуклонами, так и несовер шенство существующих математических методов расчета.

Рис. 7. Зависимость заряда $Z$ в коре нейтронной звезды от плотности, построенная по расчетам различных авторов: ВВР - [265], NV - [504], RBP - [548], ВВ - [313]

в) Плотность, больше ядерной.

Однородная ядерная материя (ЯМ), появляющаяся в результате фазового перехода, состоит из нейтронов с небольшой примесью протонов и электронов. Пока плотность ЯМ не превышает \( 2\rho _{0} \) ( \( \rho _{0}=2.8\cdot 10^{14}\mbox{~г}\cdot \mbox{см}^{-3} \) - плотность свободной ЯМ), для расчета уравнения состояния применим метод Бракнера-Бете-Голдстоуна, основанный на теории возмущений [22]. При \( \rho \gt 2\rho _{0} \) расчеты ведутся с помощью вариационного принципа, разработанного Пандарипанде [537]. При больших плотностях учитывается рождение тяжелых гиперонов [9] и возможное рождение \( \pi ^{-} \)-мезонов (пионная конденсация [158]). Расчеты уравнения состояния ядерной материи при \( \rho \geq \rho _{0} \) выполнены при различных допущениях в [275479]. В табл. 9 из [479] приводится уравнение состояния для наиболее реалистического варианта с параметрами потенциала взаимодействия, учитывающими экспериментальные данные из ядерной физики высоких энергий, а также рождение гиперонов. Возможное появление пионной конденсации слабо влияет на уравнение состояния [275]. Следует иметь в виду, что с ростом плотности все более возрастает неопределенность наших знаний о физике сильных взаимодействий и менее точными становятся методы расчета. В связи с этим даже в реалистическом варианте из табл. 9 погрешности могут достигать ~ 50%.

Как отмечалось впервые Я.Б. Зельдовичем [103110], требование принципа причинности о том, что скорость звука \( \nu _{s} \) не должна превышать скорость света \( c \), накладывает ограничение на уравнение состояния \( P\leq \epsilon =\rho c^{2} \). Важность этого ограничения связана с тем, что оно действует при сколь угодно больших плотностях, где о свойствах ядерных взаимодействий известно очень мало.


Таблица 9. Реалистические параметры вещества при больший плотности из [479]
\( n_{b},\mbox{~см}^{-3} \) \( \frac{E}{n_{b}}-m_{n}c^{2},\mbox{~МэВ} \) \( \rho ,\mbox{~г}\cdot \mbox{см}^{-3} \) \( P,\mbox{~дин}\cdot \mbox{см}^{-2} \)
1.0 (38) 12.6 1.70 (14) 1.19 (33)
1.5 (38) 16.6 2.55 (14) 2.93 (33)
2.0 (38) 21.2 3.42 (14) 6.00 (33)
2.5 (38) 26.0 4.31 (14) 1.09 (34)
3.0 (38) 32.2 5.21 (14) 1.83 (34)
4.0 (38) 46.9 7.04 (14) 4.09 (34)
5.0 (38) 64.4 8.95 (14) 7.61 (34)
6.0 (38) 83.7 1.09 (15) 1.26 (35)
7.0 (38) 109 1.31 (15) 1.99 (35)
8.0 (38) 134 1.54 (15) 2.85 (35)
9.0 (38) 160 1.76 (15) 3.71 (35)
1.0 (39) 189 2.01 (15) 4.02 (35)
1.1 (39) 215 2.26 (15) 5.02 (35)
1.25 (39) 254 2.66 (15) 6.76 (35)
1.4 (39) 295 3.08 (15) 8.81 (35)
1.5 (39) 324 3.37 (15) 1.03 (36)
1.7 (39) 383 4.00 (15) 1.38 (36)
2.0 (39) 475 5.04 (15) 2.02 (36)
2.5 (39) 639 7.02 (15) 3.40 (36)
3.0 (39) 814 9.36 (15) 5.20 (36)
При \( n_{b}\gt 1.0(39) \) имеют место асимптотические формулы
\( E=n_{b}\left( 15.05+3.03\left( \frac{n_{b}}{10^{39}}\right) ^{1.39}\right) \cdot 10^{-4} \mbox{эрг}\cdot \mbox{см}^{-3} \), \( P=403\left( n_{b}/10^{39}\right) ^{2.33}\cdot 10^{33} \mbox{дин}\cdot \mbox{см}^{-3} \),
\( \left( \rho =E/c^{2}\right) \).

г) Учет конечной температуры.

Учет температурных эффектов при плотностях, когда существенно ядерное взаимодействие, сделан для \( \rho \leq \rho _{0} \) с помощью обобщения методов, применявшихся при исследовании холодного вещества [465] (см. также [314]). В равновесии рассматривались протоны, нейтроны, ядра гелия, один тип тяжелых ядер, для которых находился минимум полной энергии. Ввиду свободного улета нейтрино и отсутствия равновесия по бета-процессам, задавалось отношение \( N_{p}/\left( N_{n}+N_{p}\right) =Y_{e} \). Выражения всех видов энергии, входящих в \( E_{tot} \) из (4.6) и из (4.8), записывались с учетом температурных поправок и к \( E_{tot} \) добавлялась энергия движения ядер, связанная с конечной температурой, а также энергии ядер гелия. К независимым величинам, которые являются аргументами функции \( E_{tot} \) в (4.6), при конечной температуре добавляются \( n_{p} \) -концентрация свободных протонов и \( n_{\alpha } \) - концентрация ядер гелия, которая определяется из условия равновесия относительно разбиения на протоны и нейтроны типа (3.3), но с учетом конечного объема (задаваемого) ядер гелия. Для определения \( n_{p} \) добавлялось условие равновесия по обмену протонами между ядрами и протонами в газе. На рис.8 из [456] представлены химические свойства при больших температурах и плотностях для \( Y_{e}=0.25 \). Интересно, что при больших плотностях \( \rho \gt 10^{14}\mbox{~г}\cdot \mbox{см}^{-3} \), когда ядра начинают занимать больше половины объема, вместо ядер, погруженных в менее плотный нуклонный газ, возникают шарики менее плотного вещества (пузыри), погруженные в более плотную ядерную материю. Важный вывод, следующий из данных расчетов, состоит в сохранении ядер до очень высоких температур \( \sim 20\mbox{~МэВ}\approx 2\cdot 10^{11}\mbox{~K} \), что является следствием учета ядерных взаимодействий и влияния окружающего газа на свойства ядер. В [456] отмечалось, что диаграмма на рис. 8 мало чувствительна к изменению \( Y_{e} \) в пределах \( 0.2\leq Y_{e}\leq 0.5 \). Учет различных типов тяжелых ядер, одновременно присутствующих в равновесии [346], также слабо влияет на полученные результаты [456]. На рис. 9 и 10 из [456] приведены изэнтропы вещества и зависимость показателя адиабаты \( \gamma _{1}\left( \rho \right) \) для этих изэнтроп. Видно, что существование ядер заметно уменьшает \( \gamma _{1} \) по сравнению с нуклонным газом при ненулевой энтропии. Как видно из рис. 8, при плотности, больше ядерной, вещество всегда однородно.

Рис. 8. Химические свойства вещества при высоких температурах и плотностях для $Y_e$ = 0,25 из [456]. Сплошная линия ограничивает область, где весовая доля ядер $X_H$ > 0,1, а в области внутри штриховой линии $X_H$ > 0,5. В области внутри штрих-пунктирной линии весовая доля гелия $X_H$ > 0,15. Тонкие сплошные линии задают массы ядер. Заштрихована область существования пузырей. Пунктирная линия определяет границу устойчивости однородной материи относительно разбиения на две фазы

Рис. 9. Адиабаты с указанными значениями безразмерной энтропии $s$ (в единицах $k/m_u$) для $Y_L$ = 0,25 (сплошные линии) и $Y_l$ = 0,35 (пунктирные линии), из [456]. Линии $X_H$ = 0,1 и $X_H$ = 0,5 имеют тот же смысл, что и на рис. 8. При $\rho$ > 1012 г/см3 вещество непрозрачно для нейтрино и сохраняется величина $Y_l$ - лептонный заряд на один барион. Указана также траектория центра звезды при коллапсе, полученная Арнеттом [255a]

Рис. 10. Зависимость показателя адиабаты $Г = \gamma_1 = (\partial\ln P/\partial\ln\rho)_S$ от плотности $\rho$ для указанных значений безразмерной энтропии $s$ (см. рис. 9) при $Y_l$ = 0,25 из [456]. В области штриховых линий кривые $\gamma(\rho)$ приведены сглаженными без некоторых несущественных деталей

Энергия Ферми идеального нейтронного газа, определяемая аналогично (2.21), есть


\begin{displaymath}
\eqalign{
& \epsilon _{Fn}=m_{n}c^{2}\left( \sqrt{1+y_{n}^{2}}-1\right)\,, \cr
& y_{n}=\frac{p_{Fn}}{m_{n}c}=\left( \frac{\rho _{n}}{6.2\cdot 10^{15}}\right)^{1/3}\,, \cr & p_{Fn}=\left( \frac{3\pi \rho _{n}}{m_{n}}\right) ^{1/3}\hbar\,. \cr }
\end{displaymath} (4.10)

При \( \rho \approx \rho _{n}\geq \rho _{0} \) имеем \( y_{n}\geq 0.36 \) и \( \epsilon _{Fn}\geq 61\mbox{~МэВ} \). Максимальные температуры, достигаемые при гравитационном коллапсе, обычно не превышают \( \sim 20\mbox{~МэВ} \), поэтому при \( \rho \geq \rho _{0} \) температурные поправки к уравнению состояния несущественны.

д) Неравновесная нейтронизация при увеличении плотности в холодном веществе.

Равновесие по ядерному составу достигается в веществе при высоких температурах, когда открыты все каналы реакций. По мере остывания, большинство каналов реакций закрывается и в холодном веществе достижение СМЭ, строго говоря, невозможно. Вещество при малых температурах всегда находится в неравновесном состоянии, однако степень не равновесности и характер ее зависят от пути, которым вещество пришло к состоянию с данными \( \rho \) и \( T \). В [6061287] рассмотрены два возможных пути: сжатие холодного вещества и остывание вещества при данной плотности с неравновесными составами, возникающими при этом.

Рассмотрим сжатие холодного вещества. Пусть при малых \( \rho \) вещество состоит из самого стабильного элемента \( ^{56}Fe \) и медленно сжимается при температуре, близкой к нулю. Когда плотность достигает величины \( 1.24\cdot 10^{9}\mbox{~г}\cdot \mbox{см}^{-3} \), \( \epsilon _{Fe}=3.81\mbox{~МэВ} \), становится энергетически выгодным захват электрона ядром \( ^{56}Fe \)7. Ввиду меньшей устойчивости нечетно-нечетных ядер, за первым захватом сразу следует второй и идет цепочка реакций

\begin{displaymath}
^{56}Fe+e^{-}\rightarrow ^{56}Mn+\nu _{e} , \quad ^{56}Mn+e^{-}\rightarrow ^{56}Cr+\nu _{e} .
\end{displaymath} (4.11)

Для второй реакции из (4.11) величина \( \epsilon _{\beta }=1.6\mbox{~МэВ} \) [99]. При \( \epsilon _{Fe}=3.81\mbox{~МэВ} \) она проходит неравновесно и сопровождается нагревом ([56], см. гл. 5), которым пренебрегаем ввиду того, что он не влияет на формирование химического состава. После образования \( ^{56}Cr \) в ходе дальнейшего сжатия становятся энергетически выгодными превращения \( ^{56}Cr\rightarrow ^{56}K\rightarrow ^{56}Ar \) и т.д., пока не образуется ядро с энергией отрыва последнего нейтрона, близкой к нулю \( Q_{n}\approx 0 \). После этого повышение плотности и захват электронов сопровождаются холодным испарением нейтронов из ядер и уменьшением \( A \) наряду с \( Z \). Ввиду различия в свойствах четных и нечетных ядер, часть испаряющихся нейтронов уносит энергию $\sim$1 МэВ/нейтрон, идущую на нагрев вещества [78].

В процессе неравновесной нейтронизации число ядер, приходящихся на один барион, не меняется. Если \( A_{0} \)- атомная масса начального ядра, то при данной плотности концентрации ядер и электронов (с учетом электронейтральности) запишутся в виде


\begin{displaymath}
n_{A,Z}=\frac{\rho }{A_{0}m_{u}} , \quad n_{e}=\frac{Z\rho }{A_{0}m_{u}} .
\end{displaymath} (4.12)

Здесь масса ядра приближенно принята равной \( Am_{u} \), а плотность \( \rho =Nm_{u} \), \( N=N_{n}+N_{p} \) полная концентрация барионов. До начала испарения нейтронов при \( A=A_{0} \) заряд ядра \( Z \) остается постоянным в интервале плотностей, где выполняется неравенство


\begin{displaymath}
\epsilon _{\beta }\left( A_{0},Z+2\right) \lt \epsilon _{Fe}\left( n_{e}\right) \lt \epsilon _{\beta }\left( A_{0},Z\right) .
\end{displaymath} (4.13)

Когда первое соотношение в (4.13) обращается в равенство, начинается очередной захват электрона при постоянных \( \epsilon _{Fe} \) и \( P_{e} \) в интервале плотностей \( \rho _{Z+2}\lt \rho \lt \rho _{Z} \):

\begin{displaymath}
\eqalign{
& \epsilon _{Fe}\left( n_{e}\right) =\epsilon _{\beta }\left( A_{0},Z+2\right) =\mbox{const}, \cr
& \frac{\rho _{Z}}{\rho _{Z+2}}=\frac{Z+2}{Z}\,. \cr }
\end{displaymath} (4.14)

Свойства вещества при холодной нейтронизации до начала испарения нейтронов приведены в табл. 10, где использовались ядерные данные из [379]. После начала испарения нейтронов концентрация их определяется соотношением


\begin{displaymath}
n_{n}=\frac{\rho }{m_{u}}\left( 1-\frac{A}{A_{0}}\right) , \quad
m_{n}\approx m_{p}\approx m_{u} .
\end{displaymath} (4.15)

Испарение нейтрона из ядра при наличии свободных нейтронов допустимо, если энергия вылетающих нейтронов превышает их энергию Ферми (4.10). Условие испарения после очередного захвата электрона в условиях (4.14) запишется в виде


\begin{displaymath}
Q_{n}\left( A,Z\right) =B_{A,Z}-B_{A-1,Z}\leq -\epsilon _{Fn}\left( n_{n}\right) .
\end{displaymath} (4.16)

Ввиду колебания функции \( Q_{n}(A,Z) \) в зависимости от четности \( A \) и \( Z \), протекает сразу цепочка реакций


\begin{displaymath}
\eqalign{
& \left( A,Z\right) +e^{-}\rightarrow \left( A,Z-1\right) +\nu _{e}\rightarrow \left( A-k_{1},Z-1\right) +k_{1}n+\nu _{e}\,, \cr & \left( A-k_{1},Z-1\right) +e^{-}\rightarrow \left( A-k_{1},Z-2\right) +\nu _{e}\rightarrow \left( A-k_{1}-k_{2},Z-2\right) +k_{2}n+\nu _{e}\,, \cr }
\end{displaymath} (4.17)

где \( k_{1}\approx k_{2}=3\div 4 \) для различных ядер [560]. Расчет неравновесной нейтронизации с испарением нейтронов по системе уравнений (4.12)-(4.16) проводился в [78] с энергией связи по формуле Вейцзеккера, а в [560] с учетом ядерного взаимодействия нуклонов и вычислении энергии связи по методике работы [265] при дополнительном учете влияния ядерных оболочек. Результаты этих расчетов приведены в табл. 10.


Таблица 10. Свойства вещества при холодной нейтронизации ядер железа
\( \rho, 10^{11}\mbox{~г}\cdot \mbox{см}^{-3} \) \( A \); \( Z \) \( \epsilon _{Fe},\mbox{~МэВ} \) \( \epsilon _{Fn},\mbox{~МэВ} \) \( P,10^{30}\mbox{~дин}\cdot \mbox{см}^{-2} \) \( P_{\mbox{п}},10^{30}\mbox{~дин}\cdot \mbox{см}^{-2} \)
0.012 56; 24 3.81 0 5.5 (-4) 0
0.14 56; 22 8.83 0 0.014 0
0.57 56; 20 14.0 0 0.080 0
1.6 56; 18 19.3 0 0.28 0
3.098 54; 18 23.84 0.069 0.569 3(-4) (6.5)
3.9 56; 16 25.2 0 0.73 0
5.01 56; 16 26.98 0 0.933  
6.17 54; 16 28.85 0.09 1.20  
6.233 48; 16 29.52 0.232 1.25 0.0059 (5.7)
7.24 46; 16 29.19 0.26 1.29  
7.664 42; 14 29.66 0.468 1.41 0.0344 (3.6)
9.689 36; 12 30.47 0.693 1.75 0.0921 (4.1)
10.1 56; 14 31.78 0 1.97 0
10.23 40; 12 31.00 0.44 1.66  
12.1 48; 12 32.06 0.43 2.07 0.029
12.88 35; 10 31.60 0.58 1.86  
14.13 32; 9 31.47 0.65 1.86  
14.88 30; 10 31.36 0.987 2.12 0.422 (4.5)
15.0 40; 10 32.35 0.79 2.25 0.13
16.8 36; 9 32.50 0.9 2.38 0.23
17.43 24; 8 32.47 1.41 2.48 0.537 (5.2)
19.2 32; 8 32.64 1.22 2.59 0.39
22.2 28; 7 32.78 1.49 2.88 0.65
25.65 18; 6 33.47 2.04 3.57 1.36 (6.0)
26.3 24; 6 32.92 1.82 3.35 1.07
44.12 12; 4 35.38 2.52 6.06 3.31 (7.8)
\( P \) - полное давление, \( P_{n} \) - давление нейтронов; слева в столбцах даны результаты из [78], справа - из [560], где $P_n$ не приводится. В середине даны результаты расчета с использованием таблиц [379135] до начала испарения нейтронов и формулы (4.20). В последнем столбце в скобках дана температура \( T(10^{9}K) \) при образовании данного состава за счет неравновесных \( \beta \)-захватов из [78].

В [560] отмечена важность пикноядерных реакций в ходе холодной не равновесной нейтронизации ввиду того, что скорость этих реакций быстро растет с ростом плотности (см. гл. 4). При \( \rho =1.4\cdot 10^{12}\mbox{~г}\cdot \mbox{см}^{-3} \) ядра с \( \left( A,Z\right) =\left( 32,9\right) \), образующиеся из ядер \( ^{56}Fe \), сливаются и после захвата двух электронов и испарения восьми нейтронов, образуются ядра с \( \left( A,Z\right) =\left( 56,16\right) \). В [560] не рассматривался тепловой эффект этой и последующих реакций слияния. Холодная неравновесная нейтронизация доведена до \( \rho =5.13\cdot 10^{13}\mbox{~г}\cdot \mbox{см}^{-3} \) с образованием ядра \( \left( 99,19\right) \) (ср. с ядром \( \left( 375,74\right) \) из табл. 8 или \( Z=30\div 40 \) из рис. 7 для равновесного состава при той же плотности). Не исключено, что выделение тепла при пикноядерных реакциях с последующим быстрым неравновесным захватом электронов и испарением нейтронов может достаточно повысить температуру для установления ядерного состава, близкого к равновесному.

Энергия отрыва последнего протона

\begin{displaymath}
Q_{p}\left( A,Z\right) =B_{A,Z}-B_{A-1,Z-1}=\epsilon _{\beta }\left( A,Z\right) +Q_{n}\left( A,Z-1\right) -\left( m_{n}-m_{p}\right) c^{2}+m_{e}c^{2} \end{displaymath} (4.18)

на границе \( Q_{n}=0 \) исследовалась П.Э. Немировским, оценки которого приведены на рис. 11. Грубо эта зависимость аппроксимируется формулой [61]
\begin{displaymath}
Q_{p}=\left( 33-\frac{Z}{7}\right) \mbox{~МэВ} , \quad A=4Z \mbox{~при~} Q_{n}=0 , \quad
Z\geq 6 \end{displaymath} (4.19)

Рис. 11. Энергия отрыва протона $Q_p$ в зависимости от $Z$ дня ядер, лежащих на границе существования с $Q_n$ = 0. Зависимость построена согласно количественным оценкам П.Э. Немировского

Пренебрегая энергией Ферми свободных нейтронов и считая, что нейтронизация с испарением нейтронов идет вдоль линии \( Q_{n}=0 \), получим из (4.14) с учетом (2.21), (4.12) и (4.19) соотношение для определения \( Z\left( \rho \right) \), \( A\left( \rho \right) \) :

\begin{displaymath}
\left[ 1+\left( \frac{Z_{\rho }}{10^{6}A_{0}}\right) ^{2/3}\right] ^{^{1/2}}=1.96\left( 33-\frac{Z}{7}\right) +2.53 ;
\quad A=4Z \end{displaymath} (4.20)

Здесь учтено, что \( m_{e}c^{2}=0.511\) МэВ, \( \left( m_{n}-m_{p}\right) c^{2}=1.293\) МэВ. Расчеты с использованием (4.20) также приведены в табл. 10, откуда видно, что различие в давлении \( P\left( \rho \right) \) для всех трех способов не превышает 20%. В то же время это давление более чем в полтора раза превышает равновесное давление при той же плотности (табл. 7).

е) Неравновесный состав при остывании горячего плотного вещества.

Когда температура в остывающем веществе станет меньше \( \left( 4-5\right) \cdot 10^{9}K \), реакции между заряженными частицами резко замедляются и концентрация ядер замораживается. В этих условиях возможно протекание реакций с нейтронами, фотоотщепление и захват нейтронов, \( e^{-} \) распады при \( \epsilon _{\beta }\gt \epsilon _{Fe} \) и \( e^{-} \)-захваты при \( \epsilon _{\beta }\lt \epsilon _{Fe} \). В условиях кинетического равновесия по бета-процессам (3.6) с \( \frac{dN_{n}}{dt}=0 \), в веществе имеется большой избыток свободных нейтронов [224]. При конечной температуре ядра могут присоединять нейтроны и отщеплять их, если

\begin{displaymath} -\epsilon _{Fn}\lt Q_{n}\lt Q_{nb}\approx 20kT-\epsilon _{Fn} .
\end{displaymath} (4.21)

В условиях избытка нейтронов \( x_{n}\gt 0.5 \) процесс формирования неравновесного химического состава при быстром остывании вещества представлен на рис. 12 из [61287]. Плоскость \( \left( A,Z\right) \) для ядер разбита на три области:

    I область с \( Q_{n}\gt Q_{nb} \),

    II область с \( -\epsilon _{Fn}\lt Q_{n}\lt Q_{nb} \), \( \epsilon _{\beta }\gt \epsilon _{Fe} \).

    III область с \( -\epsilon _{Fn}\lt Q_{n}\lt Q_{nb} \), \( \epsilon _{\beta }\lt \epsilon _{Fe} \).

Рис. 12. Образование химического состава при остывании на стадии ограниченного равновесия. Линия $Q_n = -\epsilon_{Fn}$ отделяет область существования ядер. Линия $Q_nb$, разделяет область I, где невозможно фотоотщепление нейтронов, и области II и III. Штриховые линии - уровни постоянного $\epsilon_\beta, \epsilon_\beta1 < \epsilon_\beta2 < \ldots < \epsilon_\beta,\max$. В области I $Q_n > Q_nb$, в области II $Q_n < Q_nb$, $\epsilon_Fn < \epsilon_\beta$, в области III $Q_n < Q_nb$, $\epsilon_Fn < \epsilon_\beta$. Линия со штриховкой справа отделяет область деления и альфа-распада. Заштрихованная область abed определяет границы значений (A, Z) при ограниченном равновесии с данными $Q_nb(T)$ и $\epsilon_Fn(\rho)$

При высокой концентрации нейтронов существующие ядра быстро перейдут из области I в области II и III из-за нейтронного захвата. В областях II и III имеется равновесие по отношению к захвату и отщеплению нейтронов. В области II бета-распады приводят к росту Z и ядра переходят в область выше линии \( cd \). В области III бета-захваты уменьшают \( Z \) и переводят ядра в область ниже линии \( ab \). Таким образом, в условиях ограниченного равновесия при \( \epsilon _{Fe}\gg Q_{nb}\gg kT \), когда температурные эффекты не влияют на бета-процессы, ядерный состав определяется узкой областью aecd на плоскости \( (A,Z) \) (рис. 12). Выход за пределы этой области не происходит из-за отсутствия допустимых бета-процессов и фотоотщеплений нейтронов. При \( T\leq 5\cdot 10^{8}K \) остается только одно ядро на границе (4.16).

Если пренебречь \( \epsilon _{Fn} \) и учесть (4.18) и (4.19), то ядерный состав можно приближенно найти из соотношений, аналогичных (4.20),

\begin{displaymath}
A=4Z , \quad Z=7\left\{ 33-0.511\left[ \left( \frac{\rho }{\mu _{Z}10^{6}}\right) ^{2/3}+1\right] ^{1/2}+1.293\right\} .
\end{displaymath} (4.22)

Зависимость \( Z(\rho ) \) при \( x_{n}=1/2 \) приведена на рис. 13. При \( \rho =\rho _{2}=2.24\cdot 10^{12}\mbox{~г}\cdot \mbox{см}^{-3} \) образуется ядро с \( Z=6 \), обладающее \( \epsilon _{\beta ,\max} \), а при \( \rho =\rho _{1}=1.2\cdot 10^{11}\mbox{~г}\cdot \mbox{см}^{-3} \) образуется ядро с \( Z=150 \), для которого \( Z^{2}/A=Z/4=37.5 \) и [164] время деления, зависящее от \( Z^{2}/A \), есть \( \tau _{f}=3\cdot 10^{7} \) лет. При \( \rho \lt \rho _{1} \) величина \( Z^{2}/A \) растет, а время деления уменьшается. Ядерное деление, а также альфа-распад приводят к росту числа зародышевых ядер. Таким образом, при \( \rho \gt \rho _{2} \) или \( \rho \lt \rho _{1} \) химический состав может приблизиться к равновесному в состоянии СМЭ. Существенная неравновесность при остывании возможна только при \( \rho _{1}\lt \rho \lt \rho _{2} \).

Рис. 13. Зависимость $Z(p)$ из (4.22) для неравновесного состава, образующегося при остывании плотного вещества, для $x_n$ = 1/2. Крестиками указаны примерные границы неравновесности

Отличие неравновесности, образующейся при остывании, от той, которая образуется при холодной нейтронизации, имеется, в основном, при меньших плотностях, где для горячего случая получаются очень большие \( Z \), см. табл. 10 и рис. 13. В то же время в обоих вариантах имеются принципиальные отличия от СМЭ, где \( Z(\rho ) \) увеличивается с ростом плотности, в то время как в обоих неравновесных составах \( Z(\rho ) \) быстро падает с ростом плотности.

ж) Термодинамические свойства вещества при учете кулоновского взаимодействия. Кулоновские взаимодействия в обычном веществе являются основными, определяющими агрегатное состояние, степени ионизации и поправки к термодинамическим функциям. В разреженных газах взаимодействия в основном локализованы внутри атомов и неполностью ободранных ионов. При нахождении степени ионизации по формуле Саха (1.8) происходит приближенный учет этого взаимодействия.

Основным параметром, характеризующим степень ионизации в разреженных газах, является величина \( \lambda _{ij}=I_{ij}/kT \) (см. 1). При уменьшении температуры и росте \( \lambda _{ij} \) газ становится нейтральным, а затем, когда тепловая энергия станет меньше энергии взаимодействия между атомами, превратится сначала в жидкое, а затем в твердое кристаллическое тело с упорядоченной структурой, обладающей минимумом энергии8. При сжатии вещества среднее расстояние между атомами \( l=n^{-1/3}_{A}=Z^{1/3}n^{-1/3}_{e} \) уменьшается9 и при \( \rho =\rho _{i1} \) становится порядка размера атома \( a_{Z}=Z^{-1/3}a_{0} \) в модели Томаса-Ферми (ТФ) [555], \( a_{0}=\hbar ^{2}/m_{e}e^{2} \) - Боровский радиус. При \( \rho =\rho _{i1} \) начинается ионизация давлением. Вещество становится полностью ионизованным при \( \rho =\rho _{i2} \). когда межэлектронные расстояния становятся порядка радиуса ближайшей к ядру электронной орбиты \( a_{Z0}=a_{0}Z^{-1} \) [71]. Плотность \( \rho _{i1} \) находится с учетом (2.21) из соотношения

\begin{displaymath}
\eqalign{
& \frac{l}{a_{Z}}=\frac{Z^{2/3}}{n^{1/3}_{e}a_{0}}=\left( 3\pi ^{2}\right) ^{1/3}\frac{Z^{2/3}\alpha }{y}=1\,, \cr
& n_{e}=\frac{y^{3}}{3\pi ^{2}}\frac{m^{3}_{e}c^{3}}{\hbar ^{3}}\,, \cr
& \rho _{i1}=\mu _{Z}m_{u}\frac{m^{3}_{e}c^{3}}{\hbar ^{3}}Z^{2}\alpha ^{3}=3\pi ^{2}10^{6}\mu _{Z}Z^{2}\alpha ^{3}=11.4\mu _{Z}Z^{2}\,, \cr }
\end{displaymath} (4.23)

\( \alpha =e^{2}/\hbar c=1/137 \) - постоянная тонкой структуры. Аналогично для плотности \( \rho _{i2} \) имеем
\begin{displaymath}
\eqalign{
& \frac{Z}{n^{1/3}_{e}a_{0}}=\left( 3\pi ^{2}\right) ^{1/3}\frac{Z\alpha }{y}=\theta _{i2}\,, \quad \theta _{i2}\geq 1\,, \cr
& \rho_{i2}=\frac{\mu _{Z}m_{u}}{\theta _{i2}^{3}}\frac{m^{3}_{e}c^{3}}{\hbar ^{3}}Z^{3}\alpha ^{3}=\frac{3\pi ^{2}}{\theta _{i2}^{3}}10^{6}\mu _{Z}Z^{3}\alpha ^{3}=\frac{11.4}{\theta _{i2}^{3}}\mu _{Z}Z^{3}\,. \cr }
\end{displaymath} (4.24)

Отметим, что при \( Z\lt \frac{\theta }{\alpha \left( 3\pi ^{2}\right) ^{1/2}}\approx 44\theta _{i2} \) полная ионизация давлением происходит при нерелятивистских электронах. Для железа \( ^{56}Fe \) имеем \( \rho _{i1}=1.7\cdot 10^{4}\mbox{~г}\cdot \mbox{см}^{-3} \), \( \rho _{i2}=4.4\cdot 10^{5}\theta _{i2}^{-3}\mbox{~г}\cdot \mbox{см}^{-3} \).

Рассмотрим фазовые превращения и кулоновские поправки к термодинамическим функциям в плотном, полностью ионизованном веществе. В виду малости параметра \( \beta _{Z}=\frac{Z^{2/3}}{n^{1/3}_{e}a_{0}}\left( \approx \frac{E_{c}}{Z\epsilon _{Fe}}\right) \) поправки при \( T=0 \) можно находить методом последовательных приближений [3,555]. Основную поправку к кинетической энергии электронов дает электростатическое взаимодействие ионов в решетке и свободных электронов между собой, а также обменное взаимодействие электронов. В первом приближении электроны можно считать расположенными однородно. Энергия электростатического взаимодействия \( E_{c} \) наиболее просто рассчитывается в приближении Вигнера-Зейца (\( WS \)), где учитывается только взаимодействие ионов и электронов внутри сферической ячейки радиуса $l_{WS}$

\begin{displaymath}
\frac{4\pi }{3}l^{3}_{WS}=n_{A,Z}^{-1} , \quad
l_{WS}=\left( \frac{3}{4\pi }\right) ^{1/3}l .
\end{displaymath} (4.25)

Энергия на одно ядро \( E_{c} \) в данном приближении есть [3555]


\begin{displaymath}
E_{c}=-\int\limits _{0}^{l_{WS}}\frac{Ze}{r}dq_{e}+\int\limits^{l_{WS}}_{0}\frac{q_{e}dq_{e}}{r}=-\frac{9}{10}\frac{(Ze)^{2}}{l_{WS}}=-\frac{9}{10}\left( \frac{4}{9\pi }\right) ^{1/3}Z^{5/3}\alpha m_{e}c^{2}y ,
\end{displaymath} (4.26)

где \( q_{e}=\frac{4\pi }{3}en_{e}r^{3} \) - суммарный электрический заряд электронов внутри радиуса \( r \), учтено также соотношение \( l_{WS}=\left( \frac{9\pi }{4}\right) ^{1/3}\frac{\hbar }{m_{e}c}\frac{Z^{1/3}}{y} \)10. В [555] учитывается поправка следующего порядка за счет неоднородности электронной плотности внутри ячейки в ТФ приближении. Эта поправка состоит из поправки к энергии Ферми \( E_{TF} \) и корреляционной поправки к электростатической энергии \( E_{cor} \):
\begin{displaymath}
\eqalign{
& E_{TF}=-\frac{162}{175}\left( \frac{4}{9\pi }\right) ^{2/3}\sqrt{1+y^{2}}Z^{7/3}\alpha ^{2}m_{e}c^{2}\,, \cr
& E_{cor}=\left[ 0.031\ln \left( \frac{e^{2}m_{e}}{\hbar ^{2}}l_{WS}\right) -0.048\right] Z\alpha ^{2}m_{e}c^{2}\,. \cr }
\end{displaymath} (4.27)

Последний член становится существенным при уменьшении плотности. Энергия обменного взаимодействия в общем случае есть [555]


\begin{displaymath}
E_{ex}=-\frac{3Z}{4\pi }\alpha m_{e}c^{2}y\varphi \left( y\right) ,
\end{displaymath}


\begin{displaymath}
\hbox{%
\begin{tabular}{l}
$\displaystyle{ \varphi \left( y\right) =\frac{1}{4y^{4}}\left[ \frac{9}{4}+3\left( \beta ^{2}-\frac{1}{\beta ^{2}}\right) \ln \beta -6\left( \ln \beta \right) ^{2}-\left( \beta ^{2}+\frac{1}{\beta ^{2}}\right) -\frac{1}{8}\left( \beta ^{4}+\frac{1}{\beta ^{4}}\right) \right] , }$ \\
$\displaystyle{ \beta =y+\sqrt{1+y^{2}} ; }$ \\
$\displaystyle{ \varphi \left( y\right) = \left\{ \begin{array}{cl} 1 & \quad \mbox{при~} y\ll 1 , \\ -1/2 & \quad \mbox{при~} y\gg 1 . \\ \end{array} \right. }$ \\
\end{tabular}}
\end{displaymath} (4.28)

В нерелятивистском пределе \( E_{ex} \) приведена в [3]. Таким образом, основные поправки к кинетической энергии холодного газа (2.22) за счет кулоновского взаимодействия для единицы массы равны с учетом (2.19)

\begin{displaymath}
\eqalign{ E_{q}=&\frac{1}{Am_{u}}\left( E_{c}+E_{TF}+E_{cor}+E_{ex}\right) = \cr & =\left\{ -\frac{3}{10\pi ^{2}}\left( \frac{4}{9\pi }\right) ^{1/3}\alpha Z^{2/3}y-\frac{54}{175\pi ^{2}}\left( \frac{4}{9\pi }\right) ^{2/3}\alpha ^{2}Z^{4/3}\sqrt{1+y^{2}}+ \right. \cr & \left. +\frac{0.031\alpha ^{2}}{3\pi ^{2}}\ln \left[ \left( \frac{9\pi }{4}\right) ^{1/3}\frac{\alpha Z^{1/3}}{y}\right] -\frac{0.016}{\pi ^{2}}\alpha ^{2}-\frac{\alpha }{4\pi ^{3}}y\varphi \left( y\right) \right\} \frac{m^{4}_{e}c^{5}}{\hbar ^{3}}\frac{y^{3}}{\rho }\,. \cr }
\end{displaymath} (4.29)

Из термодинамического соотношения \( P=\rho ^{2}\frac{dE}{dp} \), \( \frac{dy}{dp}=\frac{1}{3}\frac{y}{\rho } \) находим поправку к давлению

\begin{displaymath}
\eqalign{ P_{q}& =P_{c}+P_{TF}+P_{cor}+P_{ex}= \cr & =-\frac{m^{4}_{e}c^{5}}{\hbar ^{3}}y^{3}\left[ \frac{1}{10\pi ^{2}}\left( \frac{4}{9\pi }\right) ^{1/3}\alpha Z^{2/3}y+\frac{18}{175\pi ^{2}}\left( \frac{4}{9\pi }\right) ^{2/3}\alpha ^{2}Z^{4/3}\frac{y^{2}}{\sqrt{1+y^{2}}}-\frac{0.031}{9\pi ^{2}}\alpha ^{2}+\frac{\alpha }{4\pi ^{3}}\frac{\chi \left( y\right) }{y^{3}}\right]\,, \cr }
\end{displaymath} (4.30)

где

\begin{displaymath}
\hbox{%
\begin{tabular}{l}
$\displaystyle{ \chi \left( y\right) =\frac{y^{4}}{3}\frac{d}{dy}\left[ y\varphi \left( y\right) \right] =\frac{1}{32}\left( \beta ^{4}+\frac{1}{\beta ^{4}}\right) +\frac{1}{4}\left( \beta ^{2}+\frac{1}{\beta ^{2}}\right) -\frac{9}{16}-\frac{3}{4}\left( \beta ^{2}-\frac{1}{\beta ^{2}}\right) \ln \beta +\frac{3}{2}\left( \ln \beta \right) ^{2}- }$ \\
$\displaystyle{ \qquad\qquad -\frac{y}{3}\left( 1+\frac{y}{\sqrt{1+y^{2}}}\right) \left[ \frac{1}{8}\left( \beta ^{2}-\frac{1}{\beta ^{5}}\right) -\frac{1}{4}\left( \beta -\frac{1}{\beta ^{3}}\right) -\frac{3}{2}\left( \beta +\frac{1}{\beta ^{3}}\right) \ln \beta +\frac{3}{\beta }\ln \beta \right] , }$ \\[3mm]
$\displaystyle{ \chi \left( y\right) =\left\{ \begin{array}{cl} y^{4}/3 & \quad \mbox{при~} y\ll 1 , \\ -y^{4}/6 & \quad \mbox{при~} y\gg 1 . \\
\end{array}\right. }$ \\
\end{tabular}}
\end{displaymath}

Отсюда в предельных случаях имеем с учетом (2.32) и (2.33)

\begin{displaymath}
\eqalign{ \frac{P_{qr}}{P_{e}}& =-\frac{6}{5}\left( \frac{4}{9\pi }\right) ^{1/3}\alpha Z^{2/3}-\frac{216}{175}\left( \frac{4}{9\pi }\right) ^{2/3}\alpha ^{2}Z^{4/3}-\frac{0.124}{3}\frac{\alpha ^{2}}{y}+\alpha /2\pi = \cr & =-4.56\cdot 10^{-3}Z^{2/3}-1.78\cdot 10^{-5}Z^{4/3}-\frac{2.2\cdot 10^{-6}}{y}+1.16\cdot 10^{-3} \cr
}
\end{displaymath} (4.31)

для \( y\gg 1 \) [555].


\begin{displaymath}
\eqalign{ \frac{P_{qnr}}{P_{e}}& =-\frac{3}{2}\left( \frac{4}{9\pi }\right) ^{1/3}\alpha \frac{Z^{2/3}}{y}-\frac{54}{35}\left( \frac{4}{9\pi }\right) ^{2/3}\alpha ^{2}Z^{4/3}-\frac{0.155\alpha ^{2}}{3y^{2}}-\frac{5\alpha }{4\pi y}= \cr & =-5.7\cdot 10^{-3}\frac{Z^{2}}{y}-2.23\cdot 10^{-5}Z^{4/3}-\frac{2.75\cdot 10^{-6}}{y^{2}}-\frac{2.9\cdot 10^{-3}}{y} \cr }
\end{displaymath} (4.32)

для \( y\ll 1 \). Формулы (4.30), (4.32) применимы только для значений \( y \), при которых \( P_{q}/P_{e}\ll 1 \). С уменьшением плотности кулоновское взаимодействие качественно меняет вид уравнения состояния. В [647] приведено приближенное аналитическое соотношение для уравнения состояния:
\begin{displaymath}
\rho =\rho _{0}(\zeta +\varphi )^{3} , \quad \zeta ^{5}\equiv P/P_{0} ,
\end{displaymath} (4.33)

где
\begin{displaymath}
\rho _{0}=\frac{32}{3}\pi ^{-3}AZm_{u}a_{0}^{-3}=3.88ZA\cdot \mbox{г}\cdot \mbox{см}^{-3} ,
\end{displaymath} (4.34)


\begin{displaymath}
P_{0}=Z^{10/3}\frac{2^{8}\left( 2\pi \right) ^{1/3}}{15\pi ^{4}}\left( \frac{e^{2}}{\hbar c}\right) ^{2}\frac{m_{e}c^{2}}{a_{0}^{3}}=9.52\cdot 10^{13}Z^{10/3} \mbox{дин}\cdot \mbox{см}^{-2} \end{displaymath}


\begin{displaymath}
\varphi =\frac{1}{20}3^{1/3}+\frac{1}{8}\left( \frac{3}{4}\pi ^{-2}Z^{-2}\right) ^{1/3} ,
\quad a_{0}=\hbar ^{2}/m_{e}e^{2} .
\end{displaymath}

при \( \rho \gg \rho _{0} \), \( \zeta \gg \varphi \) уравнение состояния (4.33) сводится к (39.3).

В силу квантовых свойств вещества ионы в решетке при абсолютном нуле совершают колебания с частотой, определяемой взаимодействием с электронами. В \( WS \) приближении возвращающая сила, действующая на ион со стороны электронов в ячейке при отклонении от равновесия, есть

\begin{displaymath}
F\left( r\right) =-\frac{Zeq_{e}}{r^{2}}=\frac{4\pi }{3}Ze^{2}n_{e}r ,
\end{displaymath}

что приводит к частоте гармонических колебаний [125555]
\begin{displaymath}
\omega ^{2}_{i}=\frac{F(r)}{Am_{u}r}=\frac{4\pi }{3}\frac{Ze^{2}n_{e}}{Am_{u}}=\frac{4}{9\pi }\frac{Zy^{3}}{Am_{u}}\alpha \frac{m_{e}^{3}c^{4}}{\hbar ^{2}} .
\end{displaymath} (4.35)

Нулевую энергию трехмерного осцилятора \( E_{zp}=\frac{3}{2}\hbar \omega _{i} \) следует сравнивать с \( E_{c} \) из (4.26). Имеем
\begin{displaymath}
f=\frac{E_{zp}}{E_{c}}=\frac{5}{3}\left( \frac{4}{9\pi }\right) ^{1/6}\left( \frac{m_{e}}{m_{u}}\right) ^{1/2}\left( \frac{y}{\alpha AZ^{7/3}}\right) ^{1/2}=\left( 0.11\frac{y}{AZ^{7/3}}\right) ^{1/2} .
\end{displaymath} (4.36)

При \( f\geq 1 \) кулоновский кристалл разрушается уже при нулевой температуре за счет квантовых колебаний. При этом необходимы слишком большие плотности \( \rho \gt 7.5\cdot 10^{8}\mu _{Z}A^{3}Z^{7} \), так что в холодных звездах нулевые колебания не разрушают кристаллической структуры.

Более реальным является разрушение кристалла тепловыми движениями ионов [485620]. Плавление кристалла происходит при \( T=T_{m} \), когда кинетическая энергия колебаний \( kT \) составляет ~ 1/150 от кулоновской энергии \( WS \) ячейки [544245 578]11. Имеем с учетом (4.26)

\begin{displaymath}
\eqalign{
& kT_{m}=0.003\alpha m_{e}c^{2}Z^{5/3}y\,, \cr
& T_{m}=1.3\cdot 10^{5}Z^{5/3}\left( \frac{\rho }{\mu _{Z}10^{6}}\right) ^{1/3}K\,. \cr
}
\end{displaymath} (4.37)

Термодинамические свойства кристаллических тел хорошо известны [145]. При малых температурах возбуждаются только степени свободы, соответствующие длинным волнам (малым частотам). Эти моды колебаний (фононы) обладают свойствами, аналогичными фотонному газу. Напротив, при больших температурах возбуждаются все возможные моды колебаний. При этом энергия колебаний кристалла в два раза больше кинетической энергии вещества в газовом состоянии при той же температуре. В общем случае для тепловой энергии ионной решетки на единицу массы справедлива интерполяционная формула Дебая [145]

\begin{displaymath}
\eqalign{
& E_{iT}=\frac{3kT}{Am_{u}}{\cal D}\left( \frac{\theta }{T}\right)\,, \cr & {\cal D}\left( x\right) =\frac{3}{x^{3}}\int _{0}^{x}\frac{z^{3}dz}{e^{z}-1} \quad \mbox{--- функция Дебая,} \cr & \theta =0.775\frac{\hbar \omega _{i}}{k}=\frac{3.5\cdot 10^{3}\sqrt{\rho }}{\mu _{Z}}K\,, \cr }
\end{displaymath} (4.38)

\( \theta \) - дебаевская температура кулоновской решетки [145620]. Здесь \( \omega _{i} \) - частота ионных колебаний в кристалле (4.35). В предельных случаях функция \( {\cal D}\left( x\right) \) равна [145]12
\begin{displaymath}
{\cal D}\left( x\right) =\left\{ \eqalignleft{ \frac{\pi ^{4}}{5x^{3}}-3e^{-x}\,, \quad & \mbox{~при~} x\gg 1\,, \cr 1-\frac{3}{8}x+\frac{x^{2}}{20}\,, \quad & \mbox{~при~} x\ll 1\,. \cr
}\right. \end{displaymath} (4.39)

Энтропия единицы массы и давление, связанные с ионным кристаллом, равны [145]

\begin{displaymath}
S_{iT}=\frac{k}{Am_{u}}\left[ -3\ln \left( 1-e^{-\frac{\theta }{T}}\right) +4{\cal D}\left( \frac{\theta }{T}\right) \right] ,
\end{displaymath}


\begin{displaymath}
P_{iT}=\frac{3}{2}\frac{\rho kT}{Am_{u}}{\cal D}\left( \frac{\theta }{T}\right) =\frac{\rho E_{iT}}{2} .
\end{displaymath} (4.40)

В предельных случаях
\begin{displaymath}
S_{iT}=\left\{ \begin{array}{ll} \displaystyle{ \frac{3k}{Am_{u}}\left( \ln \frac{\theta }{T}+\frac{4}{3}\right) , } & \mbox{~при~} T\gg \theta , \\ \displaystyle{ \frac{4}{5}\frac{\pi ^{4}}{Am_{u}}\frac{kT^{3}}{\theta ^{3}} , } & \mbox{~при~} T\ll \theta . \\
\end{array}\right. \end{displaymath} (4.41)

В зависимости от соотношения между температурами \( T \), \( T_{m} \) и \( \theta \) имеем следующие состояния ионизованного вещества:

\( T\lt \theta \lt T_{m} \) \( \theta \lt T\lt T_{m} \) \( \theta \lt T_{m}\lt T \) \( T_{m}\lt \theta \lt T \) \( T\lt T_{m}\lt \theta \) \( T_{m}\lt T\lt \theta \)
Квантовый (вырожденный) кристалл Классический ионный кристалл Классическая жидкость (неидеальная плазма) Квантовая жидкость

Плавление ионного кристалла при \( T\gt T_{m}\gt \theta \) сопровождается выделением тепла, а при дальнейшем росте температуры происходит переход к идеальному газу с \( E_{iT}=\frac{3}{2}kT/Am_{u} \). Согласно [544245578] приближение идеального газа становится применимым при:

\begin{displaymath}
T_{g}\approx 150T_{m}\approx 2\cdot 10^{7}Z^{5/3}\left( \rho /\mu _{Z}10^{6}\right) ^{1/3}K .
\end{displaymath} (4.42)

В промежутке \( T_{m}\lt T\lt 150T_{m} \) можно использовать интерполяцию. Вопрос о теплоте плавления не вполне ясен. В [485] приводятся аргументы в пользу очень малой теплоты плавления при постоянных \( T \) и \( \rho \). В [620] из других соображений предполагается, что фазовый переход относится к первому роду и теплота плавления решетки при постоянных \( T \) и \( \rho \) есть

\begin{displaymath}
\delta U_{coul}\approx -\frac{3}{4}kT_{m} .
\end{displaymath} (4.43)

Строгого решения этой проблемы не сделано.

Наиболее сложен для количественного описания промежуточный интервал плотностей \( \rho _{i1}\lt \rho \lt \rho _{i2} \), где степень ионизации определяется как тепловой энергией, так и давлением. В ТФ приближении процесс ионизации холодного вещества давлением исследован в [126]. Оболочечные поправки приводят к скачкам давления при каждой последующей ионизации, при плотности

\begin{displaymath}
\eqalign{
& \rho _{n}=\left( \frac{k}{\pi n}\right) ^{6}\mu _{Z}m_{u}\frac{m^{3}_{e}c^{3}}{\hbar ^{3}}Z^{4}\alpha ^{3}\,, \cr n=1,2,\ldots\,, \quad k=3\left( \frac{3}{4\pi }\right) ^{1/6}\ln \frac{1+\sqrt{3}}{\sqrt{2}}\approx 1.56\,, \cr
}
\end{displaymath} (4.44)

Зависимость немонотонной поправки к давлению от плотности задается формулой [126]

\begin{displaymath}
\delta P=-12\left( \frac{9\pi }{4}\right) ^{1/6}k\frac{Z^{5/3}\upsilon ^{4/3}}{\ln \left( Z^{3}\upsilon \right) }\left[ S_{0}\right] \frac{m^{4}_{e}c^{5}}{\hbar ^{3}}\alpha ^{5} ,
\end{displaymath} (4.45)

где
\begin{displaymath}
\upsilon =\frac{Am_{u}}{\rho }\alpha \frac{m^{3}_{e}c^{4}}{\hbar ^{3}} , \quad
S_{0}=kZ^{1/2}\upsilon ^{1/6}-\pi /2 ,
\end{displaymath} (4.46)

\( \left[ f\left( S_{0}\right) \right] \) при \( -\pi /2\lt S_{0}\lt \frac{\pi }{2} \) с периодическим продолжением вне этой области. Скачки давления (4.45), возникающие при \( \rho =\rho _{n} \) и \( \upsilon =\upsilon _{n} \), соответствуют фазовым переходам первого рода. В окрестности \( \rho =\rho _{n} \) давление электронов постоянно, причем величина \( P_{n} \) находится по правилу Максвелла (равенство площадей, рис. 14). Оболочечные эффекты исчезают, согласно (4.44), при
\begin{displaymath}
\rho \gt \rho _{1}=0.015\mu _{Z}m_{u}\frac{m^{3}_{e}c^{3}}{\hbar ^{3}}\alpha ^{3}Z^{4} ,
\end{displaymath} (4.47)

когда все уровни энергии уходят в непрерывный спектр. По физическому смыслу \( \rho \leq \rho _{i2} \) из (4.24). Если принять, что при \( \rho =\rho _{i2} \) электроны становятся релятивистскими с \( y=1 \), то получаем \( \theta _{i2}=1.11 \) и условие \( Z\leq 49 \) для применимости нерелятивистского ТФ-рассмотрения процесса ионизации [126].

Рис. 14. Схематическая зависимость давления от удельного объема при учете оболочечных эффектов. Постоянные уровни давления $Р_n$ находятся из равенства заштрихованных площадей

Термодинамические свойства смеси водорода и гелия при конечной температуре с учетом взаимодействия и ионизации давлением изучались в [403]. Учитывались соединения \( H_{2} \), \( H \), \( H^{+} \), \( H^{-} \), \( He \), \( He^{+} \), \( He^{++} \). Для данной температуры \( T \) и полного числа ионов \( N_{i} \) в объеме \( V \) вычисляется свободная энергия \( F \) в виде

\begin{displaymath}
F=F_{0}+F_{c} ,
\end{displaymath} (4.48)

где \( F_{0} \) соответствует идеальному газу, a \( F_{c} \) отражает кулоновское взаимодействие. Давление \( P \), энтропия \( S \) и энергия системы \( E \) находятся из соотношений
\begin{displaymath}
\eqalign{
& P=-\left( \frac{\partial F}{\partial V}\right) _{T,N_{i}}\,, \cr
& S=-\left( \frac{\partial F}{\partial T}\right) _{V,N_{i}}\,, \cr
& E=F+TS\,. \cr
}
\end{displaymath} (4.49)

Кулоновское взаимодействие учитывается здесь также в виде поправок к идеальному газу

\begin{displaymath}
P=P_{0}+P_{c}\,, \quad S=S_{0}+S_{c}\,, \quad E=E_{0}+E_{c}\,.
\end{displaymath} (4.50)

На основе ТФ теории с дебай-хьюккелевским потенциалом вокруг ядра (ТФДХ) имеем [403]

\begin{displaymath}
\eqalign{
F_{c}=-\frac{kTV}{12\pi }k_{{\cal D}}^{3}\,, \cr
& k^{2}_{{\cal D}}=\frac{4\pi e^{2}}{kT}\left( n_{i}+n_{e}\theta _{e}\right)\,, \quad n_{i}=\frac{N_{i}}{V}\,, \quad n_{e}=\frac{N_{e}}{V}\,, \cr
& \theta _{e}=F_{1/2}^{'}\left( \eta \right) /F_{1/2}\left( \eta \right)\,. \cr }
\end{displaymath} (4.51)

где \( F_{\alpha }(\eta ) \) - функции Ферми из (2.49), \( \eta kT \) - химический потенциал свободных электронов, \( k^{-1}_{{\cal D}} \) - дебаевский радиус (см. также (8.47)), и предполагается, что все ионы имеют одинаковый средний заряд \( Z=1 \). Для рассмотренной смеси с \( x_{H}=0.739 \), \( x_{He}=0.261 \) величина \( Z \) меняется от единицы до 1.0811 при изменении состава от полностью нейтрального до полностью ионизованного, \( n_{e} \) - концентрация свободных электронов. С учетом (4.51), поправки к термодинамическим функциям равны
\begin{displaymath}
\eqalign{
& P_{c}=-F_{c}\left( \frac{1}{V}+\frac{3}{k_{\cal D}}\cdot \frac{\partial k_{\cal D}}{\partial V}\right)\,, \cr & S_{c}=-F_{c}\left( \frac{1}{T}+\frac{3}{k_{\cal D}}\cdot \frac{\partial k_{\cal D}}{\partial T}\right)\,, \cr & E_{c}=\frac{kT^{2}}{4\pi }Vk^{2}_{\cal D}\frac{\partial k_{\cal D}}{\partial T}\,. \cr
}
\end{displaymath} (4.52)

Результаты ТФДХ рассмотрения верны при

\begin{displaymath}
\Gamma_{e}=\frac{e^{2}}{kT}\left( \frac{4\pi }{3}n_{e}\right) ^{1/3}\ll 1 .
\end{displaymath} (4.53)

В [403] ТФДХ теория применяется при \( \Gamma_{e}\lt 0.1 \), при \( \Gamma_{e}\gt 1 \) используются результаты расчетов Монте-Карло, а при \( 0.1\lt \Gamma_{e}\lt 1 \) результаты обоих методов интерполируются.

Для определения степени ионизации используется уравнение Саха, в котором учитывается, в дополнение к (1.8), конечная степень вырождения свободных электронов, влияющая на их химический потенциал. Сдвиг уровней ионов и молекул при учете экранирующего ДХ потенциала приводит к уменьшению числа связанных состояний и уменьшении их глубины. Их учет в формуле Саха отражает ионизацию давлением. В целом, учет ионизации давлением сводится к сложной самосогласованной проблеме, в которой концентрация свободных электронов пе и ионов и,- определяет ДХ радиус и потенциал, а те, в свою очередь, определяют концентрацию посредством сдвига уровней, входящих в уравнение Саха. Задача осложняется отсутствием аналитических решений для уравнения Шредингера с ДХ потенциалом [240], поэтому необходимы численные расчеты. Результаты таких расчетов, выполненных в [403], представлены на рис. 15 и 16. Таблицы уравнения состояния с учетом кулоновского взаимодействия, рассчитанные аналогичным способом, даны в [358].

Рис. 15. Область на плоскости $(\lg\rho, \lg T)$, где рассчитывались термодинамические функции [403]. Слева от линии $\Gamma_e = \Gamma_1$ = 0,1 расчеты велись по ТФДХ теории, справа от линии $\Gamma_e = \Gamma_1$ = 1,0 по методу Монте-Карло (см. (4.53)). Слева от линии $I_1$ водород и гелий частично ионизованы, между линиями $I_1$ и $I_2$ водород ионизован полностью, а гелий частично, а справа от линии $I_2$ оба ионизованы полностью, штрихами указаны линии постоянного значения $F = 4F_{1/2}(\eta)/\sqrt{\pi\,}$, характеризующие вырождение, в невырожденной плазме $F \ll 1$.
Рис. 16. Отношение кулоновской поправки к давлению к полному давлению при T = const. На каждой кривой указана величина $\lg T$, для которой $\lg P_0 = 2,74118$, $\lg T = 4,58825$

Задача. Сделать гладкую интерполяцию уравнения состояния холодного вещества нейтронных звезд.

Решение [14]. Гладкая с первой производной интерполяция уравнения состояния \( P(\rho ) \) из табл. 9 [479] определяется следующими формулами (\( \rho \) - полная плотность массы-энергии с учетом взаимодействия):

\begin{displaymath}
P=\left\{ \begin{array}{ll} \displaystyle{ P^{(1)}=b_{1}\rho ^{5/3}/(1+c_{1}\rho ^{1/3}) ,\quad } & \rho \leq \rho _{1} , \\ \displaystyle{ P^{(k)}=a\cdot 10^{^{b_{k}(lg\rho -8.419)^{c_{k}}}} ,\quad } &\rho _{k-1}\leq \rho \leq \rho _{k} , \\
\end{array}\right. \end{displaymath} (1)

\( k=2,3,\ldots 6 \),
\(\displaystyle{
\begin{array}{@{}lll}
a=10^{26.1673}, & c_{1}=10^{-2.257}, & \rho_{1}=10^{9.419}, \\
b_{1}=10^{12.40483}, & c_{2}=1.1598, & \rho_{2}=10^{11.5519}, \\ b_{2}=1, & c_{3}=0.356293, & \rho_{3}=10^{12.26939}, \\
b_{3}=2.5032, & c_{4}=1.2972138, & \rho_{4}=10^{14.302}, \\
b_{4}=0.70401515, & c_{5}=2.117802, & \rho_{5}=10^{15.0388}, \\
b_{5}=0.16445926, & c_{6}=1.237985, & \rho_{6} \gg \rho_{5} \\
b_{6}=0.86746415, \\
\end{array}}\)

Непрерывность производных \( dP/d\rho \) в точках \( \rho =\rho _{k} \) из (1) достигается с помощью сглаживания в виде

\begin{displaymath}
P(\rho )=\left\{ \eqalign{ & P^{(k)},\quad \rho \in \left[ \rho _{k-1}+\xi _{k-1},\rho _{k}-\xi _{k}\right], \quad k=1,2,...6, \cr & \theta _{k}P^{(k)}+(1-\theta _{k})P^{(k+1)}, \quad \rho \in \left[ \rho _{k}-\xi _{k}\,,\right. \cr & \left.\rho _{k}+\xi _{k}\right]\,, \quad \rho _{0}+\xi _{0}=0\,,
}\right. \end{displaymath} (2)

где

\begin{displaymath}
\theta _{k}=\theta _{k}(\rho )=\frac{1}{2}-\frac{1}{2}\sin\left( \frac{\pi }{2\xi _{k}}(\rho -\rho _{k})\right) ,
\end{displaymath}


\begin{displaymath}
\xi _{k}=0.01\rho _{k} .
\end{displaymath}



<< 3. Уравнение состояния при ... | Оглавление |

Публикации с ключевыми словами: Эволюция звезд - физические процессы
Публикации со словами: Эволюция звезд - физические процессы
См. также:
Все публикации на ту же тему >>

Оценка: 2.6 [голосов: 87]
 
О рейтинге
Версия для печати Распечатать

Астрометрия - Астрономические инструменты - Астрономическое образование - Астрофизика - История астрономии - Космонавтика, исследование космоса - Любительская астрономия - Планеты и Солнечная система - Солнце


Астронет | Научная сеть | ГАИШ МГУ | Поиск по МГУ | О проекте | Авторам

Комментарии, вопросы? Пишите: info@astronet.ru или сюда

Rambler's Top100 Яндекс цитирования