Документ взят из кэша поисковой машины. Адрес оригинального документа : http://antares.sinp.msu.ru/docs/Orlov%20-%20Term%20paper.pdf
Дата изменения: Mon Dec 16 17:03:25 2013
Дата индексирования: Thu Feb 27 20:27:02 2014
Кодировка: Windows-1251
Московский Государственный Университет им. М. В. Ломоносова Физический факультет

Курсовая работа

Акустические методы детектирования нейтрино

Выполнил студент 413 группы Орлов М. В. Научные руководители к.ф-м.н. Широков Е. В. проф. Деденко Л. Г.

Москва, 2007


1 Введение
1.1 Цель работы
В данный момент в мире конструируются три больших эксперимента, IceCube, NEMO и NESTOR, направленных на поиск космических нейтрино высоких энергий (вплоть до 1020 eV). Такие нейтрино, взаимодействуя с веществом детектора, рождают электромагнитные и адронные ливни. Выделение энергии в этих ливнях приводит к нагреванию вещества, резкому изменению его плотности в области ливней и, в конечном сч?те, к образованию акустических волн. Конечной целью выполняемой работы является создание программного пакета, позволяющего с приемлимой скоростью производить расч?т акустических волн от нейтрино энергий вплоть до 1020 эВ для морского детектора с произвольным рельефом дна.

1.2 Разработанные области
Все нейтрино, детектируемые нашими установками, можно разделить на три большие группы. Во-первых, это солнечные нейтрино, то есть те, которые рождаются внутри Солнца в результате различных последовательных реакций горения. Вторую группу составляют нейтрино, рождающиеся в атмосфере Земли при рассеянии космических лучей. В третью группу мы поместим все остальные источники нейтрино, как космические, так и находящиеся в пределах Земли, и потом подойд?м к этой теме подробнее. Нейтрино обладают двумя важнейшими свойствами. Первое состоит в том, что при движении они, в отличие от космических лучей, не отклоняются электромагнитными полями. Это позволяет точно определять направление на источник. Второе свойство заключается в крайне небольшом сечении взаимодействия нейтрино с веществом. Помимо отрицательной стороны, очевидно следующих из этого проблем с их детектированием, у этого свойства есть и крайне положительная: нейтрино долетают до Земли, проходя сквозь плотные области в космическом пространстве, которые являются непроницаемыми для гамма-квантов и других частиц.

1.2.1 Солнечные нейтрино
Общеизвестно, что нейтрино бывают тр?х типов, или ?ароматов?: электронное (e ), мюонное (ч ) и таонное ( ). В реакциях на Солнце рождаются только электронные нейтрино. При этом, как показывали первые эксперименты, их поток на земной поверхности ч?тко меньше теоретически предсказанного Джоном Бакаллом. Определим для удобства величину R, равную отношению экспериментальных данных и теоретических предсказаний. Итак, R < 1. Первичной реакцией, происходящей внутри Солнца, является горение водорода с образованием дейтерия: p + p d + e+ + e (E 0.42МэВ). Получающийся позитрон аннигилирует с электроном e+ + e- 2 , а дейтерий вступает в реакцию с другим водородом: d + p 3 He + . Далее есть три возможных пути развития событий. 1


Первый состоит в реакции двух ядер гелия-3: 3 He + 3 He + 1 H + 1 H + 12.86МэВ. Он происходит в примерно 70% случаев. Второй и третий пути предвfряются реакцией гелия-3 с гелием-4: 3 He + 7 Be + . Затем второй путь состоит в захвате электрона ядром бериллия 7 Be+ e- 7 Li+ e (E = 0.86МэВ) и последующей реакции бериллия с протоном, приводящей к образованию двух -частиц. Вероятность протекания этого пути 30%. В третьем пути бериллий захватывает протон 7 Be + p 8 B + , а затем получающийся бор испускает позитрон и нейтрино: 8 B 8 Be + e+ + e (E 14.06МэВ). В конце концов, 8 Be переходит в две -частицы. Третий путь имеет вероятность в десятые доли процента. Все эти цепочки называются протон-протонными (?pp-?) циклами, в их результате восемь протонов превращаются в две -частицы, испускаются 23 нейтрино и выделяется энергия примерно в 25МэВ. Другая цепочка реакций, происходящих на Солнце, называется CNO-циклом. Она не вносит существенного вклада в выделение энергии. Из этих фактов можно заключить, что зная солнечную светимость, можно оценить поток нейтрино, который должен достигать Земли. Хорошее представление о важности реакций да?т рисунок 1.2.1. Таким образом, интерес представляют три реакции: pp d (< 0.42МэВ), поглощение электрона 7 Be (0.86МэВ) и распад 8 B (< 14.06МэВ). Для детектирования солнечных нейтрино применимы три типа взаимодействий. Так называемое взаимодействие через зарядовые токи (charged-current) происходит посредством W - -бозонов, это реакции типа e + n p + e- . Они могут быть использованы для детектирования только электронных нейтрино. Нейтрино всех ароматов подвержены упругому рассеянию x + e- x + e- , однако сечение этого процесса для электронных выше в шесть раз, чем для мюонного и таонного. Наконец, сечение взаимодействия через нейтральные токи (neutral-current), идущего с помощью Z -бозона, x + n x + n, одинаково для всех типов нейтрино. В 1968 году американским уч?ным Рэем Дэвисом в шахте Homestake был провед?н первый эксперимент по детектированию солнечных нейтрино. В н?м использовалась предложенная в 1946 году Бруно Понтекорво хлор-аргоновая реакция e + 37 Cl 37 Ar + e- . На глубине в полтора километра под земл?й располагалась цистерна с 610 тоннами C2 Cl4 . При захвате нейтрино отдельные ядра хлора превращались в ядра аргона. 37 Ar не является стабильным, он имеет период полураспада 35 дней и превращается обратно в хлор через захват электрона. Соответственно, принцип проведения эксперимента состоял в накоплении стабильного числа ядер аргона в течение какого-то времени, а затем определения их количества через детектирование распадов с помощью пропорционального сч?тчика. Средняя скорость образования аргона составляла 1 атом за двое суток. Порог применявшейся реакции составляет 0.814МэВ. Соответственно, эксперименты на е? основе чувствительны только к нейтрино от 7 Be и 8 B, а нейтрино от pp d обладают недостаточной для регистрации энергии. Дэвис получил следующий результат: R 0.3. Его значительное отклонение от единицы получило название ?проблемы солнечных нейтрино?. 2


