Skip to content

Commit d1846fb

Browse files
vobolgusclaude
andcommitted
Масса Джинса
Co-Authored-By: Claude Fable 5.1 <noreply@anthropic.com>
1 parent 58df6b2 commit d1846fb

7 files changed

Lines changed: 113 additions & 0 deletions

File tree

Lines changed: 107 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,107 @@
1+
\subsection{Масса Джинса}
2+
\label{sec:jeans-mass}
3+
Звёзды рождаются из холодных облаков межзвёздного газа. Сжимает облако его собственная гравитация, а удерживает от сжатия давление газа. При каких условиях гравитация побеждает и облако начинает сжиматься ответ даёт \term{критерий Джинса}, полученный Джеймсом Джинсом в 1902 году.
4+
5+
Рассмотрим ту же модель, что и при оценке давления в центре звезды, \lookSecRef{sec:pressure-and-temperature}: однородное шарообразное облако массы $M$, радиуса $R$, плотности $\rho$ и температуры $T$, состоящее из частиц массы $m$.
6+
7+
\paragraph{Гравитационная энергия облака} Будем собирать облако из тонких сферических слоёв, добавляя их снаружи один за другим. Когда собран шар радиуса $r$ и массы $M(r) = \frac{4}{3}\pi\rho r^3$, очередной слой толщиной $dr$ имеет массу $dm = 4\pi\rho r^2\,dr$, а его потенциальная энергия в поле уже собранного шара
8+
\begin{equation*}
9+
dU = -\frac{G M(r)\,dm}{r} = -\frac{16}{3}\pi^2 G \rho^2 r^4\,dr.
10+
\end{equation*}
11+
Проинтегрируем по всем слоям и выразим плотность через полную массу, $\rho = 3M / (4\pi R^3)$:
12+
\begin{equation}
13+
U = -\frac{16}{3}\pi^2 G \rho^2 \int\limits_0^R r^4\,dr
14+
= -\frac{16}{15}\pi^2 G \rho^2 R^5
15+
= -\frac{3}{5}\frac{GM^2}{R}.
16+
\label{eq:sphere-grav-energy}
17+
\end{equation}
18+
19+
\paragraph{Критерий Джинса} Кинетическая энергия облака~--- это энергия теплового движения его $N = M / m$ частиц, \lookSecRef{sec:mkt}. Чтобы не путать её с температурой, обозначим кинетическую энергию $K$:
20+
\begin{equation*}
21+
K = \frac{3}{2} N k T = \frac{3}{2}\frac{M}{m}kT.
22+
\end{equation*}
23+
По теореме о вириале, \lookSecRef{sec:virial}, облако находится в равновесии, если $2K = -U$. Если же тепловой энергии не хватает, $2K < -U$, гравитация побеждает и облако сжимается:
24+
\begin{gather*}
25+
3\frac{M}{m}kT < \frac{3}{5}\frac{GM^2}{R},\\
26+
M > \frac{5kT}{Gm}R.
27+
\end{gather*}
28+
В это условие радиус входит вместе с массой. Исключим его с помощью плотности, $R = \left(3M / 4\pi\rho\right)^{1/3}$:
29+
\begin{gather*}
30+
M > \frac{5kT}{Gm}\left(\frac{3M}{4\pi\rho}\right)^{1/3},\\
31+
M^{2/3} > \frac{5kT}{Gm}\left(\frac{3}{4\pi\rho}\right)^{1/3},\\
32+
M > \left(\frac{5kT}{Gm}\right)^{3/2}\left(\frac{3}{4\pi\rho}\right)^{1/2}.
33+
\end{gather*}
34+
Правая часть называется \term{массой Джинса}:
35+
\begin{equation}
36+
M_J = \left(\frac{5kT}{Gm}\right)^{3/2}\left(\frac{3}{4\pi\rho}\right)^{1/2},
37+
\label{eq:jeans-mass}
38+
\end{equation}
39+
здесь $T$ и $\rho$~--- температура и плотность облака, $m$~--- масса его частиц. Облако массы больше $M_J$ не может удерживаться давлением и сжимается под действием собственной гравитации. Масса Джинса растёт с температурой, $M_J \propto T^{3/2}$, и падает с плотностью, $M_J \propto \rho^{-1/2}$: чем холоднее и плотнее газ, тем меньше массы достаточно для сжатия.
40+
41+
То же условие можно записать и через размер облака. Подставим в $M > 5kTR / (Gm)$ массу $M = \frac{4}{3}\pi\rho R^3$:
42+
\begin{gather*}
43+
\frac{4}{3}\pi\rho R^3 > \frac{5kT}{Gm}R,\\
44+
R > \sqrt{\frac{15 kT}{4\pi G m \rho}}.
45+
\end{gather*}
46+
Правая часть называется \term{радиусом Джинса}:
47+
\begin{equation}
48+
R_J = \sqrt{\frac{15 kT}{4\pi G m \rho}}.
49+
\label{eq:jeans-radius}
50+
\end{equation}
51+
Облако данной плотности и температуры сжимается, если его размер больше $R_J$, и остаётся в равновесии, если меньше.
52+
53+
\paragraph{Смысл радиуса Джинса} Давление противостоит гравитации не мгновенно: возмущение давления распространяется по газу со скоростью звука $c_s$, которая по порядку величины равна тепловой скорости частиц, $c_s \sim \sqrt{kT / m}$. Гравитация же сжимает облако за время свободного падения $t_\text{св}$. Оценим его: частица на краю облака испытывает ускорение $g = GM / R^2 = \frac{4}{3}\pi G \rho R$ и, если считать его постоянным, проходит путь $R$ за время
54+
\begin{equation*}
55+
t_\text{св} \simeq \sqrt{\frac{2R}{g}} = \sqrt{\frac{3}{2\pi G \rho}} \sim \frac{1}{\sqrt{G\rho}}.
56+
\end{equation*}
57+
Заметим, что время свободного падения зависит только от плотности и не зависит от размера облака. Радиус Джинса \eqref{eq:jeans-radius} с точностью до числового множителя есть путь, который звук проходит за время свободного падения:
58+
\begin{equation*}
59+
R_J \sim \sqrt{\frac{kT}{m}} \cdot \frac{1}{\sqrt{G\rho}} \sim c_s t_\text{св}.
60+
\end{equation*}
61+
В облаке меньше $R_J$ давление успевает выровняться быстрее, чем гравитация заметно сожмёт облако, и сжатие останавливается. В облаке больше $R_J$ давление не успевает отреагировать, и сжатие ничто не сдерживает. Строгий расчёт, учитывающий распространение волн в газе, даёт те же зависимости $M_J$ и $R_J$ от температуры и плотности, отличаясь лишь числовыми множителями порядка единицы, поэтому формулы \eqref{eq:jeans-mass} и \eqref{eq:jeans-radius} следует понимать как оценки.
62+
63+
\begin{wrapfigure}[14]{r}{0.48\tw}
64+
\centering
65+
\vspace{-1pc}
66+
\tikzsetnextfilename{jeans-mass}
67+
\begin{tikzpicture}
68+
\footnotesize
69+
% Масса Джинса для молекулярного водорода, m = 2 m_p:
70+
% M_J / M_sun = 23.1 (T / 10 К)^{3/2} (n / 10^9 м^{-3})^{-1/2}
71+
% Степенной закон на логарифмических осях есть прямая, поэтому
72+
% каждую кривую задаём двумя крайними точками, n = 1e6 и 1e12 м^-3
73+
\begin{loglogaxis}[
74+
width = 0.43\tw,
75+
height = 5.2cm,
76+
clip = false,
77+
xmin = 1e6, xmax = 1e12,
78+
ymin = 0.1, ymax = 1e5,
79+
xlabel = {$n,~\text{м}^{-3}$},
80+
ylabel = {$M_J / M_\odot$},
81+
xtick = {1e6, 1e8, 1e10, 1e12},
82+
ytick = {0.1, 10, 1e3, 1e5},
83+
]
84+
\addplot [black] coordinates {(1e6, 730.4) (1e12, 0.7304)}
85+
node [pos = 1, anchor = west] {$10~\text{К}$};
86+
\addplot [black] coordinates {(1e6, 3795) (1e12, 3.795)}
87+
node [pos = 1, anchor = west] {$30~\text{К}$};
88+
\addplot [black] coordinates {(1e6, 23100) (1e12, 23.1)}
89+
node [pos = 1, anchor = west] {$100~\text{К}$};
90+
\addplot [dashes] coordinates {(1e6, 1) (1e12, 1)}
91+
node [pos = 0.12, anchor = south] {$M_\odot$};
92+
\addplot [only marks, mark = *, mark size = 1.5pt] coordinates {(1e9, 23.1)};
93+
\end{loglogaxis}
94+
\end{tikzpicture}
95+
\caption{Масса Джинса для молекулярного водорода при разных температурах; точка~--- облако из примера}
96+
\label{pic:jeans-mass}
97+
\end{wrapfigure}
98+
\paragraph{Пример} Звёзды образуются в молекулярных облаках: температура $T = 10~\text{К}$, концентрация молекул $n = 10^9~\text{м}^{-3}$, газ состоит в основном из молекулярного водорода, $m = 2m_p$, откуда $\rho = nm = 3.3 \cdot 10^{-18}~\text{кг}/\text{м}^3$. Подставляя в \eqref{eq:jeans-mass} и \eqref{eq:jeans-radius}, получаем, \lookPicRef{pic:jeans-mass},
99+
\begin{gather*}
100+
M_J \simeq 4.6 \cdot 10^{31}~\text{кг} \simeq 23 M_\odot,\\
101+
R_J \simeq 1.5 \cdot 10^{16}~\text{м} \simeq 0.5~\text{пк},
102+
\end{gather*}
103+
а время свободного падения составляет около $1.5$ миллиона лет. Для сравнения, в разреженном атомарном газе с $T = 100~\text{К}$, $n = 10^7~\text{м}^{-3}$ и $m = m_p$ масса Джинса равна $3 \cdot 10^4 M_\odot$, а радиус Джинса~--- $30~\text{пк}$. Именно поэтому звёзды рождаются не в тёплом межзвёздном газе, а в его самых холодных и плотных сгустках.
104+
105+
\paragraph{Фрагментация} Что происходит с облаком после того, как сжатие началось? При сжатии гравитация совершает работу, и выделяющееся тепло нагревает газ. Судьба этого тепла зависит от того, насколько облако прозрачно для собственного излучения, \lookSecRef{sec:optical-thickness}. Оптическая толщина облака по порядку величины равна $\tau \sim n \sigma R$, где $\sigma$~--- сечение взаимодействия фотона с частицей газа. Пока $\tau \ll 1$, фотон, испущенный в любой точке облака, свободно выходит наружу: облако прозрачно, тепло уносится излучением, и температура газа почти не меняется, а плотность растёт. Масса Джинса при этом падает как $\rho^{-1/2}$~--- облако сползает вправо и вниз по своей кривой на \picRef{pic:jeans-mass},~--- и части облака, которые раньше удерживались давлением, сами начинают сжиматься. Облако распадается на всё более мелкие сгустки~--- \term{фрагментирует}. Так из одного облака массой в тысячи солнечных рождается не одна гигантская звезда, а целое звёздное скопление, \lookSecRef{sec:open-clusters}.
106+
107+
При сжатии сгустка постоянной массы концентрация растёт как $n \propto R^{-3}$, поэтому его оптическая толщина увеличивается, $\tau \propto R^{-2}$. Рано или поздно $\tau$ становится порядка единицы: фотон поглощается раньше, чем успевает покинуть сгусток, и тепло остаётся внутри. С этого момента сжатие идёт с ростом температуры, а масса Джинса, пропорциональная $T^{3/2}$, снова увеличивается. Фрагментация прекращается, и такой сгусток, продолжая сжиматься уже как целое, становится \term{протозвездой}.

