Субнептуны — самый распространённый класс экзопланет в Галактике, но их внутреннее устройство остаётся предметом споров. Традиционно считалось, что под водородной атмосферой лежит чёткая граница с силикатной магмой. Однако недавние расчёты показали, что при давлениях и температурах, характерных для недр субнептунов, силикат и водород могут полностью смешиваться. Это ставит задачу описания термодинамических свойств такой сверхкритической жидкости, без чего невозможно строить корректные модели внутреннего строения планет. Работа находится в русле исследований фазовых равновесий и уравнений состояния для силикатно-водородных систем, активно развиваемых в последние годы.
Авторы провели первопринципные молекулярно-динамические расчёты в рамках теории функционала плотности с использованием программного пакета VASP и обменно-корреляционного функционала PBEsol. Моделировалась ячейка из 243 атомов (27 Mg, 27 Si, 81 O, 108 H), что соответствует составу MgSiO3H4. Сначала расплавляли кристаллический MgSiO3, затем добавляли водород. Расчёты вели при температурах 4000, 6000 и 8000 К и давлениях до ~600 ГПа. Полученные данные по давлению и внутренней энергии аппроксимировали свободной энергией Гельмгольца, что позволило вывести плотность, энтропию, теплоёмкость, тепловое расширение, модуль объёмного сжатия и параметр Грюнайзена из единой термодинамической поверхности.
Сверхкритическая жидкость MgSiO3H4 оказалась менее плотной, чем сухой силикатный расплав, при всех исследованных условиях. Отклонение от идеального смешения объёмов MgSiO3 и H2 достигает нескольких процентов и растёт с давлением. Анализ структуры показал, что с ростом давления координационное число Si–H увеличивается, а количество устойчивых связей H–H резко падает: при 4000 К доля атомов водорода, участвующих в долгоживущих контактах H–H, снижается ниже 0,1 при давлениях около 340 ГПа. Это указывает на переход от молекулярного водорода к более атомарному, встроенному в силикатный каркас. Применение уравнения состояния к моделям субнептунов массой 6 масс Земли с одинаковым общим запасом водорода 4% показало, что в зависимости от давления на границе атмосфера–недра (1 ГПа против 5 ГПа) радиус планеты меняется на 5,7%.
Работа впервые даёт самосогласованное уравнение состояния для сверхкритической силикатно-водородной жидкости, что закрывает важный пробел в моделировании субнептунов. В отличие от соседних работ, которые в основном оценивали фазовые равновесия или смешение на уровне термодинамических моделей, здесь получены прямые первопринципные данные о плотности и термодинамических производных. Это позволяет количественно учитывать растворённый водород в недрах и его влияние на радиус и тепловую эволюцию планет.
Естественное продолжение — расширение состава за счёт железа, которое может сильно влиять на плотность, окислительно-восстановительные условия и фазовые переходы. Также необходимы экспериментальные проверки при высоких давлениях и температурах, которые пока крайне редки. Совместные усилия теоретиков и экспериментаторов позволят построить уравнения состояния для полного диапазона составов, давлений и температур, релевантных субнептунам.
Результаты напрямую повлияют на модели внутреннего строения и эволюции субнептунов, а также на интерпретацию данных о массах и радиусах экзопланет.
Включить железо в расчёты и провести экспериментальные измерения плотности и сжимаемости сверхкритических силикатно-водородных жидкостей.
Работа затрагивает фундаментальную проблему планетологии: как определить состав и внутреннюю структуру экзопланет по ограниченным наблюдательным данным. Она показывает, что игнорирование неидеальности смешения силиката и водорода может приводить к систематическим ошибкам в оценках содержания водорода и теплового состояния недр.
Этот раздел адресован автору работы и его коллегам. Соседние работы нашёл смысловой поиск по нашему архиву — не по совпадению слов, а по смыслу. Мы не рецензируем: мы показываем поле вокруг.
Что сделано в работе
Авторы провели первопринципные молекулярно-динамические расчёты сверхкритической смеси MgSiO3H4 при давлениях до 600 ГПа и температурах 4000–8000 К, получили уравнение состояния и показали, что плотность смеси ниже, чем у сухого силиката, а поведение отклоняется от идеального смешения. Применив уравнение состояния к моделям субнептунов массой 6 масс Земли, они обнаружили, что распределение водорода между атмосферой и недрами влияет на радиус планеты (разница 5,7% при граничных давлениях 1 и 5 ГПа).
В чём сила
Сила работы в получении самосогласованного уравнения состояния сверхкритической силикатно-водородной жидкости из первых принципов, что закрывает пробел в моделировании субнептунов и даёт возможность учитывать реальное распределение водорода в недрах.
Куда можно двигаться
Что можно попробовать
Почему это может быть важно
Учёт сверхкритического смешения силиката и водорода может изменить оценки радиуса субнептунов на несколько процентов, что важно для интерпретации наблюдаемой зависимости радиус–масса и выявления внутренней структуры экзопланет.
Работы рядом в нашем архиве
Ближайшая работа архива — 0.65 по смысловой близости. Обычный для нашего корпуса уровень.
При давлении около 340 ГПа и температуре 4000 К лишь один из десяти атомов водорода сохраняет связь с другим водородом дольше 100 фемтосекунд — остальные уже встроены в силикатную сеть.