Солнце поглотило собственную суперземлю: как гибель планеты в пять раз массивнее Земли решила 30-летнюю проблему физики звезд

Солнце — наиболее изученная звезда во Вселенной. Астрофизики с высокой точностью знают его массу, радиус, светимость и возраст, который составляет 4,57 миллиарда лет. На основе этих фундаментальных параметров построена так называемая стандартная модель Солнца. Это система математических уравнений, описывающая физическое состояние звездного вещества от центра до внешних слоев атмосферы. Она учитывает гравитацию, давление газа, термоядерные реакции, а также перенос тепла излучением и механическим перемешиванием.

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

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

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

Кроме того, стандартная модель дает сбой еще по трем важным параметрам:

1. Положение границы конвективной зоны. Расчетный радиус основания этой зоны оказывается больше, чем показывают данные измерений, которые фиксируют эту границу строго на отметке 0,713 от полного радиуса Солнца.
2. Поверхностный химический состав. Долгое время считалось, что Солнце содержит около двух процентов элементов тяжелее водорода и гелия (в астрофизике все такие элементы называют металлами). Однако уточненный спектральный анализ, проведенный в 2000-х и 2020-х годах, снизил эту величину примерно до полутора процентов. Когда новую, пониженную металличность заложили в стандартную модель, расхождение по скорости звука увеличилось еще сильнее.
3. Дефицит лития. В наружных слоях Солнца содержание лития более чем в сто раз ниже, чем в веществе первичных метеоритов, сформировавшихся одновременно со всей планетной системой.

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

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

В исследовании, опубликованном в журнале Monthly Notices of the Royal Astronomical Society, профессор астрономии Мутлу Йылдыз из Эгейского университета доказал: все указанные несоответствия полностью разрешаются, если учесть гравитационное взаимодействие молодого Солнца с его собственным протопланетным диском и факт поглощения крупной каменистой планеты на раннем этапе эволюции.

Механика солнечных оболочек и природа расхождения

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

Энергия, рожденная в термоядерном ядре, движется к поверхности двумя путями.

От центра звезды до расстояния примерно в 71 процент от ее полного радиуса простирается лучистая зона. Здесь плотность вещества чрезвычайно высока, а плазма достаточно прозрачна для того, чтобы кванты электромагнитного излучения постепенно уносили тепло наружу. В лучистой зоне нет механических перемещений вещества: слои плазмы неподвижны относительно друг друга, а химические элементы могут лишь медленно оседать вниз под действием гравитации.

Выше отметки в 71 процент радиуса температура плазмы снижается до значений, при которых электроны начинают связываться с атомными ядрами. Газ становится менее прозрачным для излучения. Лучистый перенос тепла перестает справляться с потоком энергии из ядра. В результате вещество начинает физически перемешиваться: горячие массы плазмы устремляются вверх, отдают тепло в наружных слоях, охлаждаются и опускаются обратно вниз. Возникает конвективная зона, которая занимает внешние тридцать процентов радиуса Солнца. За счет постоянного движения вещество в этой оболочке непрерывно перемешивается до полной однородности.

Переходный слой между лучистой и конвективной зонами называется тахоклином. Именно в нем концентрируется наибольшая ошибка скорости звука стандартной модели.

Скорость распространения звука в ионизированном газе напрямую связана с его термодинамическими свойствами: она пропорциональна квадратному корню из отношения давления к плотности вещества, а также зависит от коэффициента сжатия газа при быстром изменении объема. Чем выше давление при фиксированной плотности, тем быстрее распространяются акустические колебания.

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

Увеличить скорость звука в этой области можно двумя путями: либо снизить средний молекулярный вес плазмы (например, локально уменьшив долю гелия), либо повысить температуру и давление. Однако эти величины жестко связаны со способностью плазмы задерживать тепловое излучение, то есть с ее непрозрачностью.

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

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

Однако в стабильной лучистой зоне нет внутренних физических механизмов, способных сформировать изолированный слой с высокой концентрацией металлов. Конвективное перемешивание и естественная диффузия способны лишь сглаживать химические различия, но не создавать локальные скопления элементов. Единственный физически возможный источник такого избытка тяжелых веществ — внешнее воздействие на этапе формирования звезды.

Куда исчезли внутренние планеты

Наблюдения за внесолнечными планетными системами предоставили фактический материал, которого не было у создателей первых звездных моделей.