sections/astrophysics/optical-thickness.tex

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1,4 +1,5 @@
11
\subsection{Оптическая толщина. Закон Бугера -- Ламберта -- Бера}
2+
\label{sec:optical-thickness}
23
\begin{wrapfigure}[10]{c}{0.4\tw}
34
\centering
45
\vspace{-1pc}

sections/astrophysics/pressure-and-temperature.tex

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1,4 +1,5 @@
11
\subsection{Давление и температура в центре звезды}
2+
\label{sec:pressure-and-temperature}
23
В данном разделе найдем оценку для давления и температуры в центре звезды. Для этого рассмотрим такую модель, пусть звезда~--- однородное газовое облако радиусом $R$. Пусть плотность звезды равна $\rho$, а вещество~--- водородоподобная плазма, такое допущение справедливо для звёзд на главной последовательности. Тогда средняя масса частиц газа составляется примерно $m = 0.5m_p$, потому вещество состоит из одинакового числа протонов и электронов, так как в целом оно нейтрально, а $m_e \ll m_p$.
34

45
Аналогично поиску высоты атмосферы рассмотрим тонкий слой вещества высоты $dr$ с площадью сечения $dS$ на расстоянии $r$ от центра звезды. Пусть его масса равна $dm$. На слой действуют три силы: сила давления газа снизу $F_\uparrow = p(r)\,d S$, сила давления газа сверху $F_\downarrow = -p(r + \,d r)\,d S$ и сила гравитационного притяжения $F_g = -GM(r)/r^2\,d m$ от массы, находящейся внутри сферы радиуса $r$, так как суммарное воздействие от внешней оболочка равно нулю. Стационарность звезды обеспечивается балансом этих сил, значит

