Астрономическое образование с сохранением традиций
ОБЩАЯ АСТРОФИЗИКА:

К.В.БЫЧКОВ

ТЕОРИЯ СПЕКТРАЛЬНОЙ КЛАССИФИКАЦИИ

Учебное пособие по курсу "Астрономия"



ОБЪЯСНЕНИЕ СПЕКТРАЛЬНОЙ КЛАССИФИКАЦИИ

        Слои звёзд, в которых образуются абсорбционные линии, как правило, близки к состоянию термодинамического равновесия. Во всяком случае это относится к линиям, по которым проводится спектральная классификация. Поэтому термодинамические методы активно применяются при анализе звёздных атмосфер.

        Конечно, полная задача об эквивалентной ширине спектральных линий не сводится только к термодинамике. При её решении необходимо учитывать также перенос излучения, профиль коэффициента поглощения, движение газа, вращение звезды и другие процессы. Но существует обстоятельство, позволяющее упростить полную задачу, сохраняя её физическое содержание. Дело в том, что эквивалентная ширина связана с числом поглощающих атомов, причём при заданных значениях атомных параметров эта связь однозначна и монотонна. Поэтому решение вопроса о населённости нижнего уровня перехода позволяет понять некоторые качественные аспекты спектральной классификации. Напомним, что при переходе к более горячим атмосферам определённые линии сначала усиливаются, затем проходят через максимум, а потом начинают ослабляться. Каждая линия имеет свою температуру максимума. Поставим следующие вопросы:

  1. Почему линии быстро ослабляются при переходе от максимума к более низким температурам?
  2. Почему ослабление в сторону высоких температур происходит значительно медленнее?
  3. Как связана температура максимума со структурой атома?
Ответы на них мы можем получить, не выходя за рамки термодинамики.



1  ФОРМУЛЫ САХА И БОЛЬЦМАНА.

        Этих двух формул вполне достаточно для первого знакомства с термодинамикой звёздных атмосфер. Формула Больцмана описывает населённость возбуждённых состояний, а формула Саха - состояние ионизации химического элемента. Первая выведена в предположении постоянства числа частиц, а вторая учитывает химические реакции, то есть возможное исчезновение старых частиц и появление новых.



1.1  Формула Больцмана и статистический вес.

        Формула Больцмана связывает друг с другом населённости дискретных уровней иона или атома:
 n2

n1
=  g2

g1
exp ж
и
 E1-E2

kT
ц
ш
 .
(1)
Здесь nj обозначает число атомов на j-м энергетическом уровне, Ej - его энергию, gj - статистический вес. Количество атомов рассчитываем на единицу объёма, то есть, n - плотность числа частиц.

        Статистическим весом называют число различных состояний атома, имеющих одну и ту же энергию. В случае атома водорода вес энергетического уровня с номером k равен
gk=2k2 .
(2)
         В более сложных системах различают термы и уровни1. Реально в каждом атоме существуют уровни. Вес уровня определяется его полным моментом J:
gJ=2J+1 .
(3)
        В отличие от уровня, терм - воображаемое понятие и представляет собой среднее положение нескольких близко расположенных уровней. Вес терма определяется суммарным спином электронов S и их суммарным орбитальным моментом L:
gLS=(2S+1)(2L+1) .
(4)
        Вес терма равен сумме весов всех составляющих его уровней.



1.2  Формула Саха
.
        Формула Саха описывает равновесное состояние ионизации:
 nine

na
=  2Ui

Ua
 (2p mT)3/2

h3
e-P/kT .
(5)
        Как и в случае формулы Больцмана, число частиц пропорционально статистическому весу. Нас интересует полное число ионов или атомов, поэтому вместо веса отдельного энергетического уровня фигурирует сумма по состояниям:
Ua=g0+
е
ex 
gexexp(-Eex/kT) .
(6)
        Здесь g0 - вес основного состояния, а индексом <<ex>> помечены веса и энергии возбуждённых состояний, по которым выполняется суммирование. Аналогично вычисляется сумма по состояниям иона Ui. Множитель 2 равен статистическому весу электрона - числу возможных проекций его спина.

        В приложении к звёздным атмосферам температуру удобно измерять в электрон-вольтах:
