Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Христиансен, Г. Б

.pdf
Скачиваний:
19
Добавлен:
22.10.2023
Размер:
13.7 Mб
Скачать

Не менее важен также вопрос об эффективной области реги­ страции ливней с различным значением параметра структурной функции п при заданном значении N. С другой стороны, не менее важно также измерение спектра по N с помощью одной и той же установки как в области больших, так и малых х так, чтобы про­ межуточная область приходилась на середину исследуемого ин­ тервала. Что касается области малых N, где происходит увеличе­ ние х, то там такие измерения сделаны уже в нескольких работах (например МГУ {233], Чакалтая {41]).

В

области больших

N,

где х уменьшается, фактически

этому

требованию удовлетворяют

только

измерения

на

установке [234].

В работе [338] показатель

х изменяется от 2,3 до

1,8 при

Nm\Os.

Таким образом, подтверждается уменьшение х с помощью

уста­

новки,

охватывающей

интервал

по N в

переходной области

х = 1,8+2,3.

Большой интерес представляют собой данные о спектре ливней по числу мюонов. Эти данные на обычных комплексных установ­ ках могут быть получены при достаточно точном определении на­ правления оси ливня, а следовательно, расстояния от оси до детек­ тора, при достаточной точности в определении р^ (для чего необ­ ходима большая площадь детектора) и относительно малой роли флуктуации /ц(г ) за счет соответствующего выбора регистрируе­ мого диапазона расстояний.

На установке МГУ неопределенность в нахождении

за

счет

всех этих факторов была рассчитана методом Монте-Карло

и со­

ставляла величину ~25% . Другое важное обстоятельство

пра­

вильный учет эффективности регистрации различных AV При на­ земной системе управления круг эффективной регистрации данного JVp, определяется требованием 95%-ной вероятности регистрации даже самых малых Ne с самыми большими s, соответствующими

данному

N^.

спектров по Ny, были

 

 

 

 

 

 

 

 

Исследования

проведены в-работах [227,

235, 236]. В работе [227] исследовались

мюоны с

энергией

более

10 Бэв.

N^

определялось

по усредненной

функции

fn(r),

причем г

изменялось

в интервале

1 5 < г < 5 0

м. Площадь детектора

 

мюонов

составляла

40 м2. Результаты относятся

к

интервалу

Л7^ = 3 • 103+-

-М05 . В работе [236] исследовались

потоки

мюонов с энергией бо­

лее 5 Гэв на высоте гор с помощью

детектора

мюонов площадью

30 м2. Диапазон

изменения N^ = 5- Ю 3 +5 - 10 4 . Результаты

исследо­

вания

спектров

мюонов

представлены

на

рис. 69.

Эксперимен­

тальные

данные

работ

[227, 236] приведены

к одной и той же по­

роговой

энергии Е^^Б

Гэв в соответствии

с данными

об энергети­

ческом

спектре мюонов в ш. а. л. (которые

были приведены

в § 3).

Как видно,

наблюдается

хорошее

согласие

работ

[227] и [236].

Спектр ш. а. л. по N^

так же, как и

спектр по Ne,

показывает

изменение

формы, причем

Ах,* === 0,7. Существует

и другой

незави­

симый

способ получения

спектра по А/~ц. Если

пренебречь

флуктуа-

циями

в N^

при

фиксированном Ne,

 

то

спектр

по N^, может быть

170

получен

из спектра по Ne и соотношения Л7^ — kNf.

Весь

вопрос

заключается

в том, с какой точностью

определено а

в различных

интервалах изменения

Ne.

 

 

 

 

 

 

 

Для

£ ^ ^ 1 0 Гэв

(см. § 3)

а = 0,78

во всем

исследованном ин­

тервале

значений Ne

от 2-Ю4 до 2-Ю7 . Если

рассмотреть

экспери­

ментальные

данные

отдельно

в интервале

 

jVe = 2-104 -^4-105

и

Ne = 1064-107,

где значения х е

соответственно

равны

1,5±0,08

и

