Из-за ошибки в калибровке «Джеймс Уэбб» случайно открыл быстрый способ искать экзопланеты

источник
3 фото

В наблюдательной астрономии точность фундаментальных измерений зависит от стабильности опорных источников. Чтобы переводить электрические сигналы детекторов в физические величины — потоки излучения, светимости и температуры далеких объектов, космическим обсерваториям необходимы звезды с неизменным и досконально изученным блеском. Четыре десятилетия звезда HD 101452 спектрального класса A2V в созвездии Центавра использовалась исследователями в качестве эталона для калибровки инфракрасных телескопов. В 2024 году ее включили в специализированную калибровочную программу CAL 4496 для настройки приборов космической обсерватории «Джеймс Уэбб» (JWST).

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

Проблема калибровочных стандартов и статус звезды HD 101452


Для калибровки абсолютного потока излучения исследователи чаще всего выбирают горячие белые звезды спектрального класса A. Их фотосферы имеют высокие температуры, а спектры содержат небольшое количество линий поглощения, что упрощает теоретическое моделирование их поведения. Для работы приборов телескопа «Уэбб» требовалась точность калибровки с относительной погрешностью строго менее 5%. Чтобы гарантировать такой уровень достоверности, звезда-ориентир должна быть одиночной и изолированной: рядом с ней не должно быть ни вторичных звездных компонентов, ни газопылевых дисков, вносящих дополнительное излучение.

Звездная астрономия показывает, что соблюсти эти критерии сложно. Приблизительно 25% звезд спектрального класса А входят в состав тесных кратных систем, а около 33% обладают остаточными дисками из пыли и обломков, которые генерируют избыточный поток в инфракрасном диапазоне на длинах волн около 24 микрометров.

Звезда HD 101452, расположенная на расстоянии около 165 парсек от Земли, считалась стабильным калибратором с начала 1980-х годов. Первые признаки несоответствия проявились при автоматической обработке снимков в рамках программы изучения экзопланет Rocky Worlds DDT. Аналитический программный комплекс SPARTA зафиксировал аномальный избыток инфракрасного излучения величиной 6% на длине волны 15 микрометров.

В июне 2024 года в рамках калибровочной программы CAL 4496 обсерватория провела съемку HD 101452 с помощью камеры MIRI в четырех спектральных полосах: 15, 18, 21 и 25,5 микрометра. Общее время накопления сигнала на фильтр составило от 17,3 до 41,5 секунды. На всех полученных кадрах в структуре распределения света центральной звезды отчетливо проявилось асимметричное смещение яркости.

Астрометрия и физическая природа обнаруженного спутника


Для детального анализа структуры изображения группа исследователей под руководством Рэйчел Боуэнс-Рубин применила алгоритмы разделения сигналов. Присутствие второго объекта подтвердилось во всех четырех фильтрах на одинаковом угловом удалении.

Основные измеренные параметры системы:

Угловое расстояние между компонентами составило 1,276 угловой секунды с погрешностью 0,004 угловой секунды.
С учетом тригонометрического параллакса системы это значение соответствует физическому расстоянию в проекции на небесную сферу, равному 210,9 астрономической единицы.
Позиционный угол спутника относительно основного компонента составил 84,3 градуса.
Разница в яркости между объектами в среднем инфракрасном диапазоне составила приблизительно 27 раз в пользу главной звезды.
Чтобы доказать гравитационную связь обнаруженного источника с главной звездой, астрономы сопоставили координаты MIRI с данными европейского астрометрического каталога Gaia DR3. В каталоге Gaia данный спутник был зафиксирован в оптическом диапазоне еще в 2016 году, но не был классифицирован как физический компаньон.

Сопоставление координат за восьмилетний интервал подтвердило, что оба объекта перемещаются в пространстве с одинаковой скоростью и направлением. Статистический тест подтвердил модель общего собственного движения: величина изменения функции хи-квадрат составила Δχ² = 5953 в пользу общего перемещения по сравнению с гипотезой неподвижного фонового источника. Вероятность случайного визуального наложения далекой фоновой звезды составила 0,005%.

Построение и подгонка спектрального распределения энергии по четырем точкам MIRI и фотометрии Gaia позволили рассчитать физические характеристики спутника. Его эффективная температура равна 3724 Кельвина, масса находится в пределах нижней части главной последовательности, а спектральный класс определен на границе между поздними оранжевыми карликами и ранними красными карликами (K7V-M0.5V).