TeV=  T

11605 K
 .
        Подставляя в (5) численные значения атомных констант, приходим к формуле
 nine

na
=B  Ui

Ua
TeV3/2e-P/TeV        B » 6 1021.
(7)
        Начиная с этого места условимся, что если температура измеряется в любых энергетических единицах, но не в градусах, то в этих же единицах выражены все потенциалы - ионизации и возбуждения.

        Теперь обсудим важное отличие формулы Саха от формулы Больцмана. Если (5) переписать в виде
 ni

na
=  2Ui

Ua
 (2p mT)3/2

h3ne
e-P/kT ,
то мы увидим аналогию с (1), но в правой части присутствует дополнительный безразмерный множитель
gtr=  (2pmkT)3/2

h3ne
 .
        Это не что иное, как число квантовых состояний электрона, связанное с его перемещением в пространстве. В формуле Больцмана такие множители сокращаются, так как при возбуждении число частиц остаётся прежним. А формула Саха описывает реакции с изменением числа частиц2, поэтому множитель gtr остаётся.

        В условиях звёздных атмосфер численное значение gtr лежит в диапазоне от 105 до 108. Следовательно, заметная ионизация всех химических элементов имеет место при температурах, в несколько раз меньших потенциала ионизации атома. Например, в звёздах спектрального класса А0 водород ионизован более, чем наполовину, хотя температура их атмосфер в 12-15 раз меньше потенциала ионизации атома водорода.



1.3  Проблема суммы по состояниям.

        Вернёмся к формуле (6) и рассмотрим внимательнее сумму
S=
е
ex 
gexexp(-Eex/kT)
(8)
для атома водорода. На первый взгляд может показаться, что она невелика по сравнению с g0. Например, в уже упоминавшихся звёздах класса А0 численное значение g1exp(-E1/kT) оставляет менее 0.1% от g0, а экспоненциальный множитель в следующих слагаемых ещё меньше. Однако энергия возбуждения имеет верхний предел - потенциал ионизации атома. Следовательно, для любого Eex справедливо
exp(-Eex/kT) > exp(-P/kT) .
        Заменив в (8) энергии возбуждения на потенциал ионизации, получим очевидное неравенство:
S > exp(-P/kT)×2 Ґ
е
j=1 
j2 .
        Теперь стало совершенно ясно, что такая сумма не существует: она бесконечно велика!

        Объясняется полученный парадокс просто: атом, взаимодействующий с другими частицами, имеет только ограниченное число уровней, а схема с бесконечным множеством состояний вблизи границы ионизации - это идеализация, справедливая только для уединённого атома. В атмосферах звёзд обычно реализуется не более 20-30 уровней, и вклад суммы S редко превышает слагаемое g0. Не ставя своей целью выполнение рафинированных расчётов, мы примем простое решение, положив сумму по соcтояниям равной весу нижнего уровня:
Ua=g0 .
(9)
        Аналогичная формула, конечно, справедлива и для иона.

        Всё сказанное относится не только к атому водорода, но и к любой атомной системе, так как все состояния с высокой степенью возбуждения близки к соcтояниям атома водорода или водородоподобным ионам.



2  ЛИНИИ ГЕЛИЯ.

Самой простой системой является атом водорода, и вполне логично было бы начать изложение именно с его бальмеровской серии. Однако пойти по такому пути нам мешает принятая в предыдущем разделе гипотеза о его полной ионизации. Исследуя поведение субординатных линий в области высоких температур, мы должны рассмотреть конкуренцию ионизации и возбуждения. Но формула (10) не позволяет корректно учесть уравнение Саха для водорода. Поэтому мы начнём с гелия: изучение его субординатных переходов вполне допустимо в рамках сделанных предположений. К тому же для задач спектральной классификации, как мы убедились выше, гелий значительно интереснее водорода, так как его линии более чувствительны к температуре, а их интенсивность ближе к интенсивности линий других химических элементов3.