egp(>Nj(%/iohm2сек-стер

 

^

^

 

 

 

 

* *

'

 

• Тянь-Шань

 

 

 

 

 

 

 

 

-

Сидней

 

 

 

-11

 

 

1 1

-

 

 

 

 

 

4

5

'

6

7

N 8

fgN

Рис. 69. Спектр ливней

по числу мюонов Л^, полученный

путем

пересчета

к числу мюонов

с £ m i n ^ 5

Гэв

по данным .работ: А — (227]

— 1236]: . . . — [235]

2,1 ±0,07,

то

значения а

в этих интервалах будут соответственно

0,79 ±0,05

и

0,78 ±0,03.

Отсюда значения

= — = 1,9 ±

0,15

для Л^ = 3- Ю3 -М04

и хц,= 2,6±0,15

для

= 104 + 105.

 

Что

касается влияния

флуктуации

в Np. на проведенный

таким

образом

пересчет, то оно могло бы быть

существенным при зави-

 

 

VD

(NJ

 

 

 

 

 

 

 

симости

 

 

=—=—

от jVe. Однако (см. § 3)

такая

зависимость

практически

отсутствует 1 0 3 .

 

 

 

 

 

На рис. 69 приведены

также данные о спектре ш. а. л. по числу

мюонов

с пороговой

энергией > 1 Бэв

в интервале NVI,= 106 4-2-109

[235], также приведенные к £ ^ . =5 Гэв.

 

 

 

 

 

Установка Сиднейского университета

[235] состоит

из 34

пунк­

тов регистрации, охватывающих общую площадь около 40 км2. Каждый пункт представляет собой два жидкостных сцинтиллятора площадью 6 м2 каждый, расположенных под слоем грунта и

разнесенных на расстояние 50 м. Если через

каждый из

сцинтил-

1 0 3 Строго

говоря, флуктуации не влияют при условии

постоянства х^. В

нашем

случае

изменение

приведет к некоторому искажению спектра

в

сторону

уменьшения А х а .

171

ляторов проходит три или более релятивистских частицы, то соот­ ветствующие сведения об амплитуде импульса и времени прохож­

дения записываются на

магнитную

ленту. Сопоставление лент,

относящихся к разным

пунктам,

производится с точностью

50 нсек.

 

 

Для обработки отбираются случаи срабатывания трех или бо­ лее пунктов. Угловое распределение регистрируемых ливней полу­ чается достаточно широким, и на рис. 69 отобраны случаи ливней,

близких

к вертикальным

(f> = 0-4-33°).

Для

определения

числа

авторы

используют аппроксимацию

f^C)

по

Грейзену {237] и

экспериментально показывают, как показатель /ц(г)

— зависит

от угла

Зависимость

от угла f> учитывается

далее

при обра­

ботке экспериментальных данных. Данные (235] в совокупности свидетельствуют об уменьшении Хц при iV^>105 .

Таким образом, данные как спектра по Ne, так и спектра по свидетельствуют о сложном характере формы первичного энер­

гетического спектра космических лучей сверхвысоких энергий.

Глава 5

Математические методы анализа экспериментальных данных по ш. а. л. и основные результаты исследования ядерно-каскадного процесса

§ 1. ОСНОВНЫЕ ПРЕДПОЛОЖЕНИЯ и ПРИНЦИПЫ МАТЕМАТИЧЕСКОГО АНАЛИЗА ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ

Для получения определенных количественных выводов о пара­ метрах ядерно-каскадного процесса и первичного космического излучения сверхвысоких энергий необходимы строгие математи­ ческие методы расчета различных характеристик ш. а. л. Эти ме­ тоды необходимы не только для проверки применимости различ­ ных моделей элементарного акта ядерно-каскадного процесса для

описания

ш. а. л., но и для выяснения,

какие характеристики ш. а. л.

наиболее

чувствительны к вариациям тех

или иных

параметров

элементарного

акта

и параметров 0

и А)

первичного

излучения.

В этом смысле

роль

строгих расчетов

становится стимулирующей

для эксперимента.

До обнаружения ядерно-каскадного процесса роль стимули­ рующей теории принадлежала классической теории э.-ф. лавин. Количественные расхождения в высотном ходе и пространствен­ ном распределении реальных ш. а. л. с расчетами по электромаг­ нитной каскадной теории способствовали формированию гипотезы ядерно-каскадного процесса. Аналитические методы расчета э.-ф. лавин были развиты в известных работах Ландау, Гейзенберга и др. (продольное развитие лавин) и в работах Нишимуры и Каматы (поперечное развитие лавин) (см., напр. (238]).

