CC BY
32
УДК 524.35-337
ОСОБЕННОСТИ СТРУКТУР МАГНИТНЫХ ПОЛЕЙ ХИМИЧЕСКИ
ПЕКУЛЯРНЫХ ЗВЁЗД
© 2012 Ю. В. Глаголевский
Специальная астрофизическая обсерватория РАН, Нижний Архыз, 369167 Россия Поступила в редакцию 26 января 2012 г.; принята в печать 11 апреля 2012 г.
Продолжена работа по программе исследования крупномасштабных структур магнитных полей химически пекулярных звезд. В данной работе исследовано восемь звезд, из которых три имеют структуру центрального диполя, три — структуру диполя, смещенного вдоль оси, и две — диполя, смещенного поперек оси. Для двух звезд, ИЭ 62140 и ИЭ 71866, имелись измерения с высокой точностью — а = 50 и 80 Гс. Оказалось, что модельные фазовые зависимости совпадают с измерениями в пределах точности измерений. Этот результат показывает, что предположение дипольной структуры магнитных полей СР-звезд достаточно обосновано.
Ключевые слова: звезды: химически пекулярные—звезды: магнитное поле
1. ВВЕДЕНИЕ
В этой работе мы продолжаем программу исследования структур магнитных полей химически пекулярных (СР) звезд с целью выяснения основных свойств и параметров крупномасштабных магнитных конфигураций. Эти данные нужны для понимания природы магнитных СР-звезд и особенностей их эволюции, а также для исследования влияния магнитного поля на распределение химических аномалий по поверхности звезд. Предварительные результаты статистических исследований около 50 звезд обсуждались в работах [1—3]. Ценность наших исследований структур магнитных полей состоит в том, что они выполняются единой методикой. С учетом результатов данной работы число изученных объектов приближается к 60. Так как наблюдается большое разнообразие особенностей структур магнитных полей, планируется продолжение аналогичных исследований для получения по возможности полной информации о свойствах изучаемого класса объектов. Для исследования отбираются прежде всего те звезды, для которых имеется достаточное число измерений магнитного поля и известны необходимые параметры. В работах [1—3] приводится подробное описание методики моделирования и сделан ряд интересных выводов и предположений, которые были получены на основе анализа выборки около 50 магнитных звезд. Главный результат состоит в том, что все исследованные звезды оказалось возможным промоделировать при предположении виртуального магнитного диполя. При этом для согласования модельной и наблюдаемой фазовых зависимостей изменения магнитного поля с фазой
периода вращения звезды достаточно было манипулировать только положением диполя внутри звезды и величиной магнитного момента, соответственно для подбора формы фазовой зависимости и ее амплитуды. Различие между измеренными величинами и модельной зависимостью практически всегда находится в пределах ошибок измерений (±1—2а). По схеме всех предыдущих публикаций по моделированию, рассматриваются поочередно параметры выбранных объектов и приводятся данные об особенности каждой модели. По окончании цикла работ по моделированию предполагается провести статистический анализ полученных моделей и их параметров с целью приближения к разгадке природы и эволюции магнитных полей химически пекулярных звезд, а также к уточнению связи магнитного поля и диффузии химических элементов. В случае каких-либо особенностей распределения магнитного поля по поверхности в статьях всегда приводятся карты его распределения.
Ниже приводятся результаты моделирования очередной выборки магнитных звезд, изложенные по схеме, выработанной в предыдущем цикле аналогичных работ.
2. ИЗУЧЕННЫЕ ОБЪЕКТЫ
а. НО 49976 (БгСгБи). Оценки магнитного поля взяты из работы [4], где оно измерялось фотографической техникой преимущественно по линиям Ёе и Сг. Сильная переменность отмечается только у линий Са II 3934 А и Бг II 4215 А. Линии остальных химических элементов слабо переменны, поэтому можно предположить, что магнитное
Результаты расчетов моделей
Знак Долгота Широта вр, Bs, /3, Да,
монополя монополя монополя Гс Гс град. Д*
HD 49976 А - 236 -4 ±2245 1960 86 0
+ 56 +4
В - 270 ±2 —4 ±1 +2525+ 150 1315 + 10 86 + 1 0.3
+ ЗО ±2 +4± 1 -2532+ 150
HD 51418 А + 340 + 11 + 1912 : 1126 : 79 : 0
- 160 -11
HD 54118 А (г = 7°) - 330 + 1 + 1911 : 1055 : 89 0
+ 150
А (г = 14°) - 330 + 1 +9719 : 5400 : 89 0
+ 150
HD 62140 А + 340 +2 +2382 1300 88 0
- 160 -2
В + 340 ±2 0± 1 +2019 + 200 1566 90 + 1 0.045(-) +0.005
- 160 ±2 0± 1 -2639 + 200
HD 64740 А + 180 -15 +1294 841 75 0
- 0 + 15
В + 180 ±2 + 16 ±1 + 1630+ 100 849 + 50 74 + 1 0.20(+) + 0.05
- 0 ± 2 -16 ±1 -1200+ 100
HD 71866 А + 0 +3 +5636 3323 87 0
- 180 -3
В + 355 ± 2 +2 ± 1 +5153 + 300 3470 + 170 88 + 1 0.05(+) + 0.01
- 175 ±2 —2 ± 1 -7017 + 350
HD 79158 А - 0 +8 +2998 1798 82 0
+ 180 -8
В - ЗО ±2 +2 ± 1 +3168+150 1762+ 100 88 + 1 0.1
+ 100 ±2 —2 ± 1 -3178+150
HD 108662 А - 72 ±5 +50 ± 1 + 1600+ 150 1040 + 150 40 + 1 0
+ 252 ± 5 -50 ±1
поле мало искажается вследствие неравномерного распределения элементов. В этой же работе даны необходимые для моделирования данные: v sin і = 31 км/с, начальная фаза принята в момент JD(positive crossover) = 2441298.76 + 2 d 976 E, когда Be = 0 Гс. По данным работы [5] Te = 9200 K,
абсолютная болометрическая звездная величина Мь = 0.9, откуда радиус звезды К = 2.2К©. Экваториальная скорость V = 50.6К/Р = 37 км/с, поэтому г = 56°.
Результаты расчетов моделей для этой и всех остальных звезд собраны в таблице. В первом
Я,
Рис. 1. Меркаторская карта распределения магнитного поля по поверхности ИЭ 49976.
Phase
Рис. 2. Модель смещенного поперек оси диполя
ИЭ 49976.
столбце обозначены варианты расчетов: А — центральный диполь, В — смещенный диполь. Величина смещения Да в долях радиуса звезды указана в последнем столбце. Знак указывает направление смещения: “+” — в сторону положительного монополя, “—” — в сторону отрицательного. Вр — величина поля на магнитных полюсах, В8 — среднее поверхностное поле (В8 (тах) + В8 (тт))/2, в — угол между осью диполя и осью вращения звезды.
Фазовая зависимость Ве (ф) меняется в пределах от +2000 до —1800 Гс, это определяет широту обоих магнитных полюсов 5 = ±4°, т.е. они находятся вблизи экватора вращения.
Модельная фазовая зависимость лучше соответствует наблюдениям при предположении смещенного поперек оси диполя. Меркаторская кар-
Phase
Рис. 3. Модель центрального диполя ИЭ 51418.
та распределения магнитного поля по поверхности показана на Рис. 1. Параметры такой модели приведены в таблице, а фазовая зависимость на Рис. 2. Монополи оказались смещенными от центра звезды на Да = 0.3К*, причем расстояние между монополями получилось равным I = 0.5К* (Значительная величина I наблюдается у всех звезд со смещенным поперек оси диполем). Параметр в определяется в случае модели смещенного диполя как угол между плоскостью, в которой находится диполь, и осью вращения.
Координаты магнитных полюсов в случае модели центрального или смещенного вдоль оси диполей совпадают с координатами магнитных монополей, а в случае смещения диполя поперек оси координаты магнитных полюсов другие. У звезды
Рис. 4. Модель центрального диполя HD 54118.
HD 49976 они следующие: отрицательный полюс Л = 265°, 5 = —5°; положительный полюс Л = 35°, 5 = +5°.
Точность полученных параметров зависит от точности определения угла i. Поскольку она, как правило, неизвестна, то ошибки определяются при предположении неточности, равной 5°.
b. HD 51418 (SrCrEu). Магнитное поле звезды измерялось по линиям металлов фотографическим методом в работе [6], в этой же работе сделана оценка v sin i = 15 км/с. Период вращения P = 5d4379 определен в [7]. Температура звезды Te = 9500 к, абсолютная болометрическая величина Mb = —0.4 [5]. Отсюда мы определили радиус звезды R = 4R©. Далее делаем оценку v = 37 км/с, что дает угол наклона звезды к лучу зрения i = 24°. Измерений магнитного поля мало и они неточны, поэтому рассчитана только модель центрального диполя, которую можно рассматривать как первое приближение. Параметры модели приведены в таблице. На Рис. 3 кружками показаны измеренные величины поля, а сплошной линией — модельная фазовая зависимость.
c. HD 54118 (Si). Магнитное поле измерено по линиям На [8] и Н@ [9], поэтому модель не искажена эффектом неравномерного распределения химических элементов. В этой же работе приведен момент начальной фазы JD(max Be) = 2446831.251 +3d27533E. Температура звезды Te = 10200 к [5], абсолютная болометрическая величина Mb = 0.15, откуда радиус звезды R = 2.6R©. К сожалению, неизвестно точное значение v sin i. В каталоге [10] приводится v sin i = 0 км/с, что соответствует ориентации
Рис. 5. Модель смещенного вдоль оси диполя HD 62140.
оси звезды точно на наблюдателя. Если бы это было так, то мы не наблюдали бы переменности поля. Средняя ошибка определения v sin i обычно порядка 10 км/с, следовательно можно предположить, что реальное значение v sin i находится в пределах 0—10 км/с. Наши расчеты мы сделали для двух величин v sin i = 5 и 10 км/с. Так как v = 50.6R/P = 40 км/с, получаем i = 7° и 14°. В Таблице даны параметры моделей центрального диполя для этих двух случаев.
Из полученных моделей видно, что Bp и Bs сильно различаются, однако ориентация диполя в звезде остается одинаковой. Измеренная фазовая зависимость Be(<p) обозначена на Рис. 4 точками, сплошной линией показана вычисленная зависимость для i = 7°. Из данных таблицы видно, что при малых значениях угла i ошибки приводят к неуверенным результатам.
d. HD 62140 (SrCrEu). Магнитное поле измерено LSD-методом в работе [11 ] с высокой точностью, там же даны v sin i = 23 км/с и P = 4d28. В каталоге [5] приведены Te = 10 000 K и абсолютная болометрическая звездная величина Mb = 0.4, откуда R = 2.3R©. Экваториальная скорость v = 27 км/с и i = 58°. Модель центрального и смещенного диполей получается с параметрами, приведенными в таблице. Наблюдения (точки) и модель для смещенного диполя (сплошная линия) приведены на Рис. 5. Отклонения измеренных величин от модельных оказываются в пределах ошибок измерений а = 50 Гс, что подтверждает правильность дипольного представления магнитной конфигурации.
800 -
Phase
Рис. 6. Модель смещенного вдоль оси диполя
ИБ 64740.
e. HD 64740 (He-r). Магнитное поле измерено по линии Ив в работах [12, 13] а также
по линии Не 5876 А [12]. Это позволяет считать, что модель не подвержена влиянию неравномерного распределения химических элементов. Максимальная интенсивность линий Не наблюдается в момент ЛБ 2444611.859 + 1 ?33026Е [12]. В работе [13] приведены следующие параметры звезды: V 8ш г = 160 км/с, К = 5К© , г = 57°. В работе [11] К = 3.2-6.4К©, г > 41°, в = 76° .В [14] Те = 25 000 К, К = 6.1К©.
В [5] Те = 23 850 К, Мь = —5.2, откуда К = 5.5К©, V = 50.6К/Р = 209 км/с, а г = 50°, что примерно соответствует результатам других авторов. Параметры модели центрального и смещенного вдоль оси приведены в таблице.
Много лучше модельные данные соответствуют наблюдениям при предположении смещения диполя вдоль оси в сторону положительного монополя. Поскольку диполь смещен вдоль его оси, то координаты магнитных полюсов совпадают с координатами монополей. Фазовая зависимость приведена на Рис. 6.
Зная карту распределения магнитного поля, можно сопоставить с ней фазовые изменения интенсивности линий С IV и Б1IV, приведенные в [14]. Оказалось, что интенсивности линий С IV и Б1IV достигают минимальной величины в областях обоих полюсов. Двойная волна максимумов интенсивностей показывает, что углерод и кремний концентрируются в виде кольца, совпадающего с магнитным экватором. В работах [15—17] утверждается, что Б1 должен скапливаться путем диффузии на магнитном экваторе, где магнитные силовые линии
2500 -
Phase
Рис. 7. Модель смещенного вдоль оси диполя
ИБ 71866.
горизонтальны поверхности. Это было подтверждено в [18] на примере звезды ИБ 21699. Очевидно, что С и Б1 у ИБ 64740 тоже концентрируются под влиянием лучевого давления в областях с горизонтально направленными магнитными силовыми линиями.
f. HD 71866 (SrCrEu). В данном случае мы используем высокоточные измерения поля из [11], не подверженные влиянию неравномерного распределения химических элементов. Начальная фаза определяется из эфемериды ЛБ(тахНе) = 2438297.5 + 6?80022Е. Параметры V 8ш г = 14 км/с [11], Те = 8650 Ки Мь = 0.2 из [5], откуда К = 3.5К©, V = 26 км/с и г = 32°. Две полученные модели — для центрального диполя и, более точная, для смещенного вдоль оси диполя — приведены в таблице. Фазовая зависимость Бе(ф) для второго варианта приведена на Рис. 7. Отклонения измеренных величин поля от модельных находится в пределах ошибок, что подтверждает правильность дипольного представления магнитной конфигурации.
g. HD 79158 (He-w). Измерения магнитного поля выполнены по линиям Ив в работах [19, 20], вследствие чего эффект неравномерного распределения элементов отсутствует. Величина Те = 12 800 К и абсолютная болометрическая звездная величина Мь = —1.6 взяты из [5]. Отсюда К = 3.7К©, V = 49 км/с, г = 32°. С этими данными вычисляем модель центрального и смещенного поперек оси диполя. Результаты приведены в таблице. Координаты отрицательного магнитного полюса: Л = 355°, 5 = +2°, положительного: Л = 137°, 5 = —2°, расстояние между монополями
Phase
Рис. 8. Модель смещенного вдоль оси диполя
ИЭ 79158.
I = 0.12К*. Фазовая зависимость Бе(ф) приведена на Рис. 8.
h. HD108662 (SrCrEu). Измерения выполнены фотографическим методом по линиям металлов [21 ], фазы взяты там же, они вычислены для периода вращения Р = 5?0808, определенного в [22]. Параметры Те = 9900 К, Мь = +0.6, К = 2.2К© взяты из [5]. По V 8ш г = 19 км/с [10] и V = 50.6К/Р = 22 км/с определяем г = 60°. Разброс измерений Бе велик, поэтому рассчитана только модель центрального диполя (Рис. 9).
3. ЗАКЛЮЧЕНИЕ
Характерной особенностью рассматриваемой серии работ является то, что измерения магнитных полей делались разными авторами с разной аппаратурой. Это, конечно, ухудшает качество получаемых моделей. Во многих случаях ошибки измерений столь велики, что невозможно определить тип структуры — центральный диполь, смещенный диполь или более сложная конструкция. В таких случаях модель строилась в предположении диполя в центре звезды, например, как в случае ИЭ 51428. Параметры в таком случае оказываются неточными, но они обычно не отличаются в среднем больше чем на 20% от точных результатов при смещении диполя не более 0.1К*. Этот вывод мы делаем на основе того, что для каждой звезды вычисляем модель, наилучшим образом соответствующую наблюдениям, и модель центрального диполя, и можем видеть различия.
Из восьми исследованных звезд три промоделированы с помощью модели центрального диполя,
-0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 Phase
Рис. 9. Модель центрального диполя HD 108662.
две звезды имеют смещенные поперек оси диполи, у трех звезд диполи смещены вдоль оси. Две звезды HD 62140 и HD 71866 были измерены с высокой точностью (а = 50 и 80 Гс) и оказалось, что модельные фазовые зависимости согласуются с наблюдениями в таких же пределах. Этот результат еще раз подтверждает правильность предположения о дипольном характере структур магнитного поля у CP-звезд. Это означает, что распределение магнитного поля и направления силовых линий по поверхности и в глубине звезды именно такое, какое создает гипотетический магнитный диполь. Этот вывод мы делаем также по результатам исследований более 50 магнитных звезд, выполненных нами к настоящему времени. На основании данного исследования и указанных предыдущих работ видно, что для подавляющего большинства звезд форму фазовой зависимости определяет только положение диполя внутри звезды.
Сравнение фазовых зависимостей изменения интенсивности линий CIV и Si IV, приведенных в [14] для звезды HD 64740, с распределением магнитного поля показало, что углерод и кремний концентрируются на магнитном экваторе. Это подтверждает расчеты [15—17], предполагающие, что диффузия этих химических элементов происходит под действием лучевого давления в магнитном поле.
СПИСОК ЛИТЕРАТУРЫ
1. Yu. V. Glagolevskij and E. Gerth, in Proc. of Int. Conf. “Magnetic Stars", Nizhny Arkhyz, Russia, 2010 (SAO RAS, Nizhny Arkhyz, 2011), p. 147.
2. Yu. V. Glagolevskij, Astrophysical Bulletin 66, 144 (2011).
3. Yu. V Glagolevskij, Astrophysics 54, 231 (2011).
4. C. A. Pilachowski, W. K. Bonsack, and S. C. Wolff, Astronom. and Astrophys. 37, 275 (1974).
5. Yu. V. Glagolevskij, Bull. Spec. Astrophys. Obs. 35, 33 (2002).
6. T. J. Jones, S. C. Wolff, and W. K. Bonsack, Astrophys. J. 190,579(1974).
7. B. Musielok and J. Madey, Astronom. and Astrophys. 202, 143(1988).
8. E. F. Borra and J. D. Landstreet, Publ. Astronom. Soc. Pacific 87,961 (1975).
9. D. A. Bohlender, J. D. Landstreet, and I. B. Thompson, Astronom. and Astrophys. 269, 355 (1993).
10. A. Uesugi and I. Fukuda, Contrib. from the Inst. of Astrophys. and Kwasan Obs. 189, 205 (1970).
11. G. A. Wade, J.-F. Donati, J. D. Landstreet, and S. L. S. Shorlin, Monthly Notices Roy. Astronom. Soc. 313,851 (2000).
12. D. A. Bohlender, J. D. Landstreet, D. N. Brown, and I. B. Thompson, Astrophys. J. 323, 325(1987).
13. E. F. Borra and J. D. Landstreet, Astrophys. J. 22S, 809 (І979).
14. S. N. Shore and D. N. Brown, Astrophys. J. 365, 665 (1990).
15. S. Vauclair, J. Hardorp, and D. M. Pederson, Astrophys. J. 227,526(1979).
16. G. Alecianand S. Vauclair, Astronom. and Astrophys. 101, 16(1981).
17. C. Megessier, Astronom. and Astrophys. 13S, 267 (1984).
18. Yu. V. Glagolevskij, A. V. Shavrina, J. Sylvestre, et al. Astrophysical Bulletin 64, 166(2009).
19. D. A. Bohlender, IAUS 162, 155 (1994).
20. E. F. Borra, J. D. Landstreet, and I. B. Thompson, Astrophys. J. Suppl. 53, 151 (І983).
21. Yu. S. Rustamov, A. N. Khotnyanskij, Sov. Astron. Lett. 6,202(1980).
22. G. Preston, K. Stepien and S. C. Wolff, Astrophys. J. 166,653(1969).
Specificities of Magnetic Field Structures in Chemically Peculiar Stars
Yu.V. Glagolevskij
We continue our program of the study of large-scale structures of magnetic fields in chemically peculiar stars. In this paper we analyze eight stars, out of which three stars have the structure of a central dipole, three—the structure of a dipole shifted along the axis, and two—of a dipole shifted across the axis. High-precision measurements (with а = 5O and SO G) are available for two stars (HD 62140 and HD 71866, respectively). The model phase dependences agree with the measurements within the errors. This result shows that the hypothesis about the dipole structure of the magnetic fields of CP stars is well founded.
Keywords: stars: chemically peculiar—stars: magnetic field
