Я тут порыл датасет ELODIE от 2003 для HD 23596 и есть кое-что интересное, имхо.
Модель с одной планетой даёт RMS 9.2 m/s (у меня в ASE 8.9 m/s). В публикации Перье там вообще 5 строчек об этой звезде. Упоминается о завышенном RMS, но они даже остатки не прокомментировали толком. Я сделал 100 прогонов GA/SQP в ASE без использования периодограммы и наткнулся на второй 189-дневный компонент, массой в 0.34 MJup. На периодограмме видна небольшая горка на 190 днях - выше, чем остальные пики, но не сильно. Модель с двумя планетами даёт RMS 5.9 m/s, что больше подходит для звезды главной последовательности класса F8V.
Далее - FAP тест: я взял остатки от модели с одной планетой и нагенерил 1000 синтетических сетов, перемешивая их случайно. Далее на каждый сет прогенерировал модель с двумя планетаи. Почти для каждого сета у меня Хи получается больше, чем для реальных данны. FAP получился 0.3%. Я использовал "быструю" процедуру, только минимизацию на каждый сет. Сейчас считаю более реальный FAP с использованием 50 прогонов GA/SQP на каждый сет - пока сделал 20 прогонов и FAP уже <5% получается. Нужно ещё дней 10 считать.
Ну и самое интересное. В опубликованном Перье спектре HD 23596 нет эмиссии ионизированного кальция, а это значит, что звезда малоактивна. Это только подтверждает то, что остатки должны быть похожи на остальные звёзды главное последовательности, т.е. около 5 m/s, да? В 2003 точности ELODIE не хватало, чтобы замерить метрики активности - упоминается в публикации. И ещё - в 2006 некто Саффе опубликовал их. Используя S-index от Саффе и B-V, Teff, MV от Перье и формулы Райта от 2005 можно посчитать метрик активности дельта-FCaII, который у меня получился 0.38. Так вот, в публикации Райта (2005) звёзды с такими показателями должны иметь расчётный джиттер около 5.3 м/c (80% вероятность), что близко к моей модели с двумя планетами (5.9 м/c).
Быстрый тест на орбитальную стабильность (гибридный BS и стандартный BS в MERCURY6) на 10000 лет - положительно. MEGNO индикатор буду считать как только закончу с FAP. Если FAP будет <1%, то я думаю это будет интересно.
Такие вот пироги. Поругайте кто-нибудь, пожалуйста.
