Модели долгопериодических переменных звезд шарового скопления 47 Тuc

Основная информация

• ID статьи: 2511.09065
• Название: Models of long-period variables of the globular cluster 47 Tuc
• Автор: Ю. А. Фадеев (Институт астрономии РАН)
• Классификация: astro-ph.SR (звёздная и солнечная физика)
• Дата подачи: 12 ноября 2025 г. на arXiv
• Ссылка на статью: https://arxiv.org/abs/2511.09065

Аннотация

В данном исследовании проведены эволюционные расчёты звёзд с массой на главной последовательности $M_\mathrm{ZAMS}=0.86M_\odot$ и начальной металличностью $Z=0.003$ и $Z=0.004$. Модели, отобранные из эволюционных последовательностей, использованы для расчёта радиальных колебаний на этапах RGB (красный гигант), eAGB (ранняя асимптотическая ветвь гигантов) и TP-AGB (асимптотическая ветвь гигантов с тепловыми импульсами). Исследование показывает, что пульсирующие красные гиганты в шаровом скоплении 47 Tuc не все являются переменными типа Миры: нижний предел периода колебаний на этапе TP-AGB составляет примерно 70 дней, тогда как диапазон периодов радиальных колебаний на этапе eAGB составляет примерно 5–40 дней. Периоды и светимости гидродинамических моделей пульсирующих звёзд eAGB и TP-AGB распределяются вдоль общего соотношения период–светимость. Малая масса переменных типа Миры в 47 Tuc ($0.54M_\odot\le M\le 0.70M_\odot$) является основной причиной появления нерегулярных колебаний с большой амплитудой и динамической неустойчивости внешних слоёв при периодах $\Pi > 200$ дней.

Научный контекст и мотивация

Исследуемая проблема

Данное исследование направлено на определение условий возбуждения радиальных колебаний красных гигантов в шаровом скоплении 47 Tuc на различных этапах эволюции (RGB, eAGB, TP-AGB) и уточнение эволюционного статуса этих долгопериодических переменных звёзд.

Значимость проблемы

1. Особенности наблюдений: 47 Tuc (NGC 104) — шаровое скопление с аномально высокой металличностью, превышающей металличность большинства шаровых скоплений почти на порядок величины. Это привело к обнаружению большого количества пульсирующих красных гигантов (из 117 долгопериодических переменных, обнаруженных в 102 шаровых скоплениях к 2001 году, 14 находятся в 47 Tuc)
2. Конденсация пыли: инфракрасные наблюдения показывают конденсацию пылевых частиц во внешних атмосферах долгопериодических переменных 47 Tuc, что тесно связано с периодическими ударными волнами, генерируемыми колебаниями звёзд с большой амплитудой
3. Неопределённость эволюционного статуса: существуют разногласия в литературе относительно эволюционного статуса этих пульсирующих красных гигантов, в частности, являются ли они все переменными типа Миры

Ограничения существующих методов

1. Недостаточность линейного анализа пульсаций: Лебцельтер и Вуд (2005) использовали линейные модели пульсаций для звёзд с $Z=0.004$, $M=0.9M_\odot$, но не предоставили убедительных доказательств
2. Смещение в оценке металличности: ранние исследования использовали значения металличности, значительно превышающие современные оценки наблюдений
3. Отсутствие согласованных расчётов: недостаёт согласованных расчётов, объединяющих эволюцию звёзд и нелинейные пульсации

Мотивация исследования

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

Основные вклады

1. Разработана полная система расчётов эволюция–пульсация: использован код MESA (версия r15140) для расчётов эволюции звёзд в сочетании с радиационной гидродинамикой и уравнениями нестационарной конвекции для расчёта пульсаций звёзд
2. Определены характеристики пульсаций на различных этапах эволюции:
   • Этап RGB: область неустойчивости пульсаций узкая или отсутствует (при $Z=0.003$)
   • Этап eAGB: диапазон периодов 5–40 дней, относительно регулярные амплитуды
   • Этап TP-AGB: нижний предел периода примерно 70 дней, нерегулярные колебания с большой амплитудой при длинных периодах (>200 дней)