2.1  Электронная плотность.

        В поставленной задаче населённость уровней зависит от двух параметров: температуры T и плотности числа тяжёлых частиц n. Под тяжёлыми частицами мы понимаем ионы и атомы, но не электроны. Строго говоря, для вычисления n мы должны просуммировать число частиц ионов и атомов всех химических элементов. Но если мы будем учитывать только водород, то сделаем ошибку лишь около десяти процентов, соответственно содержанию гелия - следующего по обилию химического элемента в звёздах главной последовательности. Далее, будем считать водород практически полностью ионизованным, что вполне допустимо для рассматриваемых нами звёзд горячее А0. Гелий в ОВ-звёздах может быть ионизован однократно, либо двукратно, то есть, его вклад в электронную плотность составляет от десяти до двадцати процентов. Пренебрегая этим вкладом, приходим к следующей простой формуле для вычисления электронной плотности:
ne=n ,
(10)
которой и будем пользоваться в дальнейшем.


2.2  Состояние ионизации гелия.

        Гелий может находиться в трёх состояниях ионизации: HeI, HeII и HeIII, причём только HeI и HeII имеют дискретные уровни. Напомним, что HeI обозначает атом гелия, а HeII - его первый ион. Обозначим их относительные концентрации как x1, x2 и x3. Из их определения следует
x1+x2+x3=1 .
(11)
        Формула Саха даёт два уравнения
 x2

x1
=f1 ,        x3

x2
=f2 .
(12)
        Обе функции, f1 и f2, содержат экспоненциальные множители:
fj=ajexp(-Pj/kT) .
(13)
        Здесь aj - тоже функции температуры ( µ T3/2), но они меняются значительно медленнее, чем экспонента.

        Хотя система линейных уравнений (11) и (12) решается без труда, тем не менее физическое содержание коэффициентов fj позволяет выполнить дальнейшие упрощения.

        Во-первых, в рассматриваемом диапазоне температур всегда выполнено условие:
f2 < f1 .
(14)
        Во-вторых, экспоненциальная зависимость fj позволяет рассматривать сравнительно чётко очерченные диапазоны температуры, в каждом из которых можно пренебречь тем или иным состоянием ионизации. Выпишем решение (11) и (12) в явном виде:
 1

x1
=1+f1(1+f2) ,                               
(15)
 1

x2
=1+  1

f1
+f2 ,                               
(16)
 1

x3
=1+  1

f2
ж
и
1+  1

f1
ц
ш
 .                               
(17)
        Отметим вытекающее из (12) полезное тождество:
 x3

x1
=f1f2 .
(18)
        Теперь рассмотрим различные диапазоны изменения f1 и f2.

  1. f1f2 << 1. В этом случае можно пренебречь x3 по сравнению с x1 и x2:
    x1+x2 » 1 .
    (19)   
  2. f1f2 >> 1. Здесь пренебрегаем x1 по сравнению с x2 и x3:
    x2+x3 » 1 .
    (20)   
        Каждое из уравнений (19) и (20) значительно проще, чем (11).

        Теперь мы можем непосредственно заняться интересующей нас величиной - числом атомов nl на нижнем уровне перехода. Напомним, что все наблюдаемые в оптическом диапазоне линии гелия относятся к субординатным, то есть, их нижний уровень является возбуждённым. Обозначим его энергию через El и введём безразмерную величину
y=  nl

nHe
 ,
(21)
где nHe - полное число частиц гелия, просуммированное по всем состояниям ионизации. Рассмотрим сначала переходы нейтрального гелия.



2.3  Скорость убывания и возрастания линий HeI.

        Чтобы воспользоваться формулами Саха и Больцмана, перепишем (21) для данного случая в виде
y=  nl

nHeI
 nHeI

nHe
 ,
(22)
или
y=  gl

g1
exp ж
и
-  El