Рис. 1: Кривые представляют собой потоки солнечных нейтрино от разных реакций в зависимости от энергии. Пунктиром отмечены нейтрино, рождающиеся в результате реакций CNO-цикла. Как видно, они имеют низкие энергии и одновременно потоки куда меньше, чем у нейтрино от первичной реакции pp-цикла. Самая нижняя правая кривая, помеченная hep суть поток нейтрино от реакции гелия-3 с протоном. Эта реакция предсказана только теоретически, а предполагаемый поток нейтрино столь мал, что е? можно не принимать во внимание. До этого в течение десятилетия зарождалась теория так называемых нейтринных осцилляций. Впервые идея о возможности превращений нейтрино в антинейтрино (e e ) была высказана в 1957 году Понтекорво[3]. Затем в 1962 году японские уч?ные Маки, Накагава и Саката предложили и возможность изменения аромата нейтрино. Надо заметить, что хлор-аргоновая реакция, очевидно, является примером взаимодействия через зарядовые токи, и, следовательно, с е? помощью можно регистрировать только электронные нейтрино. Таким образом, этот экспериментальный результат дал определ?нное подтверждение теории осцилляций: те семьдесят процентов электронных нейтрино, которые не были зафиксированы в эксперименте Дэвиса, вполне могли превратиться в мюонные и таонные, пока летели от Солнца до Земли. Эта теория, однако, 3


явно противоречила принятой в физике стандартной модели, которая безусловно считает нейтрино безмассовыми частицами и тем самым налагает запрет на возможность осцилляций. Соответственно, полученное Дэвисом значение R 0.3 дало мощный толчок для постройки новых детекторов с целью проверки его результатов. Следующая серия экспериментов баксанский SAGE (1990), а также GALLEX (1991) и GNO (1998) отличалась тем, что была построена на основе галлий-германиевой реакции e + 71 Ga 71 Ge + e- . В остальном принцип проведения эксперимента был таким же. Эта реакция также ид?т через зарядовые токи, но имеет предел в 0.233МэВ, и, соответственно, она чувствительна к нейтрино от солнечной реакции pp d. Результат, который дали эти три эксперимента, такой: R 0.6. Это позволило придти к выводу, что уменьшение потока электронных нейтрино не только действительно присутствует, но и зависит от энергий нейтрино. Третий тип экспериментов водные черенковские детекторы Kamiokande (1986 1995) и пришедший ему на смену SuperKamiokande (19962001). Они построены на реакции упругого рассеяния x + e- x + e- и представляют собой огромный подземный бак с водой, окруж?нный фотоумножителями. Электроны после рассеяния могут приобрести скорость, большую скорости света в воде, и таким образом начать излучать черенковский свет. Этот свет детектировался фотоумножителями. На каждый из них при этом попадало примерно по одному фотону. В отличие от хлоровых и галлиевых экспериментов, детектирование происходило в реальном времени, что позволило получить информацию об энергиях и направлениях движения каждого детектируемого нейтрино. В частности, благодаря этому было точно подтверждено, что детектируемые нейтрино прилетают из Солнца. Недостатком же реакции упругого рассеяния является предел более 6МэВ. То есть, они чувствительны только к нейтрино от распада 8 B , число которых относительно мало. Полученный результат R 0.5. Его можно объяснить следующим образом: поскольку в Камиоканде регистрировались нейтрино только очень больших энергий, то исходя из двух предыдущих результатов логично предположить, что значение R должно быть куда меньше, чем 0.3. Однако, водные эксперименты построены на реакции упругого рассеяния, а значит какая-то часть задетектированных нейтрино мюонные и таонные. Соответственно, их уч?т увеличивает реальное значение R. На СуперКамиоканде также проводились эксперименты, направленные на обнаружение осцилляций у атмосферных нейтрино. Они были основаны на предположении, что плотность мюонных нейтрино, прошедших сквозь Землю и прилетевших в детектор снизу (upward-going), будет отличаться от плотности мюонных нейтрино, прилетевших сверху (down-going). Этот эксперимент прив?л к подтверждению теории осцилляций. Тем не менее, данную методику нельзя было применить к солнечным нейтрино, поскольку диаметр Земли крайне мал по сравнению с расстоянием до Солнца. Поэтому СуперКамиоканде не дал окончательного решения проблемы солнечных нейтрино. Решила эту проблему канадская установка Sadbury Neutrino Observatory (SNO). Во многом это удалось благодаря уникальной среде детектора, чистой тяж?лой воде. Помимо упругого рассеяния на электронах, нейтрино взаимодействуют с ней двумя путя4


ми, недоступными в обычной. Первый из них зарядово-токовая реакция поглощения электронного нейтрино нейтронами дейтронов: e + 2 H p + p + e- . Порог реакции 2.23МэВ. Рождающийся электрон приобретает большую скорость и при движении созда?т детектируемый черенковский свет. Вторая новая реакция нейтрально-токовый разрыв дейтрона: x + d n + p + x . Нейтрино любого аромата освобождает нейтрон от протона. При этом не рождается никаких новых частиц, и реакция происходит с равной вероятностью для всех типов нейтрино. Е? порог такой же, как для реакции упругого рассеяния нейтрино на электронах, примерно 6МэВ. Освобожд?нный нейтрон не созда?т черенковского света, поскольку не нес?т заряда. Однако при рассеянии на ядрах он захватывается другим дейтроном, превращающимся при этом в ядро трития и испускающим фотон высокой энергии. Тот в свою очередь рассеивает электрон, который и начинает испускать черенковское излучение. Можно определить коэффициенты

RE = R

CC

/R

ES

,

RN = R

CC

/R

NC

.

Здесь RC C (RE S , RN C ) отношения между измеренным числом зарядово-токовых событий (упругих рассеяний, нейтрально-токовых событий) и их предсказанным числом в предположении отсутствия осцилляций. В результате экспериментов, провед?нных на установке SNO в 2001-м и последующих годах, когда к тяж?лой воде была добавлена соль, получились значения RE = 1, RN = 1. Поскольку SNO так же, как SuperKamiokande, состоит из большого бака, окруж?нного фотоумножителями, в ней информация так же снимается в режиме реального времени. Соответственно, принципиально возможно определить направление прил?та мюонных и таонных нейтрино. Этот факт позволяет рассматривать полученные значения как явное свидетельство существования осцилляций солнечных нейтрино и, следовательно, как подтверждение правильности стандартной солнечной модели Бакалла[7].

1.2.2 Атмосферные нейтрино
Нейтрино рождаются также и при столкновениях космических лучей с ядрами в верхних слоях земной атмосферы. Непосредственно при этих столкновениях рождаются пионы и каоны. Нейтральные пионы распадаются на два фотона ( 0 2 ), а заряженные на пару мюон плюс мюонное нейтрино: + ч+ + ч (и - . . .) Мюоны в дальнейшем претерпевают распад на электрон и два нейтрино: ч+ e+ + e + ч (и ч- . . .). Поэтому можно ожидать, что отношение количества мюонных нейтрино к электронным будет равно двум. Безусловно, это не вполне так, поскольку нейтрино также дают и каоны, а время жизни высокоэнергетичных мюонов настолько большое, что они могут не успеть распасться. Тем не менее, оценка достаточно точна. Существует два типа детекторов атмосферных нейтрино. Во-первых, это водные детекторы те же Kamiokande и SuperKamiokande, а также работавший в течение восьмидесятых годов IMB в США. Во-вторых, это железные калориметры установки 5