3. Определён основной эволюционный статус долгопериодических переменных 47 Tuc: большинство пульсирующих красных гигантов находятся на этапе eAGB, а не TP-AGB
4. Установлено соотношение период–светимость: данные период–светимость звёзд eAGB и TP-AGB распределяются вдоль общей линии регрессии
5. Объяснен механизм динамической неустойчивости: малая масса ($0.54–0.70M_\odot$) является основной причиной неустойчивости внешних слоёв при длинных периодах

Подробное описание методов

Определение задачи

Входные данные: начальные параметры звёзд (масса на главной последовательности, металличность, параметр скорости потери массы)
Выходные данные: характеристики пульсаций на различных этапах эволюции (периоды, амплитуды, устойчивость)
Ограничения: соответствие наблюдательным ограничениям для 47 Tuc (возраст 11,8–12,4 млрд лет, металличность $-0.78 \le [Fe/H] \le -0.66$)

Архитектура расчётов

1. Расчёты эволюции звёзд

• Код: MESA версия r15140
• Начальные параметры:
  • Масса на главной последовательности: $M_\mathrm{ZAMS} = 0.86M_\odot$ (соответствует возрасту на этапе AGB примерно 11,9 млрд лет)
  • Металличность: $Z = 0.003$ (соответствует $[Fe/H] = -0.67$) и $Z = 0.004$ (соответствует $[Fe/H] = -0.54$)
  • Параметр скорости потери массы (формула Реймерса): $\eta_R = 0.3$ и $\eta_R = 0.5$
• Эволюционные последовательности: всего рассчитано 4 эволюционные последовательности (2 значения металличности × 2 параметра потери массы)

2. Расчёты гидродинамических пульсаций

• Система уравнений: уравнения радиационной гидродинамики + уравнения нестационарной конвекции
• Начальные условия: модели в гидростатическом и тепловом равновесии, отобранные из эволюционных последовательностей
• Метод решения: численное решение задачи Коши
• Анализируемое содержание:
  • Идентификация мод колебаний (основная мода/обертоны)
  • Анализ удельной механической работы: $\oint PdV$ (P — полное давление, V — удельный объём)
  • Эволюция амплитуды и определение устойчивости

Технические инновации

1. Метод согласованных расчётов

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

2. Исследование многопараметрического пространства

Систематическое исследование пространства параметров металличности и скорости потери массы для оценки влияния неопределённостей на результаты.

3. Анализ физических механизмов

Посредством радиального распределения удельной механической работы $\oint PdV$ количественно анализируется вклад областей ионизации водорода и гелия в возбуждение/затухание пульсаций:

• Область частичной ионизации водорода (адиабатический индекс $\Gamma_1 < 4/3$) является основной областью возбуждения
• Положение узлов обертонов определяет моду пульсаций

4. Механизм неустойчивости при длинных периодах

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

Экспериментальная установка

Конфигурация параметров

Параметры эволюции звёзд

• Масса на главной последовательности: $0.86M_\odot$
• Металличность: $Z = 0.003$ (солнечное значение $Z_\odot = 0.014$) и $Z = 0.004$
• Скорость потери массы: формула Реймерса, $\eta_R = 0.3, 0.5$
• Этапы эволюции: вершина RGB, eAGB, TP-AGB (4 тепловых импульса)

Гидродинамические модели

• Пространственная дискретизация: примерно 600–700 сферических слоёв
• Моды пульсаций: основная мода и первый обертон
• Критерии:
  • Предельный цикл: периодические колебания с устойчивой амплитудой
  • Затухающие колебания: амплитуда постепенно уменьшается до нуля
  • Нерегулярные колебания: $\Delta R/R \approx 1$

Показатели оценки

1. Период пульсаций $\Pi$ (дни)
2. Постоянная пульсаций $Q = \Pi\sqrt{\bar{\rho}/\rho_\odot}$ (дни, $\bar{\rho}$ — средняя плотность)
3. Амплитуда относительного радиального смещения $\Delta R/R$
4. Продолжительность неустойчивости пульсаций $\Delta t_\mathrm{ev}$ (миллионы лет)
5. Светимость $\log(L/L_\odot)$

Базы для сравнения

• Результаты линейных пульсаций Лебцельтера и Вуда (2005)
• Наблюдаемый нижний предел периода переменных типа Миры в Галактике (Samus et al. 2017)
• Распределение периодов долгопериодических переменных, наблюдаемых в 47 Tuc

Результаты экспериментов