В современных астрономических каталогах преобладают так называемые горячие суперземли — каменистые тела с массами от двух до десяти масс Земли, а также близкие к звездам газовые нептуны. Большинство обнаруженных планетных систем содержат крупные тела на расстояниях менее одной десятой астрономической единицы от родительских светил (для сравнения: расстояние от Солнца до ближайшей к нему планеты, Меркурия, составляет почти четыре десятых астрономической единицы).

Внутренняя часть Солнечной системы, напротив, полностью свободна от массивных объектов.

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

В процессе роста планеты взаимодействуют с окружающим газом диска. Гравитационное сопротивление газовой среды тормозит планету, она теряет момент импульса и по спиральной траектории смещается к центральной звезде. Этот процесс называют орбитальной миграцией.

Мутлу Йылдыз рассчитал сценарий, при котором на раннем этапе эволюции Солнечной системы сформировавшаяся каменистая суперземля сошла со стабильной орбиты и упала на молодое Солнце. Автор присвоил этой гипотетической планете имя Дев Дилек (в переводе с турецкого — «Великое Желание» в исторической традиции Ближнего Востока это имя связывалось с планетой Меркурий).

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

Физика погружения: приливные силы, сжатие и абляция

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

Расчеты показали, что падение крупного каменистого тела происходит иначе, чем падение малых комет или астероидов.

1. Предел приливного разрушения

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

Если к звезде приближается газовый гигант, его протяженная и неплотная внешняя оболочка начинает разрушаться задолго до контакта со звездной поверхностью. Газ срывается приливом и равномерно распределяется по внешней атмосфере светила.

Однако плотное каменистое или железо-силикатное ядро ведет себя противоположным образом. Плотность минерального ядра суперземли достигает величин порядка 8-10 граммов на кубический сантиметр. Расчеты показали, что радиус каменистого ядра суперземли оказывается меньше критического радиуса приливного разрушения вплоть до момента его входа в фотосферу Солнца. Это означает, что твердое тело массой от пяти до десяти масс Земли пересекает видимую границу звезды целым, не распадаясь на фрагменты.

2. Сжатие веществом звезды

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

Это давление со всех сторон сжимает падающее каменистое тело. Расчет уравнения состояния силикатов и железа показал заметное механическое сокращение размеров:

• Радиус планеты массой в десять масс Земли в процессе спуска уменьшается на 26 процентов.
• Радиус планеты массой в шесть масс Земли уменьшается на 29 процентов.

3. Подавление теплового испарения

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

Скорость движения планеты в окрестностях солнечной поверхности близка ко второй космической скорости и составляет около 450 километров в секунду. На такой скорости путь сквозь всю двухсоттысячекилометровую толщу конвективной зоны занимает всего несколько часов.

За столь короткое время тепловой поток от ударной волны не успевает прогреть массивное тело насквозь. Кроме того, механическое сжатие уменьшает радиус планеты, из-за чего площадь ее поперечного сечения падает почти в два раза. Уменьшение площади сечения пропорционально снижает как аэродинамическое сопротивление встречного газа, так и передачу тепла к поверхности планеты.

В результате за время быстрого движения через конвективную зону планета теряет за счет испарения ничтожную долю своего вещества.

4. Остановка и смешивание

Падение продолжается до тех пор, пока плотность окружающей солнечной плазмы не сравнивается со средней плотностью самого сжатого планетного ядра. По расчетам, плотность плазмы достигает значений твердых минералов именно в области тахоклина — на глубине, где сосредоточено от 96 до 97,3 процента массы Солнца.

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

Тяжелые элементы, составлявшие основу планеты, оказываются заперты в стабильном слое непосредственно под подвижной зоной конвекции.

Двухэтапное формирование химического состава

Растворение суперземли обогащает тахоклин железом, кремнием и никелем. Это локально увеличивает непрозрачность плазмы и выравнивает расчетную скорость звука со значениями, измеренными гелиосейсмологами.

Однако возникает новое противоречие: если в недра Солнца попало массивное каменистое тело, почему спектральные измерения солнечной поверхности фиксируют низкую концентрацию тяжелых элементов?

Ответ дает естественная последовательность эволюции протопланетного диска, состоящая из двух разнесенных во времени этапов:

1. Ранний этап (возраст системы от 15 до 30 миллионов лет). В протопланетном диске происходит разделение пыли и газа. Пылевые частицы слипаются в крупные твердые тела, концентрируя в себе практически все тяжелые химические элементы. Образовавшаяся суперземля мигрирует к центру и падает в Солнце, доставляя тугоплавкие металлы под основание конвективной зоны.
2. Поздний этап (возраст системы от 50 до 60 миллионов лет). В окружающем пространстве остается чистый остаточный газ, лишенный пыли. Молодое Солнце притягивает этот обедненный металлами водородно-гелиевый объем на свою поверхность. Масса этого дополнительного газа составляет от одного до трех процентов общей массы Солнца.

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