NUSEX, FREJUS и SOUDAN2. Все эти детекторы измеряют величину
ч e Nobs /Nobs R= ч . e NM C /NM C

Значения, полученные в NUSEX и FREJUS, примерно равны единице, но меньше е?, а все остальные эксперименты дали значение R 0.6. Это свидетельствует о наличии осцилляций ч . Эксперименты K2K и KamLAND, проводившиеся с земными источниками нейтрино, также подтвердили наличие осцилляций нейтрино, так что сейчас сам факт их существования является непредложным. Тем не менее, вероятность осцилляций зависит от характеристик вещества, через которые движутся нейтрино (эффект МихееваСмирнова-Вольфенштейна, MSW-eect), а численные характеристики осцилляций (m и ), полученные SuperKamiokande и IMB, противоречат друг другу. Новые эксперименты, например ANTARES, в том числе нацелены на решение этой проблемы. Таким образом, формализм теории вс? ещ? не довед?н до абсолютного конца.

2

1.3 Нейтринная астрофизика
Кратко рассмотрим третью группу источников нейтрино. Отдельным особняком стоят ГЗК-нейтрино. Существует так называемый эффект Грейзена-Зацепина-Кузьмина (GZK-eect), состоящий в том, что протоны энергий выше EGZK = 4 ћ 1019 эВ взаимодействуют с реликтовым излучением, рождают мезоны и теряют при этом энергию. Длину затухания для протонов при этом можно оценить в 50 мегапарсек. На данный момент не обнаружено и не предсказано никаких источников столь энергетичных космических лучей в пределах этого расстояния. Из этого следует, что космических лучей с энергией больше EGZK не должно существовать. Это утверждение подтверждено экспериментом HiRes Fly's Eye[4], в котором измеряется флуоресцентное излучение от атмосферных ливней. Другой же эксперимент, японский AGASA, представляющий собой множество наземных сцинцилляторов, показал полное отсутствие каких-либо нам?ков на существование ГЗК-эффекта. Этой ситуации возможны два объяснения[6]: либо зарегистрированные AGASA события на самом деле соответствуют частицам более низких энергий, а используемая теоретическая модель неверна, либо в HiRes также были получены соответствующие результаты, но отброшены как шумы. Ясность в вопрос существования ГЗК-эффекта может внести ввод в строй таких установок, как The Telescope Array[9] в американском штате Юта, Pierre Auger Observatory в Аргентине и CODALEMA во Франции. Все они сейчас находятся в стадии конструирования, и их проекты объединяет использование методов детектирования космических лучей как из AGASA, так и из эксперимента HiRes. В случае, если ГЗК-эффект вс? же существует, следует ожидать наличие нейтрино, рождающихся при распаде + , образующихся при взаимодействии протонов с реликтовым излучением. Другая категория источников космических нейтрино включает пульсары, молодые остатки сверхновых и ядра активных галактик[13]. Общее у этих объектов в том, что для них характерны механизмы ускорения ионов, которые затем тормозятся в плотной 6


среде и рождают нейтральные и заряженные пионы. При распаде пионов образуется примерно равное число нейтрино и -квантов. Соответственно, поскольку от этих источников детектируется большое количество фотонов, ожидается также получать от них и такое же количество нейтрино. Существует теория появления нейтрино из распада так называемых топологических дефектов, частиц с массами более 1020 эВ, оставшихся с момента рождения вселенной. Существуют теоретические предсказания частицы под названием нейтралино, подверженной только гравитационному взаимодействию, и потому концентрирующейся в центре Земли и Солнца, а затем распадающейся с испусканием нейтрино. Общей чертой всех источников, выделенных в третью группу, является полное отсутствие даже примерной информации об ожидаемых потоках нейтрино и их энергий. Современные эксперименты до сих пор не обнаружили никаких космических нейтрино, и результатом их работы является только установка новых ограничивающих пределов на их потоки. Существует также теоретическая модель Ваксмана и Бакала (W&B, [16]), которая ставит на порядок более ж?сткий предел, чем лучшие экспериментальные.

Рис. 2: График зависимости потоков нейтрино от их энергии. На координатной плоскости отложены максимальные пределы, заданные современными экспериментами и моделью W&B, а также известные потоки атмосферных и галактических нейтрино. Последние это те нейтрино, которые испускаются небесными телами нашей галактики. Соответственно, космогенные нейтрино суть все остальные.

7


1.4 Результаты 1.5 Методы детектирования
Все методы детектирования нейтрино построены на получении информации о свойствах рождаемых ими при взаимодействии с веществом установок лептонов и адронов. При вышеописанных зарядово-токовых взаимодействиях электронных нейтрино их энергия переда?тся электрону и адрону. В лептон уходит большая е? часть, примерно 80%. Электрон при движении испускает вторичные -кванты посредством тормозного излучения (bremsstrahlung). Испущенные фотоны в свою очередь могут превращаться в e+ /e- и ч+ /ч- -пары (pair production). Далее новые лептоны повторяют последовательность, и вскоре появляется огромное количество частиц меньшей энергии. Вся их совокупность называется электромагнитным ливнем. Энергия, переданная адрону, вызывает поток других адронов, называемый адронным ливнем. Рождающиеся в ливнях частицы теряют энергию, продвигаясь в среде. Около 80% процентов энергии от нейтрино E = 1020 эВ в воде выделяется в цилиндре радиусом 20 сантиметров и длиной 10 метров. Остальные двадцать процентов уносятся -квантами низких энергий. В масштабе огромных детекторов вроде IceCube такие события являются по сути точечными, и выделение из них угловой информации крайне затруднительно. Зарядово-токовые взаимодействия ч также рождают адронные ливни. Мюон же создаст множество отдельных электромагнитных ливней и при этом будет распространяться в среде с уменьшающейся скоростью. Во взаимодействиях типа ч + d ч- + u средняя энергия ч- составляет половину энергии нейтрино, во взаимодействиях ч + u ч+ + d средняя энергия ч+ 3/4 энергии нейтрино. Таким образом может быть оценена энергия исходных нейтрино. Собственно энергия мюонов определяется: при Eч < 100ГэВ из длины их пути, а при Eч > 1ТэВ из потерь ими энергии на единицу длины. Все современные создаваемые установки, а именно ANTARES, NEMO, NESTOR и IceCube ориентированы именно на детектирование этих зарядово-токовых взаимодействий ч с помощью Черенковского излучения. Таонные нейтрино при таких взаимодействиях создают -лептоны, которые почти сразу же распадаются на мюоны и соответствующие нейтрино. Поэтому такие взаимодействия нельзя отличить от мюонных. Ими пренебрегают из-за малого числа таонных нейтрино в общем потоке. Парадигма детектирования мюонных треков состоит в помещении установки глубоко под поверхность земли. Это позволяет элиминировать вклад от атмосферных мюонов и нейтрино, а также от космических лучей, которые иначе могли бы влететь в детектор сверху. Подразумевается, что нейтрино пролетают сквозь землю и подвергаются зарядовотоковому взаимодействию с веществом уже у самого детектора. Рождается мюон, который и детектируется установкой. Эксперименты по обнаружению осцилляций детектировали атмосферные нейтрино, прилетающие как сверху, так и снизу. В нейтринных телескопах нисходящие, down-going события не учитываются. 8