Модели этапа RGB

Основные результаты (таблица 1)

Ключевые выводы:

1. Модель с $Z=0.003, \eta_R=0.3$ не пульсирует (область неустойчивости пульсаций отсутствует)
2. Увеличение металличности и параметра потери массы приводит к смещению эволюционной траектории к более низким эффективным температурам
3. Более высокая $Z$ и меньшая масса делают звезду более неустойчивой к радиальным колебаниям (более низкая плотность внешней конвективной зоны, более широкая область ионизации)
4. Продолжительность пульсаций вблизи вершины RGB очень короткая (<1 млн лет)

Физический механизм (рисунок 2)

• Колебания первого обертона (узел обертона $j \approx 318$, радиус $r \approx 0.81\bar{R}$)
• Область частичной ионизации водорода ($\Gamma_1 < 4/3$) обеспечивает основную движущую силу
• Для $Z=0.004, \eta_R=0.5$ область ионизации водорода расширяется достаточно глубоко, чтобы поддерживать колебания основной моды

Модели этапа eAGB

Характеристики эволюции (рисунки 3–4)

1. Последовательность появления колебаний: сначала появляются колебания первого обертона, затем с увеличением радиуса и светимости звезды происходит переход к основной моде
2. Диапазон периодов: 5–40 дней
3. Характеристики амплитуды: $\Delta R/R \lesssim 0.3$, колебания относительно регулярны
4. Верхний предел светимости: пульсации прекращаются при $L \gtrsim 10^3 L_\odot$
5. Продолжительность: примерно $1.2 \times 10^7$ лет (намного дольше, чем на этапе TP-AGB)

Модели этапа TP-AGB

Характеристики тепловых импульсов (рисунок 5)

• Диапазон масс: при первом тепловом импульсе $0.59–0.71M_\odot$, в начале этапа после-AGB примерно $0.54M_\odot$
• Интервал между тепловыми импульсами: в среднем $1.6 \times 10^5$ лет
• Количество тепловых импульсов: 4 (общая продолжительность этапа TP-AGB $5–6 \times 10^5$ лет)
• Условия пульсаций: пульсации появляются только при $\log(L/L_\odot) \ge 3.4$, что означает, что в течение значительной части каждого интервала между импульсами звезда не пульсирует

Периоды и устойчивость

1. Диапазон периодов: 70–240 дней
2. Порог нерегулярных колебаний: при $\Pi \approx 170$ дней происходит переход к нерегулярным колебаниям с большой амплитудой ($\Delta R/R \approx 1$)
3. Нижний предел периода: примерно 70 дней (согласуется с переменными типа Миры в Галактике)

Соотношение период–светимость (рисунок 6)

Для эволюционной последовательности с $\eta_R = 0.3$ модели eAGB и TP-AGB распределяются вдоль общей линии регрессии:

$\log(L/L_\odot) = 0.743 \log \Pi + 1.79 \quad (Z = 0.003)$

$\log(L/L_\odot) = 0.724 \log \Pi + 1.78 \quad (Z = 0.004)$

• Стандартная ошибка коэффициентов < 3%
• При $\log \Pi \gtrsim 2$ разброс увеличивается, что объясняется эволюционным ростом массы углеродного ядра и светимости звезды

Ключевые экспериментальные результаты

1. Определение эволюционного статуса: большинство пульсирующих красных гигантов в 47 Tuc находятся на этапе eAGB, а не TP-AGB
2. Граница периодов:
   • Период < 70 дней → звёзды eAGB
   • Период ≥ 70 дней → возможно звёзды TP-AGB
3. Эффект массы: малая масса ($0.54–0.70M_\odot$) является основной причиной динамической неустойчивости внешних слоёв при длинных периодах (>200 дней)
4. Механизм конденсации пыли: периодические ударные волны, генерируемые колебаниями с большой амплитудой, являются основным механизмом конденсации пылевых частиц

Связанные работы

Наблюдательные исследования

1. Clement et al. (2001): обнаружено 117 долгопериодических переменных в 102 шаровых скоплениях, 14 из них в 47 Tuc
2. Lebzelter & Wood (2005), Lebzelter et al. (2005): обнаружено больше долгопериодических переменных в 47 Tuc, общее количество удвоилось
3. Percy & Gupta (2021): продолжающиеся наблюдения с обнаружением новых переменных
4. Origlia et al. (2002), McDonald et al. (2011): инфракрасные наблюдения, обнаружены свидетельства конденсации пыли