sections/celestial-mechanics/virial.tex

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1,4 +1,5 @@
11
\subsection{Теорема о вириале}
2+
\label{sec:virial}
23
Рассмотрим скалярную функцию механической системы, состоящей из $N$ материальных точек:
34
\begin{equation}
45
V = \sum_{k=1}^N \scalarNoVec{\vec F_k}{\vec r_k},

sections/objects/open-clusters.tex

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1,4 +1,5 @@
11
\subsection{Рассеянные звёздные скопления}
2+
\label{sec:open-clusters}
23

34
\begin{wrapfigure}[9]{r}{0.3\tw}
45
\centering

sections/physics/mkt.tex

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1,4 +1,5 @@
11
\subsection{МКТ и термодинамика}
2+
\label{sec:mkt}
23
\term{Молекулярно-кинетическая теория} изучает процессы, происходящие в макроскопических телах, на основе предположения о том, что вещество состоит из частиц, движение которых подчиняется законам механики Ньютона.
34

45
\term{Уравнение Менделеева-Клайперона} является следствием этой теории и~связывает давление $(P)$, объем $(V)$ и температуру $(T)$ идеального газа:

sys/astrophysics.tex

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -19,4 +19,5 @@ \section{Астрофизика}
1919
\input{sections/astrophysics/earth-atmosphere}
2020
\input{sections/astrophysics/optical-atmosphere-effects}
2121
\input{sections/astrophysics/pressure-and-temperature}
22+
\input{sections/astrophysics/jeans-mass}
2223
\input{sections/astrophysics/gauss-theorem}

0 commit comments

Comments
 (0)