Несмотря на то что применимость квантовой электродинамики проверена экспериментально только до энергии 20 Гэв, с точки зрения современных тео­ ретических представлений (о применимости ее вплоть до расстояний порядка гравитационного радиуса электрона) экстраполяция элементарных сечений про­ цессов образования пар, тормозного излучения и пр., в область предельно вы­ соких энергий космических лучей не вызывает сомнения. Правда, в области пре-

173

 

m e c 2 X 0

 

дельно высоких энергий электронов и фотонов

Ес я — 8 я г в — •

длина /-единицы; гв — боровский радиус для атома

водорода) начинает

играть

роль эффекта подавления процессов образования пар и тормозного излучения за

счет нарушения

процессов кулоновского

рассеяния 1 0 4 .

В верхних

слоях

атмо­

сферы

при

Х0

= 3-\05 см

Ec—\0,s

эв;

в

нижних слоях,

где

Х0С^.З-\04

см,

ЕС17

эв,

у-излучение

такой

высокой

энергии

если

и возникает, то

только

в актах взаимодействия

частиц

с энергией,

по-видимому,

более

10 1 9 - М0 2 0 эв.

Поэтому учет влияния этого процесса на развитие э.-ф. лавин необходим,

может

быть, только в области предельно высоких энергий.

 

 

 

 

После

установления существенной

роли

ядерно-каскадного

процесса возник вопрос о создании адекватного задаче математи­ ческого аппарата. Экспериментальные данные о ядерных взаимо­ действиях при высоких энергиях были в то время чрезвычайно скудны. Поэтому некоторые авторы [239] пошли по пути проведе­ ния аналогии между развитием ядерной и э.-ф. лавин, предполо­ жив, что сечение ядерного взаимодействия обладает тем же свой­

ством однородности, что и сечение

элементарных

взаимодействий

для фотонов и электронов 1 0 5 . В этом случае

уравнения

для

ядер­

ного каскада можно решать теми

же аналитическими

методами,

что и для э.-ф. каскада. Однако появление теории

Ферми и затем

теории Ландау поставило вопрос

о создании более

универсаль­

ного математического аппарата (в

котором

учитывается

множе­

ственность процесса и неоднородность сечения), а также о сравне­

нии

этих теорий

с экспериментальными данными по ш. а. л.

 

Эти

вопросы

были

решены

Розенталем [240],

обобщившим

метод

последовательных

поколений,

предложенный

еще в 30-х го­

дах

в

связи с попыткой

решения некоторых задач

электромагнит­

ной каскадной теории [241]. Предложенный в [240] метод заклю­ чается в расчете ядерной лавины я.-а. частиц с энергией £ и на глубине х как суммы лавин, состоящих из я.-а. частиц, рожден­ ных в результате столкновения первичной частицы, столкновений

вторичных частиц

(возникших

 

от первичной), столкновений

частиц,

рожденных от вторичных, и т. д. до со.

 

 

 

 

 

1 0 4

Это

нарушение происходит,

если

 

многократное кулоновское рассеяние на пу­

 

ти,

равном

эффективной

длине

процесса

~

hjqu (где qu — продольная

 

компонента

разности

импульсов,

первичной

и

вторичных частиц),

создает

 

 

 

 

 

 

 

тс2

 

 

 

 

 

 

 

 

 

отклонение

на угол f> ^

•-

 

(угол

рождения

частиц). Из этого

условия

 

можно

получить

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

Х 0

 

 

 

 

 

Г

 

 

 

 

- 2 _

Y ^

2

• h/mcf-^)2

тс1 /

' Е' =

^

У

'

 

Отсюда

 

 

 

 

 

\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

Х0

1

 

 

Х0

 

 

 

 

 

 

 

Y

^ 8п

h2/e2m

8 я ' г в

 

 

1 0 5

Сечение

есть

функция

отношения

Е'/Е,

где

Е

и

Е'

— энергии взаимодей­

 

ствующей и вторичной