kT
ц
ш
x1(T,n) .
(23)
        При переходе от (22) к (23) мы учли то обстоятельство, что населённость возбуждённых уровней значительно ниже, чем основного, и, следовательно, мы можем принять число атомов гелия в основном состоянии равным полному числу атомов гелия. Действительно, температура даже самых горячих О-звёзд не превышает 4 эВ, а потенциал возбуждения даже самых низких уровней нейтрального гелия выше 19 эВ. Поэтому равновесная населённость его возбуждённых уровней содержит экспоненциальный множитель exp(-19/4) < 0.01. Но и такая оценка ещё завышена. Линии нейтрального гелия наиболее заметны у звёзд класса В, температура которых не превышает 2.5 эВ и тот же множитель меньше, чем exp(-19/2.5) < 0.001.

        Ясно, что исследование нейтрального гелия целесообразно проводить в области температур, которой отвечает уравнение (19). Подставляя его решение в (23), получим
y=  gl

g1
exp ж
и
-  El

kT
ц
ш

1+a1(n,T)exp ж
и
-  P1

kT
ц
ш
 .
(24)
        Функция y(T) не монотонна, она имеет максимум. Тем самым она отражает реальные свойства линий гелия, которые усиливаются при переходе от поздних к ранним подклассам В, а затем ослабляются у более горячих звёзд.

        Рассмотрим область настолько низких температур, что в знаменателе (24) можно пренебречь вторым слагаемым по сравнению с единицей. Физически это означает, что гелий находится в нейтральном состоянии, и населённость y(T) экспоненциально растёт с температурой по формуле Больцмана:
y µ exp(-El/kT) .
(25)

        Энергия возбуждения El оптических переходов атома гелия составляет 19-20 эВ. Температура атмосфер поздних и средних подклассов В лежит в диапазоне 1-1.5 эВ. Соответственно, показатель экспоненты (25) находится в диапазоне от 15 до 20, что означает очень быстрый рост населённости с температурой.

        В противоположном случае высоких температур мы пренебрегаем единицей в знаменателе дроби (24) и приходим к:
y »  1

a(n,T)
exp ж
и
 P1-El

kT
ц
ш
 .
(26)
        Оно описывает убывающую функцию температуры (напомним, что множитель a(n,T) пропорционален T3/2). Здесь мы видим результат влияния двух факторов. С одной стороны, увеличивается населённость возбуждённого уровня относительно основного состояния по формуле (1), но с другой - усиливается ионизация гелия. В рассматриваемом предельном случае концентрация атомов гелия x1 зависит от температуры как
x1 »  1

f1
~  1

a(n,T)
exp ж
и
 P1

kT
ц
ш
(27)
в согласии с (15) при f1 >> 1. Ионизация нейтральных атомов здесь играет большую роль, чем усиление возбуждения. Главная причина в том, что потенциал ионизации больше энергии возбуждения любого дискретного уровня.

        Формула (26) описывает значительно более слабую зависимость населённости от температуры, чем (25). В самом деле, все возбуждённым уровни атома гелия близки к границе ионизации: при P1=24.6 эВ значения El лежат в диапазоне от 19.5 до 21 эВ, поэтому разность P1-El не превышает 4 эВ.

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



2.4  Температура максимума.

        Определим температуру, отвечающую максимуму функции (24). Для удобства расчётов температуру в этом разделе измеряем в электрон-вольтах и введём новую переменную:
z=  g1

gl
 1

y
 .
        Будем искать минимум z как функции аргумента q=1/TeV. В силу монотонной зависимости как z от y, так и q от TeV, задача о минимума y(TeV) эквивалентна задаче о максимума z(q). Итак, ищем значение q, при котором выражение
z=exp(Elq)[1+A(n)q-3/2exp(-P1q)]
(28)
имеет минимум. Здесь мы выписали зависимость от температуры в явном виде и ввели множитель
A(n)=  Ui

Ua
 B

n
.
        Проинтервьюировав (28) по q и приравняй производную нулю, приходим к уравнению
exp(P1q) = A(n)q-3/2
 3

2
+(P1-El)q

Elq
 .
(29)
        Оно без труда решается численными методами, если задать конкретные значения параметров. Но наша задача заключается в другом: надо получить, хотя бы приближённо, аналитическую зависимость температуры максимума Tm от P1 и El.

        Произведение P1q является большой величиной: от десяти и выше. Это позволяет нам утверждать, что на величину Tm в гораздо большей степени влияет потенциал ионизации, чем энергия возбуждения. В справедливости сделанного утверждения убедимся сначала на модельном примере. Рассмотрим функцию