Вклад этого спутника в совокупный поток системы превышает 3,5% в зависимости от фильтра, что делает звезду непригодной для роли абсолютного фотометрического калибратора. Дополнительно астрометрические данные Gaia выявили высокий показатель приведенной погрешности единичного веса (RUWE = 2,309 при норме для одиночных звезд до 1,25). При совмещении данных миссий Hipparcos и Gaia было обнаружено вековое ускорение собственного движения со значимостью 51 сигма, соответствующее изменению тангенциальной скорости на 3300 метров в секунду. Поскольку спутник на расстоянии 211 астрономических единиц не способен вызвать столь мощное динамическое возмущение, в системе с высокой вероятностью присутствует еще один, третий компонент на более близкой орбите.

Ограничения инфракрасных детекторов и отказ от коронографа


Обнаружение спутника HD 101452 B стало важным техническим прецедентом из-за условий съемки. Традиционно регистрация тусклых объектов рядом с яркими звездами требует применения коронографа — оптического узла, блокирующего прямой свет звезды физической непрозрачной маской или фазовой пластиной.

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

Они требуют сложного центрирования звезды на микроскопической маске с точностью до долей пикселя.
Коронограф создает мертвую зону (внутренний рабочий угол), внутри которой регистрация объектов физически невозможна из-за дифракции на краях маски.
Для калибровки светового ореола требуется многократный поворот телескопа вокруг оптической оси либо наведение на отдельную опорную звезду, что увеличивает затраты времени обсерватории.
Съемка HD 101452 проводилась в базовом режиме прямого формирования изображения на субматрице детектора без коронографа. Главная трудность прямого наблюдения на малых угловых расстояниях заключается в поведении детектора инфракрасного излучения MIRI (матрицы на основе кремния, легированного мышьяком, с примесной проводимостью). При попадании высокой плотности фотонов в пиксель накапливающийся в его объеме электрический заряд начинает менять локальные электростатические границы. Вновь возникающие электроны испытывают кулоновское отталкивание и смещаются в соседние, менее насыщенные пиксели.

Этот физический процесс в полупроводниках называют эффектом расширения зарядового пятна (brighter-fatter effect). На итоговом снимке яркая звезда искусственно увеличивается в диаметре, а ее дифракционные кольца деформируются, перекрывая слабые сигналы в пространстве ближе трех дифракционных пределов оптической системы (3λ/D, где λ — длина волны излучения, а D — диаметр главного зеркала телескопа, равный 6,5 метра).

Кадры камеры MIRI до и после обработки. В верхнем ряду — исходные изображения калибровочной звезды HD 101452 на разных длинах волн инфракрасного диапазона (от 15 до 25,5 мкм). Скрытый спутник едва заметен как легкая деформация правого края звездного пятна. В нижнем ряду — результат работы алгоритма: свет главной звезды полностью удален, а ее спутник отчетливо виден в виде яркой точки на каждом из фильтров.

Авторы исследования доказали, что этот эффект можно нейтрализовать математически без применения аппаратных масок. Для этого использовали метод референтной дифференциальной визуализации (RDI). В качестве опорного источника была выбрана другая калибровочная звезда аналогичного спектрального класса — HD 2811, обладающая практически той же звездной величиной в среднем инфракрасном диапазоне (W4 = 7,0 против W4 = 6,8 у HD 101452).

Поскольку плотность фотонного потока от обеих звезд была близкой, полупроводниковая матрица MIRI сформировала одинаковые структурные искажения в обоих случаях. С помощью алгоритма анализа главных компонент (PCA), реализованного в программном пакете VIP, профиль опорной звезды был вычтен из снимка целевой звезды. Это позволило подавить остаточный звездный ореол вплоть до расстояний менее 2λ/D на длине волны 25,5 микрометра. На фильтре 15 микрометров за 17,3 секунды экспозиции был достигнут предел контраста 5,9·10⁻⁴ (разность яркостей Δmag = 8,1) на угловом расстоянии менее 2,5 угловой секунды. Спутник был выделен с отношением сигнал/шум, равным 40.

Физика атмосфер планет-гигантов в среднем инфракрасном диапазоне


Возможность получать глубокие контрасты за десятки секунд без коронографа имеет прямое отношение к поиску экзопланет. Большинство инструментов прямого наблюдения (включая наземные адаптивные системы SPHERE и GPI, а также камеру NIRCam на борту самого JWST) оптимизированы под ближний инфракрасный диапазон (длины волн от 1 до 5 микрометров). В этом диапазоне зеркало диаметром 6,5 метра обеспечивает высокое угловое разрешение, однако астрономы сталкиваются с фундаментальными ограничениями, накладываемыми физикой планетных атмосфер.

Когда температура газового гиганта опускается ниже 200-150 Кельвинов, термодинамическое равновесие газов смещается:

Водяной пар полностью конденсируется, образуя оптически плотный ярус облаков из водяного льда.
Газообразный метан накапливается в верхних слоях атмосферы, создавая мощные полосы поглощения на длинах волн от 1,6 до 3,3 микрометра.
Повышенная металличность атмосферы усиливает непрозрачность газовой оболочки.
В результате тепловой поток из внутренних слоев планеты блокируется. В диапазоне от 3,5 до 5 микрометров планеты с массами порядка массы Юпитера и температурами на уровне температур планет Солнечной системы становятся тусклыми, их светимость падает на порядки.

В среднем инфракрасном диапазоне (длины волн более 15 микрометров), где работает MIRI, ситуация меняется. На этих длинах волн планета излучает собственное тепло в виде непрерывного спектра, близкого к спектру абсолютно черного тела. Влияние облачности и молекулярного поглощения здесь минимально по сравнению с ближним диапазоном. Фотометрия на 15-25 микрометрах отражает общую светимость и температуру объекта, снижая неопределенность теоретических моделей.

Расчет чувствительности для соседних звездных систем


Опираясь на кривые контраста, полученные при анализе HD 101452, авторы работы оценили способность камеры MIRI обнаруживать реальные экзопланеты методом коротких некоронографических экспозиций. Для расчетов использовалась сетка эволюционных моделей BEX-ATMO, описывающая остывание планет с массами от долей массы Юпитера.

Пределы чувствительности камеры MIRI при экспозициях короче одной минуты. Графики показывают, планеты какой минимальной массы телескоп «Уэбб» способен различить на разных расстояниях от звезды. В фильтрах 18 и 21 мкм (F1800W и F2100W) кривые чувствительности для близких систем (например, звезды Тигардена — сплошная линия) опускаются ниже горизонтали 1 массы Юпитера на масштабах от 15 астрономических единиц.

Было рассмотрено несколько типов близких к Солнцу звезд, обладающих схожей с HD 101452 видимой яркостью в среднем инфракрасном диапазоне (величина W4 в каталоге WISE около 6,8):

Звезда Тигардена (расстояние 3,8 парсека). Старый красный карлик спектрального класса M7.5 возрастом около 8 миллиардов лет. Из-за близости к Солнечной системе угловое расстояние пересчитывается в компактные физические масштабы. За 34,6 секунды экспозиции в фильтре F2100W прибор способен зафиксировать планету с массой менее 1 массы Юпитера на удалении от 15 астрономических единиц, а предел обнаружения на больших удалениях достигает 0,75 массы Юпитера.
Молодые маломассивные звезды (расстояние 6,5-8 парсек). Группа красных карликов возрастом около 700 миллионов лет (включая GL Vir, G 13-22 и SCR J0740-4257). За счет меньшего возраста планеты в таких системах сохранили больше внутреннего тепла от гравитационного сжатия. Предел чувствительности при аналогичной экспозиции составил 0,7 массы Юпитера на орбитальных расстояниях от 8 до 10 астрономических единиц.
Молодые движущиеся группы (расстояние 15-20 парсек). Объекты ассоциации AB Дорада возрастом около 100 миллионов лет (например, звезда HIP 17695). В силу молодости даже суб-юпитерианские планеты обладают высокими температурами. В фильтрах F1500W и F1800W короткие снимки позволяют преодолеть порог массы Юпитера на удалениях от 20 до 50 астрономических единиц.
В то же время фильтр F2550W (25,5 микрометра) показал низкую пригодность для коротких экспозиций. На этой длине волны доминирует инструментальный тепловой фон оптической системы телескопа, что снижает отношение сигнал/шум при малых временах интеграции. Оптимальными для поиска холодных суб-юпитеров оказались фильтры F1500W, F1800W и F2100W.

Оптимизация рабочего графика космического телескопа


Главным фактором, ограничивающим число наблюдаемых объектов на телескопе «Уэбб», являются не технические возможности сенсоров, а непроизводительные временные затраты на переориентацию обсерватории.

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

Результаты анализа HD 101452 предлагают другое решение: попутные наблюдения. На телескопе регулярно выполняются крупные научные программы по исследованию геометрии протопланетных дисков, спектроскопии звездных атмосфер или картографированию остаточных структур с помощью ближнеинфракрасных приборов NIRCam и NIRSpec.

Если телескоп уже сориентирован на выбранную близкую звезду, смена инструмента с NIRCam на MIRI не требует механического перенаведения всей обсерватории. Переключение оптического тракта занимает около 30 минут, а смена спектральных фильтров внутри вращающегося колеса MIRI — менее 10 минут.

Добавив к уже утвержденной научной программе серию коротких снимков MIRI общей продолжительностью менее двух минут чистого времени, астрономы получают данные о наличии или отсутствии планет-гигантов в диапазоне масштабов от 5 до 50 астрономических единиц — в областях, соответствующих расположению внешних планет в Солнечной системе.

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