частиц.

 

 

 

 

 

 

 

 

174

С

помощью этого

метода [241] были рассчитаны средние числа

я.-а. частиц, электронов

и мюонов в ливнях

с различными первич­

ными

энергиями Е0.

В

работе [242] метод

[240] был усовершен­

ствован и обобщен для большого числа сортов я.-а. частиц с уче­

том

возможной зависимости сечения

взаимодействия

от

энергии,

а также с учетом ионизационных потерь я.-а. частиц.

 

 

 

Вся эта методика в дальнейших

работах была

использована

для

проверки возможности описания

ядерно-каскадного

процесса

в рамках моделей Ферми, Ландау и Гейзенберга. При этом пред­ полагалось, что элементарные взаимодействия в ядерно-каскад­ ном процессе происходят в соответствии с той или иной моделью при любых энергиях я.-а. частиц вплоть до энергий, при которых ядерно-каскадный процесс затухает 1 0 6 . Однако ни одна из моделей не выдержала проверку на такую универсальность1 0 7 . В то же время существенно менялась экспериментальная ситуация. Изуче­ ние я.-а. частиц высоких энергий в составе ш. а. л., а затем и ис­ следование потоков мюонов продемонстрировало решающую роль

флуктуации

в ядерно-каскадном процессе.

С другой

стороны,

значительно

изменились и обогатились экспериментальные

данные

о ядерных взаимодействиях при высоких энергиях.

 

 

 

Исследование флуктуации вызвало к жизни

большое

число

различных упрощенных моделей развития

ш. а. л.

В то

же

время

точный расчет флуктуации в ядерно-каскадном процессе предпо­ лагал решение значительно более сложной задачи, чем интегродифференциальные уравнения для среднего числа я.-а. частиц с энергией Е на глубине х. Строго говоря, для решений этой задачи надо было обратиться к уравнению Колмогорова — Дмитриева для стохастических процессов. Для ядерной лавины с первичной энергией £о—Ю1 5 эв точный расчет распределений числа я.-а. ча­ стиц с энергией Е на глубине х даже на современных ЭВМ про­ вести, по-видимому, невозможно [243]. Можно определить лишь

среднее число и дисперсию

распределения.

 

 

 

 

Поэтому

оставалась

единственная

возможность — расчет

ш. а. л. методом

случайных

испытаний. Однако практическая его

реализация

на

современных ЭВМ

оказалась

всецело

связанной

с

теми

приближениями, которые

возможны

при экстраполяции

на

весь

диапазон энергий,

существенный

для

развития

ядерной

лавины, основных представлений, заимствованных из области высоких энергий.

Согласно современным данным, полученным как на ускорите­ лях, так и в космических лучах, взаимодействие нуклонов в диа­

пазоне энергий

10 1 0 - М0 1 3 эв характеризуется с точностью да

10%

1 0 6 Строго

говоря,

применение модели Ландау было ограничено значением

. Е ^ Ю 1 1

эв.

В

областях £ < 1 0 и

эв при этом использовали данные

модели

Ферми.

 

 

 

 

 

№7 Этот результат

не исключает возможной применимости, например, моделей

Ландау

или

Гейзенберга только

в области сверхвысоких энергий (например,

при Е0и

эв).

 

 

17S

постоянством

сечения

неупругого взаимодействия

(пробег относи­

тельно взаимодействия

в

воздухе

к~80 г/см2),

приблизительным

постоянством

 

среднего

значения

коэффициента

неупругости

взаимодействия

(&„ = 0,5)

и определенной

формой

распределе­

ния f(kn) т .

В

процессе

взаимодействия

нуклонов

рождаются

вторичные частицы пионы, каоны и нуклон-антинуклонные пары.

Доля каонов и нуклон-антинуклонных

пар

не превосходит 20%

от

всех вторичных

частиц 1 0 9 .

Полное число

вторичных

частиц

пя

в

рассматриваемом

диапазоне

энергий

слабо зависит от энергии

( ~ 1 п £ 0 ) , составляет <Г10

и

сильно

флуктуирует.

Частицы

рас­

пределены по падающему энергетическому спектру

ф„

(Е) так,

что их средняя энергия

существенно

возрастает

с Е0. Угловое