w=exp(Elq)[1+Bexp(-P1q)]
        Она отличается от (28) тем, что зависящий от температуры множитель Aq-3/2 заменён некоторой константой C. Температура максимума w выражается аналитически:
Tmw=   P1

ln й
л
C ж
и
 P1

El
-1 ц
ш
щ
ы
 .
(30)
        Мы видим, что Tmw зависит от P1, в основном, линейно, а от El - только логарифмически. Можно надеяться, что (30) правильно передаёт и основные черты решения (29), так как множитель q3/2 меняется значительно медленнее, чем exp(P1q) и в первом приближении его можно заменить постоянной величиной.

        Ниже мы выпишем аналитическое приближение к решению уравнения (29), но сначала поставим вопрос с физической точки зрения: почему энергия возбуждения столь слабо влияет на величину Tm? Дело в том, что формула Больцмана, куда входит El, сама по себе не имеет экстремумов: чем выше температура, тем сильнее заселены возбуждённым уровни относительно основного. Появление максимума обусловлено переходом атомов в ионизованное состояние. Этот процесс описывается формулой Саха и им управляет потенциал ионизации.

        Логарифмическая зависимость Tm от отношения P1/El происходит от эффекта насыщения возбуждённого уровня. Состояния с меньшим потенциалом возбуждения заселяются при менее высоких значениях температуры и, следовательно, несколько раньше реагируют на изменение ионизации.

        Теперь модифицируем формулу (30) таким образом, чтобы она описывала максимум функции (28). Конечно, точного решения мы не получим, но нашей целью, напомним, является получение аналитического выражения, качественно правильно описывающего зависимость Tm. Прологарифмируем уравнение (29):
P1q = lnA(n)-  3

2
lnq+ln й
л
 3

2
+(P1-El)q щ
ы
-ln(Elq) .
        Среди слагаемых правой части по величине выделяется первое. Действительно, численное значение q около единицы, а сумма под знаком логарифма в третьем слагаемом и произведениеElq - порядка десяти. В то же время величина A(n) в условиях звёздных атмосфер при плотности n » 1015ё16 составляет A(n) » 105ё6. Отношение весов нижнего и верхнего уровней не превышает двойки: спин S не меняется, а орбитальное квантовое число L либо также остаётся прежним, либо меняется на единицу. Отношение gi/ga в случае гелия равно двум. Следовательно, мы не сделаем большой ошибки, пренебрегая в правой части всеми слагаемыми, кроме первого. В результате получим приближённую оценку:
q0 »  lnA(n)

P1
 .
        С помощью q0 напишем выражение для константы C:
C=A(n)q0-3/2
и подставим его в (30). В результате получим
Tm=  P1

lnA(n)+  3

2
lnP1-  3

2
lnlnA(n)+ ln ж
и
 P1

El
-1 ц
ш
 .
(31)
        Снова видна сильная зависимость от потенциала ионизации и слабая - от энергии возбуждения.



3  БАЛЬМЕРОВСКАЯ СЕРИЯ АТОМА ВОДОРОДА

        По аналогии с задачей о гелии, исследуем вопрос о населённости n2 второго уровня атома водорода - нижнего уровня всех переходов бальмеровской серии.

        Как мы убедимся ниже, решения этих задач имеют сходные черты: немонотоннная зависимость населённости возбуждённых уровней от температуры; наличие максимума, положение которого определяется, главным образом, потенциалом ионизации атома; медленное уменьшение от максимума к высоким температурам и быстрое падение в сторону холодных звёзд. Главной причиной такого поведения населённостей является сходство структуры энергетических уровней водорода и гелия:

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

Тем не менее, для водорода требуются специальные расчёты, так как здесь неприменимо предположение (10) о независимости электронной плотности ne от температуры. Водород является основным донором электронов в звёздах спектральных классов O, B и A. Поэтому в рассматриваемой задаче необходимо вычислять ne одновременно с населённостью уровней.


3.1  Состояние ионизации водорода

Введём обозначения
xi=  np

