Солнце – уникальный для науки объект, единственная звезда, поверхность который мы можем наблюдать. рис.1 Солнце 6 августа 2010г. - презентация

1 Солнце – уникальный для науки объект, единственная звезда, поверхность который мы можем наблюдать. рис.1 Солнце 6 августа 2010г

2 Для физиков – это лаборатория, где вещество хотя и нельзя потрогать, но узнать о нём можно очень много, если его наблюдать и исследовать с помощью спектрального анализа. рис.2 Фраунгоферов спектр

3 Линии поглощения, образующиеся в видимой нами плазме Солнца, зависят от: –температуры, –лучевых скоростей, –магнитных полей, –химического и ионного состава и т.д. Профили линий меняют –интенсивность –уширение –форму (симметрию) –число компонент для линий, расщепленных в магнитном поле –поляризацию –положение в шкале длин волн –соотношение друг с другом в зависимости от параметров атомных переходов рис. 3 Участок спектра 6296 – 6303 АА вблизи пятна для двух направлений круговой поляризации

4 Астрофизики получают информацию, сравнивая профили линий: наблюдаемые и полученные из расчетов по моделям. Разница между наблюдаемыми и расчетными профилями получается от несовершенства алгоритмов расчета неточного знания параметров линий, атомных параметров неточного знания моделей – т.е. математического описания состояния плазмы в наблюдаемой области. рис.4 Cравнение профилей из атласа Гуртовенко и др. с расчетными: Профили линии Ni I A для пяти положений в направлении центр - край на диске Солнца: - наблюдаемые - рассчитанные

5 Таким образом, программа, позволяющая варьировать физические условия в плазме и рассчитывать результирующие профили линий, становится таким же важным инструментом, для изучения Солнца, как и сам телескоп в паре со спектрографом. рис.5 Башенный солнечный телескоп БСТ-1 обсерватории КрАО

6 Программный комплекс SunWorld позволяет рассчитывать профили линий поглощения в присутствии магнитного поля для областей с произвольно заданным распределением с глубиной температуры и других физических параметров. Кроме того, он может восстанавливать ход с высотой газового и электронного давлений по заданному ходу температуры с помощью формул ионизационного и гидростатического равновесий. В настоящее время такой комплекс является единственным в РФ

7 Сложная структура программы Каждое физическое исследование обычно требует индивидуального подхода – нужно считать другой набор параметров, в другом порядке, по-новому организовывать сравнение результатов расчетов и наблюдений и т.д. В связи с этим программа не может быть написана раз и навсегда, она под каждый новый проект требует доработок. Кроме того, в программу включено достаточно сервисных функций: рисование профилей спектральных линий (профилей Стокса в случае присутствия магнитного поля), рисование произвольных графиков, база данных спектральных линий, работа с данными в текстовом виде, в частности, в виде текстовых таблиц, расчет мультиплетной структуры линий и т.д.

8 SunWorld умеет считать структуру расщепления линии в магнитном поле. Можно создать текст с символьной записью обозначений переходов, как на этом примере, и SunWorld посчитает расщепеления и покажет их в графическом виде: lewel_lo level_up At+Ion Lamb b 4D 1/2 z 4Po 1/2 ! FE II S 1/2 2Po 1/2 ! SI II Po 1/2 2S 1/2 ! NA I z 5P 1 e 5D 0 ! FE I z 5P 1 e 5D 1 ! FE I D 3/2 2P 1/2 ! BA II D 3/2 2P 3/2 ! BA II 5853 y 5Fo 1 f 5P 1 ! FE I y 5Fo 1 g 5D 1 ! FE I a 5F 1 z 7Fo 1 ! FE I y 5Do 1 e 5F 1 ! FE I a 4D 1/2 z 4Po 1/2 ! V I a 4D 1/2 z 4Po 3/2 ! V I a 4D 3/2 z 4Po 1/2 ! V I a 4D 3/2 z 4Po 3/2 ! V I рис.6 Структуры расщепления линий, заданных символьными обозначениями атомных переходов При обычных расчетах профилей SunWorld выбирает атомные параметры линий из базы данных и по ним рассчитывает структуры расщепления линий.

9 База данных Для пополнения базы данных - импорта данных о линиях из внешних источников SunWorld имеет специальные полуавтоматизированные процедуры. В SunWorld есть база данных. Наиболее полезная таблица в базе данных - это таблица спектральных линий. Она позволяет задавать разные критерии отбора и в диалоге подбирать линии для той или иной физической задачи. Рис.7 Выборка содержит линии ионизованного железа Fe II, для которых эквивалентная ширина в тени пятна (столбец WSW) превышает 3 мА.