Теоретические исследования

1. Lebzelter & Wood (2005): анализ линейных пульсаций, использование $Z=0.004, M=0.9M_\odot$, но выводы недостаточны
2. Willson (2000): предложен физический механизм конденсации пыли, вызванной периодическими ударными волнами
3. Fadeyev (2013, 2023, 2024): предыдущие работы автора, разработана система уравнений радиационной гидродинамики и нестационарной конвекции

Преимущества данной работы

1. Более точные параметры: использование современных оценок металличности ($Z=0.003$ соответствует $[Fe/H]=-0.67$)
2. Нелинейные расчёты: применение полных уравнений гидродинамики вместо линейного приближения
3. Согласованная система: бесшовное соединение расчётов эволюции и пульсаций
4. Систематический анализ: охват трёх этапов эволюции (RGB, eAGB, TP-AGB)

Выводы и обсуждение

Основные выводы

1. Эволюционный статус: большинство пульсирующих красных гигантов в 47 Tuc находятся на этапе eAGB, причины:
   • Продолжительность eAGB (примерно $1.2 \times 10^7$ лет) в два раза превышает продолжительность TP-AGB (примерно $6 \times 10^5$ лет)
   • На этапе TP-AGB небольшое количество тепловых импульсов (≤4)
   • В каждом интервале между импульсами звезда пульсирует только часть времени
2. Распределение периодов:
   • TP-AGB: нижний предел периода примерно 70 дней (согласуется с переменными типа Миры в Галактике)
   • eAGB: 5–40 дней
   • Вершина RGB: период ≥ 60–70 дней, но область неустойчивости узкая или отсутствует
3. Короткопериодические переменные: красные переменные 47 Tuc с периодом < 70 дней следует классифицировать как звёзды eAGB, а не RGB или TP-AGB
4. Динамическая неустойчивость: при периодах > 200 дней внешние слои динамически неустойчивы, периодические ударные волны являются основным механизмом конденсации пыли

Ограничения

1. Неопределённость потери массы: хотя рассмотрены два значения $\eta_R$, фактический механизм потери массы может быть более сложным
2. Диапазон металличности: $Z=0.004$ выходит за верхний предел наблюдаемого диапазона для 47 Tuc, используется в основном для сравнения с предыдущими работами
3. Обработка конвекции: модель нестационарной конвекции содержит теоретические неопределённости
4. Количество тепловых импульсов: смоделировано только 4 тепловых импульса, что может быть недостаточно для полного описания этапа TP-AGB

Направления будущих исследований

1. Расширение на большее количество циклов тепловых импульсов для исследования поздней эволюции TP-AGB
2. Улучшение формулы скорости потери массы с учётом усиленной потери массы, вызванной пульсациями
3. Трёхмерное моделирование гидродинамики для более точного описания конвекции и распространения ударных волн
4. Статистическое сравнение с наблюдательными данными для проверки предсказаний модели

Глубокая оценка

Достоинства

1. Инновационность методов

• Система согласованных расчётов: впервые систематически объединены расчёты эволюции MESA с нелинейными расчётами гидродинамических пульсаций, что исключает несоответствие параметров
• Анализ физических механизмов: количественный анализ областей возбуждения/затухания через $\oint PdV$ с ясной физической картиной

2. Достаточность экспериментов

• Охват параметрического пространства: 2 значения металличности × 2 параметра потери массы = 4 эволюционные последовательности
• Полнота этапов эволюции: систематическое исследование трёх этапов (RGB, eAGB, TP-AGB)
• Подробность данных: таблицы и рисунки полно представляют результаты расчётов

3. Убедительность результатов

• Разрешение наблюдательных противоречий: разумное объяснение сосуществования переменных с различными периодами в 47 Tuc
• Количественные предсказания: чёткие границы периодов (70 дней) и соотношение период–светимость
• Физическая согласованность: ясная причинно-следственная цепь от малой массы к низкой плотности к широкой области ионизации к неустойчивости

4. Ясность изложения

• Строгая структура, логичное построение
• Информативные рисунки и таблицы с чёткими обозначениями
• Глубокое и доступное физическое объяснение

