15 175 107 звёзд. Каждая получила свой возраст. Но ни одна цифра не лишена налёта неопределённости: возраст звезды — это не гвоздь, забитый на оси времени, а облачко вероятности, туман, скрывающий истинное значение. Для галактической эволюции такая временна́я метка бесценна: без неё карта Млечного Пути остаётся плоской. Ещё в начале XX века Сесилия Пейн-Гапошкина показала, что звёзды состоят в основном из водорода, а Ханс Бете объяснил, как ядерные реакции определяют их жизненный цикл. Однако, в отличие от блеска или положения на небе, возраст не измерить напрямую. Его выводят из моделей звёздной эволюции или — для избранных — из сейсмических колебаний (астросейсмологии, то есть измерения возраста по колебаниям звезды), как это делали телескоп «Кеплер» под руководством Уильяма Боруки и TESS. Но таких точных возрастов — капля в море по сравнению с миллионами звёзд, для которых есть лишь данные фотометрии (измерения яркости) и спектроскопии (разложения света по длинам волн). Миссия Gaia Европейского космического агентства предоставила низкоразрешающие спектры для 220 миллионов источников. Из них извлекают температуру, силу тяжести на поверхности, металличность (долю элементов тяжелее гелия) — параметры, несущие отпечаток возраста. Вопрос в том, как по этим косвенным признакам восстановить хронологию.
Исследователи обучили ансамбль из 135 569 звёзд из проекта APOGEE — их возрасты были известны достаточно точно. На этих эталонах натренировали NGBoost — алгоритм машинного обучения, который выдаёт не одно число, а целое распределение вероятностей. Однако стандартный NGBoost считает метки абсолютно точными, а в действительности возраста эталонов сами содержат погрешности. Поэтому функцию потерь модифицировали: полная модельная дисперсия стала учитывать квадрат ошибки обучающего возраста: \(\sigma_{\text{total}}^2 = \sigma^2 + \sigma_{A,\text{obs}}^2\). Иными словами, неуверенность модели складывается с известной неточностью меток. Но и входные данные — температура, химический состав — известны с ошибками. Чтобы оценить, как эта неопределённость просачивается в итог, применили метод Монте-Карло (моделирование с многократной случайной выборкой): для каждой звезды генерировали 200 вариантов входных параметров, слегка возмущённых в пределах погрешностей, и усредняли предсказания. Так родилась модель, которая не только называет возраст, но и честно сообщает разброс.
На тестовой выборке предсказанные возраста легли вдоль линии идеального совпадения, показав систематический сдвиг всего –0,119 миллиарда лет и разброс 1,200 миллиарда лет. Звёзды с большими неопределённостями получали и более широкие предсказанные интервалы — модель адекватно оценивала свою уверенность. Затем алгоритм применили к 15 175 107 звёздам каталога Gaia с параметрами в определённых диапазонах (температура 3700–6032 К, логарифм поверхностной гравитации 0,66–4,83, металличность от –0,76 до 0,53). Для каждой теперь известен не только возраст, но и три компонента ошибки: собственная модельная, перенесённая из входных данных и полная. Типичная полная погрешность составляет около 1,4 миллиарда лет — примерно четверть возраста Земли.
Сравнение с независимыми каталогами, полученными методами астросейсмологии и спектроскопии, подтвердило согласие: расхождения укладываются в 1,5–1,8 миллиарда лет, причём основные систематические сдвиги унаследованы от шкалы возрастов обучающей выборки. Главный результат — не просто цифры, а честная карта ошибок. Теперь астрономы могут отбирать звёзды с приемлемой точностью возраста для изучения эволюции Галактики. Вероятностный подход позволяет не маскировать неопределённость, а работать с ней: строить распределения возрастов в разных областях Млечного Пути, помня, что каждый возраст — это не гвоздь, а облачко вероятности. Это особенно ценно для истории звездообразования и химической эволюции диска. В перспективе модель сможет учитывать дополнительные индикаторы — содержание лития, вращение звёзд — и использовать данные других масштабных обзоров, таких как LAMOST или будущий спектрограф 4MOST. Уточнение шкалы возрастов с помощью новых астросейсмологических данных сузит вероятностный туман. Так, глядя сквозь него, исследователи начинают различать контуры истории Млечного Пути. Неточность возрастов остаётся главным препятствием, но честный учёт неопределённостей — это способ не заблудиться в тумане, а проложить по нему надёжный маршрут.
Возраст самой старой известной звезды в Млечном Пути (HD 140283, «Мафусаил») поначалу оценили в 16 миллиардов лет — больше возраста Вселенной! Позже ошибку исправили, но случай показывает, насколько трудно измерить возраст. Новый метод как раз и призван снизить такие казусы, честно указывая на размытость оценки.
Идея считывать возраст звезды по её спектру напоминает трикодер из «Звёздного пути», способный определить состав и возраст любого объекта за секунду. Реальность пока далека от такой точности, но каталог на 15 млн звёзд — шаг к «галактическому трикодеру».