10 Анализ моделей фотосферы SunWorld позволяет графически сравнить модели фотосферы. На рисунке представлены две модели невозмущенной фотосферы: Рис. 8. наиболее ппопулярная модель 1974 г. (HOLMU) уточняющая её модель 1999 г. (GRESA). Приведены графики температуры и газового и электронного давлений – их распределение с оптической глубиной. Как видно, различия между моделями невелики. Тем не менее, они заметно сказываются на расчетах и на параметрах среды. В частности новая модель предполагает на 40% больше атомов железа на Солнце, чем старая.

11 Табличное сравнение моделей Модели можно сравнивать не только графически, но и таблично. Программа позволяет выбрать любые параметры, которые требуется сравнить и вывести их в таблицу с шагом той модели, которая выводится первой. Табличные данные удобны для дальнейшего использования их в Excel, Matlab и т.п. Рис.9 Сравнение температуры и коэффициента поглощения в линии Fe I 6302 A для моделей HOLMU 1974 и Grevesse, Sauval (з вездочками отмечены экстраполированные значения)

12 Расчет для списка вариантов Чтобы подобрать расчетный профиль так, чтобы он соответствовал наблюдаемому, нужно варьировать входные параметры. Например, фрагмент задания на расчет может выглядеть так: ! СПИСОК ПАРАМЕТРОВ, ЗАДАННЫХ ОДНОКРАТНО: #LIN FE6302 # ! ИДЕНТИФИКАТОР ЛИНИИ #NLAM 61 # ! ЧИСЛО ТОЧЕК ПРОФИЛЯ ! ВАРЬИРУЕМЫЕ ПАРАМЕТРЫ: #A # ! ПОСТОЯННАЯ ЗАТУХАНИЯ - 3 ВАРИАНТА #ETA # ! ОТНОШ.КОЭФ-ТОВ СЕЛЕКТ.И НЕПР.ПОГЛОЩЕНИЙ 4ВАР-ТА ! ЗАПУСК АЛГОРИТМА РАСЧЕТА ПО МОДЕЛИ МИЛНА-ЭДДИНГТОНА: #RUN RACH # В данном случае мы посчитаем 3*4=12 наборов профилей Стокса по аналитическим формулам Рачковского (т.е. для модели Милна – Эддингтона). Профили сохраняться в памяти программы и при желании параллельно будут записаны в текстовый файл. Интегральные характеристики профилей можно посмотреть на специальной странице программы и распечатать. Наконец, в диалоге можно задать параметры графиков и вывести их на экран.

13 Вывод графиков на экран Рис. 10 Профили параметра стокса V (круговой поляризации) для задания с предыдущей страницы Параметр расчета ETA0 ( ), определяющий силу линии меняется слева направо. Параметр затухания A ( ), уширяющий линию меняется сверху вниз.

14 Анимация Рис. 11 Пример анимированных профилей. Изменение профилей линий Fe I 6082 и 6089 AA при перемещении по диску Солнца Профили интенсивности и круговой поляризации линий Fe I 6082 A при изменении напряженности поля от 0 до 3000 Э Если входной параметр расчета менять плавно, можно получить анимированную последовательность профилей. Анимацию можно наблюдать в программе, а можно выгрузить в gif-файл, который затем вставить в электронный документ.

15 Подбор параметров в диалоге Когда нужно делать подбор параметров для множества однотипных случаев, например для сотни линий, то удобнее использовать не матричный расчет, а диалоговый. В этом случае ползунками меняя параметры добиваемся, чтобы разность наблюдаемого и рассчитанного профилей была минимальной. Рис. 12 Две стадии подбора профиля, заданного для 5 положений на диске. Красный график – разность рассчитанного и наблюдаемого профилей. Пример Известны профили примерно ста линий для пяти положений на диске Солнца. Решаем задачу подбора параметров для модели ф/с Милна-Эддингтона

16 Расчет для участка длин волн. Обычно нас интересуют отдельные конкретные линии. Но иногда нужно сопоставлять расчет со спектром, включающем ряд бленд. На этот случай программа умеет считать профили для диапазона спектра. Например, если в задание на расчет вместо обычного задания идентификатора линии #LIN CA6102 # включить задание нижней и верхней границ диапазона длин волн: #LAM_L # #LAM_U # то программа автоматически включит в расчет все линии из базы данных: LAMB E_low lg(gf) J_l J_u g_l g_u SI I FE I CR II /2 3/ FE I SI I CA I TI II /2 9/ FE I V I /2 9/ RH I /2 11/ CA I TI I CR I FE I FE I SM I