Недостатки

1. Ограничения методов

• Единственная масса: рассчитана только $M_\mathrm{ZAMS}=0.86M_\odot$, не исследовано влияние массы на результаты
• Модель конвекции: теория нестационарной конвекции содержит неопределённости, которые могут влиять на предсказания амплитуды
• Упрощение переноса излучения: возможны упрощающие предположения в переносе излучения при расчётах гидродинамики

2. Экспериментальная установка

• Недостаточное покрытие TP-AGB: только 4 тепловых импульса, возможно упущены особенности поздней эволюции
• Отсутствие статистического сравнения: нет количественной статистической проверки соответствия наблюдаемому распределению периодов в 47 Tuc
• Обоснованность $Z=0.004$: эта металличность выходит за пределы наблюдаемого диапазона, что снижает применимость части выводов

3. Глубина анализа

• Разброс в соотношении период–светимость: хотя упомянуто, что рост массы углеродного ядра вызывает разброс, количественный анализ отсутствует
• Механизм нерегулярных колебаний: физический механизм нерегулярности при $\Pi > 200$ дней описан недостаточно подробно
• Конденсация пыли: упомянут механизм ударных волн, но конкретные условия образования пыли не рассчитаны

Влияние

1. Вклад в область

• Теоретический прорыв: уточнён эволюционный статус долгопериодических переменных 47 Tuc, разрешены давние споры
• Методологический образец: система согласованных расчётов эволюция–пульсация может быть распространена на другие звёздные скопления и типы переменных
• Руководство для наблюдений: граница периодов служит теоретической основой для будущей классификации переменных

2. Практическая ценность

• Классификация переменных: предоставляет инструмент для классификации переменных в 47 Tuc и аналогичных скоплениях с высокой металличностью
• Определение возраста: соотношение период–светимость может использоваться для оценки расстояния и возраста скопления
• Исследование образования пыли: обеспечивает физическую основу для понимания образования пыли в шаровых скоплениях

3. Воспроизводимость

• Использование открытого кода MESA (номер версии указан)
• Подробное описание параметров
• Однако код гидродинамики не открыт, что ограничивает полную воспроизводимость

Области применения

Прямое применение

1. Исследование долгопериодических переменных в скоплениях с высокой металличностью (например, $\omega$ Cen)
2. Исследование эволюции и пульсаций звёзд малой массы ($\sim 0.86M_\odot$)
3. Наблюдательное различие звёзд eAGB и TP-AGB

Требующие адаптации

1. Звёзды другой массы или металличности: необходимо пересчитать эволюционные последовательности
2. Более молодые или старые скопления: начальная масса и временные масштабы эволюции различны
3. Звёзды поля (не в скоплениях): скорость потери массы может отличаться

Неприменимые случаи

1. Скопления с низкой металличностью ($[Fe/H] < -1.5$): механизм пульсаций может быть иным
2. Звёзды большой массы ($M > 2M_\odot$): эволюционный путь совершенно иной
3. Короткопериодические переменные (например, RR Lyrae): другая полоса неустойчивости

Ключевые ссылки

1. Paxton et al. (2019): Код MESA, Astrophys. J. Suppl. Ser. 243, 10
2. Fadeyev (2013): Уравнения гидродинамики, Astron. Lett. 39, 306
3. Lebzelter & Wood (2005): Анализ линейных пульсаций, Astron. Astrophys. 441, 1117
4. McDonald et al. (2011): Инфракрасные наблюдения и пыль, Astrophys. J. Suppl. Ser. 193, 23
5. Willson (2000): Механизм ударных волн и конденсации пыли, Ann. Rev. Astron. Astrophys. 38, 573

Общая оценка: Это высокого качества теоретическое исследование в области физики звёзд, которое путём строгих численных расчётов разрешает давний спор об эволюционном статусе долгопериодических переменных 47 Tuc. Методологическая инновация (система согласованных расчётов эволюция–пульсация), чёткие выводы (большинство — звёзды eAGB) и ясная физическая картина (малая масса вызывает неустойчивость). Основные недостатки — расчёты для единственной массы и отсутствие статистического сравнения с наблюдениями. Данная работа имеет важное теоретическое значение для исследования переменных звёзд в шаровых скоплениях, а методология может быть распространена на другие астрофизические задачи.