Сечение взаимодействия нейтрино с веществом увеличивается пропорционально квадрату его энергии. Нейтрино энергий порядка 1020 эВ уже не проходят через Землю, и зафиксировать их строящимися установками нельзя. Основные методы детектирования мюонов и адронных ливней следующие.

1.5.1 Черенков
Все заряженные частицы конусообразно излучают слабый голубой свет, если движутся в веществе со скоростью, превышающей скорость света c/n в этом веществе. Этот свет называется Черенковским излучением. Угол раствора конуса связан со скоростью частицы формулой 1 cos = . n Например, для воды n 1.33, и если считать, что скорость 1, то получится значение угла около 42-х градусов. Количество фотонов, испускаемых частицей единичного заряда вдоль пути dx в области длин волн d определяется выражением

d2 N 2 = 2 sin2 d dx
Из него можно получить, что частица с |q | = 1 и 1 на сантиметр пути выделяет примерно 200 оптических фотонов. Простые характеристики Черенковского излучения позволяют регистрировать его массивом фотоумножителей и затем легко восстанавливать трек частицы и е? энергию.

1.5.2 Радио
Избыток электронов в электромагнитном каскаде, вызываемый аннигиляцией позитронов и включением в ливень комптоновских электронов приводит к когерентному излучению в радиодиапазоне. Этот метод детектирования космических лучей был предложен в 1961 году советским уч?ным Г. Аскарьяном. Сейчас он используется в экспериментах, эксплуатирующих в качестве рабочего тела соляные массивы на Земле, а также Луну.

1.5.3 Акустический
Энергия от электромагнитного и адронного каскадов выделяется за время порядка 10-8 с. Можно считать, что ограниченный объ?м вещества подвергается моментальному нагреванию. Это вызывает неравномерное изменение его плотности, которое приводит к изменению давления и развитию акустических колебаний. Колебания начинают развиваться в радиальном направлении от оси каскада, и потому по их распространению можно восстановить трек и энергию нейтрино. Акустический метод имеет важное преимущество по сравнению с черенковским: длина затухания звука в плотных средах на порядок больше, чем длина затухания света (в океане для пиковых частот: километр против 70 метров). Поэтому отдельных датчиков акустических колебаний (они называются гидрофонами) требуется меньше, чем черенковских сч?тчиков. На практике 9


Таблица 1: Акустические параметры различных сред. Данные Океан Южный полюс Температура T (C) 15 -51 Скорость звука v (м/с) 1530 3920 Расширение , (10-5 K-1 ) 25.5 12.5 Тепло?мкость CP (J/kg/K) 3900 1720 Пиковая частота fp (kHz) 7.7 20 2 Постоянная Грюнайзена = v /Cp 0.153 1.12

взяты из [14]. Соль 30 4560 11.6 839 42 2.87

это означает значительно меньшую стоимость экспериментов. На данный момент это не играет особой роли, поскольку методики детектирования черенковского света лучше изучены, а размеры установок ( 1км3 ) позволяют укладываться в разумные бюджеты. Активное изучение акустического метода производится коллаборацией SAUND на военном гидроакустическом комплексе на Гавайях, проектом AMADEUS при ANTARES'е, в итальянском NEMO, а также в IceCube, где рядом с тремя струнами черенковских сч?тчиков опущены струны с гидрофонами. Следующий этап построения нейтринных телескопов будет проходить при доминировании акустических методов.

2 Обзор установок
2.1 Современные эксперименты
Все нейтринные телескопы можно разделить на две эпохи. Это старые установки эффективных площадей порядка тысячи квадратных метров и новые проектируемые и конструируемые телескопы объ?мами в один кубический километр. Из ?старых? установок интересны прежде всего отечественные Baksan и Baikal. Первая представляет собой 3132 отдельных объ?ма с жидким сцинциллятором, наблюдаемых каждый отдельным фотоумножителем. Установкой было обнаружено более 700 мюонных нейтрино. Байкальский детектор[17] это большой водный черенковский детектор, расположенный на дне озера Байкал. На данный момент он работает в конфигурации NT200 и состоит из 192 фотоумножителей, ограничивающих площадь около 200 квадратных метров и размещ?нных на восьми вертикальных струнах. В 2005 году над этой конфигурацией была произведена надстройка под названием NT200+. На расстоянии 100 метров от центра NT200 были помещены дополнительно три струны с двенадцатью фотоумножителями каждая. В перспективе планируется создать детектор объ?мом в кубический километр (km3-scale Baikal Telescope), и три этих струны как раз являются его первой частью. Однако, проект почти не имеет финансирования и вряд ли будет реализован. Существуют теоретические наработки[15] по созданию в озере Байкал акустического нейтринного телескопа. Все они обречены на полный провал, поскольку вода в озере имеет среднюю температуру около 4 градусов и, следовательно, не расширяет10


ся при нагревании. Кроме того, половину года поверхность Байкала покрыта льдом, порождающим громкий треск при малейшем изменении температуры. Заслуживает упоминания также установка AMANDA, действовавшая на южном полюсе и ставшая помимо прочего полигоном для обкатки технологий, применяемых сейчас в IceCube. Все эти эксперименты скорее подтвердили возможность детектирования нейтрино при используемых технологиях, однако не зарегистрировали ни одну частицу от космических источников.

2.2 Перспективные проекты
2.2.1 IceCube
IceCube[5] черенковский детектор объ?мом в кубический километр, конструируемый американскими уч?ными на южном полюсе. Всего во льду планируется пробурить 80 шахт глубиной по два с половиной километра и опустить в них струны с 60-ю оптическими модулями на каждой, удал?нными друг от друга на расстояние в 17 метров. Таким образом, самые верхние оптические модули будут закреплены на глубине в 1400 метров. Струны расположены на расстоянии 125 метров друг от друга и образуют гексагональную структуру. Бурение производится с помощью горячей воды. На поверхности льда будет расположен комплекс IceTop. Он будет использоваться для отбрасывания всех нисходящих событий с энергией больше 1014 эВ а также для определения атмосферных ливней. IceCube будет способен идентифицировать мюонные треки от ч с E > 1011 эВ, электромагнитные каскады от e с Ee > 1013 эВ, а также таонные нейтрино энергий E > 1015 эВ. Постройку планируется завершить в 2010 году, и детектор сможет ловить нейтрино, прилетающие через верхнюю полусферу Земли. Как было указано выше, в три скважины помимо струн с оптическими модулями были также опущены струны с тестовыми гидрофонами. Изучение распространения звука во льду южного полюса позволит в дальнейшей перспективе построить ещ? более огромный детектор, построенный на акустическом принципе и содержащий IceCube в качестве своей центральной части.

