Санкт-Петербургский Государственный Университет
Физический факультет
Кафедра ядерно-физических методов исследования
Разработка дрейфовой камеры для измерения скорости дрейфа электронов в газовых
смесях детекторов частиц
Магистерская работа студента дневного отделения
_________________
Валияровой Адели Ринатовны
Научный руководитель:
_________________
к. ф.-м. н., Ильин Д. С.
Рецензент:
_________________
Санкт-Петербург
2016
к. ф.-м.н., с. н. с., Гаврилов Г. Е.
Оглавление
Введение......................................................................................................................................... 4
1
Большие экспериментальные установки CERN с ГРД....................................................... 7
1.1
Детекторы и газовые смеси....................................................................................... 7
1. 2
Взаимодействие заряженных частиц в газовых средах.........................................11
1. 2.1 Удельные потери энергии............................................................................... 11
1. 2.2 Ионизационные потери тяжелой заряженной частицы в среде..................11
1. 2.3 Ионизационные потери для электронов........................................................12
1. 2.4 Радиационная длина........................................................................................13
1. 2.5 Образование ионных пар................................................................................ 14
1. 3
Основные характеристики газовых смесей............................................................15
1. 3.1 Коэффициент Таунсенда.................................................................................15
1. 3.2 Скорость дрейфа..............................................................................................16
1. 3.3 Диффузия......................................................................................................... 16
1. 3.4 Диффузия ионов и электронов в нулевом магнитном поле.........................16
1. 3.5 Подвижность ионов.........................................................................................19
1. 3.6 Угол Лоренца................................................................................................... 20
1. 4
Принцип работы газоразрядного детектора с усилением – пропорционального
счетчика 20
1. 5
Метод исследования скорости дрейфа электронов................................................24
1. 5.1 Дрейфовые камеры..........................................................................................24
1. 5.2 Примеры установок измерения скорости дрейфа из литературы...............26
2
2
Разработка дрейфовой камеры............................................................................................ 30
2.1
Описание методов моделирования..........................................................................30
2.2
Моделирование детектора в программе Garfield................................................... 30
2.3
Конструкция дрейфовой камеры............................................................................. 31
2. 4
Выбор и оптимизация геометрических параметров детектора...........................32
2. 5
Вычисление скорости дрейфа..................................................................................36
2. 6
Оценка вклада факторов, влияющих на точность определения скорости дрейфа40
2. 7
Оценка толщины сцинтилляторов...........................................................................45
Заключение...................................................................................................................................47
Список литературы......................................................................................................................48
3
Введение
Актуальность. Газоразрядные детекторы ионизирующих излучений применяются в
различных областях современной физики. Они просты по конструкции, компактны, удобны в
работе, но во время работы необходимо следить за такими процессами, как взаимодействие газа
с материалами, очистка газа, загрязнение газовой смеси и т.д. [1-7].
Большая часть газовых детекторов, используемых в физике высоких энергий, работает на
основе смесей, содержащих парниковые газы, которые обеспечивают необходимо временное
разрешение и газовое усиление. Для оценки вклада этих газов в атмосферный парниковый
эффект используется потенциал глобального потепления (GWP) – параметром, численно
определяющим радиационное (разогревающее) воздействие молекулы определенного
парникового газа относительно молекулы СО 2. К примеру, GEM-детекторы [8] (газовые
электронные умножители), которые используются в экспериментальной установке CMS
(Компактный мюонный соленоид) Большого Адронного Коллайдера (БАК), работают на основе
Ar/CO2 газовой смеси [9]. Однако для достижения высокого временного разрешения к данной
смеси добавляется CF4 [10]. Нужно отметить, что потенциал глобального потепления CF4
составляет 7390 [11]. Камеры с резистивными пластинами (Resistive Plate Chamber) [12]
работают на основе R134a/C4H10/SF6 газовых смесей, средний GWP которых равен 1430. В
настоящее время возникла необходимость исследовать новые газовые смеси, характеристики
которых соответствовали бы предъявляемым требованиям. Основные газовые смеси, активно
используемые в физики высоких энергий, и их транспортные свойства при меньшем
содержании фреонов представлены в работах [13,14].
Сегодня проблема выбросов парниковых газов в атмосферу стала актуальной для всей
современной индустрии. Поэтому при разработке газоразрядных детекторов необходимо
принимать во внимание экологическое воздействие используемых газовых смесей на атмосферу.
К тому же было показано, что наиболее часто используемые в детекторах газы R134a, CF4 и SF 6
вносят наибольший вклад в парниковый эффект. В настоящее время основным решением
данной проблемы является поиск новых экологически безопасных газовых смесей.
Целью данной работы является разработка и оптимизация параметров дрейфовой
камеры для исследования характеристик новых газовых смесей в детекторах частиц.
Основные положения, выносимые на защиту:
4
1. Для создания модуля дрейфовой камеры для исследования скорости дрейфа
электронов в газовых смесях при Е=500÷5000 В/см разработана модель детектора в
программном пакете GARFIELD.
2. Проведена оптимизация электрического поля, и определены диапазоны рабочих
положений источника 90Sr, в которых обеспечивается требуемая однородность
электрического поля.
3. Построены временные зависимости для газовой смеси 60%Ar + 30%CO2 + 10%CF4 во
всем рабочем диапазоне напряженностей электрического поля.
4.
Проведена успешная проверка алгоритма определения скорости дрейфа для
указанной газовой смеси, и определены вклады основных ошибок.
Научная новизна. Особенностью установки является возможность перемещения
источника 90Sr вдоль модуля с малым шагом 200 мкм ÷ 1 мм, т.е. получения ряда
экспериментальных точек и соответственно более точной аппроксимации данных, по сравнению
с аналогами. Кроме этого, компактность детектора позволяет поместить его в соленоид
небольшого размера для измерения скорости дрейфа от напряженности магнитного поля.
Научная и практическая значимость. Разработанный дрейфовый детектор позволяет
исследовать характеристики новых рабочих газовых смесей, в частности измерения скорости
дрейфа электронов при напряженностях полей Е=500÷5000 В/см с высокой точностью.
Апробация работы. Основные результаты работы были представлены на Школе ПИЯФ
по физике конденсированного состояния (ФКС-2016).
Публикации. По результатам конференции опубликованы:
1. Валиярова, А. Р. Разработка дрейфовой камеры для измерения скорости дрейфа
электронов в газовых смесях детекторов частиц / А.Р. Валиярова, Г.Е. Гаврилов, Д. С.
Ильин // Сборник тезисов и список участников: 5 0 - я Школа ПИЯФ по физике
конденсированного состояния (ФКС-2016), 14-19 марта 2016 г., Санкт-Петербург / ФГБУ
«ПИЯФ» НИЦ «Курчатовский институт». – Гатчина, 2016. – 216 с.
2. Валиярова А. Р, Гаврилов Г.Е., Ильин Д. С. Модель дрейфовой камеры для измерения
скорости дрейфа электронов в газовых смесях детекторов частиц / Валиярова А.Р.,
Гаврилов Г.Е., Ильин Д. С. – Гатчина, 2016. – 20 с. – (Препринт / «ПИЯФ» НИЦ
«Курчатовский институт»; №29914)
5
Структура и объем. Диссертация состоит из введения, двух глав и заключения. Полный
объем диссертации составляет 51 страницу, 27 рисунков и 8 таблиц.
6
1
1.1
Большие экспериментальные установки CERN с ГРД
Детекторы и газовые смеси
Проблема загрязнения атмосферы – одна из наиболее актуальных проблем, с которой
столкнулось человечество. Опасность промышленного загрязнения атмосферы заключается не
только в том, что в воздух попадают губительные для всего живого вещества, но и в
формировании условий, стимулирующих развитие парникового эффекта, вызывающего
глобальные изменения климата Земли.
Основным источником загрязнения атмосферы в результате деятельности человека, как
показывают последние исследования [15], являются промышленные выбросы углекислого газа.
За последние 200 лет концентрация CO2 в атмосфере Земли выросла почти на 30%. Процесс
загрязнения настолько масштабен, что привел к наблюдаемым глобальным экологическим
явлениям, а именно к потеплению климата и сопровождающим его явлениям. Источниками
загрязнения атмосферы могут быть не только промышленные предприятия, но и масштабные
экспериментальные физические установки, такие как БАК.
На сегодняшний день на Большом Адронном Коллайдере в CERN (European Organization
for Nuclear Research) работают 4 основные экспериментальные установки – детекторы ATLAS,
ALICE, LHCb
и CMS.
Каждая из этих гигантских установок на 70-80% состоит из
газоразрядных субдетекторов пропорциональных камер, дрейфовых камер и резистивных
плоско-параллельных камер. Как показано на Рис. 1.1, все эти гигантские детекторы
расположены на кольце БАК вокруг точек столкновения пучков. Основным назначением
детекторов ATLAS и CMS является поиск бозона Хиггса и открытие «нестандартной физики»,
в частности тёмной материи. Детектор ALICE предназначен для изучения кварк-глюонной
плазмы в столкновениях тяжёлых ионов свинца, а LHCb – для исследования физики b-кварков,
что позволит лучше понять различия между материей и антиматерией.
7
Рисунок 1.1: Экспериментальные установки БАК [16].
В состав вышеназванных экспериментальных установок БАК входит 10 различных
газоразрядных детекторов, которые снабжаются рабочими газовыми смесями, состоящими из 10
газов: Ar, CO2, CF4, C4F10, iC4H10, SF6, nC5H12, H2, O2, CH4, N2 (в зависимости от типа детектора).
Для обеспечения детекторов рабочими газовыми смесями на БАК используется 28 различных
газовых систем. Из них:
•
14 систем замкнутой рециркуляции газовой смеси;
•
11 систем открытого продува газовой смеси;
•
3 системы подачи сжиженного газа.
Система газового снабжения каждого эксперимента представляет собой комплекс
сложных установок и аппаратуры, занимающих более нескольких сотен метров, и обеспечивает
чрезвычайно высокую надежность с точки зрения устойчивости
качества газовой смеси,
поступаемой в детекторы. Подача рабочих газовых смесей осуществляется из газосмесительных
установок, находящихся на поверхности. Под высоким давлением (24 бар) газ поступает на
распределительную станцию в шахте на глубине до 235 м. Пройдя через детекторы
экспериментальной установки, газовая смесь собирается из распределительных линий в
циркуляционную систему низкого давления (
1.0 бар) с эффективностью циркуляции до 90%.
В Таблице 1.1 представлены объемы и виды газовых смесей для каждого из
субдетекторов, а также объем используемого газа и вклад каждого из них в формирование
парникового эффекта, при скорости подачи рабочей газовой смеси со скоростью 1 объем в
8
сутки.
Таблица 1.1: Газовые смеси, используемые в экспериментальных установках БАК.
Субдетекторы
Экспериментальная
установка
Газовая смесь
Объем газовой
смеси
Эффективность
циркуляции газов
Вклад в парниковый
эффект
RPC (Resistive Plate
CSC (Cathode Strip
GEM (Gas Electron
Chambers)
ALICE, ATLAS,
Chambers)
Multiplier)
CMS
C2H2F4, iC4H10, SF6
CMS
LHCb
Ar, CO2, CF4
Ar, CO2, CF4
95%, 4.5%, 0.3%
40%, 50%, 10%
45%, 15%, 40%
15 м3 (каждый)
90 м3
< 1 м3
90%
(ALICE открытый
90%
открытый цикл
цикл)
75%
10%
6%
Рассмотрим с экологической точки зрения работу одного из субдетекторов CSC (Cathode
strip chamber – катодные стриповые камеры) установки CMS, для того чтобы оценить
воздействие газа CF4 на окружающую среду. Катодные многопроволочные стриповые камеры
являются самым большим по объёму газа субдетектором
эксперимента СМS, при этом
используется газовая смесь 40%Ar + 50%CO2 + 10%CF4. На Рис. 1.2 представлен поперечный
срез детектора CMS.
Рисунок 1.2: Сектор детектора CMS [17].
9
Общее количество катодных стриповых камер в эксперименте CMS составляет 540 штук.
Каждый
модуль CSC состоит из
пространственное
шести слоёв камер, каждая из которых обеспечивает
разрешение мюонов 150 микрон. Необходимость перекрытия больших
телесных углов для регистрации мюонов, вылетающих из точки взаимодействия, обусловила
значительные
размеры всех модулей CSC, в результате общий газовый объем детекторов
составляет 90 м3. С помощью системы газового снабжения в модули катодных стриповых камер
установки CMS рабочая газовая смесь подается со скоростью 2 объема в сутки. Таким образом,
общий объем газа, поступающий в камеры CSC, равен 180 м3. Газовая система обеспечения CSC
работает по схеме рециркуляции с эффективностью 90%, поэтому при каждой смене объема
около 10% смеси заменяется в системе на смесь из внешней линии снабжения. Для извлечения
СF4 из отработавшей газовой смеси 40%Ar + 50%CO2 + 10%CF4 в CMS создана специальная
система рекуперации CF4. Эффективность этой системы составляет менее 75%. Следовательно,
за сутки в атмосферу из газовой системы CSC выбрасывается 0.45 м3 газа СF4, что, учитывая его
плотность СF4 = 3.65 г/л, составляет 1.6 кг. Как уже было сказано, парниковый эффект для газов
нормируется в единицах воздействия на атмосферу одного килограмма углекислого газа, и в
единицах GWP СF4 имеет коэффициент 7390 кг СО 2 [18]. Поэтому эквивалентом суточного
выброса камер CSC является 1.6 7390 12 тонн СО2. Годовой выброс СF4 в единицах GWP
составляет тогда около 4380 тонн СO2. Необходимо отметить, что время жизни СО 2 составаляет
50-200 лет, а время жизни CF 4 в атмосфере составляет 50000 лет. Учитывая, что одна тепловая
электростанция, работающая на природном газе, выбрасывает в год около 11000 тонн С O2,
можно сказать, что работа только одних камер СSС является мощным источником загрязнения
атмосферы Земли. Суммарный выброс всех установок БАК сопоставим в единицах GWP с
выбросом десятка тепловых электростанций. Очевидно, что такая система газового снабжения
нуждается в оптимизации состава рабочей газовой смеси.
Потенциалы глобального потепления некоторых парниковых газов перечислены в
Таблице 1.2. В качестве эталонного газа взят CO 2.
10
Таблица 1.2: Газы-кандидаты и парниковые газы, используемые в качестве рабочих смесей [15].
ODP
Chemical
Refrigerant
formula
identifier
Chloropentafluoroethane
C2ClF3
R115
7370
0.44
Hexafluoroethane
C2F6
R116
-
-
Propane
C3H8
R290
3
0
Isobutane
C4H10
R600a
3
0
Carbon Dioxide
CO2
R744
1
0
Tetrafluoromethane
CF3I
R1311
0
0
Trifluoromethane
CHF3
R23
0
0
Sulfur Hexafluoride
SF6
R7146
23000
0.04
1, 1, 1, 2 – Tetrafluoroethane
CH2FCF3
R134a
1430
0
Molecular name
GWP
(Ozone Depletion
Potential)
Для замены CF4 и других газов на более экологичные газовая смесь должна обеспечивать
все требуемые раюочие характеристики детектора. Чтобы разработать рабочую газовую смесь
для того или иного типа детектора, необходимо изучить, как различные газовые компоненты и
состав смеси влияют на характеристики детектора.
1. 2
Взаимодействие заряженных частиц в газовых средах
1. 2.1
Удельные потери энергии
Как известно, при прохождении частицы через вещество часть энергии частицы
передается атомам вещества. Энергия, теряемая заряженной частицей в слое данного вещества
единичной толщины, называется тормозной способностью и может быть охарактеризована
величиной удельных потерь
⟨ ⟩
dE
,
dx
где E – энергия частицы, x – длина отрезка, прошедшего
частицей. Сама величина существенно зависит от скорости частицы.
1. 2.2
Ионизационные потери тяжелой заряженной частицы в среде
Потери энергии заряженной частицы в среде главным образом вызваны
взаимодействиями с атомными оболочками. Среднее значение потерь элементарной заряженной
11
частицы массой m >> me в результате взаимодействий с атомными оболочками при прохождении
пути длиной dx описывается формулой Бете-Блоха:
( )
dE
dx
=
4 π Na α
ion
me β
2
2
ρ
Z
( C+ln β2 γ2 ) (1.1)
A
2
me – масса электрона, α=
частицы),
γ=
e
1
≈
hc 137
(постоянная тонкой структуры), β
= v/c (скорость
1
, Z – заряд среды, A – массовое число среды, NA = 6,02∙1023 моль, ρ –
√1− β2
плотность среды, C ≈ 8 (в зависимости от материала).
Потери энергии не зависят от массы проходящей через вещество частицы, но
существенно зависят от скорости частицы. Кривая Бете-Блоха (Рис. 1.3) после резкого спада
достигает своего минимума при β γ ≈ 4. При βγ → 1 dE/dx постепенно возрастает, такой
"релятивистский подъем" вызван логарифмическим членом в Ур. 1.1.
Рисунок 1.3: Зависимость ионизационных потерь энергии от энергии тяжелых частиц [].
1. 2.3
Ионизационные потери для электронов
В области низких энергий электронов (E < 1 МэВ) определяющий вклад в потери
энергии дают неупругие ионизационные процессы взаимодействия с атомными электронами,
включающие ионизацию атомов. Передаваемая в одном столкновении энергия в среднем очень
мала и при движении в веществе потери складываются из очень большого числа таких малых
потерь.
Электроны с большей энергией теряют большую часть своей энергии из-за тормозного
излучения. Для релятивистских электронов формула Бете-Блоха приобретает вид
12
2
( )
dE
dx
(
)
4π NAα Z 1
=
ρ
ln γ 2 + D D ≈ 7(1.2)
m
A
2
ion , e
e
Ионизационные потери энергии ограничены величиной, называемой критической
энергией [19]:
EC =
710 МэВ
(1.3)
Z +0.92
Для аргона (Z = 18) EC = 37,5 МэВ.
1. 2.4
Радиационная длина
Радиационная длина – это средняя толщина вещества, на которой энергия электрона
уменьшается в e раз. Радиационные длины для различных веществ приведены в Таблице 1.3.
Радиационная длина атома может быть найдена [20]
( )
( )
1
ħ 2
183
=4
nZ (Z +1)α 3 ln 1
X0
mc
Z3
(1.4)
Таблица 1.3: Значения X0 для различных веществ.
Вещество
H (жидкость)
Воздух
Ar (жидкость)
C
Al
Fe
Cu
Pb
1. 2.5
X0
г/см2
см
63.04
36.62
19.55
42.70
24.01
13.84
12.86
7.43
890.4
3.04·104
14.00
19.32
8.897
1.757
1.436
0.5612
Образование ионных пар
При прохождении налетающей частицы через вещество она взаимодействует с
веществом и теряет часть своей энергии на ионизацию атомов. В результате этого процессах
образуется ионная пара. Число ионизаций, образованных налетающей частицей на 1 см пути,
выражается через NP в единицах см-1. Если каждая образовавшаяся ионная пара обладает
13
достаточной кинетической энергией, то она способна также произвести ионную пару, и этот
процесс будет называться вторичной ионизацией. Число, равное сумме первичных и вторичных
ионов на единицу длины, обозначается как NT, и главным образом зависит от вещества, энергии
и массы налетающей частицы. Величина этого параметра определяет количество и размеры
лавин, образующихся заряженной частицей в детекторе, газовая смесь которого находится под
действием электрического поля.
Связь между первичной ионизацией и атомным номером может быть найдена с помощью
экспериментальных данных [21]:
N P=3.996
Zm
Z
0.4
−0.025
( )
Zm
Z
0.4
см−1 ,(1.5)
которая выполняется при нормальных условиях (1 атм и 20°С). Эта формула дает хорошее
согласие с экспериментальными данными для углеводородов, но почти не подходит
для
молекул, содержащих главным образом фтор, к примеру, для CF4. Отличие результатов,
полученных экспериментально и по этой формуле, составляет до 30%.
Довольно трудно оценить общее число ионов. Несмотря на то, что не существует общей
формулы, есть самый простой способ – это использование сечения рассеяния для расчета
электронов первичной ионизации и использование моделирования Монте-Карло для
прослеживания вторичных электронов. Число пар ионов не зависит от энергии и массы
налетающей частицы. Для налетающей частицы величина W = ∆E/N определяет среднюю
энергию, необходимую для образования одной пары ионов. Величина W представляет собой
функцию, слабо меняющуюся с изменением энергии частицы [22], и таким образом, ее можно
рассматривать постоянной для некоторого интервала значений энергии. Общее число ионов на
единицу длины можно найти, используя
N T =ρ
dE 1
(1.6)
dx W
Среднее значение W для газовой смеси может быть найдено, если известны значения W
для газовых компонент [21]
W =∑ [ f ( nm ) Z ( nm ) W ( nm ) ] / ∑ [ f ( nm ) Z ( nm ) ] ,(1.7)
m
m
гд е nm указывает на индекс молекулы, а f(nm) – относительное число молекул данного вида
газовой смеси.
14
Тяжело предсказать значение W, и не существует прямого способа получить точную
оценку, основанную только на экспериментальных данных.
В Таблице 1.4 приведены параметры взаимодействия заряженных частиц с газовыми
смесями. Приведённые данные для газовых смесей вычислены на основе известных рабочих
параметров для составляющих газовых компонент [15].
Таблица 1.4. Некоторые параметры газовых смесей [15].
−
Газовая смесь
1. 3
⟨ ⟩
dE
dx
, МэВ∙см2/г
min
L ,
г/см
2
N P , см-1
40%Ar+50%CO2+10%CF4
1.67
29.21
36.72
40%Ar+55%CO2+5%CF3I
1.667
28.31
45.97
40%Ar+59%CO2+1%CF3I
1.67
29.29
36.43
40%Ar+50%CO2+10%CF3I
1.64
27.07
57.56
Основные характеристики газовых смесей
1. 3.1
Коэффициент Таунсенда
Среднее расстояние, которое прошел электрон между двумя последовательными
столкновениями в газе, называется средней длиной свободного пробега. Обратная величина,
определяющая число актов ионизаций на пути в один сантиметр, называется первым
коэффициентов Таунсенда. Этот параметр определяет газовое усиление для различных газов.
Если n0 – число первичных электронов в электрическом поле в отсутствии усиления, n – число
электронов в лавине, пройденной расстояние x, которое определяется n = n0eαx, а газовое
усиление G дается выражением G ≡ n/n0 = eαx. Первый коэффициент Таунсенда зависит от
природы газа, электрического поля и давления. Второй коэффициент Таунсенда η, равный числу
электронов, выбиваемых с поверхности одним положительным ионом. С учетом эмиссии
электронов газовое усиление примет вид
αx
G≡
1. 3.2
e
(1.8)
1−η(e αx−1)
Скорость дрейфа
15
В отсутствии электрического поля движение электронов в газе хаотичное, при этом они
обладают тепловой энергией
двигаться в направлении
3
kT . Под действием электрического поля электроны начинают
2
приложенного поля со средней скоростью vd (среднее расстояние,
проходимое дрейфующим электроном, в единицу времени). Скорость дрейфа зависит от
давления, температуры, а также наличие загрязняющих веществ, таких как вода и кислород,
может повлиять на величину скорости дрейфа. Компьютерная программа MAGBOLTZ
позволяет вычислять скорости дрейфа электронов с использованием метода Монте-Карло путем
отслеживания многочисленных столкновений электронов с молекулами газа на
микроскопическом уровне.
1. 3.3
Диффузия
Электроны и ионы в газах под действием электрического поля двигаются вдоль
направления поля. Но присутствует часть электронов, которая вследствие рассеяния на
молекулах газа отклоняются от этого направления. Рассеяние, приводящее к изменению
скорости дрейфа, называется продольной диффузией, а рассеяние, приводящее к смещениям,
называется поперечной диффузией. В холодных газах, таких как углекислый газ, диффузия
мала, а величина скорости дрейфа не велика. Это приводит к нелинейной X-T-зависимости.
Такие газы, как аргон, обладают более высокой диффузией; но при смешивании с
многоатомными (органическими) газами в большинстве случаев диффузия уменьшается, а
скорость дрейфа возрастает.
1. 3.4
Диффузия ионов и электронов в нулевом магнитном поле
Заряженные частицы, образованные при ионизации, довольно быстро теряют свою
энергию из-за многократных столкновений с молекулами газа, поэтому необходимо учитывать
распределение средней тепловой энергии газа. Среднее значение тепловой энергии в
кинетической теории газов дается выражением сТ = 3/2kT = 0,04 эВ при нормальных условиях, а
функция распределение Максвелла имеет вид
−(
F ( c )=C √ ϵ e
ϵ
)
kT
.(1.9)
При отсутствии других эффектов распределение заряженных частиц подчиняется закону
Гаусса:
16
dN
1
−( x / 4 Dt )
=
e
dx ,(1.10)
N √ 4 πDt
2
где dN/N
– часть заряженных частиц, находящихся в элементе dx на расстоянии x от
начала координат и через некоторое время t; D обозначает коэффициент диффузии.
Среднеквадратичное распределения, или стандартное отклонение, определяется соответственно
σ X =√ 2 Dt или σ V = √ 6 Dt ,(1.11)
для линейной и объемной диффузии. В качестве примера на Рис. 1 . 4
изображено
пространственное распределение ионов, образованных в воздухе при нормальных условиях и
различных временных интервалов; значения коэффициентов диффузии нескольких ионов в
соответствующих газах приведены в Таблице 1.5.
Рисунок 1.4: Пространственное распределение ионов, образованных в воздухе, при нормальных условиях
и при различных временных интервалах [23].
17
Таблица 1.5: Значения длин свободного пробега, скоростей дрейфа, коэффициентов диффузии и
подвижности для различных газов при нормальных условиях [23-26].
Газ
λ, см
v, см/с
D+, см2/с
μ+, см2 с-1 В-1
H2
He
Ar
O2
H2O
1.8 10-5
2.8 10-5
1.0 10-5
1.0 10-5
1.0 10-5
2 105
1.4 105
4.4 104
5.0 104
7.1 104
0.34
0.26
0.04
0.06
0.02
13.0
10.2
1.7
2.2
0.7
В процессе диффузии происходит столкновение ионов с молекулами газа; средняя длина
свободного пробега, средняя скорость, и время между столкновениями также приведены в
таблице. Электроны движутся гораздо быстрее, чем ионы (из-за их небольшой массы), их
средняя тепловая скорость составляет примерно 10 7 см/сек; Кроме того, из-за их небольшого
размера, их средняя длина свободного пробега обычно в четыре раза больше, чем средняя длина
свободного пробега ионов в том же газе.
В газе положительный ион может рекомбинировать с отрицательно заряженным ионом
или с электроном. Процесс передачи заряда также возможен с молекулой собственного газа или
с молекулами другого рода, обладающими более низкими потенциалами ионизации.
Электроны же могут быть нейтрализованы ионом, могут быть присоединены к молекуле,
имеющей электронное сродство (электроотрицательное), или могут быть поглощены на стенках.
Вероятность прилипания h практически равна нулю для водорода и всех благородных газов, но
для других газов это величина имеет вполне определенные значения, которые приведены в
Таблице 1.6. В таблице также представлено среднее время прилипания t = (hN)-1, где N – число
столкновений за единицу времени. Можно заметить, что, например, в кислороде среднее время,
за которое произойдет прилипание электрона, составляет 140 нс. Коэффициент прилипания
сильно зависит от электрического поля.
18
Таблица 1.6: Значения коэффициентов прилипания, чисел столкновений и средних времен прилипания
для электронов в различных газах при нормальных условиях [27, 23, 26].
Газ
h
N, с-1
t, с
CO2
O2
H2O
Cl
6.2 10-9
2.5 10-5
2.5 10-5
4.8 10-4
2.2 1011
2.1 1011
2.8 1011
4.5 1011
0.71 10-7
1.9 10-7
1.4 10-7
4.7 10-9
1. 3.5
Подвижность ионов
При приложении электрического поля к объему газа происходит движение ионов вдоль
направления поля. Средняя скорость этого медленного движения (не путать с непосредственной
скоростью ионов) называется дрейфовой скоростью w+, и она оказывается прямо
пропорциональной уменьшающемуся полю Е/p до очень высоких полей, p – давление газа.
Поэтому величину μ+, называемую подвижностью, удобно определить
+¿
w
.(1.12)
E
+¿=¿
¿
μ
Значение подвижности для каждого иона, движущегося в данном газе, вполне
конкретное. Постоянная подвижность является прямым следствием того факта, что вплоть до
очень высоких полей средняя энергия ионов практически не изменена; позже будет видно, что
это для электронов это совсем не так.
Связь подвижности с коэффициентом диффузии имеет вид
+¿
μ =
kT
.(1.13)
e
+¿
D
¿
¿
Значения подвижности и коэффициента диффузии для ионов, движущихся в газе,
приведены в Таблице 1.5, а Таблица 1.7 дает подвижность некоторых ионов, дрейфующих в
газах, которые обычно используются в пропорциональных и дрейфовых камерах.
19
Таблица 1.7. Экспериментальные значения подвижности для нескольких ионов в различных газах при
нормальных условиях [25].
Газ
Ar
IsoC4H10
(OCH3)2CH2
Ar
IsoC4H10
Ar
CH4
Ar
CO2
Ионы
(OCH3)2C +¿¿
H2
(OCH3)2C +¿¿
H2
(OCH3)2C +¿¿
H2
IsoC4 +¿¿
H 10
IsoC4 +¿¿
H 10
+¿
C H¿
4
+¿
C H¿
4
+¿
C O¿
2
+¿
C O¿
2
Подвижность, см2 с-1 В-1
1.51
0.55
0.26
1.56
0.61
1.87
2.26
1.72
1.09
В многокомпонентной смеси G1, G2, … , Gn подвижность μi+ иона Gi+ определяется
соотношением (закон Бланка):
pj
+¿ ,(1.14)
μij
n
μi =∑ ¿
+¿
j=1
1
¿
где рj – это объемная концентрация газа j в смеси, μij+ – подвижность иона Gi+ в газе Gj [28].
1. 3.6
Угол Лоренца
В скрещенных электрическом и магнитном полях сила Лоренца заставляет электроны
дрейфовать под углом к силовым линиям электрического поля, называемым углом Лоренца и
определяемым произведением величины магнитного поля на время между взаимодействиями.
Таким образом, угол Лоренца тем больше, чем меньше сечение взаимодействия электрона с
молекулами газа. Также угол Лоренца зависит от величин магнитного и электрического полей.
Обычно угол Лоренца велик при малых электрических полях и падает с увеличением
электрического поля.
20
1. 4
Принцип работы газоразрядного детектора с усилением – пропорционального
счетчика
Рассмотрим тонкий слой газа, например, 1 см слой газовой смеси 70% Ar + 30% Изобутан
при нормальных условиях между двумя плоскими электродами. Минимально ионизирующая
заряженная частица образуют в результате ионизации в этом слое около 120 пар ионов; если
этот заряд собирается на одном электроде, соответствующий сигнал будет равен
V=
ne
,(1.15)
C
где при n = 120 и типичной емкости С = 10 пФ, V
~ 2 мкВ. Такой сигнал невозможно
зарегистрировать, поскольку он будет соразмерен с шумами предусилителей.
Самыми распространенными детекторами ионизирующих излучений являются
газонаполненные детекторы. Конструктивно газонаполненные детекторы представляют собой
замкнутые газовые объемы с двумя электродами, к которым приложена разность потенциалов.
Если между электродами приложено сильное электрическое поле, то происходит лавинное
усиление, что приводит к повышению амплитуды сигнала на несколько порядков. Но такая
геометрия детектора обладает несколькими серьезными ограничениями. К примеру, величина
зарегистрированного сигнала не пропорциональна энергии, образованной в результате
ионизации. Коаксиальная геометрия электродов позволяет преодолеть это ограничение (Рис.
1.5).
Рисунок 1.5: Пропорциональный счетчик и форма электрического поля вблизи анодной проволочки.
Вблизи анода поле возрастает настолько, что становится возможным размножение электронов [28].
Тонкая металлическая проволока натягивается вдоль оси цилиндра и изолируется от него
для того, чтобы между электродами можно было приложить разность потенциалов. Полярность
выбирается таким образом, чтобы центральная проволочка была положительна относительно
21
внешнего цилиндра. Электрическое поле максимально на поверхности анодной проволоки и
быстро спадает как r-1 в направлении катода; при использовании тонких проволочек можно
получить очень высокие значения полей вблизи анода. Работа счетчика может быть объяснена
следующим образом. Электроны, образовавшиеся в результате первичной ионизации, дрейфуют
к аноду, а положительные ионы соответственно к катоду. Но вблизи анода, как правило, на
расстоянии нескольких радиусов проволочки, где поле становится настолько сильным,
происходит усиление; начинает развиваться каплеобразная лавина, в передней части которой
располагаются электроны, а в задней части – ионы.
Рисунок 1.6: Развитие лавины в пропорциональном счетчике [28].
За счет диффузии и малого радиуса анода эта лавина окружает проволочку, как показано
на Рис. 1.6. Положительные ионы начинают дрейфовать к катоду.
Зависимость развития различных мод газового усиления в газоразрядном детекторе от
напряжения приведена на Рис. 1.7. Здесь показана вольт-амперная характеристика
газоразрядного промежутка, где красная кривая соответствует α-частицам, обладающим
высокой удельной ионизацией, а синяя – β-частицам.
Область рекомбинации
В этой области напряжение мало для создания электрического поля, достаточного для
сбора всех образованных в результате первичной ионизации частиц. Вместо этого происходят
два процесса: медленный дрейф ионов к электродам и их рекомбинация. Очевидно, что с
увеличением напряжения число собранных ионов возрастает, но рост происходит до тех пор,
пока не будет достигнуто значение насыщенности. По этой причине газообразные
ионизационные камеры не работают при этих напряжениях.
22
Рисунок 1.7: Вольт-амперная характеристика газоразрядного промежутка. Красная линия
соответствует α-частицам, синяя – β-частицам [29].
Область насыщения
При значениях напряжения, равных около 100 В, достигается насыщение: этого
напряжения недостаточно, чтобы обеспечить сбор всех ионных пар, образованных заряженной
частицей. Но в то же время напряжение еще не столь велико для создания вторичной
ионизации. Соответственно, газовое усиление в этой области отсутствует.
Пропорциональная область
Напряжение в этой области настолько велико, что возникает большой градиент
потенциала вблизи анода. Оказывается возможным газовое усиление, так как скорости
первичных электронов достаточно для создания вторичной ионизации, лавины Таунсенда. В
этой области амплитуда импульса напряжения пропорциональна количеству образованных
ионов.
Область ограниченной пропорциональности
По мере увеличения напряжения электроны вылетают с поверхности катода под
действием электрического поля, что приводит к образованию лавины Таунсенда,
распространяющейся вдоль анода, в то время как положительные ионы стремятся остаться там,
23
где они были образованы. На этом участке уже сказывается влияние пространственного заряда
положительных ионов и соседних лавин, мешающих взаимному развитию. В этом случае
величина начальной ионизации влияет на коэффициент газового усиления таким образом, что
его возможный рост ограничивается при сильной первичной ионизации по сравнению с
наиболее быстрым увеличением при слабой ионизации.
Область Гейгера-Мюллера
На этом участке газовое усиление возрастает настолько, что собираемый заряд не
зависит от первичной ионизации. Это так называемая область Гейгера-Мюллера. Однако разряд,
как и в предыдущих областях, остается вынужденным, то есть начинается после прохождения
ионизирующей частицы. Это область работы счётчиков Гейгера-Мюллера.
Область непрерывного разряда
И, наконец, напряжение увеличивается настолько, что происходит непрерывный разряд.
Плато Гейгера – это последняя область, используемая для непрерывно приложенных
напряжений. Области, находящиеся выше, не используются для регистрации частиц.
1. 5
Метод исследования скорости дрейфа электронов
1. 5.1
Дрейфовые камеры
Одной из основных характеристик газовой смеси является скорость дрейфа электронов,
поскольку эта характеристика влияет на временное разрешение детектора и его быстродействие.
Дрейфовые камеры (Рис. 1.8) – это газонаполненные детекторы, предназначенные для
регистрации заряженных частиц. Заряженная частица, проходя через газовый объем в момент
времени t = 0 (этот момент регистрируется триггерной системой – обычно органический
сцинтиллятор), ионизирует атомы газа, что приводит к образованию электрон-ионных пар. В
электрическом поле, которое прикладывается к газовому объему, образовавшиеся ионы и
электроны начинают дрейфовать в противоположных направлениях: электроны двигаются по
направлению к аноду, а положительно заряженные ионы – к катоду. В качестве анода обычно
используется тонкая проволочка диаметром 20 мкм - 50 мкм для получения электрического поля
высокого напряжения вблизи анода. В окрестности анодной проволочки электроны ускоряются
этим полем настолько, что они становятся способными ионизировать другие молекулы газа.
Поскольку этот процесс повторяется несколько раз, то электроны образуют размножающуюся
24
лавину, которая регистрируется как отрицательный электрический импульс. Процесс
размножения электронов в газе называется газовым усилением.
Рисунок 1.8: Принцип работы дрейфовой камеры.
Обычно для измерения скорости дрейфа применяется следующий подход. Измеряется время
дрейфа электронов ионизации от трека частицы с измеренной координатой - x, до анодной
проволочки в однородном электрическом поле. Скорость дрейфа определяется из измеренной
зависимости времени сбора электронов на аноде при заданных положениях источника
электронов по формуле
v=
dx
,(1.16)
dt
где dt = t2 - t1 – разница между временами прихода электронов на анод, соответствующим двум
различным положениям источника x2 и x1, соответственно, dx = x2 - x1.
На графике ниже (Рис. 1.9) представлена зависимость скоростей дрейфа для большинства
широко используемых газовых смесей от приложенного электрического поля.
Рисунок 1.9: Зависимости скоростей дрейфа электрона от электрического поля для нескольких
стандартных газовых смесей при нормальных условиях [28].
25
1. 5.2
Примеры установок измерения скорости дрейфа из литературы
В литературе представлено несколько методик по определению скоростей дрейфа в газах.
Они отличаются применяемой техникой образования электронов. В некоторых
экспериментальных установках электроны образовывались под действием УФ-лазера, либо
радиоактивным источником β-частиц, например 90Sr.
Экспериментальная установка, которая использовалась в работе [30], схематически
представлена на Рис. 1.10. Дрейфовая камера содержит 40 электродов, выполненных в виде
колец из нержавеющей стали, расположенных на расстоянии 10 мм ± 0.1 мм с использованием
тефлоновой прокладки. Электроды соединены между собой посредством цепи резисторов для
создания однородного поля. Самый верхний электрод подсоединен к источнику высокого
напряжения, который обеспечивает выходное напряжение вплоть до 50 кВ, самый нижний
электрод находится под нулевым потенциалом. В каждом электроде просверлена дырка
диаметром 25 мм. Электроны в однородном поле дрейфуют сквозь эти дырки к сигнальной
проволочке, расположенной ниже самого последнего электрода. Проволочка частично окружена
u-образным катодом, находящимся под потенциалом земли. Для ограждения сигнальной
проволочки от наведенного заряда отверстие в нижнем электроде покрыто никелевой сеткой.
Газовая смесь в камере ионизируется под действием лазера на азоте с длиной волны 337 нм,
длительностью импульса 500 пс. Лазерный пучок разделяется на два отдельных, которые
фокусируются вдоль оси дрейфовой камеры. Ультрафиолетовый лазер образует электронные
кластеры в основном путем ионизации примесей в газовой смеси. Резервуар высокого давления
оборудован четырьмя парами кварцевых окошек, расположенных напротив друг друга, которые
используются в качестве входа и выхода лазерного пучка. Измеряя ∆t, разность времен прихода
между двумя отдельными сигналами, соответствующие одному лазерному импульсу, мы можем
вычислить v = ∆x/∆t, компонента скорости дрейфа, параллельная электрическому полю E.
26
Рисунок 1.10: Схематичное изображение экспериментальной установки [30].
В начале измерений в резервуаре высокого давления откачивается воздух, и несколько
раз он продувается аргоном. Измерения выполняются вскоре после того, как камера полностью
наполнится газом. Ошибка в определении дрейфовой скорости составляет 0.5%.
В работе [31] представлена другая методика измерения дрейфовой скорости электрона в
дрейфовой камере. Дрейфовая камера (Рис. 1.11) наполнялась газовой смесью аргона и метана
(Ar + 40%CH4) под давлением 1 атм. Дырка диаметром 2 мм, сделанная в катоде, была закрыта 5
мкм медной фольгой. Радиоактивный источник 210Po был размещен на фольге. Кремниевый
детектор, диаметр которого составлял 15 мм, регистрировал α-частицы и был расположен на
расстоянии 6 мм от источника α-частиц. Сигнал от кремниевого детектора использовался в
качестве "старта" для измерений времени дрейфа. Ядра отдачи, образованные в результате αраспада, испускаются в объеме дрейфовой камеры и образовывают около 10 3 электронов вблизи
(≤ 50 мкм) катода.
27
Рисунок 1.11: Схематическое изображение дрейфовой ячейки, используемой для измерения скоростей
дрейфа электронов. S – сигнальная проволочка диаметром 30 мкм, P – полеформирующая проволочка
диаметром 75 мкм [31].
В работе [32] дрейфовая ячейка была разделена на две части (Рис. 1.12). Область дрейфа
электронов сделана из латунных и тефлоновых пластинок. Латунные пластинки, которые
подключены к делителю напряжения, определяют форму дрейфового поля. Сигналы
формируются на анодной проволочке диаметра 20 мкм, расположенной внутри малой камеры,
которая ограничивает область дрейфа. Для проведения измерений дрейфовой скорости самое
большое значение электрического поля составило 2.6 кВ/см.
Рисунок 1.12: Вид дрейфовой ячейки в разрезе [32].
Рисунок 1.13: Вид экспериментальной
установки [32].
Экспериментальная установка изображена на Рис. 1.13. Дрейфовая ячейка установлена
на столике, который перемещается с высокой точностью. Источник электронов 90Sr помещен над
дрейфовой камерой. В самой ячейке имеются четыре окна шириной 0.38 мм, два окна
расположены сверху, два оставшихся — снизу, которые покрыты тонкой пленкой из пластика.
Через эти окнам электроны от источника 90Sr попадают в дрейфовую область. Помимо
свинцовых коллиматоров,
латунные пластины осуществляют окончательную коллимацию
28
электронного пучка. На небольшом пластиковом сцинтилляторе генерируется начальный сигнал
от электронов, проникающих внутрь камеры, который далее попадает на преобразователь
время-код. При этом и источник электронов, и фотоэлектронный умножитель остаются
неподвижными. Информацию о дрейфовой скорости получают из величины дрейфовой области
(расстояние между верхним окном и нижним) и разницей между измеренными временами
прихода сигнала. Расстояние между окнами может быть измерено с точностью до ±6 мкм.
Типичное распределение времен прихода показано на Рис. 1.14. Ширина полученного пика
составляет 6-7 нс. Статистическая ошибка в определении времен прихода сигнала равна ±0.2 нс,
которая соответствует ошибки в определении расстояния примерно ±10 мкм.
Рисунок 1.14: Распределение времен прихода электронов [32].
К недостаткам перечисленных методов относится постоянное положение источников
относительно анода, т.е. время дрейфа определяется только по двум экспериментальным точкам.
В основе данной работы лежит новый принцип с использованием системы передвижных
коллиматоров. При таком подходе время прихода электронов от различных положений
источника описывается набором экспериментальных точек, что в перспективе дает возможность
измерения скорости дрейфа с более высокой точностью.
В связи с идеологией измерений особые требования должны предъявляться к
механической точности системы перемещения, а также получения однородного электрического
поля в детекторе в рабочей области облучений.
29
2 Разработка дрейфовой камеры
Для изучения характеристик новых рабочих газовых смесей, в частности измерения
скорости дрейфа электронов при низких напряженностей полей Е = 500÷5000 В/см (область
работы большинства детекторов), был разработан компактный дрейфовый детектор. За основу
конструкции детектора взят модуль дрейфового детектора FTC [33].
2.1
Описание методов моделирования
Моделирование дрейфовой камеры осуществлялось в программном пакете GARFIELD
версии 7.45 – это программа с открытым кодом, разработанная Др.Робом Веенхофом из CERN ,
предназначена для расчета плоских и трехмерных дрейфовых камер. Данный программный
пакет позволяет получать распределение полей (как электрического, так и магнитного) внутри
дрейфовых камер, строить линии дрейфа электронов и ионов. Также эта программа может
помочь в поиске оптимальных значений полей для конкретной геометрии [34].
Для расчета характеристик газовых смесей, таких как коэффициента Таунсенда,
продольной и поперечной диффузии, скорости дрейфа электронов, используется программа
MAGBOLTZ [35]. Программный пакет MAGBOLTZ предназначен для численного решения
кинетического уравнения Больцмана для электронов в газовых смесях под воздействием
электрического и магнитного полей.
Для расчета ионизационных потерь энергии высокоэнергичной частицы в газах с учетом
дельта-электронов используется программа HEED [36]. Также существует возможность учета
многократного рассеяния влетающей частицы. HEED может воспроизводить процесс
поглощения фотонов в газовых детекторах при фотоионизации. Данная программа, как и
MAGBOLTZ, в основном использовалась (и продолжает использоваться) как компонент
упомянутой выше программы GARFIELD.
2.2
Моделирование детектора в программе Garfield
Принципиальная схема установки представлена на Рис. 2.1. Дрейфовая камера состоит из
двух областей: рабочей области детектора с однородным электрическим полем и
регистрирующей области (анодная плоскость).
Формирование узкого пучка электронов от источника 90Sr осуществляется системой из
двух соосно связанных коллиматоров (7), которые перемещаются вместе с источником, в
плоскости перпендикулярной плоскости анода. При прохождении электрона через систему
30
коллиматоров сцинтилляционные счетчики S1 и S2, предназначенные для регистрации времени
пролета электронов от источника 90Sr, дают быстрые «стартовые» сигналы. При этом электроны
от источника, проходя через рабочий объем камеры, ионизируют молекулы газа. Электроны
первичной ионизации под действием электрического поля дрейфуют к сигнальной проволочке
(1), инициируя «стоповые» сигналы.
Модуль «И»
ПУ
Рисунок 2.1: Блок-схема установки для измерения скорости дрейфа электронов: 1 – сигнальная
проволочка, 2 – полеформирующие проволочки, 3 – дрейфовый катод, 4, 5 – катоды, 6 – краевые
проволочки, 7 – система связанных соосных коллиматоров.
Таким образом, событием является наличие сигналов на преобразователе «время-код»
(ПВК): «СТАРТ» сигнал, формируемый логикой «И» при наличии двух сигналов со сборки
сцинтилляторов (и фотоэлектронных умножителей – ФЭУ) S1 и S2, который свидетельствует о
прохождении прямолинейного трека электронов от источника и «СТОП» - цифровой сигнал с
анодного предусилителя (ПУ), фиксирующий время развития лавины на аноде, Рис. 2.1.
Во временном спектре преобразователя «время-код» (ПВК) для каждого положения
источника наблюдаются пики, соответствующие времени собирания электронов на аноде.
Диапазон рабочих положений источника соответствует области однородного электрического
поля для каждого рабочего напряжения на электроде (3) и рассчитывается с использованием
GARFIELD. Измерение разности между пиками на спектре ПВК для различного положения
источника позволяет определить скорость дрейфа электронов. При этом вклад области
нелинейного поля вблизи анодной нити учитывается как постоянная составляющая, не
зависящая от положения источника.
2.3
Конструкция дрейфовой камеры
31
Регистрирующая плоскость дрейфовой камеры, Рис. 2.1, состоит из анодных и
полеформирующих проволочек. Центральная анодная проволочка (1) является сигнальной, а
полефомирующие проволочки необходимы для сглаживания поля между катодами (4) и (5). Эти
катоды расположены на расстоянии 2.5 мм от анодной плоскости.
Для минимизации разброса измеряемого времени дрейфа электронов осуществляется
сбор первичных электронов на аноде лишь с малого участка трека. Для этого в верхнем катоде
(4) имеется узкая щель шириной 1 мм. Все остальные электроны из трека собираются без
газового усиления на полеформирующих
анодных проволочках, пройдя через проволочный
катод (4).
Выбор напряженности электрического поля в дрейфовом промежутке осуществляется
установкой напряжения на дрейфовом электроде (3), расположенный на расстоянии 21 мм от
анода. Формирование однородного электрического поля обеспечивается системой электродов:
катодные плоскости (3) и (4) и боковыми электродами (6) с распределенным потенциалом.
Двойной ряд проволочек электродов (6) позволяет обеспечить более однородное электрическое
поле внутри дрейфового объема и эффективное сглаживание краевого эффекта [37-38].
2. 4
Выбор и оптимизация геометрических параметров детектора
Однородность электрического поля внутри дрейфового объема является одним из
ключевых параметров, определяющих точность измерений данной установки. Этот параметр
оценивался величиной ΔΕ/E0, где E0 = E (x = 0, y = 1.2 см) – напряженность электрического поля
в центре области дрейфа (Рис. 2.2).
32
Область дрейфа
(3)
Рисунок 2.2: Конфигурация электрического поля в дрейфовой камере. Напряжение на катоде (3) HVcath
= 8000 В/см.
Однородность поля в диапазоне рабочих напряженностей 500÷5000 В/см достигалась
подбором потенциалов рабочих напряжений на электродах (2), (4), (5). Напряжение на
сигнальной проволочке составляет 0 В.
На Рис. 2.3 приведен график зависимости напряженности электрического поля E = E (x
= 0 , y) при различных напряжениях на дрейфовом катоде (3). Значения напряжений на
дрейфовом катоде HVcath и соответствующие им значения напряженностей электрического
поля в центре области дрейфа E0 представлены на Рис. 2.4.
33
Object 67
Рисунок 2.3: График зависимости напряженности E = E (x = 0, y) при различных напряжениях на
дрейфовом катоде HVcath.
На Рис. 2.3 следует отметить, что при HVcath > 6000 В требуется дополнительно
оптимизировать напряжение на катодах (4) и (5), поскольку наблюдается частичная экранировка
электронов на катоде (4) и снижение эффективности сбора электронов на аноде.
Object 69
Рисунок 2.4: График зависимости напряженности электрического поля в центре области дрейфа E0 от
напряжения на дрейфовом катоде HVcath. Ошибки соответствуют однородности электрического поля
ΔΕ/E0 = 3%.
Оценка рабочего диапазона перемещений источника 90Sr
Одной из задач моделирования было определение допустимого диапазона перемещений
источника 90Sr относительно анода, в котором обеспечивается заданная однородность
электрического поля Е. Результаты оценки, сделанной по данным Рис. 2.3, приведены на Рис.
34
2.5, где для параметра ΔΕ/E0 установлены значения 1 и 3%. Численные значения приведены в
Таблице 2.1.
Object 72
Рисунок 2.5: Диапазоны рабочих положений источника 90Sr при различных напряженностях
электрического поля E0 в дрейфовом промежутке. Показаны верхние и нижние границы диапазонов для
электрических полей с однородностью на уровне ΔΕ/E0 = 1% (▲) и ΔΕ/E0 = 3% (●). Пунктирная линия
соответствует центру дрейфового промежутка y = 12 мм. Пунктирная линия соответствует центру
дрейфового промежутка y = 12 мм, сплошные линии – катоды (3) y = 21 мм и (4) y = 2.5 мм.
На Рис. 2.5 видно, что наибольший диапазон однородного электрического поля в
дрейфовом промежутке (y = 0.5÷2.0 см при ΔΕ/E0 = 3%) достигается при напряженностях поля
E = 1000-1500 В/см, а при E > 1500 В/см рабочий диапазон сужается до y = 0.8÷1.6 см (ΔΕ/E0 =
3%), но сохраняется постоянным вплоть до E = 5000 В/см, что вполне достаточно для
проведения измерений скорости дрейфа.
35
Таблица 2.1: Рабочий диапазон положений источника быстрых электронов 90Sr относительно анода (x =
0, y = 0) и соответствующая однородность электрического поля ΔΕ/E0, где E0=E (x = 0, y = 1.2 см).
HVcath, В
2800
3500
4000
5000
6000
7000
8000
9000
10000
11000
E0, В/см
515.6
894.2
1164.0
1705.0
2246.0
2787.0
3328.0
3869.0
4410.0
4950.0
ymin ÷ ymax, мм
∆E/E0 = 1%
11.2 ÷ 12.9
4.5 ÷ 19.7
4.7 ÷ 20.0
9.8 ÷ 14.8
10.3 ÷ 13.7
10.6 ÷ 13.5
10.5 ÷ 13.3
10.7 ÷ 13.3
10.7 ÷ 13.3
10.8 ÷ 13.2
∆E/E0 = 3%
9.4 ÷ 14.6
4.0 ÷ 20.0*
4.0 ÷ 20.0*
4.0 ÷ 20.0*
6.5 ÷ 19.7
7.3 ÷ 17.1
7.8 ÷ 16.7
8.2 ÷ 16.3
8.5 ÷ 15.9
8.6 ÷ 15.9
∆E/E0, %
(4 мм ÷ 21 мм)
5.9
1.6
1.7
2.8
3.3
3.8
4.1
4.4
4.6
4.7
*ΔΕ/E0 < 3%
2. 5
Вычисление скорости дрейфа
Скорость дрейфа электронов
В качестве тестовой газовой смеси для симуляций взята рабочая газовая смесь
Ar/CO2/CF4, которая используется в катодных стриповых камерах CSC установки CMS (БАК).
На Рис. 2.6 приведен график зависимости скорости дрейфа электронов от напряженности
электрического поля, вычисленный для газовой смеси CSC программой MAGBOLTZ, входящей
в программный пакет GARFIELD. Данный график вычисляется стандартным образом для
любых газов, содержащихся в библиотеке программы. Приведенная зависимость v(E) имеет
широкое плато в области Е = 3000÷20000 В/см, что ослабляет влияние неоднородности поля,
вызванной погрешностью изготовления детектора, на измеряемую скорость дрейфа. Однако
существует большое количество рабочих газовых смесей, где такого плато нет вообще.
36
Рисунок 2.6: Скорость дрейфа электронов для газовой смеси 60%Ar + 30%CO2 + 10%CF4, вычисленная
в программе MAGBOLTZ.
Поэтому задача симуляции включала в себя оценку влияния параметров конструкции
детектора и газовой смеси на точность определяемой скорости дрейфа по алгоритму,
описанному в 1.5.
Симуляция прохождения частицы от источника 90Sr
Образование первичных кластеров ионизации при прохождении быстрых электронов от
источника 90Sr вычислялось программой HEED, также включенной в программный пакет
GARFIELD. Результат вычисления в программе приведен на Рис.2.7.
37
Линии дрейфа ионов
Трек электрона
Линии дрейфа электронов
Рисунок 2.7: Симуляция прохождения электрона от источника 90Sr через объем дрейфовой камеры.
Зеленым цветом выделены кластеры ионизации, красные линии соответствуют линиям дрейфа ионов,
желтые – линиям дрейфа электронов. Напряжение на катоде HVcath = 9100 В.
Для образовавшихся первичных электронов с использованием метода Монте-Карло
строились линии дрейфа для электронов, и вычислялось среднее время дрейфа электронов до
анода (Рис.2.8).
Рисунок 2.8: Временное распределение электронов. Напряженность электрического поля 3912 В/см.
Положение источника 1.2 см, коллиматор 1 мм. Газовая смесь 60%Ar + 30%CO2 + 10%CF4.
38
Зависимость времени дрейфа от напряженности поля
На Рис.2.9. приведен график зависимости определенного в GARFIELD времени дрейфа
электронов до анода от положения источника 90Sr.
1164
V/cm
=E0
Linear
Fit
0.003*t
0.1554
+= x
cmx,
6.1
2.1
8.0
002
003
004
t, ns
005
Рисунок 2.9: Среднее время дрейфа электронов в детекторе в газовой смеси 60%Ar + 30%CO2 +
10%CF4 при напряженности в дрейфовой области E0 = 1164 В/см. Время дрейфа рассчитано для
положений источника x1 = 0.9 см и x2 = 1.6 см ΔΕ/E0 = 3%.
Из Ур. (1.16) скорость дрейфа для данной напряженности поля E0 составит:
¿ t 2 >−¿t 1>¿ ,
y 2− y1
v ( E 0 )=
¿
(2.1)
где s = y2 - y1 – расстояние между положениями источника y2 и y1 для рабочего диапазона
источника (Рис. 2.9) при ΔΕ/E0 = 3%, а <ti> - среднее время дрейфа электронов для i-го
положения источника. Видно, что при удалении от анода погрешность величины t
увеличивается, что связано с увеличением вклада диффузии электронов в газе (см. 2.6).
39
Object 76
Рисунок 2.10: График зависимости скорости дрейфа электронов v(E) в газовой смеси 60%Ar + 30%CO2
+ 10%CF4: сплошная линия – данные MAGBOLTZ (Рис. 2.6), точки – расчетные значения скорости по
Ур.(2). Ошибки по х соответствуют ΔE/E=3%, по у – суммарные погрешности1 для двух положений
источника, с учетом ширины коллиматора 1 мм.
Результат определения скорости дрейфа для напряженностей поля в диапазоне E =
500÷5000 В/см приведен на Рис. 2.10. Видно, что значения скорости дрейфа, определенные при
симуляции детектора по описанному выше алгоритму, совпадают с исходными данными
MAGBOLTZ:
∆ v v изм−v Garf
=
v
v Garf
< 2.5%.
Это означает, что предложенная конструкция
дрейфового детектора позволит правильно интерпретировать экспериментальные данные при
работе реального детектора.
2. 6
Оценка вклада факторов, влияющих на точность определения скорости дрейфа
Основные факторы, которые вносят вклад в погрешность измерения скорости дрейфа,
можно разделить на две группы:
•
систематические погрешности: неоднородность электрического поля, диффузия
•
электронов, ширина коллиматора и др.;
статистические погрешности: статистика набора.
Неоднородность электрического поля
Основным фактором, влияющего на точность определения величины скорости дрейфа,
является неоднородность электрического поля в дрейфовом промежутке детектора.
1 Для наглядности приведены средние квадратические отклонения σ(v). Погрешности, вычисленные из Ур. (2.6), составляют <
1.5% и лежат в пределах размера точек.
40
Из Рис. 2.6 видно, что заданная величина неоднородности электрического поля в
различных областях зависимости v(E) для газовой смеси 60%Ar + 30%CO2 + 10%CF4 дает
существенно отличающиеся результаты. Пусть Δv/v – стабильность скорости, равная
∆ v /v=k Δ E 0 (X%)/v( E 0 )
(2.2)
гд е k – коэффициент, характеризующий динамику зависимости v(E) на данном участке E,
аппроксимированного прямой линией (Рис.2.11), ΔE0(Х%) – абсолютное значение флуктуации
напряженности поля E0, соответствующее ΔE/E = Х%, v(E0) – скорость дрейфа электронов при E
= E0.
Object 84
Рисунок 2.11: График зависимости скорости дрейфа электронов от напряженности электрического
поля в дрейфовом объеме для газовой смеси 60%Ar + 30%CO2 + 10%CF4, построенный по данным
MAGBOLTZ. ki – наклон прямой аппроксимации.
Тогда имеем, для двух областей E
1)
Область E < 2500 В/см, k1 = 29E-4(см/мкс)/(В/см)
E0 = 1164 В/см (HVcath = 4000 В), v = 2.93 см/мкс:
∆ v=k 1 Δ E 0 (1 )=0.034 см/ мкс
=>
∆ v /v ≈ 1.16
∆ v=k 1 Δ E 0 (3 )=0.101 см / мкс
=>
∆ v /v ≈ 3.45
2)
Область E > 2500 В/см, k2 = 2E-4 (см/мкс)/(В/см)
E0 = 3869 В/см (HVcath = 9000 В), v0 = 8.08 см/мкс:
∆ v=k 2 Δ E 0 ( 1 )=0.008 см / мкс
=>
∆ v=k 2 Δ E 0 ( 3 )=0.021 см / мкс
=>
∆ v /v (E 0 )≈ 0.09
∆ v /v (E 0 )≈ 0.29
41
Таким образом, наибольшая ошибка измерения скорости дрейфа электронов в газовой
смеси 60%Ar + 30%CO2 + 10%CF4, обусловленная неоднородностью напряженности поля в
дрейфовой области детектора, будет наблюдаться в области низких напряженной поля E < 2500
В/см, и не будет превышать
∆ v /v
< 3.5% в рабочей области у = 4.1 ÷ 21.0 мм (ΔE/E = 3%).
Для 2500 В/см < E < 5000 В/см ошибка не превышает
∆ v/v
< 0.3% в указанной области у.
Диффузия электронов
Среднеквадратичное отклонение, обусловленное диффузией электронов при дрейфе к
аноду, определяется:
σ D=
√
2 D(E)∙ l
,
v
(2.3)
где D(E) – коэффициент диффузии электронов в электрическом поле с напряженностью E, l и v
– длина и скорость дрейфа электронов соответственно.
Вследствие того, что над дрейфующими электронами в электрическом поле совершается
дополнительная работа, коэффициент диффузии D(E) обладает изотропными свойства.
Соответственно вводят коэффициенты продольной DL(E) и поперечной диффузии DT(E).
При оценке вклада диффузии электронов в ошибку времени прихода электронов на анод
учитывалась только продольная диффузия (вдоль направления дрейфа). График зависимости
σD(E) вычислялся в программе MAGBOLTZ (Рис.2.12).
42
Рисунок 2.12: Зависимости диффузии электронов σD(E) от напряженности электрического поля для
газовой смеси 60%Ar + 30%CO2 + 10%CF4, вычисленные в MAGBOLTZ: σL – продольная диффузия, σT –
поперечная диффузия.
Видно, что в рабочем диапазоне E = 500÷5000 В/см при положении источника 90Sr на
расстоянии y = 1.2 см от анода вклад диффузии электронов не превысит
Для E = 1164 В/см (HVcath = 4000 В, vD = 2.93 см/мкс) σ D ( t )=
σD( y)
=¿
v
σ D ( y)
√l = 153 мкм.
5.2 нс, что составляет
1.5% от среднего времени дрейфа электронов <T> ≈ 345 нс.
Неопределенность положения начальной ионизации (электронов) относительно анода
Первый тип погрешности σy(t) и наибольший определяется шириной коллиматора Dcoll
D co≪¿
v (E)
∆t =¿
,
σ y ( t )=
∆t
√ 12
(2.4)
где v(E) – скорость дрейфа из Рис. 2.9.
Второй тип погрешности σx(t) определяется шириной отверстия dcath = 1 мм в катоде (4),
при этом σx(t) << σy(t), и ею можно пренебречь.
43
Суммарная ошибка при учете диффузии электронов, как независимого случайного
фактора, составляет
σ (t )= √ σ 2y (t )+σ 2D (t)
(2.5).
Тогда при напряженности поля E = 1164 В/см для двух вариантов коллиматоров имеем2:
Dcoll = 0.5 мм: σ(t) = 7.1 нс (σGarf(t) = 6.4 нс);
Dcoll = 1.0 мм: σ(t) = 10.9 нс (σGarf(t) = 10.4 нс).
Оценки, получаемые из Ур. (2.5), хорошо согласуются с результатами моделирования в
GARFIELD. Несколько завышенное значение σ(t) возникает из-за учета максимальной ошибки,
вносимой диффузией электронов.
Следует отметить, что финальный выбор ширины коллиматоров следует определять
исходя из приемлемого соотношения времени набора, интенсивности источника и требуемого
временного разрешения.
Статистические ошибки
Для вычисления скорости дрейфа используется Ур. (2.1). Выражение для погрешности
имеет вид
∆ v=
где
∆ μ=
σ (μ)
−¿
√N
√(
)(
)(
2
2
1
s
s
∆s +
∆
t
1 +
2
2 ∆ t2
t 2−t 1
( t 2−t 1 )
( t 2−t 1 )
)
2
,
(2.6)
погрешность среднего, соответствующей величины, N – размер выборки, s
= x2-x1. Если принять, что Δs = 0 (рассматривается только статистическая погрешность) и Δ t1 =
Δt2 = Δt, то получим
∆ v=√ 2
s ∆t
2 ∆t
= √2 ∆ v
. (2.7)
2
s
( t 2−t 1 )
Из Ур. (2.7) следует, что погрешность измеряемой скорости дрейфа
∆ v D обратно
пропорциональна расстоянию между положениями источника s и пропорциональна
2
vD .
2 Соответствующие значения σGarf(t) в скобках получены в GARFIELD с учетом диффузии электронов.
44
Оценим необходимый размер выборки N, при котором
∆ v сист = 0.5%. Используя Ур. (2.7) и
∆ μ=
∆ v стат<∆ vсист , пологая
σ (μ)
, получаем
√N
4
N=
Тогда для v = 2.93 см/мкс, σ(t) = 11 нс, s = 1 мм, то
2∆v 2
σ (t ) .
s2 ∆ v 2
(2.8)
∆ v стат < 0.5% достигается при N >
∆ v стат <0.1 , необходимо N > 1.7·105 событий.
7·103 событий, а для уровня
Другие ошибки
Следующие факторы должны быть приняты во внимание при анализе реальных данных с
детектора:
1.
Стабильность высокого напряжения на катоде (3): ΔHVcath ~ ΔE;
2.
Допуски номиналов сопротивлений на электродах с распределенным потенциалом (6), Рис.
2.1;
3.
Временное разрешение ПВК. Из Ур. (2.6) имеем
∆ v=
s
∆t
t2
или
∆ v ∆ t ∆t
= = v .
v
t
s
Если Δt = 100 пкс, s = 1 мм, vD = 3 см/мкс, то Δv/v = 0.3%;
4.
Точность позиционирования коллиматоров Δv ~ Δs (Ур. (2.6), Δt1 = Δt2 = 0);
5.
Разрешение триггерной системы: флуктуации числа фотонов в сцинтилляторах при
регистрации электронов, временное разрешение ФЭУ и электроники;
6.
Стабильность температуры и давления;
7.
Стабильность состава и продува газовой смеси.
2. 7
Оценка толщины сцинтилляторов
В качестве радиоактивного источника электронов взят 90Sr, который распадается
распадается следующим образом:
90
❑
90
90
Sr 0.54 МэВ ❑Y 2.27 МэВ ❑ Zr
→
→
45
В процессе β-распада 90Y образуются электроны с энергией 2.27 МэВ, достаточной для
регистрации этих электронов вторым сцинтиллятором.
В качестве сцинтилляционного материала детекторов рассмотрен вариант пластикового
сцинтиллятора BC-400 [39] с быстрым временем высвечивания 2.4 нс, что является ключевым
параметром для данной установки, поскольку определяет точность определения времени
сигнала «СТАРТ». Этот вид сцинтиллятора производится на основе поливинилтолуола
плотностью ρ = 1.032 г/см3.
Электроны, пройдя первый сцинтиллятор S1, должны достичь второго сцинтиллятора S 2.
Наличие двух сигналов от этих сцинтилляторов подтвердит факт пролета электрона.
Средняя энергия электронов равна
⟨ E ⟩=¿
933.5 кэВ. Произведем оценку толщины
сцинтиллятора, исходя из того, что в первом сцинтилляторе должна поглотиться половина
средней энергии, то есть
⟨ E ⟩ /2 = 466.75 кэВ.
E (E0 = 1 МэВ), МэВ
2.50
2.00
1.50
1.00
0.50
0.00
0.0
0.2
0.4
0.6
0.8
1.0
1.2
1.4
Пробег, г/см2
Рисунок 2.13: График зависимости потерь энергии от пробега электронов. Красная кривая
соответствует электронам с начальной энергией E0 = 935.5 кэВ, синяя кривая – электронам с энергий
E0 = 2,27 МэВ.
Из Рис. 2.13 (красная кривая) видно, что потери составят 466.75 кэВ при толщине
сцинтиллятора R1(см) = R1(г/см2)/ρ(г/см3)
= 0.16 см. Это и будет толщина первого слоя
сцинтиллятора. Толщина второго сцинтиллятора должна соответствовать оставшейся от полной
энергии электронов (синяя кривая). Таким образом, получаем
R2(см) = R(см) - R1(см) = 0.95 см,
где R(см) – пробег электронов с начальной энергией E0 = 2.27 МэВ.
46
Заключение
Газоразрядные детекторы частиц широко используются в самых различных физических
экспериментах, в том числе и на Большом Адронном Коллайдере
(БАК) в CERN. Состав
рабочей газовой смеси для детекторов обычно выбирается исходя из ряда взаимосвязанных
факторов: конструкции, планируемой радиационной загрузки, набора требуемых характеристик
детектора. На сегодняшний день интерес к разработке новых газовых смесей обусловлен ещё и
экологическими требованиями.
Для изучения характеристик новых рабочих газовых смесей, в частности измерения
скорости дрейфа электронов при низких напряженностей полей Е
= 500 ÷ 5000 В/см,
коэффициента газового усиления была разработана компактная дрейфовая в программном
пакете GARFIELD.
Особенностью установки является возможность перемещения источника 90Sr вдоль
модуля с малым шагом 200 мкм ÷ 1 мм, т.е. получения ряда экспериментальных точек и
соответственно более точного фитирование данных, по сравнению с аналогами. Кроме этого,
компактность детектора позволяет поместить его в соленоид небольшого размера для измерения
скорости дрейфа от напряженности магнитного поля.
Проведена оптимизация электрического поля, и определены диапазоны рабочих
положений источника 90Sr, при котором обеспечивается заданная однородность электрического
поля. Построены временные зависимости для газовой смеси 60%Ar + 30%CO2 + 10%CF4 во
всем рабочем диапазоне напряженностей электрического поля. Проведена успешная проверка
алгоритма определения скорости дрейфа для указанной газовой смеси, и определены вклады
основных ошибок.
Проведенная работа будет использована в дальнейшем при разработке прототипа
дрейфового детектора для установки по измерению скорости дрейфа.
Благодарности. Выражаю особую благодарность научному руководителю к. ф.-м. н. Д.
С. Ильину за помощь на всех этапах выполнения диссертации и помощь при оформлении
диссертации, а также к. ф.-м. н. Г.Е. Гаврилову за участие в обсуждении полученных
результатов.
48
Список литературы
1.
M. Abbrescia, A. Colaleo, R. Guida, G. Iaselli, R. Liuzzi, F. Loddo, M. Maggi and B. Marangelli
et al., Gas analysis and monitoring systems for the RPC detector of CMS at LHC,
physics/0701014.
2.
S. Colafranceschi, L. Benussi, S. Bianco, L. Fabbri, M. Giardoni, B. Ortenzi, A. Paolozzi and L.
Passamonti et al., Operational experience of the gas gain monitoring system of the CMS RPC
muon detectors, Nucl. Instrum. Meth. A 617 (2010) 146.
3.
L. Benussi, S. Bianco, S. Colafranceschi, D. Colonna, L. Daniello, F. L. Fabbri, M. Giardoni and
B. Ortenzi et al., The CMS RPC gas gain monitoring system: An Overview and preliminary
results, Nucl. Instrum. Meth. A 602 (2009) 805 [arXiv:0812.1108 [physics.ins-det]].
4.
L. Benussi, S. Bianco, S. Colafranceschi, F. L. Fabbri, M. Giardoni, B. Ortenzi, A. Paolozzi and
L. Passamonti et al., Sensitivity and environmental response of the CMS RPC gas gain
monitoring system, JINST 4 (2009) P08006 [arXiv:0812.1710 [physics.ins-det]].
5.
S. Colafranceschi, L. Benussi, S. Bianco, L. Passamonti, D. Piccolo, D. Pierluigi, A. Russo and
G. Saviano et al., Performance of the Gas Gain Monitoring system of the CMS RPC muon
detector and effective working point fine tuning, JINST 7 (2012) P12004 [PoS RPC 2012 (2012)
046] [arXiv:1209.3893 [hep-ex]].
6.
S. Colafranceschi, R. Aurilio, L. Benussi, S. Bianco, L. Passamonti, D. Piccolo, D. Pierluigi and
A. Russo et al., A study of gas contaminants and interaction with materials in RPC closed loop
systems, PoS RPC 2012 (2012) 056 [JINST 8 (2013) T03008] [arXiv:1210.1819 [physics.insdet]].
7.
L. Benussi, S. Bianco, S. Colafranceschi, F. L. Fabbri, F. Felli, M. Ferrini, M. Giardoni and T.
Greci et al., Study of gas purifiers for the CMS RPC detector, Nucl. Instrum. Meth. A 661 (2012)
S241 [arXiv:1012.5511 [physics.ins-det]].
8.
F. Sauli and A. Sharma, Micropattern gaseous detectors, Ann. Rev. Nucl. Part. Sci. 49 (1999) 341
9.
A. Sharma, Muon tracking and triggering with gaseous detectors and some applications, Nucl.
Instrum. Meth. A 666 (2012) 98.
10.
M. Alfonsi, G. Bencivenni, P. De Simone, F. Murtas, M. Poli Lener, W. Bonivento, A. Cardini
and D. Raspino et al., Aging measurements on triple-GEM detectors operated with CF 4 based gas
mixtures, Nucl. Phys. Proc. Suppl. 150 (2006) 159.
49
11.
I n t e r n a t i o n a l P a n e l o n c l i m a t e c h a n g e s . http://www.ipcc.ch/publications and
data/ar4/wg1/en/ch2s2-10-2.html
12.
R. Santonico and R. Cardarelli, Development of Resistive Plate Counters, Nucl. Instrum. Meth.
187 (1981) 377.
13.
Y. Assran, PhD Thesis 2013, University of Cairo, Egypt.
14.
Y. Assran and A. Sharma, Transport Properties of operational gas mixtures used at LHC,
arXiv:1110.6761 [physics.ins-det].
15.
Properties of potential eco-friendly gas replacements for particle detectors in high-energy
physics. L. Benussi, S. Bianco, D. Piccolo, G. Saviano, S. Colafranceschi, J. Kj_lbro, A. Sharma,
D. Yang, G. Chen, Y. Ban, Q. Li. May 2, 2015.
16.
F. Jensen: Analysis, Characterization and Optimization of an Analogue Optical Link for Data
Readout from a Particle Detector. Ph.D. thesis; 1999
17.
Турчихин, Семен. Мюонные детекторы CMS: дрейфовые трубки, катодные полосковые
камеры, камеры с резистивными пластинами [Электронный ресурс] / С. Турчихин. –
Э л е к т р о н . т е к с т о в ы е д а н . – М о с к в а : М Г У. – Р е ж и м д о с т у п а :
http://nuclphys.sinp.msu.ru/students/md_dt/index.html, свободный.
18.
GHG
Protocol
Global
Wa r m i n g
Potential
Va l u e s ,
2015.
http://www.ghgprotocol.org/files/ghgp/tools/Global-Warming-Potential-Values.pdf.
19.
Eidelman S. et al., (Particle Data Group) Review of Particle Physics: Experimental methods and
colliders // Physics Letters B. 2004. V. 592 (1-4), 1-1109. January 10, 2006. http://pdg.ldl.gov/.
20.
Y. S. Tsai, Pair Production and Bremsstrahlung of Charged Leptons, Rev. Mod. Phys. 46 (1974)
815 [Erratum-ibid. 49 (1977) 521].
21.
I. B. Smirnov, “Modeling of ionization produced by fast charged particles in gases,” Nucl.
Instrum. Meth. A 554 (2005) 474.
22.
Average energy required to produce an ion pair. Technical report, ICRU Report 31, 1979.
23.
L.B. Loeb, Basic processes of gaseous electronics (University of California Press, Berkeley,
1961).
50
24.
S.C. Curran and J.D. Craggs, Counting tubes, theory and applications (Butterworths, London,
1949).
25.
G. Schultz, G. Charpak and F. Sauli, Mobilities of positive ions in some gas mixtures used in
proportional and drift chambers, Rev. Physique Appliquge 12, 67 (1977).
26.
J. Townsend, Electrons in gases (Hutchinson, London, 1947).
27.
H.W. Fulbright, Ionization chambers in nuclear physics, in Encyclopedia of Physics (Ed. S.
Fliigge) (Springer-Verlag, Berlin, 1958), p. 1
28.
F. Sauli, Principles of Operation of Multiwire Proportional and Drift Chambers, CERN-77-09
(1977). [
29.
William Frass C4: Particle Physics Major Option. Particle Detectors. Michaelmas, 2009.
30.
Th. Berghöfer, J. Blümer and J. R. Hörandel, A measurement of drift velocities of electrons in
xenon-methane mixtures, Nucl. Instr. Meth. A 525 (2004) 544.
31.
A.A. Vorobyov, V.S. Guselnikov, A.Yu. Maiorov, Monitoring of electron drift velocity in drift
chambers with alpha-decay recoils, Nucl. Instr. Meth. A 301 (1991) 167.
32.
T. Zhao, Y.Chen, S. Han, J. Hersh, A study of electron drift velocity in Ar–CO2 and Ar–CO2–CF4
gas mixtures, Nucl. Instr. Meth. A 340 (1994) 485.
33.
G. Alkhazov et. al. Forward tracking chamber (FTC) in L3 detector. Part1. General description.
PNPI preprint EP-23-1996, 2115.
34.
R. Veenhof. Garfield - simulation of gaseous detectors. http://garfield.web.cern.ch/garfield/
35.
S. Biagi. Magboltz - transport of electrons in gas mixtures.
http://magboltz.web.cern.ch/magboltz/
36.
I. B. Smirnov. HEED – Interaction of particles with gases. http://heed.web.cern.ch/heed/
37.
G. Gavrilov et. al. Nucl. Instr. Meth. A 356 (1995) 189
38.
V. Andreev et. al. PNPI preprint 2247, 1998
39.
O rg a n i c S c i n t i l l a t i o n M a t e r i a l s . S a i n t - G o b a i n C e r a m i c s & P l a s t i c s , I n c .
http://www.crystals.saint-gobain.com/
51
Отзывы:
Авторизуйтесь, чтобы оставить отзыв