17 Результат расчета для участка длин волн Рис. 13 Верхняя половина – участок спектра вблизи линии Ca I 6102 A для невозмщенной фотосферы (вверху) и тени пятна (внизу) Нижняя половина – рассчитанные профили интенсивности для фотосферы и интенсивности и круговой поляризации для тени модель невозмущенной фотосферы – HOLMU 1974 модель тени – Stellmacher-Wiehr 1975 (SW75), напряженность поля 3000 Э HOLMU RI SW75 RI SW75 RV ТЕНЬ невозмущенная фотосфера Ca I 6102 A

18 Текстовый редактор и подсистема графики Встроенный текстовый редактор SunWorld имеет две функции, которые отличают его от обычного блокнота, это: функция гипертекста возможность редактирования таблиц. Встроенная графика делится на несколько подсистем: ряд простых, для рисования векторов моделей фотосферы, структур зеемановского расщепления и две более развитых: система рисования профилей Стокса и графическая система общего назначения

19 Функция гипертекста текстового редактора Благодаря функции гипертекста редактор позволяет организовывать огромное множество текстовых данных в единое целое. Произвольный текст Имена файлов

20 Редактирование таблиц в текстовом редакторе Любой участок текста, который содержит колонки цифр или просто слов, разделенных пробелами, считается таблицей. Для таблиц возможны следующие операции: Выравнивание – все слова выстраиваются в колонки, слова из букв равняются по первому символу, числа по десятичной точке: Сортировка. Сортируется вся таблица именно по той колонке, в которой находится курсор. Если в колонке есть слова, состоящие из букв, происходит сортировка по алфавиту, если колонка содержит числа – то по значениям. Удаление колонки. Текущую колонку можно удалить целиком (вырезать). Запоминание и вспоминание колонки. Колонку можно запомнить и затем вспомнить (получить её дубль) в новой позиции. Запоминание и вспоминание в формате Excel. Таблица целиком может быть запомнена в виде, пригодном для вставки её в Excel. Так же может быть вставлена в текст таблица, ранее запомненная в программе Excel.

21 Подсистема рисования посчитанных профилей Стокса В типовом расчете SunWorld задается несколько входных параметров, которые меняются. Например мы можем посчитать профили для трех линий, для трех значений угла наклона поля и для 20 значений напряженности поля: ( ФРАГМЕНТ ВХОДНОГО ФАЙЛА – ЗАДАНИЯ НА РАСЧЕТ) #LIN FE6089 FE6093 FE6094 # ! три линии #GM # ! три значения угла наклона поля #$H # ! значения поля от 1600 до 3500 с шагом 100 Э Задача подсистемы рисования состоит в том, чтобы распределить графики по желаемым осям X, Y, T(время), выбрать параметр или параметры Стокса (RI, RV и т.д.), задать одинаковый для всех масштаб по длинам волн и по шкале интенсивностей FE6089FE6093FE6094 GM=30 GM=45 GM=60 Напряженность поля H распределена вдоль оси времени T

22 Подсистема графики общего назначения Подсистема графики общего назначения работает с данными сохраненными в виде текстовых таблиц. В специальном текстовом файле мы описываем график, а потом интерпретатор "проигрывает" это описание и получаются графики. Примеры таких графиков приведены ниже: Эта система скромнее, чем например, графика системы MatLab, зато позволяет получать графики оперативно и хранить их описания вместе с файлами данных проекта

23 Подсистема управления проектами В SunWorld могут выполняться десятки самых разнообразных проектов. Для этой цели создана система управления проектами. Данные разных проектов абсолютно независимы, их можно архивировать или удалять с компьютера без ущерба для других проектов. В рамках одного проекта может быть несколько разных расчетов – и однотипных и разнотипных. Профили этих разных расчетов можно сравнивать, они активны в SunWorld одновременно. Все файлы данных в проектах SunWorld (кроме готовых графиков) – это текстовые файлы. Данные одного проекта подразделяются на: Входные данные расчета Файлы, которые содержат результаты расчета – профили Стокса Файлы, которые описывают графики и готовые файлы с графиками

24 Ряд работ, статей сделанных с помощью SunWorld Об отсутствии нетепловой скорости в холодном компоненте тени крупного пятна Труды XIII консультативного совещания по физике Солнца. Одесса. The Upper Temperature Limit of Bright Points in Large Sunspot Umbrae IAU Symposium 138, Kiev May Зависимость профилей Стокса от температуры модели фотосферы 1993г Астрономический журнал (всего 8 работ того периода) Конкурсный проект СПЕКТР – один из первых грантовых проектов Грант – первый грант УАФО, который получил финансирование. Полностью основан на идее использования SunWorld Идея магнитографа для пятен, измеряющего абсолютные значения поля Рукопись Метод "пробного слоя" в вычислении функций вклада и отклика Готовится к публикации Аномальное расщепление линий ванадия V I 6058 и 6111 AA Готовится к публикации Порядка 20 работ А.В.Баранова в соавторстве за последние 3 года.