2.2.2 NEMO
NEMO это итальянский проект по конструированию похожего на IceCube детектора, но в Средиземном море. В 1998 году было выбрано место для будущего строительства в 80 километрах к юго-востоку от Сицилии. Глубина моря там составляет три с половиной километра. В течение девяти лет на месте производятся измерения скорости течения, акустического шума и прочих параметров. На данный момент не существует финального проекта установки. Более того, не решено, использовать ли гидрофоны, или ограничиться черенковскими сч?тчиками. Предполагается, что в перспективе NEMO будет работать в партн?рстве с IceCube, поскольку видные им области небесной 11


сферы значительно перекрываются, а значит получаемые данные можно будет сверять и калибровать. На данный момент одновременно производятся два эксперимента. Nemo Phase-1/2 Pro ject направлен на подробное изучение места будущего расположения эксперимента и стресс-тестирование оборудования. Он будет заверш?н к концу 2007 года. Ocean Noise Detection Experiment запущен в январе 2005 года. Он представляет собой систему из нескольких гидрофонов, установленную на дне и замеряющую акустические параметры моря[12].

2.2.3 NESTOR, ANTARES, KM3NeT
Два других аналогичных NEMO проекта греческий NESTOR и французский ANTARES[11]. Первый планируется завершить в течение тр?х лет. Уже готов и установлен один модуль, с помощью которого набирается статистика по космическим ливням. Во французском проекте разрабатываются концепции замены отдельных оптических модулей гидрофонами, но при этом решительно неясны конкретные планы по вводу установки в строй. KM3NeT общеевропейский проект по разработке и постройке третьего аналога NEMO. Он поддерживается европейскими бюрократическими институтами и, судя по всему, рассматривается ими как единственный из всех четыр?х проектов, который будет действительно заверш?н. При этом в течение двух лет планируется завершить собственно разработку проекта, а затем ещ? пять лет ожидается потратить на его материализацию.

3 Работа
3.1 Концепция
Итак, поставленная задача состоит в том, чтобы уметь быстро рассчитывать создание и развитие акустического поля, создаваемого в воде нейтрино энергии E0 1020 эВ. Задача делится на два этапа: на получение начального распределения энергии и на последующий расч?т акустического поля на его основе. Второй этап не представляет собой принципиальных трудностей, поскольку сводится к решению уравнения гиперболического типа. Первый же этап требует более изощр?нного подхода, чем тот, который можно предложить в первую очередь, то есть прямое моделирование последовательного рождения частиц и потерь ими энергии. Оно невозможно по причине крайне большого числа рождающихся в каскадах частиц и, соответственно, нерациональных требуемых компьютерных ресурсов.

3.2 Гибридный метод
Предложенный[8] Леонидом Григорьевичем Деденко так называемый ?гибридный? метод решения первого этапа задачи состоит в разбиении энергетического диапазона на 12


три части и на применение в каждом из них своего метода. При энергиях от 0 до Ethr 1010 эВ (1-й диапазон) и от ELPM 1014 эВ до E0 (3-й диапазон) используется метод Монте-Карло, а в среднем (2-м) диапазоне, при энергиях от Ethr до ELPM метод транспортных уравнений. Этот метод имеет реализованные аналоги для атмосферных ливней в воздухе, но не реализован для морской среды. Из компонент адронного ливня нам интересны только нейтральные пионы, распадающиеся с испусканием -квантов. Можно положить, что нейтральные пионы стабильны при энергиях выше 1018 эВ, а остальные адроны при всех энергиях вообще. В результате почти вся энергия рожд?нных адронов перевед?тся в высокоэнергетические -кванты. Таким образом, во втором и третьем диапазонах энергии (соответственно через транспортные уравнения и с помощью метода Монте-Карло) отдельно моделируется распространение адронов и лептонов. Выбор значения Ethr обусловлен тем, что при более низких энергиях в электромагнитном ливне становится значительным число частиц, удал?нных от его оси. Кроме того, при соответствующих энергиях уже не играют роли отдельные распады пионов, так что можно считать, что ниже Ethr существуют только электромагнитные каскады. Величина ELPM характерный предел по энергии, выше которого начинает проявляться так называемый эффект Ландау-Померанчука-Мигдала (LPM-eect). Он заключается в резком уменьшении сечений процессов рождения пар и тормозного излучения. В адронную составляющую переда?тся примерно 20% энергии нейтрино, и можно считать, что эффект ЛПМ не будет играть особой роли в этом случае. Напротив, случай, когда большая часть энергии переда?тся электрону, достаточно сложен. Если энергия электрона будет больше 1019 эВ, тогда из-за эффекта ЛПМ длина каскада значительно возраст?т. Например, длину каскада от электрона энергии 1018 эВ можно оценить в 50 метров, электрон же энергии 1021 эВ даст каскад длиной уже порядка километра. Также возникнет эффект, связанный с тем, что длины свободного пробега для процессов рождения пар и тормозного излучения станут намного больше, чем длины отдельных субкаскадов, которые уже не подвержены эффекту ЛПМ. Поэтому моделирование развития отдельных каскадов возможно только с помощью метода Монте-Карло. Затем при увеличении числа вторичных частиц флуктуации придут к определ?нному балансу, поскольку энергия этих частиц уже будет ниже характерной энергии ЛПМ ELPM 1014 эВ. Поперечное распределение частиц станет значительным намного позднее, при энергиях ниже Ethr . Поэтому в этом втором диапазоне энергий можно использовать значительно более экономный по сравнению с методом Монте-Карло метод транспортных уравнений. Итак, как было указано выше, развитие электромагнитных каскадов в первом диапазоне энергии моделируется методом Монте-Карло. Исходными данными при этом служит информация, полученная в результате моделирований во втором диапазоне энергии. Эта информация представляет собой плотности распределения фотонов и электронов разной энергии в зависимости от глубины (поперечная координата отсутствует, по13


скольку мы приближ?нно считаем каскад при энергиях E Ethr сконцентрированным около оси): Le (E , z ) и L (E , z ). Каждая из этих двух функций называется ?функцией источника?. Принципиальным моментом является то, что на этом этапе собственно моделирование с помощью метода Монте-Карло выполняется не каждый раз, а только единожды. Оно состоит в расч?те пространственных распределений энергии Se (Ei , z , r) и S (Ei , z , r) от электронов и гамма-квантов отдельных энергий Ei Ethr . Для каждого Ei проводится множество моделирований, а затем полученные данные усредняются и сохраняются для последующего использования. После этого искомое распределение энергии Q(z , r) может быть легко оценено:
z Ethr