распределение вторичных частиц определяется эмпирически уста­

новленным в экспериментах

на

ускорителях

 

и в

космических

лучах распределением1 1 0

поперечных

импульсов

р^е~Р±/р'>dp±

при р 0 = 1,6- Ю8

эв/с.

 

 

 

 

 

 

 

 

 

Взаимодействие

вторичных

частиц

пионов

происходит

иначе.

Пробег

относительно взаимодействия

в

1,5

раза

больше —

120 г/см2.

Если

не

делать

различие

между пионом, вызывающим

взаимодействие, и вторичными пионами, то kn

— \

и энергетический

спектр

вторичных

частиц

получается

более

 

жесткий, чем при

взаимодействии

нуклонов.

Средние

энергии

 

вторичных

частиц

отличаются в два раза. Множественность вторичных частиц не­ сколько больше (в 1,5 раза), а ее зависимость от Е0, возможно, более сильная, чем при взаимодействии нуклонов, и не противо­ речит закону EQI*.

Если экстраполировать эти представления на область сверх­ высоких энергий, то характерная особенность ядерно-каскадного процесса, создаваемого первичным нуклоном, будет заключаться в том, что его флуктуации в основном будут связаны с конкрет­ ным поведением первичного нуклона. Это объясняется тем, что

энергия, которая сохраняется на нуклоне

среднем),

значитель­

но превосходит

энергию, приходящуюся

на

каждую из

вторичных

частиц,

так как

вторичных

частиц

рождается достаточно много.

Таким

образом,

первичный

нуклон

здесь играет роль

лидирую­

щей или «ведущей» ливень частицы, во всяком случае до тех пор, пока его энергия превосходит энергию других я.-а. частиц на дан­ ной глубине в атмосфере.

1 0 8

Вопрос

о слабом росте сечения неупругого взаимодействия

нуклонов в

диа­

 

пазоне

энергий

101 0 -н101 3

эв все еще дискутируется в литературе.

 

1 0 9

Согласно

новым

данным,

полученным

на ускорителях, доля

нуклон-антинук­

 

лонных

 

пар возрастает

в интервале

ЗОН-1000 Гэв

от

долей

процента

до

 

- 5 % .

 

 

 

 

 

 

 

 

 

 

1 1 0

Вопрос

о

справедливости

такого распределения поперечных импульсов при

 

энергиях более нескольких единиц на

101 2 эв остается открытым. Приведен­

 

ное распределение pj_ ~

справедливо

до pj_ ~ 2

Гэв/с.

При

больших

 

начинается более медленный спад функции распределения.

176

Идея лидирующего нуклона была предложена С. Н. Верновым и впервые рассчитана с учетом флуктуации в развитии ядерной лавины в работе [244]. Расчет проводился методом случайных ис­ пытаний в соответствии со схемой развития ливня (рис. 70). В этой схеме лидирующий нуклон испытывает i-тое взаимодейст­

вие на глубине xt

и передает

вторичным

частицам

долю

от той

энергии,

которую он

сохранил

после

(i—1)-го

взаимодействия.

Место

 

взаимодействия

разыгрывается

с

помощью

таблицы

слу­

чайных

чисел

в

соответствии

 

с

распределением

глубин

взаимо­

действия

по

закону

 

-xixn

_dx_

где

х-

глубина,

отсчитываемая

 

 

 

 

 

 

от

места

предыдущего

 

взаимодействия

Коэффициент

неупруго-

сти k разыгрывается в соответствии с

 

 

 

 

 

функцией

/«_(&),

 

принятой

в

 

[244]

в

виде'

 

 

 

 

 

f(k)=8(k—k).

 

В

каждой точке я,- энергия

 

 

 

 

 

kiEi-x

(Ei-i — энергия

нуклона в г-том

вза­

 

 

 

 

 

имодействии)

передается

вторичным

части­

 

 

 

 

 

цам. Лавины я.-а.

частиц,

электронов и

 

 

 

 

 

мюонов

от этих

 

вторичных

частиц

могут

 

 

 

 

 

быть рассчитаны в среднем, например, ме­

 

 

 

 

 

тодом

 

последовательных

поколений. Если

 

 

 

 

 

