КОСМИЧЕСКИЕ ИССЛЕДОВАНИЯ, 2012, том 50, № 6, с. 452-461
УДК 523.26
О НАДЕЖНОСТИ ОПРЕДЕЛЕНИЯ ОРБИТАЛЬНЫХ ПАРАМЕТРОВ ЭКЗОПЛАНЕТ ДОПЛЕРОВСКИМ МЕТОДОМ © 2012 г. Б. Р. Мушаилов, В.С. Теплицкая
Государственный астрономический институт им. П.К. Штернберга МГУ им. М.В. Ломоносова Поступила в редакцию 23.03.2011 г.
Рассмотрено гипотетическое влияние непостоянства скорости света, обусловленного параметрами движения источника излучения, на результаты спектрометрических измерений звезд в рамках программы поиска экзопланет. Проведен учет ускорений звезд относительно барицентра системы звезда—планета (планеты).
Показано, что корректный учет компоненты барицентрического лучевого ускорения звезд приводит к существенному увеличению (вплоть до порядка величины) больших полуосей орбит обнаруживаемых кандидатов в экзопланеты, что согласуется с современной теорией образования и динамической эволюции планетных систем.
ВВЕДЕНИЕ
В настоящее время (2011 г.) около 70% кандидатов в экзопланеты определены на основании "метода лучевых скоростей" — спектроскопического измерения радиальных компонент скоростей звезд [20].
Планета (или планеты) обнаруживается по наблюдаемой с Земли вариации лучевой скорости V, своей звезды, которая в нерелятивистском случае, когда V, существенно меньше скорости с регистрируемого излучения звезды, определяется доплеровской формулой
V, (!) = -ет)//о. (1)
Здесь ДД0 = /(?) — /0, ДО — регистрируемая частота излучения от звезды избранной спектральной линии, соответствующей "лабораторной" (при V, = 0) частоте/0. Значение лучевой скорости неотрицательно в случае, когда звезда не приближается к наблюдателю, а при сближении звезды с наблюдателем V, < 0.
Регистрируемая со временем вариация лучевой скорости звезды, в случае пренебрежения взаимного гравитационного возмущения экзо-планет, интерполируется суперпозицией кепле-ровских орбитальных компонент, соответствующих N гипотетическим планетам, обращающимся относительно звезды (строго говоря, центра масс экзопланетной системы) [11]:
N
(t) = Y/i (cos[®i + е» P)] +
i=i
(2)
+ e i cos щ) + Y ck
k=1
Искомыми величинами, определяемыми по критерию минимума среднеквадратичной разности модельной лучевой скорости (2) и полученной из наблюдений, являются эксцентриситеты е, средние долготы X,, орбитальные периоды Р, экзопланет, через которые могут быть выражены истинные аномалии 9,- экзопланет, а также величины аргументов перицентров орбит планет ю,-, коэффициенты Ск полиномиального тренда лучевой скорости звезды, обусловленного ограниченностью продолжительности наблюдений, амплитудные коэффициенты V, связанные с минимальными массами экзопланет.
Несмотря на достигнутый за прошедшие 20 лет с момента результативного применения допле-ровского метода к обнаружению экзопланет определенный прогресс в точности измерений лучевой скорости звезд, алгоритмов статистической обработки и планирования наблюдений, большинство экзопланет официально именуются всего лишь кандидатами в экзопланеты. Данный факт обусловлен не только косвенной природой метода лучевых скоростей. Наблюдаемая вариация лучевой скорости звезды может быть вызвана не гравитационным влиянием планет, а, например различными эффектами активности в звездной атмосфере [6]. Кроме того, применяемые с целью поиска периодических колебаний, потенциально свидетельствующих о наличии у звезды невидимых планет, периодограммы типа Ломба— Скаргла или различные ее модификации [21] непосредственно не позволяют оценивать статистическую значимость выявленных гармоник, т.е. отличить полезные сигналы (искомые гармоники) от шумовых эффектов. Полная дисперсия ошибок временных рядов лучевых скоростей звезд, наряду с соответствующей компонентой непо-
средственно инструментальной погрешности измерений, содержит и значительный добавочный компонент, вызванный нерадиальными звездными осцилляциями [17], а также условиями и методикой проведения и обработки результатов наблюдений.
Наличие в измерительных ошибках не исключенных систематических погрешностей увеличивает эффективную величину осцилляций лучевой скорости, что, в свою очередь, изменяет соотношение статистических весов наблюдений и, в конечном итоге, влияет на получаемые значения оценок орбитальных параметров и масс экзопла-нет. Если "астрофизическая компонента" осцил-ляций лучевой скорости априорно оценивается на основе эмпирических корреляций со скоростью вращения и уровнем активности звезды, измеряемым по К и Н линиям кальция [23], то влияние характеристик спектрографа, условий и методики проведения наблюдений для различных звезд и обсерваторий на высокоточные измерения лучевых скоростей в рамках программы поиска экзопланет не исследовано.
Наличие систематических ошибок при измерении лучевых скоростей приводит к тому, что получаемые оценки орбитальных параметров и масс экзопланет оказываются существенно искаженными и ненадежными.
Получаемые в настоящее время методом лучевых скоростей, основанном на эффекте Доплера, орбитальные параметры и массы экзопланет приводят к парадоксальным орбитальным конфигурациям предвычисляемых экзопланет (рис. 1). Наличие "горячих Юпитеров" — малых по величине орбитальных периодов и больших полуосей предвычисляемых экзопланет не согласуется с теоретическими воззрениями на происхождение и эволюцию орбит планет, в том числе, и в Солнечной системе [13, 35, 42].
Более того, столь тесное орбитальное сосуществование звезд и горячих Юпитеров не реализуемо на космогонических интервалах времени. После значительного прогресса в точности регистрации спектров излучений звезд стремление в приоритетном порядке обнаружить как можно больше кандидатов в экзопланетные системы не способствовало корректности интерпретации получаемых результатов [33].
УЧЕТ НЕИНЕРЦИАЛЬНОСТИ ДВИЖЕНИЙ
Как следует из фундаментального принципа наименьшего действия состояние механической системы полностью определяется заданием координат и скоростей компонент системы. Однако силовое взаимодействие (уравнения движения) в механической системе определяются производными от скоростей, то есть ускорениями исследу-
140
100
60
20 0
i
i — i_i_i_i_i_i
0123456789 10
Рис. 1. Гистограмма распределения экзопланет по большим полуосям орбит в абсолютных величинах (светлые столбцы) и в самосогласованной системе единиц (темные) [33].
емых компонент (уравнения Лагранжа или Ньютона не содержат производные от ускорений) [30].
Эффект Доплера, выраженный в форме (1), не учитывает силовое (гравитационное) воздействие, приводимое к центру масс динамической системы "звезда—планеты" и характеризуемое ускорением звезды, спектр которой измеряется.
Не умаляя общности, будем считать, что вариации лучевой скорости звезды Б обусловлены ее обращением вокруг силового центра (центра масс О гипотетической экзопланетной системы) с постоянной угловой скоростью ^ = v/r0 по круговой орбите с радиусом г0 (рис. 2). Будем считать, что наблюдатель находится в плоскости орбиты звезды Б на расстоянии г > г0 от центра О, который относительно наблюдателя имеет компоненту лучевой скорости Предполагая, что и V ^ с, показания часов в системе координат, связанной со звездой Б, будем обозначать через а в случае наблюдателя — ? (при V ^ с часы считаются синхронизированными). Ввиду значительного расстояния г от звезды до наблюдателя звезда Б считается точечным источником излучения.
Излучение от звезды Б, сгенерированное в момент времени 1Б0, с учетом принципа относительности Галилея (нерелятивистский случай; обоснование приводится ниже), достигнет наблюдателя в момент времени
to =■
r(tso) - rocos^t,
so
С - v(tso)sin Qtso - wr(tso)
(3)
Излучение от звезды, сгенерированное в 1Б1 через период излучения ТБ, достигнет земного наблюдателя в
tn = Ts +-
r(tsi) - ro cos QlS]
(4)
с - v(tsi)sin Qtsi - wr(tsi) Так как tS1 = tS0 + TS, aS = v2/r0, то, вводя обозначения ar(tS0) = aScos(fitS0), vr(t) = v(t)sinfít и
вой скоростью wr(t).
учитывая, что Q,TS < 1, wr = const за период Ts, в рассматриваемом случае jv | = const, будем иметь
(5)
V, (tsl) = V. ^8о)со8(О.Т) + V 8Ш(ОТ5 )С08(П^о)
« V.о) + Т8йг($5о), = о) + .
Следовательно, период Т = ?П — ^ регистрируемого наблюдателем излучения, периодичностью Т, звезды Б согласно (3)—(5), с учетом того, что г^^) = г > г0, с рассматриваемой точностью будет равен:
T = T,
+
w ,
1 + —
с - w, - Tsar (tsо)
as cos(Qtso)r(tso)
(6)
(с - W.)(с - . - Твйго))у
где Wг = ^ + V. (Б)).
И, наконец, учитывая, что в рассматриваемом случае |wг| <§ с, а8Т8 = &?г0Т8 < с, получим:
T = 1+w,+ra, тs
(7)
S
с с
так что в случае соответствующих частот будем иметь
bflfs = -w rfc - raj с2, где hf=f(t) -f.
(8)
При аг = 0 (в случае равномерного движения источника излучения) выражение (8) будет эквивалентно общеизвестной нерелятивистской формуле, полученной в 1842 г. К. Доплером применительно к звуковым явлениям и распространенной на электромагнитные волны А. Физо в 1848 г.
Однако в общем случае, при неравномерном движении источника излучения изменение регистрируемого частотного спектра обусловлено также наличием радиальной компоненты ускорения излучающего объекта (звезды), так что изменение частоты регистрируемого от звезды излучения происходит и при нулевой величине лучевой скорости звезды. При этом значимость второго слагаемого в (8) как раз и существенна для удаленных источников, когда г/с > 1, что и реализуется, в отличие от Солнечной системы, для гипотетических экзопланетных систем.
Формально появление второго слагаемого в (8) может быть интерпретировано в рамках Эйнштейновской общей теории относительности как результат, обусловленный гравитационным смещением в эквивалентной гравитационному полю ускоренно движущейся системе. В поле тяготения с разностью потенциалов ДФ относительная девиация частот имеет вид [41]:
А/// = ДФ/ с2,
(9)
где ДФ = Ф2 — Ф1, гравитационный потенциал Ф1 отвечает кванту излучения, а Ф2 — поглощения.
+
В случае, когда источник излучения и наблюдатель, расположенные на расстоянии г друг от друга, сближаются с ускорением аг вдоль оси, проходящей через их центры, ДФ = гаг,
Для дальнейшего прочтения статьи необходимо приобрести полный текст. Статьи высылаются в формате PDF на указанную при оплате почту. Время доставки составляет менее 10 минут. Стоимость одной статьи — 150 рублей.