Q(z , r) =
0

d
0

dE L (E , ) S (E , z - , r) + Le (E , ) Se (E , z - , r)

Значения S -функций в этой формуле усредняются из сохран?нных S -данных.

3.2.1 Расч?т акустического поля
На основе получаемого распределения энергии Q(z , r) легко расчитать развитие в воде акустических волн. Неоднородное волновое уравнение для акустического давления p имеет следующий вид: 1 2p Q p - 2 2 = - c t Cp t В н?м c скорость звука, постоянная теплового расширения воды, а Cp тепло?мкость воды при заданном давлении. Плотность поглощ?нной энергии Q можно связать с темпом е? поглощения q : Q = q , и кроме того можно ввести постоянную Грюнайзена t

=

c2 . Cp

Тогда неоднородное волновое уравнение можно переписать в виде

p -

1 2p 1 2 (q ) =- 2 c2 t2 c t2

Начальные условия для уравнения имеют вид

p(x, y , z , t = 0) = 0
Энергия выделяется в среде за время c 7 ћ 10-8 с, оно определяется скоростью распространения в среде ультрарелятивистских частиц. При этом характерное время распространения звука в поперечном направлении от оси каскада определяется скоростью

14


звука и равно по порядку величины s sin 7 ћ 10-5 с. Поэтому можно использовать приближение, когда энергия в каскаде выделяется практически мгновенно, то есть перейти от неоднородного уравнения к однородному

p -
с новыми начальными условиями

1 2p =0 c2 t 2

p(x, y , z , t = 0) = q
Его решение можно определить по формуле Пуассона:

p(x, y , z , t = 0) =

4

R

q dS , R

где интеграл бер?тся по сферической поверхности R = ct с центром в точке (x,y ,z ).

4 Нынешняя ситуация
Целостный комплекс приложений разработан Л. Деденко на языке FORTRAN. Его расч?ты, однако, носят скорее оценочный характер, и в частности, моделирование методом Монте-Карло в первом диапазоне энергий требуется более точно произвести с использованием пакета Geant4[1]. Именно эта часть вышеизложенного гибридного метода была реализована в работе.

4.1 Программа
При выполнении работы было реализовано три консольных программы, предназначенных для моделирования методом Монте-Карло пространственных распределений энергии Se (Ei , z , r) и S (Ei , z , r).

4.1.1 Библиотека данных
Получение распределений Se (Ei , z , r) и S (Ei , z , r) производится с помощью программы, использующей созданный в ЦЕРН'е пакет моделирования прохождения частиц через вещество под названием Geant4. Geant4 представляет собой не готовый программный продукт, подобный пакетам Mathematica или Matlab, а распространяющуюся в исходных кодах библиотеку, предназначенную для написания собственных программ с е? использованием. Предыдущая версия Geant'а, GEANT3, имевшая распространение в середине 90-х годов, была написана на формульном языке FORTRAN. Это ограничивало практически до нуля возможности по какой-либо серь?зной модификации пакета под собственные нужды, а также очень сильно затрудняло его использование, поскольку предлагало крайне ж?сткий путь включения вычислений в собственные программы. 15


Geant4 появился в самом конце 90-х годов, и уже к 2003 году после исправления всех серь?зных ошибок стал полнофункциональным и над?жным инструментом. В частности, он использовался в самом ЦЕРНе при моделировании параметров входящего в строй огромного ускорителя LHC. Разработка четв?ртой версии пакета производилась с нуля на языке C++ в парадигме объектно-ориентированного программирования. Это позволило перейти от устаревшей концепции вызывающих друг друга функций, использующих сторонние данные, к концепции независимых классов объектов, заключающих в себе необходимую для выражаемых ими операций логику. При этом, если логика различных операций однотипна, то для не? существует так называемый ?абстрактный? класс, а классы, реализующие эти конкретные операции, проявляют по отношению к соответствующим абстрактным классам операцию ?наследования?. Другими словами, объекты этих классов являются также и объектами абстрактных классов. Например, объекты классов G4GammaConversion и G4eBremsstrahlung, моделирующих процессы рождения пар и тормозного излучения, являются объектами G4VProcess. Пакет смоделирован таким образом, что управляющие классы используют в своей работе указатели на объекты именно абстрактных классов. То есть, существует возможность написать собственный класс для моделирования неучт?нного в стандартном пакете Geant4 физического процесса и наравне с другими использовать его при моделировании. Тем не менее, несмотря на безупречную как с точки зрения логики построения, так и с точки зрения производительности, структуру пакета, работа с ним возможна только под операционной системой Linux. Дело в том, что компиляция программ, использующих Geant4, очень здорово завязана на разные параметры конкретной установки пакета, такие как параметры визуализации, расположение и структура файлов и прочее. Компонента, который бы учитывал эти параметры при компиляции программы из-под Microsoft Visual Studio, стандартной среды для разработки программ под Windows, не существует. Он отсутствует по той банальной причине, что в научных институтах, где профессионально занимаются Geant4, применяется Linux, и такая возможность не представляет необходимости. Под Windows возможно установить программу Cygwin, являющуюся по сути эмулятором Linux'а, и компилировать программы в ней. Тем не менее, этот подход бессмысленен, поскольку не позволяет пользоваться богатыми возможностями Visual Studio по отладке программ, и, кроме того, весьма неудобен. Таким образом, для написания программы под Geant4 нужно использовать операционную систему Linux. Принцип написания такой программы следующий. Во-первых, нужно создать собственные реализации абстрактных классов G4VUserDetectorConstruction, G4VUserPrimaryGeneratorAction и G4VUserPhysicsList. Они предназначены соответственно для определения геометрии системы (?детектора?), для задания начальных частиц, прохождение которых через систему будет моделироваться, а также для выбора тех процессов, которые будут использоваться при этом моделировании. Во-вторых, нужно создать объект класса G4RunManager и передать ему с помощью его специальных методов указатели на объекты тр?х реализованных классов. Затем достаточно вызвать его функцию BeamOn(), и начн?тся собственно моделирование. Когда 16