обозначить число частиц в такой

парциаль­

 

 

 

 

 

ной лавине через

N{(xXi),

то полное

чис­

 

 

 

 

 

ло

 

частиц

 

данного

 

сорта

 

дается

 

 

 

 

 

суммированием

 

по

всем

 

 

значениям

 

 

 

 

 

 

Таким образом, в этом подходе предпо­

 

 

 

 

 

лагается,

что случайная

реализация

функ­

 

 

 

 

 

ции N(x)

определяется

в основном

случай­

 

 

 

 

 

ной реализацией величин х\, ku

а

флуктуа­

 

 

 

 

 

ции в развитии лавин от вторичных частиц

 

 

 

 

 

не играют

существенной

роли.

 

 

Справедли­

 

 

 

 

 

вость этого утверждения для довольно ши­

 

 

 

 

 

рокого класса моделей с лидирующим нук­

 

 

 

 

 

лоном была математически доказана в про­

 

 

 

 

 

цессе

расчетов

флуктуации ядерного каска­

Рис.

70.

Ядерно-кас­

да

по

методу

 

Колмогорова

— Дмит­

риева

[245].

 

 

 

 

 

 

 

 

 

 

 

 

 

 

кадный

процесс

по

 

Одна из целей исследования ш. а. л. —

модели

лидирующего

 

нуклона

С. Н. Верно-

это получение информации о ядерных

взаи­

ва:

х — место

вза­

модействиях

при

 

сверхвысоких

 

энергиях

имодействия, k — ко­

101 4 -102 0

эв.

Естественно,

что

эксперимен­

эффициент

неупруго­

 

 

сти

 

 

тально

 

наблюдаются

лавины

 

вторичных

 

 

 

 

 

 

 

 

 

 

 

частиц

с

энергией,

во

всяком

 

случае

не

превышающей,

как мы

видели,

101 2 4-101 3

эв. Чтобы

на

основании

их изучения можно

12 г. Б. Христиансен

177

было получить информацию о ядерно-каскадном процессе при бо­ лее высоких энергиях, необходимо иметь тщательные и по возмож­

ности

однозначные

данные

о

взаимодействиях

при

энергиях

< Ю 1 3 эв. Эти данные должны быть заложены в расчет в

качестве

фундамента для

любого

реального

анализа

эксперимента по

ш. а. л.

Наиболее

очевидным

первым

шагом

в таком

анализе

должна

быть экстраполяция

этих данных на область сверхвысоких

ипредельно-высоких энергий.

Вработах [244, 246] описанным выше методом были рассчи­ таны распределения чисел электронов на различных глубинах в

атмосфере

при

различных

значениях

первичной

энергии

Е0= 101 3 Ч-

101 5 эв.

В последующих работах

[247, 227,

245], были

проведены расчеты флуктуации других экспериментально наблю­

даемых

параметров ш. а. л. (числа

я.-а.

частиц,

числа мюонов,

потоков

энергии Фэ.-а. и парметра

s),

а также

определена кор­

реляция между этими параметрами. В многочисленных экспери­ ментальных работах, как мы видели выше, анализ данных прово­

дился для ливней с фиксированными Ne

или

Np. Поэтому

с

точки

зрения

эксперимента

интерес представляли распределения

раз­

личных параметров не при фиксированном Е0,

а, например,

при

фиксированном Ne.

s, Фя.-а.)

 

 

 

 

 

 

 

 

 

