Приборы и техника эксперимента, 2022, № 6, стр. 19-28
ИССЛЕДОВАНИЕ КОЛЛЕКТИВНЫХ ЭФФЕКТОВ ПУЧКА НАКОПИТЕЛЯ-ОХЛАДИТЕЛЯ ИНЖЕКЦИОННОГО КОМПЛЕКСА ВЭПП-5
В. В. Балакин a, b, *, Д. Е. Беркаев a, Ф. А. Еманов a, c
a Институт ядерной физики им. Г.И. Будкера СО РАН
630090 Новосибирск, просп. Академика Лаврентьева, 11, Россия
b Новосибирский государственный технический университет
630073 Новосибирск, просп. Карла Маркса, 20, Россия
c Новосибирский государственный университет
630090 Новосибирск, ул. Пирогова, 1, Россия
* E-mail: balakinvitalyv@gmail.com
Поступила в редакцию 11.03.2022
После доработки 10.06.2022
Принята к публикации 17.06.2022
- EDN: GSTJVK
- DOI: 10.31857/S0032816222060015
Аннотация
Представлены результаты исследования коллективных эффектов пучка накопителя-охладителя инжекционного комплекса ВЭПП-5. Созданные пучком заряженных частиц wake-поля искажают потенциальную яму ускоряющей высокочастотной системы, что приводит к удлинению пучка и искажению формы его продольного распределения. Приведены результаты измерений продольного профиля пучка диссектором и стрик-камерой и теоретически обоснован процесс искажения потенциальной ямы высокочастотной системы. Построена модель импеданса связи в виде эквивалентной RLC-цепи и оценены ее параметры путем сравнения данных моделирования с экспериментальными. Полученная модель импеданса связи была использована для прогнозирования поведения пучка накопителя-охладителя при модификации вакуумной системы ускорителя.
ВВЕДЕНИЕ
Инжекционный комплекс ВЭПП-5 [1, 2] с 2015 года является источником высокоэнергетических электронов и позитронов при проведении фундаментальных экспериментов по физике высоких энергий двумя коллайдерами Института ядерной физики (ИЯФ СО РАН): ВЭПП-2000 [3] и ВЭПП-4М [4]. На данный момент инжекционный комплекс полностью покрывает потребности коллайдеров в пучках электронов и позитронов, однако в дальнейшем планируется его использовать как источник для еще одного, третьего коллайдера, проект которого сейчас разрабатывается в ИЯФ – μ+μ−-трон (Мюмютрон) [5]. Этот комплекс предназначен для исследования димюония − связанного состояния μ+μ− и потребует 0.5−1 ⋅ 1010 позитронов в секунду. От производительности и стабильности работы инжекционного комплекса напрямую зависит эффективность работы коллайдеров, потому разработки, направленные на повышение этих характеристик, являются приоритетными.
ИСКАЖЕНИЯ ПОТЕНЦИАЛЬНОЙ ЯМЫ ВЧ-СИСТЕМЫ
В циклических ускорителях пучок большой интенсивности, взаимодействуя с окружением электромагнитно, создает новые электромагнитные поля, известные как wake-поля. Они действуют на сам пучок, приводя при неблагоприятных обстоятельствах к увеличению амплитуды wake-полей и возникновению коллективных неустойчивостей с последующей потерей частиц.
Wake-поле $W(s)$ определяется как нормированный интеграл силы Лоренца, действующей на пробный заряд, находящийся на расстоянии $s$ за источником этого поля (точечным зарядом $q$):
(1)
$W(s) = - \frac{1}{q}~\mathop \smallint \limits_{ - \infty }^\infty ({\mathbf{E}} + {\mathbf{v}} \times {\mathbf{B}})ds.$Его фурье-образ называется импедансом вакуумной камеры:
(2)
$Z({{\omega }}) = \frac{1}{c}\mathop \smallint \limits_{ - \infty }^\infty W(s)~{{e}^{{i{{\omega }}s/c}}}ds$Для большинства элементов вакуумной камеры, вносящих вклад в импеданс, таких как пикапы, фланцы, элементы, согласующие вакуумные камеры различного поперечного профиля, сильфоны и т.д., продольный импеданс на низких частотах можно приближенно полагать индуктивным, потому удобнее использовать нормализованный продольный импеданс ${{Z}_{{||}}}{\text{/}}n$, где n − номер гармоники частоты обращения.
Также импеданс можно рассматривать в виде эквивалента RLC-цепи [6]:
(3)
$Z({{\omega }}) = \frac{{{{R}_{s}}}}{{1 + iQ\left( {\frac{{{\omega }}}{{{{{{\omega }}}_{r}}}} - \frac{{{{{{\omega }}}_{r}}}}{{{\omega }}}} \right)}},$Такому представлению импеданса соответствует следующая wake-функция:
(4)
$W(s) = \left\{ \begin{gathered} 2\alpha {{R}_{s}}{{e}^{{{{\alpha }}s/c}}}\left( {{\text{cos}}\frac{{{{{{\omega }}}_{p}}s}}{c} + \frac{\alpha }{{{{{{\omega }}}_{p}}}}{\text{sin}}\frac{{{{{{\omega }}}_{p}}s}}{c}} \right),~~s < 0; \hfill \\ \alpha {{R}_{s}},~~s = 0; \hfill \\ \hfill \\ 0,~~s > 0, \hfill \\ \end{gathered} \right.$Электрическое поле, возбуждаемое пучком, с характерным распределением индуктивности L(s) вакуумной камеры вдоль кольца будет иметь вид ($c$ – скорость света) [7]:
(5)
$E(s) = \frac{{dL}}{{ds}}c\frac{{dI}}{{ds}} \approx ~\frac{{dL}}{{ds}}c\frac{{{{I}_{p}}}}{{{{\sigma }}_{z}^{2}}}s.$Интегрируя E(s) по всей длине кольца, можно получить выражение, соответствующее добавочному напряжению ${{V}_{{in}}}$ к ускоряющей высокочастотной (ВЧ) системе:
где $\frac{{{{Z}_{{||}}}}}{n} = {{{{\omega }}}_{0}}L$ – нормализованный продольный импеданс вакуумной камеры.Раскладывая до первого порядка малости выражение для напряжения ВЧ-системы, получим $V({{\varphi }}) = {{V}_{0}}({\text{sin}}({{{{\varphi }}}_{s}}) + {\text{cos}}({{{{\varphi }}}_{s}}){{\varphi }})$, где ${{{{\varphi }}}_{s}}$ – равновесная синхротронная фаза. Учитывая связь продольной координаты s и фазы φ ВЧ-колебаний $\left( {s = - \frac{R}{h}{{\varphi }}} \right)$, получаем полное значение напряжения, действующего на пучок:
(7)
$\begin{gathered} V({{\varphi }}) = {{V}_{0}}{\text{sin}}{{{{\varphi }}}_{s}} + {{V}_{0}}{\text{cos}}{{{{\varphi }}}_{s}} \times \\ \times \left[ {1 + ~\frac{{\sqrt {2\pi } {{I}_{b}}}}{{h{{V}_{0}}{\text{cos}}{{{{\varphi }}}_{s}}}}{\text{Im}}\left( {\frac{{{{Z}_{{||}}}}}{n}} \right){{{\left( {\frac{R}{{{{\sigma }_{z}}}}} \right)}}^{3}}} \right]{{\varphi }}. \\ \end{gathered} $Вышеописанный эффект ведет к уменьшению глубины ВЧ-ямы, в которой находится пучок, что впоследствии приводит также к некогерентному сдвигу синхротронной частоты и росту продольного размера пучка. Учитывая тот факт, что синхротронная частота ${{\omega }}_{s}^{2} \propto \left| {\frac{{dV}}{{d{{\varphi }}}}} \right|$, получаем:
(8)
$\begin{gathered} {{\left( {\frac{{{{{{\omega }}}_{s}}}}{{{{{{\omega }}}_{{s0}}}}}} \right)}^{2}} = {{\left( {\frac{{{{\sigma }_{{z0}}}}}{{{{\sigma }_{z}}}}} \right)}^{2}} = ~1 + ~\frac{{\sqrt {2{{\pi }}} {{I}_{b}}}}{{h{{V}_{0}}{\text{cos}}{{{{\varphi }}}_{s}}}}{\text{Im}}\left( {\frac{{{{Z}_{{||}}}}}{n}} \right){{\left( {\frac{R}{{{{\sigma }_{z}}}}} \right)}^{3}} = \\ = 1 - \frac{e}{{\sqrt {2\pi } }}\frac{{{{\eta }}{{I}_{b}}}}{{EQ_{{s0}}^{2}}}{\text{Im}}\left( {\frac{{{{Z}_{{||}}}}}{n}} \right){{\left( {\frac{R}{{{{\sigma }_{z}}}}} \right)}^{3}}. \\ \end{gathered} $Отсюда получаем кубическое уравнение, описывающее удлинение пучка:
(9)
${{\left( {\frac{{{{\sigma }_{z}}}}{{{{\sigma }_{{z0}}}}}} \right)}^{3}} - \left( {\frac{{{{\sigma }_{z}}}}{{{{\sigma }_{{z0}}}}}} \right) - ~~\frac{e}{{\sqrt {2\pi } }}\frac{{{{\eta }}{{I}_{b}}}}{{EQ_{{s0}}^{2}}}~{\text{Im}}\left( {\frac{{{{Z}_{{||}}}}}{n}} \right){{\left( {\frac{R}{{{{\sigma }_{{z0}}}}}} \right)}^{3}} = 0.$При этом происходит также смещение синхронной фазы пучка на величину [8]
(10)
$\Delta {{\varphi }} = \frac{{2\pi {{I}_{b}}}}{{{{V}_{0}}{\text{cos}}({{{{\varphi }}}_{s}})}}\frac{{\sqrt {2\pi } }}{{{{{{\omega }}}_{0}}\sigma }}{\text{Re}}\left( {\frac{{{{Z}_{{||}}}}}{n}} \right).$Следует отметить, что данное рассмотрение эффекта искажения потенциальной ямы неприменимо в случае очень коротких сгустков (${{\sigma }_{z}} \ll b$, где b – поперечный размер вакуумной камеры), когда вклад высокочастотных резонансных составляющих импеданса становится преобладающим.
Уравнение Власова описывает коллективное поведение системы, состоящей из множества заряженных частиц, под действием электромагнитных полей. Пусть уравнение движения частицы задается обобщенными координатой q и импульсом p:
(11)
$\frac{{\partial \psi }}{{\partial t}} + \frac{{dq}}{{dt}}\frac{{\partial \psi }}{{\partial q}} + \frac{{dp}}{{dt}}\frac{{\partial \psi }}{{\partial p}} = 0.$Стационарное решение задается выражением $\psi = F(H(q,~p,t~))$. Здесь $H$ – гамильтониан системы:
где δ – отклонение энергии частицы от равновесной; ${{\eta }} = {{\alpha }_{p}} - 1{\text{/}}{{\gamma }^{2}}$, ${{\alpha }_{p}}$ – так называемый коэффициент пространственного уплотнения орбит; $z$ – продольная координата частицы; $C$ – периметр синхротрона; γ – гамма-фактор; ${{r}_{0}}$ – классический радиус электрона; $\rho (z)$ – продольное распределение частиц пучка; $W{\kern 1pt} '$ − продольная wake-функция.В случае резистивного импеданса (который определяется в основном ВЧ-резонатором и другими большими полостями) $W{\text{'}}(z) = S\delta (z)$, и тогда решение уравнение Власова дает следующее выражение для продольного распределения частиц пучка:
(13)
$\rho (z) = \frac{{\sqrt {2{\text{/}}\pi } ~{{e}^{{ - \frac{{{{z}^{2}}}}{{2\sigma _{{\text{z}}}^{2}}}}}}}}{{\theta {{\sigma }_{z}}\left[ {\coth \left( {\frac{{{{\theta }}N}}{2}} \right) - {\text{erf}}\left( {\frac{z}{{\sqrt 2 {{\sigma }_{z}}}}} \right)} \right]}},$ИЗМЕРЕНИЕ ПАРАМЕТРОВ ПУЧКА НАКОПИТЕЛЯ-ОХЛАДИТЕЛЯ
С использованием установленного на накопителе-охладителе диссектора (рис. 1) был проведен ряд измерений длины пучка в зависимости от его тока и напряжения ВЧ-резонатора. Для выявления мнимой части импеданса была проанализирована зависимость длины пучка от его тока (рис. 2). Параметры диссектора приведены ниже:
− Напряжение фотокатода, кВ | 12 |
− Напряжение сетки, кВ | 10 |
− Напряжение фокусирующего электрода, кВ | 9 |
− Напряжение вторичного электронного умножителя, кВ | 2 |
− Частота ВЧ-развертки, МГц | 78.6 |
− Разрешение диссектора, пс | ~20 |
Используя измеренные данные и проведя их аппроксимацию выражением (9), можно получить значение мнимой части импеданса связи накопителя, равное ${\text{Im}}\left( {\frac{{{{Z}_{{||}}}}}{n}} \right) = 5.71~\, \pm \,{\text{0}}{\text{.5}}$ Ом.
По проекту накопителя-охладителя начальная оценка этого импеданса должна была составлять не более 1 Ом [9]. Однако этот расчет был проведен со старым медным резонатором 64-й гармоники. В новых же реалиях комплекс работает с ферритовым резонатором 1-й гармоники (с физическими габаритами, сопоставимыми со старым резонатором) [2], что, по-видимому, и дает значительно большую величину импеданса.
Была также измерена зависимость сдвига синхронной фазы пучка от его тока (рис. 3), из которой, исходя из формулы (10), было получено значение действительной части импеданса связи ${\text{Re}}\left( {\frac{{{{Z}_{{||}}}}}{n}} \right)~ = {\text{15}}{\text{.68\;}} \pm 0.47~$ Ом.
При относительно низких токах продольное распределение частиц хорошо описывается нормальным распределением (рис. 4а), среднеквадратичное отклонение которого задает характерный размер сгустка. При токе порядка 50 мА начинает происходить искажение формы пучка, и его аппроксимация функцией Гаусса становится некорректной (рис. 4б, 4в). Для аппроксимации данных при искажении профиля пучка необходимо использовать выражение (13). Программное обеспечение диссектора позволяет это реализовать (рис. 4г).
Аналогичные измерения были проведены стрик-камерой PS-1/S1, ее параметры приведены ниже:
− Ширина аппаратной функции стрик-камеры, пс | 10 |
− Спектральный диапазон фотокатода, нм | 350−1300 |
− Диапазоны длительности развертки (на экран 25 мм), нс | 0.25, 1, 5, 20, 100, 200 |
− Динамическое пространственное раз-решение (вдоль щели), пар линий/мм | 10 |
− Временное разрешение для трубки с фотокатодом, пс | <1 |
− Максимальный динамический диапазон регистрации | 100 |
− Временная нелинейность разверток, % | ≤10 |
− Задержка срабатывания, нс | ≤15 |
− Нестабильность срабатывания, нс | 0.12 (развертка 10 нс) |
− Запускающий импульс (положительный), В | 1−5 |
− Частота запуска, кГц | 1 |
− Питание от сети | 100−240 В; 50−60 Гц |
Результаты измерений профиля пучка в точках, соответствующих току в диапазоне 2.2−45 мА, приведены на рис. 5. Как было сказано выше, с учетом искажения потенциальной ямы это распределение должно аппроксимироваться моделью (9). По результатам аппроксимации значение мнимой части импеданса составило ${\text{Im}}\left( {\frac{{{{Z}_{{||}}}}}{n}} \right)$ = 6.15 ± 0.18 Ом.
Результаты измерений с помощью стрик-камеры также отчетливо демонстрируют изменение продольного профиля пучка, вызванное искажением потенциальной ямы, наведенным пучком wake-полей. Измерения также находятся в согласии с ранее предложенной моделью (13). Демонстрация этого эффекта показана на рис. 6.
Чтобы исследовать процесс захвата сгустков частиц, вышедших из линейного ускорителя инжекционного комплекса в ВЧ-систему накопителя-охладителя, был проведен ряд измерений стрик-камерой продольного профиля пучка на впуске. Инжектируемый из линейного ускорителя пучок содержит 16 сгустков (рис. 7), сформировавшихся во время ускорения. Попадая в потенциальную яму, образованную ускоряющим ВЧ-резонатором первой гармоники частоты обращения, пучок начинает перегруппировываться, образуя один сгусток. Как показал эксперимент, характерное время перегруппировки составляет около 270 мкс (≈3000 оборотов), что значительно меньше времени радиационного затухания пучка, равного 18 мс.
Для изучения влияния импеданса связи на продольную перегруппировку пучка в течение первой тысячи оборотов после инжекции было выполнено моделирование данного процесса. Программный код был реализован на языке программирования Python 3. В данной модели продольный импеданс связи накопителя-охладителя задавался как эквивалентная RLC-цепь, согласно (3), с соответствующей импедансу wake-функцией (4). Для оценки параметров эквивалентной цепи проводилась минимизация целевой функции φ, задаваемой как:
(14)
${{\varphi }} = \mathop \sum \limits_{}^N {{({{I}_{{model}}}(N,{\text{\;}}{{R}_{s}},{\text{\;}}Q,{\text{\;}}{{{{\omega }}}_{r}}) - {{I}_{{mes}}}(N))}^{2}},$Рис. 8 иллюстрирует соответствие измеренного стрик-камерой профиля пучка результату моделирования при наличии импеданса связи. Очевидно, что импеданс кольца и связанные с ним коллективные эффекты оказывают заметное влияние на продольную динамику пучка в кольце. Причина такого быстрого процесса перегруппировки пучка заключается в возникновении микроволновой неустойчивости, которая в фазовом пространстве смешивает между собой отдельные сгустки. Проявление микроволновой неустойчивости особенно четко демонстрирует рис. 9. Ее характерной особенностью является появление на профиле пучка модуляции плотности, которая уже никак не связана с начальным продольным распределением частиц, полученным при инжекции из линейного ускорителя.
Оцененный по параметрам накопителя-охладителя порог микроволновой неустойчивости составляет 1.62 А (средний ток пучка 20 мА при длине пучка при нулевом токе 13.5 см). Физически в накопителе пиковый ток охлажденного пучка составляет 1.15 А при среднем значении тока пучка 20 мА (ввиду того что длина пучка 19 см при заданном токе). При этом порог не достигается, и микроволновая неустойчивость себя не проявляет. Пиковый ток сгустков пучка линейного ускорителя превышает значение пикового тока порога микроволновой неустойчивости (уже при среднем токе одного сгустка 0.44 мА и длине пучка 10 пс), поэтому динамика продольного профиля пучка после инжекции определяется этой неустойчивостью.
МОДЕЛИРОВАНИЕ ИЗМЕНЕНИЯ ИМПЕДАНСА СВЯЗИ НАКОПИТЕЛЯ-ОХЛАДИТЕЛЯ
Для исследования димюония (связанного состояния μ+μ−) [5] в Институте ядерной физики разрабатывается проект коллайдера μ+μ−-трон (Мюмютрон). Источником частиц для него станет инжекционный комплекс.
Однако накопитель-охладитель работает с резонатором 1-й гармоники (10.94 МГц) и значительно более длинным пучком, чем предполагается использовать в μ+μ−-троне. Как итог, перепуск такого пучка в коллайдер с высокой частотой 338.98 МГц [4] будет приводить к значительным потерям частиц в момент захвата пучка коллайдером.
Один из вариантов решения данной проблемы – установка второго резонатора в накопитель для предварительного сжатия (перегруппировки) пучка в самом накопителе-охладителе перед его транспортировкой в коллайдер.
Подходящим для установки в накопитель является резонатор, аналогичный тому, что разрабатывается для проекта сибирского кольцевого источника фотонов (СКИФ), и имеющий следующие характеристики: частота 357 МГц, добротность 17000, шунтовое сопротивление 2.55 МОм, ускоряющее напряжение 400 кВ. Для исследования влияния на динамику пучка при инжекции дополнительного резонатора его параметры были добавлены в уже построенную модель импеданса связи.
Согласно приведенным ниже результатам, вставка дополнительного резонатора практически не меняет динамику пучка при инжекции – все так же возникает микроволновая неустойчивость на ранних оборотах, которая перегруппировывает пучок из 16-ти сгустков в один длинный.
На рис. 10 представлены характерные профили тока пучка и его распределения в фазовом пространстве для случая без импеданса дополнительного резонатора и с ним на 200-м и 3400-м оборотах.
Рис. 11 демонстрирует дальнейшее затухание пучка с дополнительным резонатором на 40 000-м и 200 000-м оборотах. Из этих распределений видно, что за время, соответствующее характерному времени радиационного затухания пучка, происходит его группировка в имеющейся сепаратрисе.
Далее, затухший пучок необходимо перегруппировать путем отключения резонатора 10.94 МГц, в который он инжектировался, и включением нового резонатора 357 МГц. Если резко выключить первую ВЧ-станцию и тут же запустить вторую, то перегруппировка пучка будет сопровождаться потерями. Чтобы избежать этого, было предложено постепенно увеличивать ускоряющее напряжение в резонаторе 357 МГц до 400 кВ и после достижения этой величины начать уменьшать поле в резонаторе 10.94 МГц до нуля. Таким способом удается осуществить медленную перегруппировку пучка без возникновения каких-либо существенных потерь, связанных с резким исчезновением потенциальной ямы, удерживающей накопленный пучок, и появлением новой от резонатора 357 МГц. Процесс такой перегруппировки пучка иллюстрирует рис. 12.
ЗАКЛЮЧЕНИЕ
В рамках данной работы был изучен феномен удлинения пучка накопителя-охладителя и искажения его формы. Показано, что причиной этого является искажение потенциальной ямы ВЧ-системы накопителя-охладителя. По результатам обработки измерений продольного профиля пучка стрик-камерой и диссектором была определена мнимая часть импеданса связи накопителя-охладителя, а с помощью измерений сдвига синхронной фазы пучка диссектором − действительная его часть.
При проведении измерения продольного профиля пучка при его инжекции из линейного ускорителя в накопитель было обнаружено, что слияние 16-ти сгустков, созданных ВЧ-системой линейного ускорителя, происходит значительно быстрее характерного времени радиационного затухания пучка. Моделирование захвата пучка позволило не только определить, что причиной тому является возникновение микроволновой неустойчивости, но и оценить параметры эквивалентной RLC-цепи для представления импеданса накопителя-охладителя.
В дальнейшем построенная модель была использована для исследования поведения пучка во время инжекции при модификации вакуумной камеры накопителя путем вставки дополнительного резонатора. Было показано, что дополнительный элемент вакуумной системы не приводит к возникновению неустойчивостей при инжекции и охлаждении пучка. Также моделирование подтвердило гипотезу сохранения частиц пучка при его перегруппировке из основного резонатора в дополнительный путем медленного отключения первой ВЧ-станции и включения второй.
Список литературы
Astrelina K.V., Blinov M.F., Vsevolozhskaya T.A., Dikanskii N.S., Emanov F.A., Lapik R.M., Logachev P.V., Martyshkin P.V., Petrenko A.V., Rybitskaya T.V., Skrinskii A.N., Shiyankov S.V., Yaskina T.A. // Journal of Experimental and Theoretical Physics. 2008. V. 106. № 1. P. 77. https://doi.org/10.1134/S1063776108010068
Maltseva Yu., Andrianov A., Astrelina K., Balakin V., Batrakov A., Belikov O., Berkaev D., Blinov M., Bolkhovityanov D., Butakov A., Bykov E., Dikansky N., Emanov F., Frolov A., Gambaryan V. et al. // Proc. IPAC’2018. Vancouver, Canada, 2018. Paper MOPMK011. https://doi.org/10.18429/JACoW-IPAC2018-MOPMK011
Shatunov P., Belikov O., Berkaev D., Gorchakov K., Kasaev A., Kirpotin A., Koop I., Krasnov A., Lysenko A., Motygin S., Prosvetov V., Rabusov D., Semenov A., Shatunov Yu., Timoshenko M. et al. // Proc. IPAC’17, Copenhagen, Denmark, 2017. Paper WEPIK02. https://doi.org/10.18429/JACoW-IPAC2017-WEPIK029
Журавлев А.Н., Карнаев С.Е., Левичев Е.Б., Мешков О.И., Никитин С.А., Николаев И.Б., Пиминов П.А. // Письма в ЭЧАЯ. 2020. Т. 17. № 7. С. 876.
Bogomyagkov A., Druzhinin V., Levichev E., Milstein A., Sinyatkin S. // EPJ Web of Conferences. 2018. V. 181. P. 01032.https://doi.org/10.1051/epjconf/201818101032
Chao A. Physics of Collective Beam Instabilities in High Energy Accelerators. Dallas: Wiley-Interscience Publication “John Wiley & Sons, INC”, 1993.
Hofmann A., Maidment J. Current dependent phenomena in LEP. LEP note 168. (28.6.79)
Chao A. Handbook of Accelerator Physics and Engineering. New Jersey: World Scientific Publishing Co. Pre. Ltd., 2013.
Анашин В. Препринт ИЯФ 92-44. Новосибирск, 1992.
Дополнительные материалы отсутствуют.
Инструменты
Приборы и техника эксперимента