все частицы либо вылетят за пределы детектора, либо для стабильных потеряют всю кинетическую энергию, моделирование будет завершено, и функция верн?т управление. Извлекать информацию о протекании моделирования, либо влиять на него, например, уничтожая частицы, можно, реализовав дополнительно собственные классы, наследующие от G4UserSteppingAction, G4UserTrackingAction, G4UserEventAction и также передав указатели на их объекты объекту G4RunManager. У каждого из этих классов должны быть определены специальные функции. У первого она будет вызываться для каждой частицы при каждом шаге моделирования, у второго в начале и в конце моделирования для каждой частицы, а у третьего в начале и в конце работы функции BeamOn(). Также извлекать информацию о ходе моделирования можно реализовав класс G4VSensitiveDetector и прикрепив его объект к тому объ?му в составе детектора, который непосредственно интересен. Именно последний подход был реализован в программе. Модуль под названием w01 представляет собой собственно моделирующую программу, использующую библиотеку Geant4. Он был написан в Linux'е на языке C++ с использованием среды KDevelop. Два других модуля, combine.exe и makeimage.exe, были написаны в Windows, в среде Microsoft Visual Studio, на языке C# под платформу .net Framework 2.0. Первая из этих программ предназначена для усреднения отдельных распределений Si , создаваемых w01, а вторая для быстрого построения графических изображений этих данных. Использование возможностей C# позволило легко написать и отладить два этих модуля. Несмотря на то, что они были откомпилированы в среде Windows, их работа под Linux'ом также возможна благодаря существующей под эту операционную систему реализации платформы .net под названием Mono[2]. При этом не происходит какого-то падения производительности. Дело в том, что бинарные файлы этих программ не являются исполняемыми сами по себе и не привязаны к конкретной компьютерной системе. Их содержимое максимально абстрагированный код, доводящийся до машинного .net-платформой при запуске.

4.1.2 Модуль w01
Реализованная программа w01 помимо вышеописанной структуры стандартной Geant4программы включает в себя ещ? два структурных компонента, позволяющих автоматизировать е? запуск с помощью возможностей скриптовых файлов операционной системы Linux. Эти компоненты реализованы в классах с названиями Config и Energizer. Config представляет собой мощнейший модуль по работе с настройками. В рабочей директории программы лежит специальный текстовый файл под названием app.config, на каждой строке содержащий пару ключ-значение. Его типичный пример таков:

TUBUS_HALFLENGTH_METERS=4 TUBUS_RADIUS_CENTIMETERS=15 ENERGY_EV=1e10 ENERGY_MIN_EV=0
17


GUN_PARTICLE=gamma SESSION=1 BATCH_VIS=vis.mac QUIET=1 LOG_DATA=1 LOG_BINARY=0 LOG_STEPS_R=300 LOG_STEPS_Z=400 LOG_FOLDER=/mnt/fat32/works/geant/photono/logs RUN_COUNT=1
Кроме того, те же параметры можно задавать при запуске программы из командной строки. Например, вызов

w01 ENERGY_EV=1e8 RUN_COUNT=5000 GUN_PARTICLE=eпривед?т к повторению моделирования пять тысяч раз, прич?м использоваться будет не десятигэвный -квант, а стомэвный электрон. Чтение настроек производится в программе с помощью тр?х доступных функций: Config::GetString(), Config::GetInt() и Config::GetDouble(). Как правило, они используются следующим образом:

G4double halflen = Config::GetDouble("TUBUS_HALFLENGTH_METERS", DEFAULT_TUBUS_HALFLENGTH_METERS);
В этом примере, взятом из класса конструирования геометрии детектора, узна?тся значение длины водяного цилиндра, в котором происходит рождение электромагнитного ливня. Первый параметр функции название переменной, второй значение по-умолчанию. Оно будет возвращено функцией GetDouble(), если значение переменной не задано, или не является числом. Константа DEFAULT_TUBUS_HALFLENGTH_METERS ж?стко определена в заголовочном файле:

#define DEFAULT_TUBUS_HALFLENGTH_METERS 10.
Такая конструкция позволила добиться тр?х вещей:

ћ Программу можно запускать с разными параметрами без перекомпилирования кода. ћ Возможность задавать параметры через командную строку позволила использовать богатые возможности скриптовых файлов Linux для автоматизированного запуска программы. ћ Система со значениями по умолчанию позволяет не указывать ненужные параметры в файле app.config и тем самым помогает избегать ненужных ошибок.

18


Кроме того следует отметить, что также реализованы и Set-функции, позволяющие записывать настройки в конфигурационный файл. Весь модуль Config, таким образом, предоставляет огромные возможности по использованию настроек в программе. Модуль Energizer содержит весь код, необходимый для сохранения вовне программы извлеч?нных из моделирования данных. Как было указано выше, выделение энергии от нейтрино энергии 1020 эВ в морской воде происходит внутри цилиндра длиной 15 метров и радиусом 20 сантиметров. В программе в качестве объ?ма моделирования используется водяной цилиндр произвольно задаваемых размеров (z , r). При каждом шаге моделирования реализацией G4VSensitiveDetector вызывается функция

Energizer::AddEnergy(double energy, G4ThreeVector& vec)
Е? параметры полное выделение энергии частицей за шаг моделирования и е? координата. Внутри себя Energizer содержит массив чисел размерами [LOG_STEPS_Z, LOG_STEPS_R], являющийся квантованным представлением водяного цилиндра. Ячейка с индексами [i, p] отвечает части цилиндра в форме прямоугольного тороида с координатами

i i z; z + dz , LOG_STEPS_Z LOG_STEPS_Z
где

p p r; r + dr , LOG_STEPS_R LOG_STEPS_R dr = r/LOG_STEPS_R.

dz = z /LOG_STEPS_Z,

При каждом вызове AddEnergy() определяется ячейка массива, соответствующая той точке цилиндра, в которой произошло выделение энергии, и к е? значению добавляется количество этой выделившейся энергии. После завершения моделирования массив пробегается в цикле, и значения его ячеек делятся на объ?мы соответствующих колец. Далее массив сохраняется в текстовый файл, и числа в н?м представляют распределение плотности выделенной энергии в цилиндрических координатах. Задание параметров LOG_STEPS_Z и LOG_STEPS_R, таким образом, зада?т точность получаемого распределения. Типичное название текстового файла следующее:

PgammaE1e10R150Z8000S300_400T1179352011.plg
После первой буквы P ид?т название частицы, затем после E е? энергия. После R и Z радиус и длина цилиндра. Два числа после S, раздел?нные знаком подч?ркивания используемые размеры массива. Наконец, последнее число выражает время, когда было произведено моделирование. Если программа была запущена два раза в течение одной секунды, то перед точкой в имя файла добавляется восклицательный знак и последовательные числа. Таким образом, любые два моделирования гарантированно имеют разные имена файлов. Содержимое файлов раздел?нные пробелами числа. Для экономии места на диске существует также возможность записи чисел в бинарном формате. 19


4.1.3

combine.exe

Если w01 запускается с параметром RUN_COUNT, не равным единице, то производится соответствующее количество моделирований, и результаты каждого записываются в отдельный файл. Программа combine.exe предназначена для усреднения распределений энергии, сохран?нных во множество таких файлов. В качестве своего параметра она принимает строку, с которой сверяет имена файлов в рабочей директории до первого совпадения. Затем она выбирает все файлы с той же частью названия от буквы P до буквы T, производит усреднение, и записывает результат в том же формате в файл с именем вроде

