Проанализируйте цветную диаграмму Герцшпрунга–Рассела для трёх кластеров разного возраста (молодой рассеянный кластер, открытый кластер среднего возраста и старый шаровой кластер): как по положению поворотной точки и последовательности главной ветви оценить возраст, металличность и первоначальную массу звёзд в кластере; приведите модельный расчёт для выбранного примера
Кратко — метод и признаки на CMD, затем модельный расчёт для примера (открытый кластер среднего возраста). 1) Как по CMD оценивают возраст, металличность и начальные массы: - Возраст: определяется положением поворотной точки (turn‑off, TO) главной ветви (MS). Чем старше кластер, тем краснее (более низкая TeffT_{\rm eff}Teff) и слабее (меньшая LLL) TO; соответственно масса звёзд на TO меньше. Для предварительной оценки используют аппроксимацию продолжительности главной последовательности: tMS≃1010 yr(MM⊙)−p,
t_{\rm MS}\simeq 10^{10}\,\mathrm{yr}\left(\frac{M}{M_\odot}\right)^{-p}, tMS≃1010yr(M⊙M)−p,
где \(p\sim 3\mbox{--}3.5\) (в зависимости от диапазона масс). Обращая её: MTO≃M⊙(1010 yrtcluster)1/p.
M_{\rm TO}\simeq M_\odot\left(\frac{10^{10}\,\mathrm{yr}}{t_{\rm cluster}}\right)^{1/p}. MTO≃M⊙(tcluster1010yr)1/p.
- Металличность ZZZ (или [Fe/H][{\rm Fe/H}][Fe/H]): сдвигает изохроны в цвете и яркости. При прочих равных: - более металличные звёзды — краснее и немного менее яркие на MS/TO и имеют более «крутую» и красную красную ветвь гигантов (RGB); - металленно бедные — TO и RGB смещены в синие цвета, горизонтальная ветвь (HB) у шаровых может быть сильно синей. Сопоставляют наблюдаемый цвет/яркость TO и RGB с модельными изохронами разной ZZZ. - Первоначальные массы звёзд: TO даёт текущую максимальную массу ещё живых MS‑звёзд MTOM_{\rm TO}MTO — это оценка массы звезды, которая только что покидает MS. Диапазон первичных масс в кластере: от минимальной (примерно Mmin∼0.08 M⊙M_{\min}\sim 0.08\,M_\odotMmin∼0.08M⊙) до первоначально наиболее массивных звёзд (до \(\sim 50\mbox{–}100\,M_\odot\) в очень молодых кластерах). Для оценки суммарной первоначальной массы кластера используют IMF (например, Сальпетер): ξ(M)∝M−2.35,
\xi(M)\propto M^{-2.35}, ξ(M)∝M−2.35,
и интегрирование по диапазону масс: Mtot=∫MminMmaxM ξ(M) dM.
M_{\rm tot}=\int_{M_{\min}}^{M_{\max}} M\,\xi(M)\,dM. Mtot=∫MminMmaxMξ(M)dM. 2) Ожидаемые CMD‑признаки трёх типов кластеров - Молодой рассеянный кластер (возраст \(\sim 1\mbox{–}100\) Myr): TO очень яркий и синий (O–B типа), присутствуют массивные голубые звёзды и возможно эмиссионные линии/продолжительные пред‑MS ветви. Выводы: возраст маленький, в кластере были звёзды с большими начальными массами (несколько — десятки M⊙M_\odotM⊙), металличность чаще близка к солнечной (если диск). - Открытый кластер среднего возраста (возраст \(\sim 0.1\mbox{–}2\) Gyr): TO находится у A–F типов, MS короче по цвету, появляются субгиганты/молодые красные гиганты. Из положения TO оценивают возраст порядка \(10^8\mbox{–}10^9\) лет; MTOM_{\rm TO}MTO несколько M⊙M_\odotM⊙. Металличность обычно около солнечной или чуть ниже; RGB не очень выражена. - Старый шаровой кластер (возраст \(\sim 10\mbox{–}13\) Gyr): TO слабый и красноватый ( \(M_{\rm TO}\sim 0.8\mbox{–}1.0\,M_\odot\)), яркая и длинная RGB и характерная HB (цвет HB сильно зависит от ZZZ и второго параметра). Металличность низкая ([Fe/H]≲−1[{\rm Fe/H}]\lesssim -1[Fe/H]≲−1 для многих шаровых). Большинство массивных звёзд уже ушли; текущие живые MS‑звезды — маломассивные. 3) Модельный расчёт для примера: «открытый кластер среднего возраста», возьмём tcluster=1t_{\rm cluster}=1tcluster=1 Gyr и предположим Z≈Z⊙Z\approx Z_\odotZ≈Z⊙. - Оценка массы TO (используем p=3.5p=3.5p=3.5): MTO≃M⊙(1010 yr109 yr)1/3.5=M⊙⋅101/3.5≃1.93 M⊙.
M_{\rm TO}\simeq M_\odot\left(\frac{10^{10}\,\mathrm{yr}}{10^{9}\,\mathrm{yr}}\right)^{1/3.5} =M_\odot\cdot 10^{1/3.5}\simeq 1.93\,M_\odot. MTO≃M⊙(109yr1010yr)1/3.5=M⊙⋅101/3.5≃1.93M⊙.
С учётом неопределённости в ppp диапазон будет примерно \(M_{\rm TO}\simeq 1.8\mbox{–}2.2\,M_\odot\). - Значение светимости/температуры: для MS приближённо L∝MαL\propto M^\alphaL∝Mα с \(\alpha\sim 3\mbox{–}4\). Для M≃1.93 M⊙M\simeq 1.93\,M_\odotM≃1.93M⊙ и α=3.5\alpha=3.5α=3.5: \[ \frac{L}{L_\odot}\simeq (1.93)^{3.5}\simeq 10\mbox{–}15, \] что соответствует типичному абсолютному MVM_VMV для TO порядка \(+1\mbox{–}+2\) (зависит от конверсии L→MVL\to M_VL→MV и поглощения). - Оценка доли первоначальной массы в звёздах массой выше MTOM_{\rm TO}MTO: используем Сальпетер (интегралы даю только конечный результат). Доля массы в звёздах с M>MTO=2 M⊙M>M_{\rm TO}=2\,M_\odotM>MTO=2M⊙ при Mmin=0.08 M⊙M_{\min}=0.08\,M_\odotMmin=0.08M⊙ и Mmax=100 M⊙M_{\max}=100\,M_\odotMmax=100M⊙ равна fM=∫2100M−1.35dM∫0.08100M−1.35dM≃0.26.
f_M=\frac{\int_{2}^{100} M^{-1.35}dM}{\int_{0.08}^{100} M^{-1.35}dM}\simeq 0.26. fM=∫0.08100M−1.35dM∫2100M−1.35dM≃0.26.
То есть примерно \(25\mbox{–}30\%\) первоначальной массы приходилось на звёзды более массивные >2 M⊙>2\,M_\odot>2M⊙. Это даёт оценку потерь массы от взрывов сверхновых и ветров: если ныне в кластере видна малая часть этих массивных звёзд, значительная доля массы уже ушла. - Итак, для нашего примера: возраст ∼1\sim 1∼1 Gyr, MTO≃1.9 M⊙M_{\rm TO}\simeq 1.9\,M_\odotMTO≃1.9M⊙, типичная начальная массовая функция — Сальпетер, ~25% первоначальной массы были в звёздах >2 M⊙>2\,M_\odot>2M⊙. Металличность Z≈Z⊙Z\approx Z_\odotZ≈Z⊙ будет давать TO и RGB цвета, совпадающие с солнечными изохронами; при пониженной ZZZ TO был бы чуть светлее и синее при том же возрасте. 4) Практика и уточнения: - Точные оценки возраста и ZZZ получают наложением модельных изохрон (например, PARSEC, MIST) на наблюдаемый CMD; простой t↔MTOt\leftrightarrow M_{\rm TO}t↔MTO даёт порядка и масштаб, но требует уточнения модели (конвекция, повороты, бинарные звёзды, дифференциальное поглощение). - Для суммарной первоначальной массы кластера нужны дополнительные данные: число наблюдаемых звёзд, поле зрения, потери из-за динамики; тогда интегрируют IMF и нормируют на число/массу наблюдаемых звёзд. Если нужно, могу: а) подобрать конкретные изохроны (MIST/PARSEC) и показать наложение на примерный CMD; б) выполнить интегрирование IMF и дать числовую оценку полной первоначальной массы при заданном числе звёзд сейчас.
1) Как по CMD оценивают возраст, металличность и начальные массы:
- Возраст: определяется положением поворотной точки (turn‑off, TO) главной ветви (MS). Чем старше кластер, тем краснее (более низкая TeffT_{\rm eff}Teff ) и слабее (меньшая LLL) TO; соответственно масса звёзд на TO меньше. Для предварительной оценки используют аппроксимацию продолжительности главной последовательности:
tMS≃1010 yr(MM⊙)−p, t_{\rm MS}\simeq 10^{10}\,\mathrm{yr}\left(\frac{M}{M_\odot}\right)^{-p},
tMS ≃1010yr(M⊙ M )−p, где \(p\sim 3\mbox{--}3.5\) (в зависимости от диапазона масс). Обращая её:
MTO≃M⊙(1010 yrtcluster)1/p. M_{\rm TO}\simeq M_\odot\left(\frac{10^{10}\,\mathrm{yr}}{t_{\rm cluster}}\right)^{1/p}.
MTO ≃M⊙ (tcluster 1010yr )1/p. - Металличность ZZZ (или [Fe/H][{\rm Fe/H}][Fe/H]): сдвигает изохроны в цвете и яркости. При прочих равных:
- более металличные звёзды — краснее и немного менее яркие на MS/TO и имеют более «крутую» и красную красную ветвь гигантов (RGB);
- металленно бедные — TO и RGB смещены в синие цвета, горизонтальная ветвь (HB) у шаровых может быть сильно синей.
Сопоставляют наблюдаемый цвет/яркость TO и RGB с модельными изохронами разной ZZZ.
- Первоначальные массы звёзд: TO даёт текущую максимальную массу ещё живых MS‑звёзд MTOM_{\rm TO}MTO — это оценка массы звезды, которая только что покидает MS. Диапазон первичных масс в кластере: от минимальной (примерно Mmin∼0.08 M⊙M_{\min}\sim 0.08\,M_\odotMmin ∼0.08M⊙ ) до первоначально наиболее массивных звёзд (до \(\sim 50\mbox{–}100\,M_\odot\) в очень молодых кластерах). Для оценки суммарной первоначальной массы кластера используют IMF (например, Сальпетер):
ξ(M)∝M−2.35, \xi(M)\propto M^{-2.35},
ξ(M)∝M−2.35, и интегрирование по диапазону масс:
Mtot=∫MminMmaxM ξ(M) dM. M_{\rm tot}=\int_{M_{\min}}^{M_{\max}} M\,\xi(M)\,dM.
Mtot =∫Mmin Mmax Mξ(M)dM.
2) Ожидаемые CMD‑признаки трёх типов кластеров
- Молодой рассеянный кластер (возраст \(\sim 1\mbox{–}100\) Myr): TO очень яркий и синий (O–B типа), присутствуют массивные голубые звёзды и возможно эмиссионные линии/продолжительные пред‑MS ветви. Выводы: возраст маленький, в кластере были звёзды с большими начальными массами (несколько — десятки M⊙M_\odotM⊙ ), металличность чаще близка к солнечной (если диск).
- Открытый кластер среднего возраста (возраст \(\sim 0.1\mbox{–}2\) Gyr): TO находится у A–F типов, MS короче по цвету, появляются субгиганты/молодые красные гиганты. Из положения TO оценивают возраст порядка \(10^8\mbox{–}10^9\) лет; MTOM_{\rm TO}MTO несколько M⊙M_\odotM⊙ . Металличность обычно около солнечной или чуть ниже; RGB не очень выражена.
- Старый шаровой кластер (возраст \(\sim 10\mbox{–}13\) Gyr): TO слабый и красноватый ( \(M_{\rm TO}\sim 0.8\mbox{–}1.0\,M_\odot\)), яркая и длинная RGB и характерная HB (цвет HB сильно зависит от ZZZ и второго параметра). Металличность низкая ([Fe/H]≲−1[{\rm Fe/H}]\lesssim -1[Fe/H]≲−1 для многих шаровых). Большинство массивных звёзд уже ушли; текущие живые MS‑звезды — маломассивные.
3) Модельный расчёт для примера: «открытый кластер среднего возраста», возьмём tcluster=1t_{\rm cluster}=1tcluster =1 Gyr и предположим Z≈Z⊙Z\approx Z_\odotZ≈Z⊙ .
- Оценка массы TO (используем p=3.5p=3.5p=3.5):
MTO≃M⊙(1010 yr109 yr)1/3.5=M⊙⋅101/3.5≃1.93 M⊙. M_{\rm TO}\simeq M_\odot\left(\frac{10^{10}\,\mathrm{yr}}{10^{9}\,\mathrm{yr}}\right)^{1/3.5}
=M_\odot\cdot 10^{1/3.5}\simeq 1.93\,M_\odot.
MTO ≃M⊙ (109yr1010yr )1/3.5=M⊙ ⋅101/3.5≃1.93M⊙ . С учётом неопределённости в ppp диапазон будет примерно \(M_{\rm TO}\simeq 1.8\mbox{–}2.2\,M_\odot\).
- Значение светимости/температуры: для MS приближённо L∝MαL\propto M^\alphaL∝Mα с \(\alpha\sim 3\mbox{–}4\). Для M≃1.93 M⊙M\simeq 1.93\,M_\odotM≃1.93M⊙ и α=3.5\alpha=3.5α=3.5:
\[
\frac{L}{L_\odot}\simeq (1.93)^{3.5}\simeq 10\mbox{–}15,
\]
что соответствует типичному абсолютному MVM_VMV для TO порядка \(+1\mbox{–}+2\) (зависит от конверсии L→MVL\to M_VL→MV и поглощения).
- Оценка доли первоначальной массы в звёздах массой выше MTOM_{\rm TO}MTO : используем Сальпетер (интегралы даю только конечный результат). Доля массы в звёздах с M>MTO=2 M⊙M>M_{\rm TO}=2\,M_\odotM>MTO =2M⊙ при Mmin=0.08 M⊙M_{\min}=0.08\,M_\odotMmin =0.08M⊙ и Mmax=100 M⊙M_{\max}=100\,M_\odotMmax =100M⊙ равна
fM=∫2100M−1.35dM∫0.08100M−1.35dM≃0.26. f_M=\frac{\int_{2}^{100} M^{-1.35}dM}{\int_{0.08}^{100} M^{-1.35}dM}\simeq 0.26.
fM =∫0.08100 M−1.35dM∫2100 M−1.35dM ≃0.26. То есть примерно \(25\mbox{–}30\%\) первоначальной массы приходилось на звёзды более массивные >2 M⊙>2\,M_\odot>2M⊙ . Это даёт оценку потерь массы от взрывов сверхновых и ветров: если ныне в кластере видна малая часть этих массивных звёзд, значительная доля массы уже ушла.
- Итак, для нашего примера: возраст ∼1\sim 1∼1 Gyr, MTO≃1.9 M⊙M_{\rm TO}\simeq 1.9\,M_\odotMTO ≃1.9M⊙ , типичная начальная массовая функция — Сальпетер, ~25% первоначальной массы были в звёздах >2 M⊙>2\,M_\odot>2M⊙ . Металличность Z≈Z⊙Z\approx Z_\odotZ≈Z⊙ будет давать TO и RGB цвета, совпадающие с солнечными изохронами; при пониженной ZZZ TO был бы чуть светлее и синее при том же возрасте.
4) Практика и уточнения:
- Точные оценки возраста и ZZZ получают наложением модельных изохрон (например, PARSEC, MIST) на наблюдаемый CMD; простой t↔MTOt\leftrightarrow M_{\rm TO}t↔MTO даёт порядка и масштаб, но требует уточнения модели (конвекция, повороты, бинарные звёзды, дифференциальное поглощение).
- Для суммарной первоначальной массы кластера нужны дополнительные данные: число наблюдаемых звёзд, поле зрения, потери из-за динамики; тогда интегрируют IMF и нормируют на число/массу наблюдаемых звёзд.
Если нужно, могу: а) подобрать конкретные изохроны (MIST/PARSEC) и показать наложение на примерный CMD; б) выполнить интегрирование IMF и дать числовую оценку полной первоначальной массы при заданном числе звёзд сейчас.