n
,       xa=  na

n
,
(32)
соответственно, для относительной концентрации ионов и нейтральных атомов водорода. Здесь np - число протонов, а na - число нейтральных атомов водорода в единице объёма. Введём параметр n, в этом разделе он равен плотности водорода по числу частиц:
n=np+na.
(33)
        Принятая нами величина примерно на 10% отличается от той, что приводится в справочниках, так как содержание водорода в звёздах главной последовательности составляет около 90%. В проводимых здесь расчётах столь небольшим различием можно перенебречь.

        Мы предполагаем, что все электроны происходят только за счёт ионизации водорода. Следовательно, их число равно числу протонов:
ne=np.
(34)
Это предположение вносит в расчёты ошибку, заведомо меньшую 10%, так как гелий ионизован значительно слабее водорода.

        Формула Саха с учётом (34) даёт квадратное уравнение для xa:
 (1-xa)2

xa
=  B

2n
 TeV3/2  exp ж
и
-  P

TeV
ц
ш
.
(35)
Здесь учтено, что вес основного состояния атома водорода, согласно (2), равен двум:
ga=2,
(36)
а состояние иона не имеет вырождения:
gi=1.
(37)
Напомним, что величина B известна из (7).



3.2  Зависимость от плотности и от температуры

        Полную плотность водорода n в наших расчётах будем считать постоянным параметром. В действительности плотность атмосферы разная у звёзд разных спектральных классов:

Sp O5B0 B5 A0 A5 F0 F5
lgn 15.0 15.0 15.0 15.2 15.6 16.1 16.6

        Как видно из таблицы, в диапазоне от B5 до F0 величина n меняется примерно в 40 раз. Тем не менее, мы будем пренебрегать вариациями плотности по сравнению с влиянием температуры на населённость уровней. В самом деле, при переходе от класса B5 к F0 температура атмосферы уменьшается примерно вдвое: от 15500 K до 7400 K. При этом экспоненты exp(-E2/TeV) и exp(-P/TeV) падают, соответственно, в 4·103 и 7·104 раз. Поэтому для звёзд главной последовательности в первом приближении эффект плотности можно не рассматривать. Для определённости параметр n примем равным плотности атмосферы звезды класса А0:
lgn=15.2.
(38)
Именно у этого класса звёзд главной последовательности наиболее сильна серия Бальмера.



3.3  Населённость второго уровня

        Преобразуем уравнение (35) к виду
 (1-xa)2

xa
=f,
(39)
где
f=a (-P/ TeV), a=C TeV3/2, C=B/(2n).
        Одно уравнение (39) заменяет целую систему уравнений задачи о гелии (11), (12), но оно является нелинейным. Так проявляется зависимость электронной плотности от температуры.

Введём безразмерную величину
y=  1

4
 n2

n
,
(40)
характеризующую населённость второго уровня атома водорода. Множитель 1/4 в правой части введён для удобства последующего использования формулы Больцмана. Он равен отношению статистических весов g1/g2. Далее, воспользовавшись формулами Больцмана (1) и Саха (39), выражаем y как функцию температуры и плотности:
y=exp ж
и
-  E2

TeV
ц
ш
· xa(TeV,n),
(41)
в которой xa(TeV,n) является решением уравнения (39).



3.4  Высокие и низкие температуры

        Снова, как и в случае гелия, убеждаемся, что функция xa(T) не монотонна. В области высоких температур, когда правая часть (39) мала, то есть выполнено условие
f >> 1,
(42)
для относительной концентрации нейтрального водорода справедлива простая формула:
xa »  1

f
,
(43)
откуда следует
y=  2n

B
  TeV-3/2 exp ж
и
 P-E2

TeV
ц
ш
.
(44)
Мы получили сравнительно медленное уменьшение населённости второго уровня с ростом температуры. Как и в случае гелия, это обусловлено близостью возбуждённого состояния к границе ионизации.

В области предельно низких температур (f << 1) ионизационное равновесие сдвинуто в сторону атомов:
xa » 1.
(45)
В этом случае из (41) следует
y » exp ж
и
-  10.2