AV#100PgammaE1e10R150Z8000S300_400T11793673634925320.plg
Число после AV# количество файлов, среди которых было произведено усреднение. Естественно, возможно производить усреднение и среди таких файлов. Тогда префикс AV# будет добавлен вторично.

4.1.4

makeimage.exe

Создание makeimage.exe имеет простую цель быстрой и л?гкой визуализации получаемых данных. Она принимает в качестве параметров имена файлов, читает их содержимое и сохраняет в виде графического файла в формате png. Командная строка Linux имеет интересную возможность, отличающую е? от таковой в Windows: при передаче параметров программе можно использовать подстановочные знаки вроде ? и *. Тогда параметр будет рассмотрен как шаблон, и программе на самом деле будут переданы все имена файлов, удовлетворяющие этому шаблону. Например, вызов

./makeimage.exe AV*E1e10*.plg
позволит создать графические файлы для всех усредн?нных данных для энергий 1010 эВ. Существует возможность выводить не только пространственное распределение энергии, но и обратно умножать числа в массиве на объ?мы колец детектора с тем, чтобы строить первоначальные накопленные данные.

4.2 Результаты
На данный момент программы полностью отлажены и с их помощью проверяются и изучаются разные особенности Geant4, которые могут повлиять на получаемый результат. В первую очередь это выбор набора физических процессов. Дело в том, что Geant4 не является ч?рным ящиком, который умеет делать вс? и идеальным образом. Как используемые в моделировании частицы, так и применяемые к ним процессы должны задаваться явным образом. При этом помимо обычных электромагнитных процессов доступны также и специальные наборы процессов для низких энергий, LowEnergy и Penelope. Они предназначены для использования в диапазоне энергий от 250эВ до 100ГэВ, то есть используемый нами первый диапазон покрывается с лихвой. 20


Рис. 3: Примерный вид графического представления получаемого распределения, создаваемый модулем makeimage.exe. Программа w01 включает в себя возможность использования пяти наборов процессов. Это, во-первых, стандартный набор QGSP, предназначенный для универсального моделирования адронных процессов и, само собой, включающий в себя все стандартные электромагнитные процессы. Другие четыре набора были реализованы самостоятельно. Второй набор включает только электроны, позитроны и фотоны. Для фотонов используются стандартные процессы фотоэффекта, комптоновского рассеяния и рождения пар. Для лептонов процессы ионизации, тормозного излучения, и для позитрона процесс его распада. Третий набор отличается от второго добавлением мюонов и соответствующих им процессов. Четв?ртый и пятый наборы похожи на второй, но стандартные процессы заменены на свои LowEnergy- и Penelope-аналоги. С помощью скриптового файла были автоматически получены распределения энергии для гамма-кванта с энергиями 105 , 106 , . . ., 1011 электронвольт для каждого из пяти наборов процессов, прич?м для каждого случая была накоплена статистика из пяти ты21


сяч событий. Обзор графических результатов и просмотр технического вывода модуля w01 позволили сделать следующие первые выводы:

ћ Данные Penelope-процессов сильно разнятся при энергиях более 109 эВ с данными других наборов. При этом характер получаемых результатов весьма странен, и судя по всему, эти процессы нельзя применять при таких энергиях несмотря на официальную документацию. ћ QGSP, а также наборы с обычными электромагнитными процессами дают слабо различающиеся результаты. Следует понять причину этих вс? же наблюдающихся различий и выбрать один из тр?х пакетов с уч?том такого фактора, как время работы программы с каждым из них. ћ LowEnergy-процессы дают более мощное выделение энергии у центра пучка и более узкий пучок в целом. Следует изучить параметры и настройки LowEnergyпроцессов. ћ Вс? вышесказанное верно для графиков ненормированных значений выделенной энергии в ячейках массива. Графики плотности выделения энергии в цилинлрических координатах практически неотличимы друг от друга. Все они представляют собой несколько пиковых точек при минимальном значении радиуса около нуля по оси z . Во всех остальных точках выделение на порядки меньше. Нужно понять, действительно ли это может быть так.

4.3 Перспективы и планы
В работе было описано общее состояние современной физики по детектированию нейтрино, дан обзор предложенного Л.Г. Деденко гибридного метода по расч?ту акустического поля, создаваемого в воде нейтрино больших энергий, а также подробно описано устройство созданной и применяемой программы на основе пакета Geant4. Дальнейшие планы состоят в реализации в качестве программ всех остальных частей гибридного метода, их сопряжении и последующем изучении первых полученных результатов. Одновременно нужно разработать стратегию анализа данных, получаемых с помощью w01, и выбрать наилучший способ использования программы.

Список литературы
[1] http://geant4.cern.ch/. [2] http://www.mono- project.com/Main_Page. [3] E. Kh. Akhmedov. Neutrino oscillations: theory and phenomenology, hep-ph/0610064v2. [4] D.R. Bergman. Observation of the gzk cuto using the hires detector. Nuclear Physics B Proceedings Supplements, 165:1926. 22


[5] IceCube Collaboration. First year performance of the icecube neutrino telescope, astroph/0604450v2. [6] N. W. Evans, F. Ferrer, and S. Sarkar. Clustering of ultrahigh energy cosmic rays and their sources. Phys. Rev. D, 67(10):103005, May 2003. [7] Yuval Grossman. Tasi 2002 lectures on neutrinos, hep-ph/0305245v1. [8] Dedenko et al. L.G. Hybrid scheme of simulation of electron-positron and electronhadron cascades in dense medium at ultra-high energies. In Nahnhauser [10], pages 4549. [9] Kai Martens. The telescope array and its low energy extension. Nuclear Physics B Proceedings Supplements, 165:3336. [10] Rolf Nahnhauser, editor. Proceedings of the International Workshop (ARENA 2005), 2005. [11] Lahmann et al. R. Integration of acoustic equipment into antares. In Nahnhauser [10], pages 227231. [12] Piera Sapienza. Status of the nemo pro ject, astro-ph/0611105v1. [13] Andrew L. Stachyra. A Search for Astrophysical Point Sources of Neutrinos with SuperKamiokande. PhD thesis, University of Washington, 2002. [14] Justin Vandenbroucke. ph/0611503v1. Summary: Acoustic detection of ehe neutrinos, astro-

[15] Aynutdinov et al. V.M. A device for detection of acoustic signals from super high energy neutrinos. In Nahnhauser [10], pages 117121. [16] Eli Waxman and John Bahcall. High energy neutrinos from astrophysical sources: An upper bound. Phys. Rev. D, 59(2):023002, Dec 1998. [17] R. Wischnewski. ph/0609743v1. The baikal neutrino experiment: From nt200 to nt200+, astro-

23