Если WEa{Ne,

Ny,,

 

многомерное

распределение

на

некоторой

глубине

в

атмосфере

по параметрам Ne,

Np,

s, Фя .-а.

при

фиксированной

первичной энергии

Е0,

то

интересующее

нас

распределение при фиксированном Ne можно

получить,

если

из­

вестен первичный

энергетический

спектр F

(E0)dE0.

 

 

 

час­

Действительно, число

ливней,

происходящих

от первичных

тиц

с энергией Е0,

 

Е0

+ dE0

и

имеющих

параметры

Л^,

 

Ne

и т.

д.,

есть

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

AE0-<y+»dE0WEt(Ne>

 

A V

ф я - а . ,

s)dNlidOa..a.dsdNe.

 

 

 

Тогда распределение,

например, по

параметру

s

при

фиксированном

Ne дается

выражением

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Щ И ^ М » ,

s, Фя.-а.,

 

 

N^dO^.dNpAE-iv+VdEods.

 

 

В работе [245а] даны простые и практически полезные форму­

лы,

позволяющие

вычислять

любые распределения

параметров

ш. а. л. при фиксировании

любых

других,

если

Wf0

представляет­

ся в

виде

многомерной

гауссовой

функции,

зависящей как

от

дисперсий

параметров,

так и

от

коэффициентов их

корреляции.

Несмотря на всю мощь используемых в настоящее время экс­

периментальных установок, сведения о таких важных

компонен­

тах ш. а. л., как я.-а. и мюонная, экспериментаторы

получают

обычно на определенном расстоянии от оси индивидуального ре­ гистрируемого ливня. Целый ряд интересных явлений, наблюдае­

мых в

ш. а. л. (многоствольные ливни,

пучки

мюонов и пр.),

нельзя

анализировать, ограничиваясь

решением

одномерной за-

178

дачи. Наконец, само по себе пространственное распределение лив­

невых частиц

интересно

с

точки зрения

анализа распределения

поперечных

импульсов

в

ядерных

взаимодействиях.

По этим

причинам уже давно назревала необходимость расчета

трехмер­

ной картины

развития ш. а. л. С другой стороны, развитие техники

ЭВМ создавало для этого

благоприятную

возможность.

Расчет

был начат независимо в работах [248,

249].

 

 

В работе [249] был предложен новый подход к расчету трех­

мерной картины, позволяющий существенно сэкономить

машин­

ное время и решать задачу для любых отношений Е0

0

энергия первичной частицы, Е — энергия вторичных частиц). Весь расчет проводится методом случайных испытаний. Сущность под­

хода заключается в том, что на первом этапе проводится

расчет

ядерной лавины частиц высоких энергий,

например

 

лавины

частиц с энергией Е, для которых

Е/Ео^Ю-^.

 

 

 

 

 

Для вычисления лавины частиц любых малых энергий е, про­

исходящих

от первичной

частицы Е0,

проводится

расчет

лавин

частиц

малых

энергий

е,

создаваемых

 

частицами

с

энер­

гией

Е^10~4

EQ

С точками зарождения

на

различных

глубинах

х0

в

атмосфере. Расчет

лавин

частиц с

энергией

Е и

более от

первичной

частицы с энергией

Е0

проводится

со статистикой, тре­

буемой

задачей. Расчет лавин

частиц

малых

энергий е от

частиц

с

энергией

Е проводится

с

ограниченной

статистикой — значи­

тельно меньшей, чем та, которая получилась бы, если бы каждая лавина от Е0 разыгрывалась до энергии е. Окончательный резуль­ тат получается объединением розыгрыша до энергий Е со слу­ чайной выборкой из розыгрышей лавин малых энергий е, зарож­ дающихся на разных глубинах Хо. Этот метод целесообразно назвать методом случайных испытаний с ограниченной выборкой. Точность его определяется выбором интервалов АЕ и Ах, для ко­

торых составляется ограниченная

выборка.

Как

показано

в ре­

зультате

практических расчетов,

точность

метода может

быть

сделана достаточно высокой, не хуже 10% [250].

 

 

 

Проводя такие расчеты для относительно малых Е0 при раз­

личных

точках

зарождения лавины х0, можно

затем

перейти к

еще большим

Е0 и дойти таким

образом до

предельно

высоких

первичных энергий. Именно такой подход был использован в ра­ боте [251] для расчета лавин мюонов в случае предельно высоких энергий 10 1 7 - М0 2 0 эв.

§ 2. МЕТОДЫ РАСЧЕТА СРЕДНИХ ХАРАКТЕРИСТИК ЯДЕРНО-АКТИВНОЙ, ЭЛЕКТРОННО-ФОТОННОЙ И МЮОННОЙ ЛАВИН

Рассмотрим простейший вариант расчета методом последова­ тельных поколений. Чтобы подчеркнуть физический смысл мате­ матических операций, введем некоторые упрощающие предполо­ жения:

12*

179

Соседние файлы в папке книги из ГПНТБ