TeV
ц
ш
.
(46)
Здесь относительно большая величина потенциала возбуждения приводит к экспоненциально быстрому падению населённости второго уровня по мере уменьшения температуры.

Итак, функция y(T) возрастает при переходе к самым холодным звёздам и убывает в области предельно высоких температур. Следовательно, она имеет максимум. Перейдём к определению температуры максимума.



3.5  Температура максимума населённости

        Выполним замены переменных, которые нам несколько упростят вычисления. Максимум функции y(T) приходится на ту же температуру, что и минимум обратной ей величины
z=  1

y
(47)
как функции обратной температуры
q =  1

TeV
.
(48)
Эквивалентность обеих задач вытекает из монотонности функций и . Итак, преобразуется в
z=  exp(E2 q)

xa(q,n)
.
(49)
Минимум функции z(q) ищем в области небольших значений xa. А именно, будем считать, что можно пренебречь квадратом концентрации:
xa2 << 1.
(50)
В этом случае левая часть становится линейной:
(1-xa)2 » 1-2xa.
(51)
Подставим в и решим полученное линейное уравнение относительно xa:
xa »  1

2+f
.
(52)
Ясно, что последняя формула справедлива только при выполнении условия . В этой области рассмотрим функцию
w(q)=exp(E2 q)(2+f).
(53)
Оно получается из , если заменить точное значение xa на его приближение . Мы исследуем на минимум функцию , а затем проверим выполнение условия .

Сочетание экспоненциального и степенного множителей в выражениях для f в приводит к тому, что уравнение
 dw

dq
=0,
(54)
выражающее требование минимума w(q), не имеет аналитического решения. Для дальнейших упрощений воспользуемся известным из наблюдений фактом, что в звёздах класса А0, где наиболее сильны линии бальмеровской серии, температура атмосферы значительно, в 10-15 раз меньше потенциалов ионизации и возбуждения водорода. Как мы уже писали в разделе , экспоненциальные множители при таких значениях аргумента очень сильно меняются с температурой. Их относительные изменения велики по сравнению с изменением множителя TeV3/2. Следовательно, в первом приближении температуру в степенном множителе можно положить равной некоторой величине
TeV=T0,
(55)
которую можно подобрать, исходя из особенностей задачи. В данном случае разумно принять для T0 значение, близкое к температуре звезды класса А0:
T0=1 eV.
(56)
С учётом сказанного продифференцируем функцию w(q), удерживая только быстро меняющиеся слагаемые:
 dw

dq
=exp(E2q)·[E2(2+f)-aPexp(-Pq)].
(57)
и приравняем нулю производную dw/dq. В результате получим алгебраическое уравнение, решением которого является температура
Tmax =  P

ln[0.85·106·T03/2 (P/E2-1)]
.
(58)
Подставляя сюда параметры атома водорода:
 P

E2
=4/3,
с учётом получим
Tmax=1.07  eV
(59)
в хорошем согласии с наблюдениями.

Теперь проверим выполнимость условия . Подставив в , получим
f(Tmax)=6.35.
Соответствующее решение квадратного уравнения получается равным
xa(Tmax)=0.12.
Таким образом, величина xa2 действительно даёт малый вклад в числитель левой части и приближение вполне оправдано.

Температура максимума населённости n2, как следует из , определяется, главным образом, потенциалом ионизации
Tmax µ P,
в то время как от потенциала возбуждения она зависит только логарифмически. Причина этого рассмотрена в разделе 2.4 предыдущей главы.



4  ЗАКЛЮЧЕНИЕ

        Итак, на примере линий водорода и гелия в звёздах классов O, B и A мы убедились, что методы термодинамики, выражающиеся в простых формулах Саха и Больцмана, позволяют описать различие спектров у звёзд с разной температурой. Анализ эффектов давления, то есть различия спектров звёзд разной светимости при одной и той же температуре, требует привлечения сведений из атомной физики и теории переноса излучения в линиях и континууме.





Footnotes:

1Их следует отличать от энергетических уровней

2Один атом превращается в пару ион-электрон

3Линии водорода являются самыми сильными во всех звёздах горячее F0.