В результате формируется устойчивое слоистое распределение:

• Поверхность Солнца и вся его внешняя конвективная оболочка состоят из вещества с умеренной металличностью (около 1,65 процента), что точно соответствует современным спектроскопическим измерениям фотосферы.
• Под конвективной зоной, в тахоклине, сохраняется повышенная концентрация металлов, оставшаяся от разрушения планеты. Этот слой обеспечивает требуемый температурный градиент и полностью устраняет аномалию скорости звука.
• Расчетная граница конвекции сдвигается точно на измеренную отметку 0,713 радиуса Солнца.

Эта же схема объясняет парадокс лития.

Ядра лития термодинамически нестабильны при температурах выше двух с половиной миллионов градусов: они разрушаются при столкновении с протонами, превращаясь в гелий. Внутри солнечной конвективной зоны температура недостаточно высока для этой реакции. Чтобы концентрация лития упала до современных значений, вещество из зоны конвекции должно медленно переноситься турбулентными потоками чуть глубже — в верхние горячие слои лучистой зоны.

В модели Йылдыза глубина турбулентного перемешивания строго ограничена: она составляет около 3,6 процента радиуса Солнца ниже конвективной границы. Этого расстояния достаточно, чтобы перемешивать поверхностный литий в зону его выгорания, но недостаточно, чтобы разрушить слой тяжелых элементов, сформированный погибшей суперземлей.

Математическое подтверждение сценария

Чтобы исключить вероятность случайного совпадения параметров, автор применил к расчетам статистический метод — Байесовский информационный критерий (BIC), предложенный математиком Гидеоном Шварцем.

Этот метод предназначен для выбора лучшей физической теории среди нескольких конкурирующих гипотез. Критерий рассчитывается как сумма двух слагаемых:

• Первое слагаемое — это математическая сумма квадратов отклонений теоретических расчетов от реальных экспериментальных измерений во всех контролируемых точках. Чем точнее расчет воспроизводит наблюдения, тем меньше эта сумма.
• Второе слагаемое — это штрафной коэффициент за математическую сложность модели. Он вычисляется как количество настраиваемых параметров модели, умноженное на натуральный логарифм от числа использованных экспериментальных измерений.

Если исследователь вводит в модель лишние свободные параметры ради искусственной подгонки результата, второе слагаемое резко возрастает. Предпочтительной признается та теория, которая дает наименьшее итоговое значение информационного критерия.

В исследовании сравнивались 74 независимых измерения (профиль скорости звука на 71 фиксированной глубине внутри звезды, а также радиус конвективной зоны, доля поверхностного гелия и поверхностная металличность):

1. Модель TurMixaV — классическая модель без падения планеты, где авторы пытались устранить расхождения исключительно оптимизацией параметров конвекции и турбулентности (использовано 6 свободных параметров).
2. Модель PureENGF — расчет, учитывающий только падение планеты, но без дополнительного турбулентного перемешивания (5 свободных параметров).
3. Модель DD1020 — комплексный расчет, объединяющий раннее падение планеты с калиброванным турбулентным переносом вещества в тахоклине (10 свободных параметров).

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

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

Для модели DD1020 масса растворившегося твердого вещества составила 5,6 массы Земли (около трети массы Нептуна). С учетом всех погрешностей допустимый диапазон массы упавшей суперземли составляет от 4,6 до 5,8 массы Земли.

Значение для звездной астрофизики

Выводы исследования меняют подход к моделированию звезд главной последовательности.

В звездной астрономии долгое время действовал принцип, согласно которому химический состав внешних слоев звезды считался прямым отражением химического состава исходного облака межзвездного газа. Работа Мутлу Йылдыза доказывает обратное: рождение планетной системы меняет химическое распределение элементов внутри родительской звезды.

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

Это объясняет парадокс, обнаруженный при спектральном анализе звезд в двойных системах. Звезды-близнецы в широких парах формируются из одного газопылевого облака и должны обладать идентичным элементным составом. Тем не менее телескопы регулярно выявляют пары (например, солнцеподобную звезду TOI-1173 A и ее компаньона), где один из компонентов содержит заметно больше тяжелых элементов, чем второй. Такая разница прямо указывает на поглощение каменистых тел массой до десятков масс Земли.

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

Источник: Monthly Notices of the Royal Astronomical Society