УДК 524.827-77
ОСЕВЫЕ СИММЕТРИИ В ДАННЫХ WMAP ^
© 2010 В. С. Беркутов1*, Я. В. Найден1**, О. В. Верходанов2***
1Санкт-Петербургский государственный университет, Петродворец, Университетский проспект, 28,
Санкт-Петербург, 148504 Россия 2Специальная астрофизическая обсерватория, Нижний Архыз, 369167 Россия Поступила в редакцию 5 октября 2009 г.; принята в печать 12 декабря 2009 г.
Предлагается тест для проверки гауссовости фонового сигнала в зонах, симметрично расположенных относительно экватора. Тест основан на исследовании доминирующих гармоник в симметричных зонах и распределения корреляционных коффициентов в зависимости от номера выделенного мультиполя.
В галактической и эклиптической системах координат выбраны симметрично расположенные зоны и исследованы на наличие антисимметричного сигнала. Для нескольких зон в распределении корреляционных коэффициентов обнаружены негауссовы свойства.
Ключевые слова: радиоисточники
1. ВВЕДЕНИЕ
Космическая миссия WMAP1 (Wilkinson Microwave Anisotropy Probe) посвящена исследованию реликтового излучения (CMB — Cosmic Microwave Background). Спутник обращается вокруг точки Лагранжа L2 системы Солнце — Земля и измерил распределение температуры микроволнового фона на небесной сфере за 7 лет. Команда WMAP предоставила астрономическому сообществу исходные данные и результаты обработки годичного, трехлетнего, пятилетнего и семилетнего циклов измерения температуры CMB [1—7]. Для восстановления сигнала CMB из многочастотных наблюдений в этих работах применялся метод внутренней линейной комбинации (ILC — Internal Linear Combination) фоновых компонент, в результате которого была получена карта реликтового излучения, также называемая картой ILC и используемая для анализа низких гармоник фоновых компонент с номерами мультиполей £ < 100. Для построения карты ILC брались данные в пяти наблюдательных каналах: 23 ГГц (полоса К), 33 ГГц (полоса Ka), 41 ГГц (полоса Q), 61 ГГц (полоса V) и 94 ГГц (полоса W).
Полная карта микроволнового фона неба на сравнительно низких мультиполях (£ < 100) чрезвычайно интересна для исследования гауссовости распределения фонового сигнала. Предоставление
E-mail: dair-targ@inbox.ru E-mail: yarnaid@mail.ru E-mail: vo@sao.ru 'httpi/Aambda.gsfc.nasa.gov
данных в общее пользование привело к большому числу работ по исследованию статистики сигнала в различных дипазонах мультиполей. Как было отмечено в ряде работ [8—10], имеются серьезные аргументы, говорящие о том, что в данной карте на различных угловых масштабах присутствует остаточный вклад фоновых компонент, который, в частности, и дает обнаруживаемую негауссовость. Этот вклад может проявляться также в обнаруженной ранее связи в квадруполе между очищенной картой микроволнового фона и галактическими компонентами излучения [11 — 13]. Негауссовость сигнала в данных WMAP, кроме фазового анализа [8,9, 14], была также обнаружена другими методами на различных угловых масштабах: с помощью функционалов Минковского [15], вейвлет-анализа [16], корреляционных функций [17], как асимметрия в распределении сигнала [18], а также другими методами.
В некоторых работах [19, 20], посвященных проблемам негауссовости, исследуется симметрия в данных карты 1ЬС. Так, в работе [19] анализируется хиральность крупномасштабной анизотропии СМВ, включающей четность и нечетность карты, обусловленной возможным подавлением четных гармоник. В ней обсуждается “видимость” нечетности карты и делается вывод, что обнаруживаемая ось симметрии обусловлена положительным (четным) зеркальным отражением. Авторы [20] исследовали две схожие области на карте, положение которых почти центрально симметрично и определяется пятнами октуполя. Подобная симметрия негауссова и авторы обсуждают возможные
объяснения эффекта с помощью нетривиальной пространственной топологии.
Ранее мы также проверяли достоверность визуальных эффектов в карте 1ЬС [21], исследуя проявление симметричных структур в зависимости от уровня шума, добавляемого к карте нечетных гармоник. Как уже упоминалось в работе [19], канал
V в данных WMAP содержит наиболее сильный сигнал реликтового фона относительно фоновых компонент излучения Галактики. При разделении гармоник на четные и нечетные (четность определяется четностью суммы номера мультиполя и его моды: £ + т) большая часть галактического фона попадает в четные гармоники в силу четности ди-польной компоненты нашей Галактики. В нечетных гармониках прослеживается фоновое излучение. Сказанное демонстрируется Рис. 1.
Кроме того, были смоделированны шумовые карты с целью определения уровня шума, при котором симметрия перестает выделяться визуально. Как было ранее продемонстрировано [21], добавление шумового сигнала на уровне 70% от основного, определяемого случайными гауссовыми возмущениями в ЛCDM-космологии, к чистому
2
антисимметричному сигналу , задаваемому только нечетными гармониками, “замывает” визуально выделяемый сигнал (Рис. 2), в то время как при более низком уровне добавляемого шума антисимметрии обнаруживаются. Уровень корреляции карт в зависимости от вклада шума показан на Рис. 3.
2. АНТИСИММЕТРИЯ ОТНОСИТЕЛЬНО ЭКВАТОРА
Напомним, что для исследования карт используется их разложение по сферическим гармоникам (мультиполям):
^ Ш=£
АТ (9,ф) = ЕЕ а£т'^£т(^>ф') ■ (1)
1=2 т=-1
Здесь АТ (в, ф) — вариации температуры фона в полярных координатах, £ — номер мультиполя, т — номер моды мультиполя.
где Pm(x) — присоединенные полиномы Лежандра. Для непрерывной функции AT(x, ф) коэффициенты разложения am выражаются как
YU0, Ф) = y(2^1>(< + m)!PrWe‘m>- |2)
x — cos в ,
г2ж
a£m — J dx
AT(x, ф)Уе*т(х, ф)йф, (3)
Мы называем сигнал антисимметричным, когда он при отображении север-юг относительно экватора меняет знак на противоположный.
где У1т обозначет комплексное сопряжение Уцт.
Используя карту нечетных мультиполей, можно выделить антисимметричные области, т.е. такие, которые расположены симметрично относительно галактического экватора, но имеют противоположное поведение сигнала. Наличие подобных зон в карте 1ЬС обсуждалось в [22, 23]. Здесь мы проводим статистический анализ сигнала на выбранных площадках. Сначала обратим внимание, что визуальное сравнение распределения сигнала нечетных гармоник и карты 1ЬС позволяет обнаружить на этой карте такие же зоны. На Рис. 4 указаны отобранные площадки. Необходимо заметить, что, кроме галактической системы координат, другой выделенной сеткой является эклиптическая система координат, которая проявляется в данных WMAP, в частности, в распределении шума. В этой системе координат (СК) мы также визуально выделили несколько зон (Рис. 4), характеристики которых тоже будут проверены. Сразу следует сказать, что выделение антисимметричных зон достаточно условно, так как их форма и положение определяются в первую очередь визуально. Тем не менее, проверка их существования — один из тестов гауссовости данных реликтового излучения, чем бы их происхождение не определялось. Итак, на карте микроволнового фона мы выбрали визуально в галактической системе координат 7 пар зон, а в эклиптической — 11 пар (Таблицы 1 и 2).
В данной работе мы проводим статистическое исследование гармоник, которые определяют свойства выбранных зон. Предлагаемая методика может быть применена также для анализа гауссово-сти карт реликтового излучения на полной сфере.
2.1. Статистический анализ гармоник в выделенных зонах CMB
В стандартном разложении (1) сигнал в каждом пикселе карты определяется суммой всех гармоник. Выделяя зоны фиксированного размера, мы ограничиваем вклад низких гармоник и усиливаем вклад высоких коррелированных мультиполей. Вычитая определенные гармоники, удается выделить основные мультиполи разложения и определить гармонические корреляционные свойства симметричных зон. Используя данный подход, можно исследовать статистические свойства пятен
1
0
Correlation coefficient
Рис. 1. Слева направо сверху вниз: данные V-канала, карта четных гармоник в этом канале, карта нечетных гармоник, карта 1ЬС. Карты построены для диапазона мультиполей £ < 100. Уровни яркости соответствуют диапазону изменения температуры фона: —0.2, 0.2 мК. В карте нечетных гармоник для увеличения контраста уровень задан как —0.1, 0.1 мК.
Noise coefficient
Рис. 3. Зависимость коэффициента корреляции от уровня шума. Ломаными линиями показаны корреляционные зависимости для различных генераций шумовых карт. Штрихами обозначена средняя аппроксима-ционная зависимость в виде 1.047 exp(—0.421x), где x — относительный уровень добавляемого случайного шума.
и уровень значимости их корреляционных (анти-корреляционных) свойств. На Рис. 5 (для галактической СК) и 6 (для эклиптической) демонстрируется поведение корреляционного коэффициент Кг (р), рассчитанного для р-ой пары симметрично расположенных областей карты, построненной для суммы гармоник с исключеним мультиполя £ = в:
Kr(p) =
1
Е (тг(пр„)-тг(пр„)) E(2V(fip.)-Tr(fip.))
N 2 Np
(4)
где &рп и 0,рв — зоны с номером р в северном и южном полушариях соответственно, N — количество пикселов в этих зонах, апр — дисперсия сигнала в р-ой площадке, Тг (Пр) — распределение температуры в зоне с вычтенным вкладом мульполя г:
Tr(в, ф) — Е Е a£mY£ m (в, ф) - Е armYrm(в, ф) •
£=2 m=-£ m=-r
Ps
m
m=r
Рис. 4. Визуально отобранные зоны на карте в галактической (слева) и эклиптической (справа) системах координат, в которых флуктуации сигнала могут иметь противоположный знак. Верхний ряд: только нечетные гармоники, нижний ряд — данные из полного набора мультиполей карты 1ЬС.
Таблица 1. Номера пар зон в галактических координатах, их форма и координаты: “с” - круг (параметры lat/—lat, lon, radius соответственно), “r” — прямоугольник (параметры lat 1 /—lat 1, lon1, lat2/—lat2, lon2 соответственно), где (lat, lon) — галактические широта и долгота и “—lat” — широта симметричной зоны.
Несколько площадок под одним номером обозначают комбинированную зону.
Пара Форма Координаты
1 с 15°/-15°, 160°, 17°
2 с 23°/-23°, 52°, 12°
3 г 75°/-75°, 60°, 60°/-60°, — 15d
г 65°/-65°, 15°, 50°/—50°, —50d
4 с 30°/-30°,-25°, 15°
с 20°/-20°, -40°, 15°
с 13°/-13°, -52°, 12°
5 с 60°/-60°, -90°, 15°
6 с 27°/—27°, —141°, 8°
7 с 15°/—15°, —179°, 4°
В данной работе мы ограничивались разрешением карты 1ЬС на уровне £тах = 100. Таким образом, пробегая весь цикл г = 2,£тах, мы рассчитываем соответствующую корреляцию.
Таблица 2. Номера пар зон в эклиптической системе координат. Обозначения такие же как в Таблице 1. Несколько площадок под одним номером обозначают комбинированную зону.
Пара Форма Координаты
1 с 32°/-32°, 160°, 6°
2 с 37°/-37°, 97°, 12°
3 г 65°/-65°, 78°, 45°/—45°, 37°
4 г 89°/-89°, 135°, 75°/-75°, 40°
5 с 24°/—24°, 46°, 6°
6 с 8°/—8°, 0°, 8°
с 8°/—8°, —15°, 8°
7 с 56°/-56°, -10°, 12°
8 с 31°/-31°, -45°, 14°
9 с 9°/—9°, —40°, 8°
с 10°/-10°, -50°, 9°
с 8°/—8°, -60°, 7°
10 с 6°/-6°, -115°, 5°
с 8°/—8°, -120°, 7°
с 7°/—7°, -129°, 6°
11 с 26°/—26°, —149°, 8°
с 12°/—12°, —149°, 9°
Multipole
Multipole
-1---------------------г~
кЛМг-
_i___________________i_
1O 2O 3O
4O 5O 6O
Multipole
7O 8O 9O 1OO
Multipole
Multipole
Multipole
Рис. 5. Корреляционные коэффициенты пар симметричных зон в галактической системе координат (в данном примере — с номерами 1, 3,4 - сверху вниз), рассчитанные для карт(1) с одним удаленным мультиполем (левый столбец), для карт (2) только с одним мультиполем (правый столбец). Номера соответствующих мультиполей отложены по оси абсцисс. Штриховыми линиями показаны среднее и границы вариаций коэффициентов на уровне 1а для аналогичных моделированных данных в ЛCDM-космологии.
Для оценки допустимого интервала разброса корреляционных коэффициентов, мы сгенерировали 100 случайных карт, соответствующих гауссовым возмущениям плотности в ЛCDM-космологии используя процедуры пакета GLESP [24, 25]. На Рис.5 и 6, 1 ст-уровень отклонений от среднего, рассчитанный для случайных флуктуаций, показан пунктирной линией. Сплошная линия демонстрирует зависимость коэффициента корреляции от номера вычитаемого мультиполя для реальных
данных. Кроме того, на тех же рисунках показано поведение корреляции зон, расчитанной для одиночного мультиполя. При исследовании полученных зависимостей среди выделенных площадок обращает на себя внимание зона 1, отмеченная в галактической системе координат. Ее особенность состоит в том, что при удалении из карты определенных мультиполей (£ = 5,6, 7), антикорреляция симметричных площадок возрастает и достигает уровня 2а, что уже является
10 20 30 40 50 60 70 80 90 100
Multipole
-1----------------------1—
_i___________________і____
10 20 30 40 50 60 70 80 90 100
Multipole
Multipole
10 20 30 40 50 60 70 80 90 100
Multipole
Multipole
Рис. 6. Корреляционные коэффициенты пар симметричных зон в эклиптической системе координат (в данном примере с номерами 5, 6, 7 — сверху вниз), рассчитанные для карт (1) с одним удаленным мультиполем (левый столбец), карт (2) только с одним мультиполем (правый столбец). Номера соответствующих мультиполей отложены по оси абсцисс. Штриховыми линиями показаны среднее и границы вариаций коэффициентов на уровне 1а для аналогичных моделированных данных в ЛCDM-космологии.
0.5
существенным показателем антисимметричности (и, как следствие, негауссовости) этих зон. Эти гармоники, подавляющие корреляции, сформированы сигналом из эклиптической плоскости (Рис. 7). Корреляции отмеченных зон в отдельных гармониках позволяют выявить угловые масштабы симметричного сигнала, дающего максимальный вклад в исследуемые площадки (правая колонка на Рис. 5). Интересно, что в симметричных зонах
наблюдается вклад низких (£ < 20) нечетных гармоник.
В площадках, выделенных аналогичным образом в эклиптической системе координат, обнаруживаются особенности, не наблюдаемые в симметричных зонах галактической системы. Так, для некоторых областей в пятнах присутствует слабый коррелированный сигнал, почти не зависящий от удаляемых мультиполей (см. Рис. 6, зоны 6 и 7). Несмотря на противоположный знак в зонах 7 в
Рис. 7. Карта мультиполей £=5, 6 и 7 из данных CMB ^С в галактической СК. На изображение наложена координатная сетка эклиптической системы.
эклиптической системе, коррелированный сигнал является положительным, что говорит о разных уровнях фона у этой пары областей карты.
Развитием предложенного метода является переход от выделения одного мультиполя і = г к диапазону і є [г0, г0 + Дг]:
п = 1
sq + Ar m = r
Tro (в,Ф)= £ £ a£rnY£ m (@, Ф) - £ £ armYrm (
1=2 m=-t rQ m = -r
КФ).
(6)
Здесь г0 и Дг — соответственно скользящий начальный мультиполь (например, зависящий от г) и величина диапазона. На Рис. 8 приведены результаты вычисления коэффициента корреляции симметричных зон в зависимости от удаленных из карт мультиполей для величины диапазона Дг = 10. Гладкими линиями показаны 1ст границы допустимых вариаций корреляции в ЛCDM-модели, полученные по результатам 100 случайных реализаций гауссова процесса.
Данный подход позволяет исследовать поведение комбинации гармоник и их вклада в корреляционные свойства. Отметим сравнительно высокое отклонение статистики сигнала от ожидаемых гауссовых свойств в зонах 1, 2, 5 и 6. В зоне 1 при выделении набора гармоник основная доля корреляционного сигнала при £ < 60 находится вне ожидаемого диапазона гауссовости.
3. ЗАКЛЮЧЕНИЕ
Предложен тест для проверки гауссовости фонового сигнала в зонах, симметрично расположенных относительно экватора. В основе подхода лежит проверка распределения корреляционных коэффициентов в зависимости от выделяемого мультиполя или набора мультиполей.
В галактической и эклиптической системах координат выделены симметричные относительно экватора зоны и исследованы на наличие в них антисимметричного сигнала. В нескольких зонах карты ILC в распределении корреляционных коэффициентов обнаружены негауссовы свойства, выраженные в сущестовании корреляции площадок для ряда гармоник. Приведены координаты этих зон. Симметричные свойства зон относительно экватора предполагают их объяснение возможными эффектами разделения компонент в процедуре чистки CMB и построения карты ILC. Также представляет интерес наличие фоновой корреляции в ряде зон, выделенных в эклиптической системе координат.
Предложенная методика может быть использована для исследования статистических свойств карт и поиска негауссовых зон в данных миссии Planck Европейского космического агенства.
БЛАГОДАРНОСТИ
Мы выражаем благодарность за возможность использования доступного архива NASA (NASA Legacy Archive), откуда были взяты
Multipole
Multipole
Multipole
Multipole
Multipole
Multipole
Multipole
Рис. 8. Корреляционные коэффициенты для семи пар зон (нумерация справа налево и сверху вниз) в галактической системе координат. Зоны выделены в картах, построенных по данным 1ЬС с удалением диапазона мультиполей Дг = 10. Положение начала диапазона задается номером по оси абсцисс.
данные WMAP Мы также признательны за использование пакета HEALPix3 [26], с помощью которого были преобразованы карты WMAP в коэффициенты a,£m. В работе использован пакет GLESP4 второй версии [27] для дальнейшего анализа данных CMB на сфере. Работа частично поддержана грантом РФФИ №09-07-00159. О.В.В. благодарит за частичную поддержку Фонд Содействия Отечественной Науке (программа “Молодые доктора РАН”) и Фонд “Династия”.
СПИСОК ЛИТЕРАТУРЫ
1. C. L. Bennett, M. Halpern, G. Hinshaw, et al., Astrophys. J. Suppl. 148, 1 (2003), astroph/0302207.
2. C. L. Bennett, R. S. Hill, G. Hinshaw, et al., Astrophys. J. Suppl. 148, 97 (2003), astroph/0302208.
3. D. N. Spergel, L. Verde, H. V. Peiris, et al., Astrophys. J. Suppl. 148, 175(2003), astro-ph/0302209.
4. G. Hinshaw, D. N. Spergel, L. Verde, et
al., Astrophys. J. Suppl. 170, 288 (2007),
astro-ph/0603451.
5. D. N. Spergel et al., Astrophys. J. Suppl. 170, 377
(2007), astro-ph/0603449.
6. G. Hinshaw, J. L. Weiland, R. S. Hill, et
al., Astrophys. J. Suppl. 180, 225 (2009),
arXiv:0803.0732.
7. E. Komatsu, J. Dunkley, M. R. Nolta, et
al., Astrophys. J. Suppl. 180, 330 (2009),
arXiv:0803.0547.
8. L.-Y. Chiang, P. D. Naselsky, O. V. Verkhodanov, and M. J. Way, Astrophys. J. 590, L65 (2003), astroph/0303643.
9. P. D. Naselsky, A. G. Doroshkevich, and
O. V. Verkhodanov, Astrophys. J. 599, L53 (2003), astro-ph/0310542.
10. P. D. Naselsky, A. G. Doroshkevich, and
O. V. Verkhodanov, Monthly Notices Roy. Astronom. Soc. 349, 695 (2004), astro-ph/0310601.
11. P. D. Naselsky and O. V. Verkhodanov, Int. J. Mod. Phys. D. 17, 179 (2008), astro-ph/0609409.
12. P. D. Naselsky and O. V. Verkhodanov, Astrophysical Bulletin 62, 203 (2007).
3http://www.eso.org/science/healpix/
4http://www.glesp.nbi.dk
13. P. D. Naselsky, O. V. Verkhodanov, and M. T. B. Nielsen, Astrophysical Bulletin 63, 216
(2008), arXiv:0707.1484.
14. P Coles, P. Dineen, J. Earl, and D. Wright, Monthly Notices Roy. Astronom. Soc.350, 989 (2004).
15. H. K. Eriksen, D. I. Novikov, P. B. Lilje, et al. Astrophys. J. 612, 64 (2004).
16. J. D. McEwen, M. P. Hobson, A. N. Lasenby, and D. J. Mortlock, Monthly Notices Roy. Astronom. Soc. 359, 1583 (2005).
17. R. Tojeiro, P. G. Castro, A. F. Heavens, and S. Gupta, Monthly Notices Roy. Astronom. Soc. 365, 265 (2006).
18. F. K. Hansen, P Cabella, D. Marinucci, and N. Vittorio, Astrophys. J. 607, L67 (2004).
19. K. Land and J. Magueijo, Phys. Rev. D 72, 101302 (2005), arXiv:astro-ph/0507289.
20. V. G. Gurzadyan, A. A. Starobinsky, A. L. Kashin, et al., Mod. Phys. Lett. A 22, 2955 (2007), arXiv:0709.0886.
21. Ya. V. Naiden, V. S. Berkutov, and O. V. Verkhodanov, in Actual problems of extragalactic astronomy, XXVI Conf. Abstr. book, (Puschino, PRAO AstroSpace Center Phys. Lebedev Inst., 2009) p. 20.
22. O. V Verkhodanov, P. D. Naselsky, L.-Y. Chiang, et al., in Particles and Cosmology, Proc. 14th Internat. School, Ed. by S. V. Demidov, V A. Matveev, and V A. Rubakov, (Inst. Nucl. Phys., 2008) p. 119.
23. O. V. Verkhodanov, in Practical Cosmology,
V.II, Proc. Internat. Conf. “Problems of Practical Cosmology”, Ed. by Yu. Baryshev, I. N. Taganov, and P Teerikorpi, (St.Petersburg, Russian Geograph. Soc., 2008), p. 206.
24. A. G. Doroshkevich, P. D. Naselsky, O. V Verkhodanov, et al., Int. J. Mod. Phys. D 14, 275 (2003), astro-ph/0305537.
25. O. V. Verkhodanov, A. G. Doroshkevich, P D. Naselsky, et al., Bull. Spec. Astrophys. Obs. 58, 40 (2005).
26. K. Gorski, E. Hivon, A. J. Banday, et al., Astrophys. J. 622, 759 (2005).
27. A. G. Doroshkevich, O. B. Verkhodanov, O. P. Naselsky, et al., arXiv0904.2517 (2009).
AXIAL SYMMETRIES IN WMAP ILC DATA
V.S. Berkutov, Ya.V. Naiden, O.V. Verkhodanov
A test of the Gaussianity of the background signal in the zones located symmetrically with respect to the Galactic equator is proposed. The test is based on the analysis of the dominating harmonics in symmetric zones and the distribution of the correlation coefficients as a function of the number of selected multipole. Zones located symmetrically in the Galactic and ecliptical coordinate systems are chosen and analyzed for the presence of an antisymmetric signal. Nongaussian properties are found in the distribution of the correlation coefficients for several zones.
Key words: radio sources АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 65 №2 2010