Физики провели 259 трехмерных симуляций столкновения пары твердых планетоподобных тел околоземной массы и получили функцию, которая приближенно связывает параметры таких соударений с долей теряемой атмосферы. По словам авторов, этот результат можно будет встраивать в модели образования планет и использовать для анализа сценариев ударного формирования Луны на основе данных об истории атмосферы Земли. Работа опубликована в The Astrophysical Journal Letters.
Считается, что землеподобные планеты образуются из десятков планетезималей, которые рождаются из протопланетного диска и сталкиваются друг с другом. Одно из таких столкновений, в частности, могло привести к образованию Луны из молодой Земли и гипотетической планеты Тейя.
Среди последствий соударений можно выделить потерю газовой оболочки (полную или частичную) небесного тела — это событие может повлиять на итоговый состав атмосферы планеты и, предположительно, частично объяснить наблюдаемые различия в составах атмосфер похожих планет в одной и той же системе. Тем не менее характер такого влияния окончательно не изучен — до недавнего времени исследования в основном останавливались на одномерных моделях, которые, кроме того, ограничивались лобовым столкновениям или рассматривали слишком малое число сценариев для построения прогнозов в широком диапазоне параметров.
Ученые из Великобритании и США под руководством Джейкоба Кегеррейса (Jacob Kegerreis) из английского Института вычислительной космологии смоделировали столкновения, подобные тому, в котором могла образоваться Луна, и проследили за влиянием таких ударов на потерю атмосферы. В качестве цели и налетающего тела выступали шаровидные объекты, содержащие железное ядро (на которое приходилось 30 процентов массы) и мантию из горных пород (несущую остальные 70 процентов массы). В некоторых симуляциях также участвовали однородные тела (полностью железные или каменистые).
Объект-цель (а в ряде случаев и налетающее тело) дополнительно обладал атмосферой, масса которой составляла сотую долю от общей массы ядра и мантии. Для анализа поведения атмосферы в высоком разрешении ученые использовали трехмерные симуляции в высоком разрешении на основе уравнений гравитации и гидродинамики. В них участвовало по 10 миллионов атмосферных частиц в расчете на одну массу Земли в случаях с крупными массами сталкивающихся тел и на 10–0,5 (около 0,316) масс Земли — в случаях с менее массивными участниками соударения.
Всего исследователи провели 259 симуляций, которые разделили на основные три группы. В первом наборе сценариев масса цели и налетающего объекта менялись в диапазоне 10–1,25–100,25 (то есть 0,056–1,778) масс Земли. По скорости случаи делились на медленные — в них относительная скорость налетающего тела была равна взаимной второй космической скорости системы (то есть той, которая необходима каждому из тел, чтобы преодолеть тяготение другого), и быстрые — в них налетающее тело двигалось втрое быстрее. Кроме того, ученые моделировали ситуации лобового столкновения — когда налетающее тело двигалось по направлению на центр цели, и задевающего, в котором прицельный параметр составлял долю в 0,7 от суммы радиусов тел.
Во второй группе симуляций для подмножества пар цель–налетающее тело физики подробнее исследовали зависимость результата от относительной скорости налетающего тела (были добавлены сценарии с промежуточным значением, в которых она вдвое превышала взаимную вторую космическую скорость) и от прицельного параметра (дополнительно к 0 и 70 процентам исследовались значения в 30, 50 и 90 процентов от суммы радиусов).
В третьем наборе сценариев при фиксированных массе цели (одна масса Земли), прицельном параметре (70 процентов суммы радиусов) и скорости (утроенная взаимная вторая космическая) авторы заменяли один или оба объекта на однородные (полностью каменные или полностью железные) тела такой же массы.
По итогам каждой из симуляций авторы вычисляли долю атмосферы, которую небесное тело теряет в результате столкновения, а затем исследовали зависимость этого числа от менявшихся параметров модели — прицельного параметра, скорости налетающего тела, масс и плотностей цели и налетающего тела.
В результате ученые смогли приближенно представить потерянную долю атмосферы в виде степенной функции, которая зависит от отношения скорости налетающего тела ко второй космической скорости, отношения массы налетающего тела к полной массе системы, отношения плотностей тел и доли взаимодействующей (непосредственно участвующей в соударении) массы, зависящей от прицельного параметра. По словам авторов, в будущих исследованиях обнаруженные закономерности можно будет встраивать в уже существующие модели формирования планет.
Кроме того, физики отметили, что в наиболее близких к различным версиям образования Луны сценариях протоземля могла потерять порядка 10–60 процентов атмосферы. В дальнейшем с помощью этих прогнозов, вероятно, можно будет уточнить модели формирования спутника на основе данных об истории земной атмосферы.
Ранее мы рассказывали о том, как компьютерное моделирование «взрывного» рождения Луны объяснило схожесть ее изотопного состава с земным и как ученые смогли воспроизвести наблюдаемое наклонение орбиты спутника к экватору Земли в таком сценарии.
Николай Мартыненко